Sebastian Schunert; Yousry Y. Azmy; Damien Fournier
2011-05-01
We present a comprehensive error estimation of four spatial discretization schemes of the two-dimensional Discrete Ordinates (SN) equations on Cartesian grids utilizing a Method of Manufactured Solution (MMS) benchmark suite based on variants of Larsen’s benchmark featuring different orders of smoothness of the underlying exact solution. The considered spatial discretization schemes include the arbitrarily high order transport methods of the nodal (AHOTN) and characteristic (AHOTC) types, the discontinuous Galerkin Finite Element method (DGFEM) and the recently proposed higher order diamond difference method (HODD) of spatial expansion orders 0 through 3. While AHOTN and AHOTC rely on approximate analytical solutions of the transport equation within a mesh cell, DGFEM and HODD utilize a polynomial expansion to mimick the angular flux profile across each mesh cell. Intuitively, due to the higher degree of analyticity, we expect AHOTN and AHOTC to feature superior accuracy compared with DGFEM and HODD, but at the price of potentially longer grind times and numerical instabilities. The latter disadvantages can result from the presence of exponential terms evaluated at the cell optical thickness that arise from the semianalytical solution process. This work quantifies the order of accuracy and the magnitude of the error of all four discretization methods for different optical thicknesses, scattering ratios and degrees of smoothness of the underlying exact solutions in order to verify or contradict the aforementioned intuitive expectation.
NASA Astrophysics Data System (ADS)
Maginot, Peter G.; Morel, Jim E.; Ragusa, Jean C.
2012-08-01
We present a new nonlinear spatial finite-element method for the linearized Boltzmann transport equation with Sn angular discretization in 1-D and 2-D Cartesian geometries. This method has two central characteristics. First, it is equivalent to the linear-discontinuous (LD) Galerkin method whenever that method yields a strictly non-negative solution. Second, it always satisfies both the zeroth and first spatial moment equations. Because it yields the LD solution when that solution is non-negative, one might interpret our method as a classical fix-up to the LD scheme. However, fix-up schemes for the LD equations derived in the past have given up solution of the first moment equations when the LD solution is negative in order to satisfy positivity in a simple manner. We present computational results comparing our method in 1-D to the strictly non-negative linear exponential-discontinuous method and to the LD method. We present computational results in 2-D comparing our method to a recently developed LD fix-up scheme and to the LD scheme. It is demonstrated that our method is a valuable alternative to existing methods.
Triangle geometry processing for surface modeling and cartesian grid generation
Aftosmis, Michael J [San Mateo, CA; Melton, John E [Hollister, CA; Berger, Marsha J [New York, NY
2002-09-03
Cartesian mesh generation is accomplished for component based geometries, by intersecting components subject to mesh generation to extract wetted surfaces with a geometry engine using adaptive precision arithmetic in a system which automatically breaks ties with respect to geometric degeneracies. During volume mesh generation, intersected surface triangulations are received to enable mesh generation with cell division of an initially coarse grid. The hexagonal cells are resolved, preserving the ability to directionally divide cells which are locally well aligned.
Triangle Geometry Processing for Surface Modeling and Cartesian Grid Generation
NASA Technical Reports Server (NTRS)
Aftosmis, Michael J. (Inventor); Melton, John E. (Inventor); Berger, Marsha J. (Inventor)
2002-01-01
Cartesian mesh generation is accomplished for component based geometries, by intersecting components subject to mesh generation to extract wetted surfaces with a geometry engine using adaptive precision arithmetic in a system which automatically breaks ties with respect to geometric degeneracies. During volume mesh generation, intersected surface triangulations are received to enable mesh generation with cell division of an initially coarse grid. The hexagonal cells are resolved, preserving the ability to directionally divide cells which are locally well aligned.
Onset of buoyancy-driven convection in Cartesian and cylindrical geometries
NASA Astrophysics Data System (ADS)
Myint, Philip C.; Firoozabadi, Abbas
2013-04-01
We perform a linear stability analysis to examine the onset of buoyancy-driven convection relevant to subsurface carbon dioxide sequestration in confined, porous Cartesian and cylindrical domains. Our work amends the analysis in an earlier study on cylindrical geometries. We consider Cartesian geometries where the aspect ratio between the two horizontal dimensions is not necessarily equal to one. Two key elements of the stability analysis are: (1) the critical time and (2) the critical wavenumber. Lateral boundaries have a much greater influence on the critical wavenumber than on the critical time. The confinement due to these boundaries impedes the onset of convection to the extent that convection cannot even occur in domains that are smaller than a certain size. Large aspect ratios can significantly reduce boundary effects. Patterns of the earliest-growing perturbation mode in the horizontal plane reveal many interesting dynamics which have not been examined in previous stability analyses. We illustrate several differences between patterns in Cartesian geometries and patterns in cylindrical geometries. Based on observations from earlier papers, we hypothesize that the contrasts between the Cartesian and cylindrical patterns may lead to significantly different behavior in the two geometries after the onset of convection. Our results may guide future numerical studies that can investigate this hypothesis and may help with understanding the onset of buoyancy-driven convection in real systems where lateral boundary effects are significant.
Cartesian grid simulations of gas-solids flow systems with complex geometry
Dietiker, Jean-Francois; Li, Tingwen; Garg, Rahul; Shahnam, Mehrdad
2013-02-01
Complex geometries encountered in many applications of gas–solids flow need special treatment in most legacy multiphase flow solvers with Cartesian numerical grid. This paper briefly outlines the implementation of a cut cell technique in the open-source multiphase flow solver—MFIX for accurate representation of complex geometries. Specifically, applications of the Cartesian cut cell method to different gas–solids fluidization systems including a small scale bubbling fluidized bed with submerged tube bundle and a complete pilot-scale circulating fluidized bed will be presented. In addition to qualitative predictions on the general flow behaviors inside each system, quantitative comparison with the available experimental data will be presented. Furthermore, some results on extending the current cut-cell technique to Lagrangian–Eulerian simulations will be presented.
Aerodynamic Design of Complex Configurations Using Cartesian Methods and CAD Geometry
NASA Technical Reports Server (NTRS)
Nemec, Marian; Aftosmis, Michael J.; Pulliam, Thomas H.
2003-01-01
The objective for this paper is to present the development of an optimization capability for the Cartesian inviscid-flow analysis package of Aftosmis et al. We evaluate and characterize the following modules within the new optimization framework: (1) A component-based geometry parameterization approach using a CAD solid representation and the CAPRI interface. (2) The use of Cartesian methods in the development Optimization techniques using a genetic algorithm. The discussion and investigations focus on several real world problems of the optimization process. We examine the architectural issues associated with the deployment of a CAD-based design approach in a heterogeneous parallel computing environment that contains both CAD workstations and dedicated compute nodes. In addition, we study the influence of noise on the performance of optimization techniques, and the overall efficiency of the optimization process for aerodynamic design of complex three-dimensional configurations. of automated optimization tools. rithm and a gradient-based algorithm.
Geometry optimization for peptides and proteins: comparison of Cartesian and internal coordinates.
Koslover, Elena F; Wales, David J
2007-12-21
We present the results of several benchmarks comparing the relative efficiency of different coordinate systems in optimizing polypeptide geometries. Cartesian, natural internal, and primitive internal coordinates are employed in quasi-Newton minimization for a variety of biomolecules. The peptides and proteins used in these benchmarks range in size from 16 to 999 residues. They vary in complexity from polyalanine helices to a beta-barrel enzyme. We find that the relative performance of the different coordinate systems depends on the parameters of the optimization method, the starting point for the optimization, and the size of the system studied. In general, internal coordinates were found to be advantageous for small peptides. For larger structures, Cartesians appear to be more efficient for empirical potentials where the energy and gradient can be evaluated relatively quickly compared to the cost of the coordinate transformations. PMID:18154373
Geometry optimization for peptides and proteins: Comparison of Cartesian and internal coordinates
NASA Astrophysics Data System (ADS)
Koslover, Elena F.; Wales, David J.
2007-12-01
We present the results of several benchmarks comparing the relative efficiency of different coordinate systems in optimizing polypeptide geometries. Cartesian, natural internal, and primitive internal coordinates are employed in quasi-Newton minimization for a variety of biomolecules. The peptides and proteins used in these benchmarks range in size from 16 to 999 residues. They vary in complexity from polyalanine helices to a β-barrel enzyme. We find that the relative performance of the different coordinate systems depends on the parameters of the optimization method, the starting point for the optimization, and the size of the system studied. In general, internal coordinates were found to be advantageous for small peptides. For larger structures, Cartesians appear to be more efficient for empirical potentials where the energy and gradient can be evaluated relatively quickly compared to the cost of the coordinate transformations.
GBL-2D Version 1.0: a 2D geometry boolean library.
McBride, Cory L. (Elemental Technologies, American Fort, UT); Schmidt, Rodney Cannon; Yarberry, Victor R.; Meyers, Ray J.
2006-11-01
This report describes version 1.0 of GBL-2D, a geometric Boolean library for 2D objects. The library is written in C++ and consists of a set of classes and routines. The classes primarily represent geometric data and relationships. Classes are provided for 2D points, lines, arcs, edge uses, loops, surfaces and mask sets. The routines contain algorithms for geometric Boolean operations and utility functions. Routines are provided that incorporate the Boolean operations: Union(OR), XOR, Intersection and Difference. A variety of additional analytical geometry routines and routines for importing and exporting the data in various file formats are also provided. The GBL-2D library was originally developed as a geometric modeling engine for use with a separate software tool, called SummitView [1], that manipulates the 2D mask sets created by designers of Micro-Electro-Mechanical Systems (MEMS). However, many other practical applications for this type of software can be envisioned because the need to perform 2D Boolean operations can arise in many contexts.
Alj, Domenico; Caputo, Roberto; Umeton, Cesare
2014-11-01
We report on the realization of a liquid crystal (LC)-based optical diffraction grating showing a polar symmetry of the director alignment. This has been obtained as a natural evolution of the POLICRYPS technique, which enables the realization of highly efficient, switchable, planar diffraction gratings. Performances exhibited in the Cartesian geometry are extended to the polar one by exploiting the spherical aberration produced by simple optical elements. This enables producing the required highly stable polar pattern that allows fabricating a circular optical diffraction grating. Results are promising for their possible application in fields in which a rotational structure of the optical beam is needed. PMID:25361314
Adjoint Sensitivity Computations for an Embedded-Boundary Cartesian Mesh Method and CAD Geometry
NASA Technical Reports Server (NTRS)
Nemec, Marian; Aftosmis,Michael J.
2006-01-01
Cartesian-mesh methods are perhaps the most promising approach for addressing the issues of flow solution automation for aerodynamic design problems. In these methods, the discretization of the wetted surface is decoupled from that of the volume mesh. This not only enables fast and robust mesh generation for geometry of arbitrary complexity, but also facilitates access to geometry modeling and manipulation using parametric Computer-Aided Design (CAD) tools. Our goal is to combine the automation capabilities of Cartesian methods with an eficient computation of design sensitivities. We address this issue using the adjoint method, where the computational cost of the design sensitivities, or objective function gradients, is esseutially indepeudent of the number of design variables. In previous work, we presented an accurate and efficient algorithm for the solution of the adjoint Euler equations discretized on Cartesian meshes with embedded, cut-cell boundaries. Novel aspects of the algorithm included the computation of surface shape sensitivities for triangulations based on parametric-CAD models and the linearization of the coupling between the surface triangulation and the cut-cells. The objective of the present work is to extend our adjoint formulation to problems involving general shape changes. Central to this development is the computation of volume-mesh sensitivities to obtain a reliable approximation of the objective finction gradient. Motivated by the success of mesh-perturbation schemes commonly used in body-fitted unstructured formulations, we propose an approach based on a local linearization of a mesh-perturbation scheme similar to the spring analogy. This approach circumvents most of the difficulties that arise due to non-smooth changes in the cut-cell layer as the boundary shape evolves and provides a consistent approximation tot he exact gradient of the discretized abjective function. A detailed gradient accurace study is presented to verify our approach
Facial biometrics based on 2D vector geometry
NASA Astrophysics Data System (ADS)
Malek, Obaidul; Venetsanopoulos, Anastasios; Androutsos, Dimitrios
2014-05-01
The main challenge of facial biometrics is its robustness and ability to adapt to changes in position orientation, facial expression, and illumination effects. This research addresses the predominant deficiencies in this regard and systematically investigates a facial authentication system in the Euclidean domain. In the proposed method, Euclidean geometry in 2D vector space is being constructed for features extraction and the authentication method. In particular, each assigned point of the candidates' biometric features is considered to be a 2D geometrical coordinate in the Euclidean vector space. Algebraic shapes of the extracted candidate features are also computed and compared. The proposed authentication method is being tested on images from the public "Put Face Database". The performance of the proposed method is evaluated based on Correct Recognition (CRR), False Acceptance (FAR), and False Rejection (FRR) rates. The theoretical foundation of the proposed method along with the experimental results are also presented in this paper. The experimental results demonstrate the effectiveness of the proposed method.
2D Colloidal Wigner crystals in confined geometries
NASA Astrophysics Data System (ADS)
Higler, Ruben; Sprakel, Joris
2015-03-01
Crystallization of bulk systems has been widely studied using colloids as a model system. However, study into colloidal crystallization in confined geometries has been sparse and little is known about the effects of strong confinement on the dynamics of colloidal crystal. In our research we prepare 2D crystals from charged colloids in an apolar solvent to study crystal dynamics, formation, and structure in circular confinements. These confining geometries are made using softlithography techniques from SU-8. In order to broaden the parameter space we can reach in experiments we employ brownian dynamics simulations to supplement our experimental results. Using single-particle tracking we have subpixel resolution positional information of every particle in the system. We study the vibrational modes of our confined crystals and find well defined modes unique to confined systems, such as a radially symmetric compression (or breathing) mode, a collective rotation mode, and distinct resonance modes. Furthermore, due to the circular nature of our constrictions, defectless crystals are impossible, we find, for sufficiently high area fractions, that the defects order at well defined points at the edge. The effect of this ``defect-localization'' has a clear influence on the vibrational modes.
A Parallel Cartesian Approach for External Aerodynamics of Vehicles with Complex Geometry
NASA Technical Reports Server (NTRS)
Aftosmis, M. J.; Berger, M. J.; Adomavicius, G.
2001-01-01
This workshop paper presents the current status in the development of a new approach for the solution of the Euler equations on Cartesian meshes with embedded boundaries in three dimensions on distributed and shared memory architectures. The approach uses adaptively refined Cartesian hexahedra to fill the computational domain. Where these cells intersect the geometry, they are cut by the boundary into arbitrarily shaped polyhedra which receive special treatment by the solver. The presentation documents a newly developed multilevel upwind solver based on a flexible domain-decomposition strategy. One novel aspect of the work is its use of space-filling curves (SFC) for memory efficient on-the-fly parallelization, dynamic re-partitioning and automatic coarse mesh generation. Within each subdomain the approach employs a variety reordering techniques so that relevant data are on the same page in memory permitting high-performance on cache-based processors. Details of the on-the-fly SFC based partitioning are presented as are construction rules for the automatic coarse mesh generation. After describing the approach, the paper uses model problems and 3- D configurations to both verify and validate the solver. The model problems demonstrate that second-order accuracy is maintained despite the presence of the irregular cut-cells in the mesh. In addition, it examines both parallel efficiency and convergence behavior. These investigations demonstrate a parallel speed-up in excess of 28 on 32 processors of an SGI Origin 2000 system and confirm that mesh partitioning has no effect on convergence behavior.
Aerodynamic Optimization of Rocket Control Surface Geometry Using Cartesian Methods and CAD Geometry
NASA Technical Reports Server (NTRS)
Nelson, Andrea; Aftosmis, Michael J.; Nemec, Marian; Pulliam, Thomas H.
2004-01-01
Aerodynamic design is an iterative process involving geometry manipulation and complex computational analysis subject to physical constraints and aerodynamic objectives. A design cycle consists of first establishing the performance of a baseline design, which is usually created with low-fidelity engineering tools, and then progressively optimizing the design to maximize its performance. Optimization techniques have evolved from relying exclusively on designer intuition and insight in traditional trial and error methods, to sophisticated local and global search methods. Recent attempts at automating the search through a large design space with formal optimization methods include both database driven and direct evaluation schemes. Databases are being used in conjunction with surrogate and neural network models as a basis on which to run optimization algorithms. Optimization algorithms are also being driven by the direct evaluation of objectives and constraints using high-fidelity simulations. Surrogate methods use data points obtained from simulations, and possibly gradients evaluated at the data points, to create mathematical approximations of a database. Neural network models work in a similar fashion, using a number of high-fidelity database calculations as training iterations to create a database model. Optimal designs are obtained by coupling an optimization algorithm to the database model. Evaluation of the current best design then gives either a new local optima and/or increases the fidelity of the approximation model for the next iteration. Surrogate methods have also been developed that iterate on the selection of data points to decrease the uncertainty of the approximation model prior to searching for an optimal design. The database approximation models for each of these cases, however, become computationally expensive with increase in dimensionality. Thus the method of using optimization algorithms to search a database model becomes problematic as the
Reentry-Vehicle Shape Optimization Using a Cartesian Adjoint Method and CAD Geometry
NASA Technical Reports Server (NTRS)
Nemec, Marian; Aftosmis, Michael J.
2006-01-01
A DJOINT solutions of the governing flow equations are becoming increasingly important for the development of efficient analysis and optimization algorithms. A well-known use of the adjoint method is gradient-based shape. Given an objective function that defines some measure of performance, such as the lift and drag functionals, its gradient is computed at a cost that is essentially independent of the number of design variables (e.g., geometric parameters that control the shape). Classic aerodynamic applications of gradient-based optimization include the design of cruise configurations for transonic and supersonic flow, as well as the design of high-lift systems. are perhaps the most promising approach for addressing the issues of flow solution automation for aerodynamic design problems. In these methods, the discretization of the wetted surface is decoupled from that of the volume mesh. This not only enables fast and robust mesh generation for geometry of arbitrary complexity, but also facilitates access to geometry modeling and manipulation using parametric computer-aided design (CAD). In previous work on Cartesian adjoint solvers, Melvin et al. developed an adjoint formulation for the TRANAIR code, which is based on the full-potential equation with viscous corrections. More recently, Dadone and Grossman presented an adjoint formulation for the two-dimensional Euler equations using a ghost-cell method to enforce the wall boundary conditions. In Refs. 18 and 19, we presented an accurate and efficient algorithm for the solution of the adjoint Euler equations discretized on Cartesian meshes with embedded, cut-cell boundaries. Novel aspects of the algorithm were the computation of surface shape sensitivities for triangulations based on parametric-CAD models and the linearization of the coupling between the surface triangulation and the cut-cells. The accuracy of the gradient computation was verified using several three-dimensional test cases, which included design
Development of a new two-dimensional Cartesian geometry nodal multigroup discrete-ordinates method
Pevey, R.E.
1982-07-01
The purpose of this work is the development and testing of a new family of methods for calculating the spatial dependence of the neutron density in nuclear systems described in two-dimensional Cartesian geometry. The energy and angular dependence of the neutron density is approximated using the multigroup and discrete ordinates techniques, respectively. The resulting FORTRAN computer code is designed to handle an arbitrary number of spatial, energy, and angle subdivisions. Any degree of scattering anisotropy can be handled by the code for either external source or fission systems. The basic approach is to (1) approximate the spatial variation of the neutron source across each spatial subdivision as an expansion in terms of a user-supplied set of exponential basis functions; (2) solve analytically for the resulting neutron density inside each region; and (3) approximate this density in the basis function space in order to calculate the next iteration flux-dependent source terms. In the general case the calculation is iterative due to neutron sources which depend on the neutron density itself, such as scattering interactions.
Quantum geometry of 2D gravity coupled to unitary matter
NASA Astrophysics Data System (ADS)
Ambjørn, J.; Anagnostopoulos, K. N.
1997-02-01
We show that there exists a divergent correlation length in 2D quantum gravity for the matter fields close to the critical point provided one uses the invariant geodesic distance as the measure of distance. The corresponding reparameterization invariant two-point functions satisfy all scaling relations known from the ordinary theory of critical phenomena and the KPZ exponents are determined by the power-like fall-off of these two-point functions. The only difference compared to flat space is the appearance of a dynamically generated fractal dimension d h in the scaling relations. We analyze numerically the fractal properties of space-time for the Ising and three-states Potts model coupled to two-dimensional quantum gravity using finite size scaling as well as small distance scaling of invariant correlation functions. Our data are consistent with dh = 4, but we cannot rule out completely the conjecture dH = -2 α1/ α-1, where α- n is the gravitational dressing exponent of a spinless primary field of conformal weight ( n + 1, n + 1). We compute the moments < L> and the loop-length distribution function and show that the fractal properties associated with these observables are identical, with good accuracy, to the pure gravity case.
Boundary treatments for 2D elliptic mesh generation in complex geometries
Technology Transfer Automated Retrieval System (TEKTRAN)
This paper presents a boundary treatment method for 2D elliptic mesh generation in complex geometries. Corresponding to Neumann- Dirichlet boundary conditions (sliding boundary conditions), the proposed method aims at achieving orthogonal and smooth nodal distribution along irregular boundaries. In ...
NASA Astrophysics Data System (ADS)
Polkowski, Marcin
2016-04-01
Seismic wave travel time calculation is the most common numerical operation in seismology. The most efficient is travel time calculation in 1D velocity model - for given source, receiver depths and angular distance time is calculated within fraction of a second. Unfortunately, in most cases 1D is not enough to encounter differentiating local and regional structures. Whenever possible travel time through 3D velocity model has to be calculated. It can be achieved using ray calculation or time propagation in space. While single ray path calculation is quick it is complicated to find the ray path that connects source with the receiver. Time propagation in space using Fast Marching Method seems more efficient in most cases, especially when there are multiple receivers. In this presentation a Python module pySeismicFMM is presented - simple and very efficient tool for calculating travel time from sources to receivers. Calculation requires regular 2D or 3D velocity grid either in Cartesian or geographic coordinates. On desktop class computer calculation speed is 200k grid cells per second. Calculation has to be performed once for every source location and provides travel time to all receivers. pySeismicFMM is free and open source. Development of this tool is a part of authors PhD thesis. National Science Centre Poland provided financial support for this work via NCN grant DEC-2011/02/A/ST10/00284.
NASA Astrophysics Data System (ADS)
Cai, Tao
2016-04-01
In this paper, we have described a 'stratified' semi-implicit spectral method to study compressible convection in Cartesian geometry. The full set of compressible hydrodynamic equations are solved in conservative forms. The numerical scheme is accurate and efficient, based on fast Fourier/sin/cos spectral transforms in the horizontal directions, Chebyshev spectral transform or second-order finite difference scheme in the vertical direction, and second order semi-implicit scheme in time marching of linear terms. We have checked the validity of both the fully pseudo-spectral scheme and the mixed finite-difference pseudo-spectral scheme by studying the onset of compressible convection. The difference of the critical Rayleigh number between our numerical result and the linear stability analysis is within two percent. Besides, we have computed the Mach numbers with different Rayleigh numbers in compressible convection. It shows good agreement with the numerical results of finite difference methods and finite volume method. This model has wide application in studying laminar and turbulent flow. Illustrative examples of application on horizontal convection, gravity waves, and long-lived vortex are given in this paper.
Westerly, David C.; Mo, Xiaohu; Tomé, Wolfgang A.; Mackie, Thomas R.; DeLuca, Paul M.
2013-01-01
Purpose: Pencil beam algorithms are commonly used for proton therapy dose calculations. Szymanowski and Oelfke [“Two-dimensional pencil beam scaling: An improved proton dose algorithm for heterogeneous media,” Phys. Med. Biol. 47, 3313–3330 (2002)10.1088/0031-9155/47/18/304] developed a two-dimensional (2D) scaling algorithm which accurately models the radial pencil beam width as a function of depth in heterogeneous slab geometries using a scaled expression for the radial kernel width in water as a function of depth and kinetic energy. However, an assumption made in the derivation of the technique limits its range of validity to cases where the input expression for the radial kernel width in water is derived from a local scattering power model. The goal of this work is to derive a generalized form of 2D pencil beam scaling that is independent of the scattering power model and appropriate for use with any expression for the radial kernel width in water as a function of depth. Methods: Using Fermi-Eyges transport theory, the authors derive an expression for the radial pencil beam width in heterogeneous slab geometries which is independent of the proton scattering power and related quantities. The authors then perform test calculations in homogeneous and heterogeneous slab phantoms using both the original 2D scaling model and the new model with expressions for the radial kernel width in water computed from both local and nonlocal scattering power models, as well as a nonlocal parameterization of Molière scattering theory. In addition to kernel width calculations, dose calculations are also performed for a narrow Gaussian proton beam. Results: Pencil beam width calculations indicate that both 2D scaling formalisms perform well when the radial kernel width in water is derived from a local scattering power model. Computing the radial kernel width from a nonlocal scattering model results in the local 2D scaling formula under-predicting the pencil beam width by as
Probability distribution of the index in gauge theory on 2d non-commutative geometry
NASA Astrophysics Data System (ADS)
Aoki, Hajime; Nishimura, Jun; Susaki, Yoshiaki
2007-10-01
We investigate the effects of non-commutative geometry on the topological aspects of gauge theory using a non-perturbative formulation based on the twisted reduced model. The configuration space is decomposed into topological sectors labeled by the index ν of the overlap Dirac operator satisfying the Ginsparg-Wilson relation. We study the probability distribution of ν by Monte Carlo simulation of the U(1) gauge theory on 2d non-commutative space with periodic boundary conditions. In general the distribution is asymmetric under ν mapsto -ν, reflecting the parity violation due to non-commutative geometry. In the continuum and infinite-volume limits, however, the distribution turns out to be dominated by the topologically trivial sector. This conclusion is consistent with the instanton calculus in the continuum theory. However, it is in striking contrast to the known results in the commutative case obtained from lattice simulation, where the distribution is Gaussian in a finite volume, but the width diverges in the infinite-volume limit. We also calculate the average action in each topological sector, and provide deeper understanding of the observed phenomenon.
Amundsen, Morten; Linder, Jacob
2016-01-01
An extension of quasiclassical Keldysh-Usadel theory to higher spatial dimensions than one is crucial in order to describe physical phenomena like charge/spin Hall effects and topological excitations like vortices and skyrmions, none of which are captured in one-dimensional models. We here present a numerical finite element method which solves the non-linearized 2D and 3D quasiclassical Usadel equation relevant for the diffusive regime. We show the application of this on three model systems with non-trivial geometries: (i) a bottlenecked Josephson junction with external flux, (ii) a nanodisk ferromagnet deposited on top of a superconductor and (iii) superconducting islands in contact with a ferromagnet. In case (i), we demonstrate that one may control externally not only the geometrical array in which superconducting vortices arrange themselves, but also to cause coalescence and tune the number of vortices. In case (iii), we show that the supercurrent path can be tailored by incorporating magnetic elements in planar Josephson junctions which also lead to a strong modulation of the density of states. The finite element method presented herein paves the way for gaining insight in physical phenomena which have remained largely unexplored due to the complexity of solving the full quasiclassical equations in higher dimensions. PMID:26961921
Amundsen, Morten; Linder, Jacob
2016-01-01
An extension of quasiclassical Keldysh-Usadel theory to higher spatial dimensions than one is crucial in order to describe physical phenomena like charge/spin Hall effects and topological excitations like vortices and skyrmions, none of which are captured in one-dimensional models. We here present a numerical finite element method which solves the non-linearized 2D and 3D quasiclassical Usadel equation relevant for the diffusive regime. We show the application of this on three model systems with non-trivial geometries: (i) a bottlenecked Josephson junction with external flux, (ii) a nanodisk ferromagnet deposited on top of a superconductor and (iii) superconducting islands in contact with a ferromagnet. In case (i), we demonstrate that one may control externally not only the geometrical array in which superconducting vortices arrange themselves, but also to cause coalescence and tune the number of vortices. In case (iii), we show that the supercurrent path can be tailored by incorporating magnetic elements in planar Josephson junctions which also lead to a strong modulation of the density of states. The finite element method presented herein paves the way for gaining insight in physical phenomena which have remained largely unexplored due to the complexity of solving the full quasiclassical equations in higher dimensions. PMID:26961921
E2GPR - Edit your geometry, Execute GprMax2D and Plot the Results!
NASA Astrophysics Data System (ADS)
Pirrone, Daniele; Pajewski, Lara
2015-04-01
In order to predict correctly the Ground Penetrating Radar (GPR) response from a particular scenario, Maxwell's equations have to be solved, subject to the physical and geometrical properties of the considered problem and to its initial conditions. Several techniques have been developed in computational electromagnetics, for the solution of Maxwell's equations. These methods can be classified into two main categories: differential and integral equation solvers, which can be implemented in the time or spectral domain. All of the different methods present compromises between computational efficiency, stability, and the ability to model complex geometries. The Finite-Difference Time-Domain (FDTD) technique has several advantages over alternative approaches: it has inherent simplicity, efficiency and conditional stability; it is suitable to treat impulsive behavior of the electromagnetic field and can provide either ultra-wideband temporal waveforms or the sinusoidal steady-state response at any frequency within the excitation spectrum; it is accurate and highly versatile; and it has become a mature and well-researched technique. Moreover, the FDTD technique is suitable to be executed on parallel-processing CPU-based computers and to exploit the modern computer visualisation capabilities. GprMax [1] is a very well-known and largely validated FDTD software tool, implemented by A. Giannopoulos and available for free public download on www.gprmax.com, together with examples and a detailled user guide. The tool includes two electromagnetic wave simulators, GprMax2D and GprMax3D, for the full-wave simulation of two-dimensional and three-dimensional GPR models. In GprMax, everything can be done with the aid of simple commands that are used to define the model parameters and results to be calculated. These commands need to be entered in a simple ASCII text file. GprMax output files can be stored in ASCII or binary format. The software is provided with MATLAB functions, which
E2GPR - Edit your geometry, Execute GprMax2D and Plot the Results!
NASA Astrophysics Data System (ADS)
Pirrone, Daniele; Pajewski, Lara
2015-04-01
In order to predict correctly the Ground Penetrating Radar (GPR) response from a particular scenario, Maxwell's equations have to be solved, subject to the physical and geometrical properties of the considered problem and to its initial conditions. Several techniques have been developed in computational electromagnetics, for the solution of Maxwell's equations. These methods can be classified into two main categories: differential and integral equation solvers, which can be implemented in the time or spectral domain. All of the different methods present compromises between computational efficiency, stability, and the ability to model complex geometries. The Finite-Difference Time-Domain (FDTD) technique has several advantages over alternative approaches: it has inherent simplicity, efficiency and conditional stability; it is suitable to treat impulsive behavior of the electromagnetic field and can provide either ultra-wideband temporal waveforms or the sinusoidal steady-state response at any frequency within the excitation spectrum; it is accurate and highly versatile; and it has become a mature and well-researched technique. Moreover, the FDTD technique is suitable to be executed on parallel-processing CPU-based computers and to exploit the modern computer visualisation capabilities. GprMax [1] is a very well-known and largely validated FDTD software tool, implemented by A. Giannopoulos and available for free public download on www.gprmax.com, together with examples and a detailled user guide. The tool includes two electromagnetic wave simulators, GprMax2D and GprMax3D, for the full-wave simulation of two-dimensional and three-dimensional GPR models. In GprMax, everything can be done with the aid of simple commands that are used to define the model parameters and results to be calculated. These commands need to be entered in a simple ASCII text file. GprMax output files can be stored in ASCII or binary format. The software is provided with MATLAB functions, which
Geometry and dynamics of a coupled 4 D-2 D quantum field theory
NASA Astrophysics Data System (ADS)
Bolognesi, Stefano; Chatterjee, Chandrasekhar; Evslin, Jarah; Konishi, Kenichi; Ohashi, Keisuke; Seveso, Luigi
2016-01-01
Geometric and dynamical aspects of a coupled 4 D-2 D interacting quantum field theory — the gauged nonAbelian vortex — are investigated. The fluctuations of the internal 2 D nonAbelian vortex zeromodes excite the massless 4 D Yang-Mills modes and in general give rise to divergent energies. This means that the well-known 2 D C{P}^{N-1} zeromodes associated with a nonAbelian vortex become nonnormalizable.
Self-calibration of cone-beam CT geometry using 3D–2D image registration
Ouadah, S; Stayman, J W; Gang, G J; Ehtiati, T; Siewerdsen, J H
2016-01-01
Robotic C-arms are capable of complex orbits that can increase field of view, reduce artifacts, improve image quality, and/or reduce dose; however, it can be challenging to obtain accurate, reproducible geometric calibration required for image reconstruction for such complex orbits. This work presents a method for geometric calibration for an arbitrary source-detector orbit by registering 2D projection data to a previously acquired 3D image. It also yields a method by which calibration of simple circular orbits can be improved. The registration uses a normalized gradient information similarity metric and the covariance matrix adaptation-evolution strategy optimizer for robustness against local minima and changes in image content. The resulting transformation provides a ‘self-calibration’ of system geometry. The algorithm was tested in phantom studies using both a cone-beam CT (CBCT) test-bench and a robotic C-arm (Artis Zeego, Siemens Healthcare) for circular and non-circular orbits. Self-calibration performance was evaluated in terms of the full-width at half-maximum (FWHM) of the point spread function in CBCT reconstructions, the reprojection error (RPE) of steel ball bearings placed on each phantom, and the overall quality and presence of artifacts in CBCT images. In all cases, self-calibration improved the FWHM—e.g. on the CBCT bench, FWHM = 0.86 mm for conventional calibration compared to 0.65 mm for self-calibration (p < 0.001). Similar improvements were measured in RPE—e.g. on the robotic C-arm, RPE = 0.73 mm for conventional calibration compared to 0.55 mm for self-calibration (p < 0.001). Visible improvement was evident in CBCT reconstructions using self-calibration, particularly about high-contrast, high-frequency objects (e.g. temporal bone air cells and a surgical needle). The results indicate that self-calibration can improve even upon systems with presumably accurate geometric calibration and is applicable to situations where conventional
Self-calibration of cone-beam CT geometry using 3D-2D image registration.
Ouadah, S; Stayman, J W; Gang, G J; Ehtiati, T; Siewerdsen, J H
2016-04-01
Robotic C-arms are capable of complex orbits that can increase field of view, reduce artifacts, improve image quality, and/or reduce dose; however, it can be challenging to obtain accurate, reproducible geometric calibration required for image reconstruction for such complex orbits. This work presents a method for geometric calibration for an arbitrary source-detector orbit by registering 2D projection data to a previously acquired 3D image. It also yields a method by which calibration of simple circular orbits can be improved. The registration uses a normalized gradient information similarity metric and the covariance matrix adaptation-evolution strategy optimizer for robustness against local minima and changes in image content. The resulting transformation provides a 'self-calibration' of system geometry. The algorithm was tested in phantom studies using both a cone-beam CT (CBCT) test-bench and a robotic C-arm (Artis Zeego, Siemens Healthcare) for circular and non-circular orbits. Self-calibration performance was evaluated in terms of the full-width at half-maximum (FWHM) of the point spread function in CBCT reconstructions, the reprojection error (RPE) of steel ball bearings placed on each phantom, and the overall quality and presence of artifacts in CBCT images. In all cases, self-calibration improved the FWHM-e.g. on the CBCT bench, FWHM = 0.86 mm for conventional calibration compared to 0.65 mm for self-calibration (p < 0.001). Similar improvements were measured in RPE-e.g. on the robotic C-arm, RPE = 0.73 mm for conventional calibration compared to 0.55 mm for self-calibration (p < 0.001). Visible improvement was evident in CBCT reconstructions using self-calibration, particularly about high-contrast, high-frequency objects (e.g. temporal bone air cells and a surgical needle). The results indicate that self-calibration can improve even upon systems with presumably accurate geometric calibration and is
Self-calibration of cone-beam CT geometry using 3D-2D image registration
NASA Astrophysics Data System (ADS)
Ouadah, S.; Stayman, J. W.; Gang, G. J.; Ehtiati, T.; Siewerdsen, J. H.
2016-04-01
Robotic C-arms are capable of complex orbits that can increase field of view, reduce artifacts, improve image quality, and/or reduce dose; however, it can be challenging to obtain accurate, reproducible geometric calibration required for image reconstruction for such complex orbits. This work presents a method for geometric calibration for an arbitrary source-detector orbit by registering 2D projection data to a previously acquired 3D image. It also yields a method by which calibration of simple circular orbits can be improved. The registration uses a normalized gradient information similarity metric and the covariance matrix adaptation-evolution strategy optimizer for robustness against local minima and changes in image content. The resulting transformation provides a ‘self-calibration’ of system geometry. The algorithm was tested in phantom studies using both a cone-beam CT (CBCT) test-bench and a robotic C-arm (Artis Zeego, Siemens Healthcare) for circular and non-circular orbits. Self-calibration performance was evaluated in terms of the full-width at half-maximum (FWHM) of the point spread function in CBCT reconstructions, the reprojection error (RPE) of steel ball bearings placed on each phantom, and the overall quality and presence of artifacts in CBCT images. In all cases, self-calibration improved the FWHM—e.g. on the CBCT bench, FWHM = 0.86 mm for conventional calibration compared to 0.65 mm for self-calibration (p < 0.001). Similar improvements were measured in RPE—e.g. on the robotic C-arm, RPE = 0.73 mm for conventional calibration compared to 0.55 mm for self-calibration (p < 0.001). Visible improvement was evident in CBCT reconstructions using self-calibration, particularly about high-contrast, high-frequency objects (e.g. temporal bone air cells and a surgical needle). The results indicate that self-calibration can improve even upon systems with presumably accurate geometric calibration and is
Effects of Training Method and Gender on Learning 2D/3D Geometry
ERIC Educational Resources Information Center
Khairulanuar, Samsudin; Nazre, Abd Rashid; Jamilah, H.; Sairabanu, Omar Khan; Norasikin, Fabil
2010-01-01
This article reports the findings of an experimental study involving 36 primary school students (16 girls, 20 boys, Mean age = 9.5 years, age range: 8-10 years) in geometrical understanding of 2D and 3D objects. Students were assigned into two experimental groups and one control group based on a stratified random sampling procedure. The first…
Evaluation of 2D shallow-water model for spillway flow with a complex geometry
Technology Transfer Automated Retrieval System (TEKTRAN)
Although the two-dimensional (2D) shallow water model is formulated based on several assumptions such as hydrostatic pressure distribution and vertical velocity is negligible, as a simple alternative to the complex 3D model, it has been used to compute water flows in which these assumptions may be ...
World-sheet geometry and baby universes in 2D quantum gravity
NASA Astrophysics Data System (ADS)
Jain, Sanjay; Mathur, Samir D.
1992-07-01
We show that the surface roughness for c<1 matter theories coupled to 2D quantum gravity is described by a self-similar structure of baby universes. There exist baby universes whose neck thickness is of the order of the ultraviolet cutoff, the largest of these having a macroscopic area ~A1/(1-γ), where A is the total area and γ the string susceptibility exponent.
Determining Transition State Geometries in Liquids Using 2D-IR
Harris, Charles; Cahoon, James F.; Sawyer, Karma R.; Schlegel, Jacob P.; Harris, Charles B.
2007-12-11
Many properties of chemical reactions are determined by the transition state connecting reactant and product, yet it is difficult to directly obtain any information about these short-lived structures in liquids. We show that two-dimensional infrared (2D-IR) spectroscopy can provide direct information about transition states by tracking the transformation of vibrational modes as a molecule crossed a transition state. We successfully monitored a simple chemical reaction, the fluxional rearrangement of Fe(CO)5, in which the exchange of axial and equatorial CO ligands causes an exchange of vibrational energy between the normal modes of the molecule. This energy transfer provides direct evidence regarding the time scale, transition state, and mechanism of the reaction.
Boundary to Constructive Solid Geometry Mappings: a Focus on 2-D Issues
NASA Technical Reports Server (NTRS)
Peterson, D. P.
1985-01-01
The problem of converting boundary representation (B-rep) to constructive solid geometry representation (CSG-rep) (and vice versa) is discussed in two phases. The first phase entails finding a CSG-rep that defines the region bounded by a polygonal profile curve. The second phase utilizes the results of the first phase to find a CSG-rep for many non-polygonal profile curves. A mathematically concise representation of a region bounded by a polygonal is presented. Namely, any polygonal region bounded by an n sided polygon may be represented by a binary tree which has at most n planar halfspaces as leaves. A structure for this representation and an algorithm for calculating is discussed.
NASA Astrophysics Data System (ADS)
Moustafa, Salli; Févotte, François; Lathuilière, Bruno; Plagne, Laurent
2014-06-01
The past few years have been marked by a noticeable increase in the interest in 3D whole-core heterogeneous deterministic neutron transport solvers for reference calculations. Due to the extremely large problem sizes tackled by such solvers, they need to use adapted numerical methods and need to be efficiently implemented to take advantage of the full computing power of modern systems. As for numerical methods, one possible approach consists in iterating over resolutions of 2D and 1D MOC problems by taking advantage of prismatic geometries. The MICADO solver, developed at EDF R&D, is a parallel implementation of such a method in distributed and shared memory systems. However it is currently unable to use SIMD vectorization to leverage the full computing power of modern CPUs. In this paper, we describe our first effort to support vectorization in MICADO, typically targeting Intel© SSE CPUs. Both the 2D and 1D algorithms are vectorized, allowing for high expected speedups for the whole spatial solver. We present benchmark computations, which show nearly optimal speedups for our vectorized implementation on the TAKEDA case.
NASA Astrophysics Data System (ADS)
Ghazanfarian, J.; Abbassi, A.
2012-03-01
Analytical and numerical solutions of the 2D transient dual-phase-lag (DPL) heat conduction equation are presented in this article. The geometry is that of a simplified metal oxide semiconductor field effect transistor with a heater placed on it. A temperature jump boundary condition is used on all boundaries in order to consider boundary phonon scattering at the micro- and nanoscale. A combination of a Laplace transformation technique and separation of variables is used to solve governing equations analytically, and a three-level finite difference scheme is employed to generate numerical results. The results are illustrated for three Knudsen numbers of 0.1, 1, and 10 at different instants of time. It is seen that the wave characteristic of the DPL model is strengthened by increasing the Knudsen number. It is found that the combination of the DPL model with the proposed mixed-type temperature boundary condition has the potential to accurately predict a 2D temperature distribution not only within the transistor itself but also in the near-boundary region.
geomIO: A tool for geodynamicists to turn 2D cross-sections into 3D geometries
NASA Astrophysics Data System (ADS)
Baumann, Tobias; Bauville, Arthur
2016-04-01
In numerical deformation models, material properties are usually defined on elements (e.g., in body-fitted finite elements), or on a set of Lagrangian markers (Eulerian, ALE or mesh-free methods). In any case, geometrical constraints are needed to assign different material properties to the model domain. Whereas simple geometries such as spheres, layers or cuboids can easily be programmed, it quickly gets complex and time-consuming to create more complicated geometries for numerical model setups, especially in three dimensions. geomIO (geometry I/O, http://geomio.bitbucket.org/) is a MATLAB-based library that has two main functionalities. First, it can be used to create 3D volumes based on series of 2D vector drawings similar to a CAD program; and second, it uses these 3D volumes to assign material properties to the numerical model domain. The drawings can conveniently be created using the open-source vector graphics software Inkscape. Adobe Illustrator is also partially supported. The drawings represent a series of cross-sections in the 3D model domain, for example, cross-sectional interpretations of seismic tomography. geomIO is then used to read the drawings and to create 3D volumes by interpolating between the cross-sections. In the second part, the volumes are used to assign material phases to markers inside the volumes. Multiple volumes can be created at the same time and, depending on the order of assignment, unions or intersections can be built to assign additional material phases. geomIO also offers the possibility to create 3D temperature structures for geodynamic models based on depth dependent parameterisations, for example the half space cooling model. In particular, this can be applied to geometries of subducting slabs of arbitrary shape. Yet, geomIO is held very general, and can be used for a variety of applications. We present examples of setup generation from pictures of micro-scale tectonics and lithospheric scale setups of 3D present-day model
Ramanathan, Muruganathan; Lokitz, Bradley S.; Messman, Jamie M.; Stafford, Christopher M.; Kilbey II, S. Michael
2013-01-01
We report a simple, one step process for developing wrinkling patterns in azlactone-based polymer thin films and brushes in 2D and 3D surfaces. The polymer used in this work wrinkles spontaneously upon deposition and solidification on a substrate without applying any external strain to the substrate, with the mode of deposition defining the direction of the wrinkles. Wrinkle formation is shown to occur on a variety of substrates over large areas. We also find that a very thin brush-like layer of an azlactone-containing block copolymer also exhibits wrinkled topology. Given the spontaneity and versatility of wrinkle formation, we further demonstrate two proofs-of-concept, i) that these periodic wrinkled structures are not limited to planar surfaces, but are also developed in complex geometries including tubes, cones and other 3D structures; and ii) that this one-step wrinkling process can be used to guide the deposition of metal nanoparticles and quantum dots, creating a periodic, nanopatterned film.
Stable boundary conditions for Cartesian grid calculations
NASA Technical Reports Server (NTRS)
Berger, M. J.; Leveque, R. J.
1990-01-01
The inviscid Euler equations in complicated geometries are solved using a Cartesian grid. This requires solid wall boundary conditions in the irregular grid cells near the boundary. Since these cells may be orders of magnitude smaller than the regular grid cells, stability is a primary concern. An approach to this problem is presented and its use is illustrated.
NASA Astrophysics Data System (ADS)
Rapoport, Diego L.
2011-01-01
In this transdisciplinary article which stems from philosophical considerations (that depart from phenomenology—after Merleau-Ponty, Heidegger and Rosen—and Hegelian dialectics), we develop a conception based on topological (the Moebius surface and the Klein bottle) and geometrical considerations (based on torsion and non-orientability of manifolds), and multivalued logics which we develop into a unified world conception that surmounts the Cartesian cut and Aristotelian logic. The role of torsion appears in a self-referential construction of space and time, which will be further related to the commutator of the True and False operators of matrix logic, still with a quantum superposed state related to a Moebius surface, and as the physical field at the basis of Spencer-Brown's primitive distinction in the protologic of the calculus of distinction. In this setting, paradox, self-reference, depth, time and space, higher-order non-dual logic, perception, spin and a time operator, the Klein bottle, hypernumbers due to Musès which include non-trivial square roots of ±1 and in particular non-trivial nilpotents, quantum field operators, the transformation of cognition to spin for two-state quantum systems, are found to be keenly interwoven in a world conception compatible with the philosophical approach taken for basis of this article. The Klein bottle is found not only to be the topological in-formation for self-reference and paradox whose logical counterpart in the calculus of indications are the paradoxical imaginary time waves, but also a classical-quantum transformer (Hadamard's gate in quantum computation) which is indispensable to be able to obtain a complete multivalued logical system, and still to generate the matrix extension of classical connective Boolean logic. We further find that the multivalued logic that stems from considering the paradoxical equation in the calculus of distinctions, and in particular, the imaginary solutions to this equation
A Computer Code for 2-D Transport Calculations in x-y Geometry Using the Interface Current Method.
Energy Science and Technology Software Center (ESTSC)
1990-12-01
Version 00 RICANT performs 2-dimensional neutron transport calculations in x-y geometry using the interface current method. In the interface current method, the angular neutron currents crossing region surfaces are expanded in terms of the Legendre polynomials in the two half-spaces made by the region surfaces.
Shapiro, A.B.
1983-08-01
The computer code FACET calculates the radiation geometric view factor (alternatively called shape factor, angle factor, or configuration factor) between surfaces for axisymmetric, two-dimensional planar and three-dimensional geometries with interposed third surface obstructions. FACET was developed to calculate view factors for input to finite-element heat-transfer analysis codes. The first section of this report is a brief review of previous radiation-view-factor computer codes. The second section presents the defining integral equation for the geometric view factor between two surfaces and the assumptions made in its derivation. Also in this section are the numerical algorithms used to integrate this equation for the various geometries. The third section presents the algorithms used to detect self-shadowing and third-surface shadowing between the two surfaces for which a view factor is being calculated. The fourth section provides a user's input guide followed by several example problems.
NASA Astrophysics Data System (ADS)
Chaudhuri, A.; Sekhar, M.; Descloitres, M.; Godderis, Y.; Ruiz, L.; Braun, J. J.
2013-11-01
Stochastic modelling is a useful way of simulating complex hard-rock aquifers as hydrological properties (permeability, porosity etc.) can be described using random variables with known statistics. However, very few studies have assessed the influence of topological uncertainty (i.e. the variability of thickness of conductive zones in the aquifer), probably because it is not easy to retrieve accurate statistics of the aquifer geometry, especially in hard rock context. In this paper, we assessed the potential of using geophysical surveys to describe the geometry of a hard rock-aquifer in a stochastic modelling framework. The study site was a small experimental watershed in South India, where the aquifer consisted of a clayey to loamy-sandy zone (regolith) underlain by a conductive fissured rock layer (protolith) and the unweathered gneiss (bedrock) at the bottom. The spatial variability of the thickness of the regolith and fissured layers was estimated by electrical resistivity tomography (ERT) profiles, which were performed along a few cross sections in the watershed. For stochastic analysis using Monte Carlo simulation, the generated random layer thickness was made conditional to the available data from the geophysics. In order to simulate steady state flow in the irregular domain with variable geometry, we used an isoparametric finite element method to discretize the flow equation over an unstructured grid with irregular hexahedral elements. The results indicated that the spatial variability of the layer thickness had a significant effect on reducing the simulated effective steady seepage flux and that using the conditional simulations reduced the uncertainty of the simulated seepage flux. As a conclusion, combining information on the aquifer geometry obtained from geophysical surveys with stochastic modelling is a promising methodology to improve the simulation of groundwater flow in complex hard-rock aquifers.
NASA Astrophysics Data System (ADS)
Huerta-Viga, Adriana; Amirjalayer, Saeed; Domingos, Sérgio R.; Meuzelaar, Heleen; Rupenyan, Alisa; Woutersen, Sander
2015-06-01
Salt bridges play an important role in protein folding and in supramolecular chemistry, but they are difficult to detect and characterize in solution. Here, we investigate salt bridges between glutamate (Glu-) and arginine (Arg+) using two-dimensional infrared (2D-IR) spectroscopy. The 2D-IR spectrum of a salt-bridged dimer shows cross peaks between the vibrational modes of Glu- and Arg+, which provide a sensitive structural probe of Glu-⋯Arg+ salt bridges. We use this probe to investigate a β-turn locked by a salt bridge, an α-helical peptide whose structure is stabilized by salt bridges, and a coiled coil that is stabilized by intra- and intermolecular salt bridges. We detect a bidentate salt bridge in the β-turn, a monodentate one in the α-helical peptide, and both salt-bridge geometries in the coiled coil. To our knowledge, this is the first time 2D-IR has been used to probe tertiary side chain interactions in peptides, and our results show that 2D-IR spectroscopy is a powerful method for investigating salt bridges in solution.
Bailey, T S; Adams, M L; Chang, J H
2008-10-01
We present a new spatial discretization of the discrete-ordinates transport equation in two-dimensional cylindrical (RZ) geometry for arbitrary polygonal meshes. This discretization is a discontinuous finite element method that utilizes the piecewise linear basis functions developed by Stone and Adams. We describe an asymptotic analysis that shows this method to be accurate for many problems in the thick diffusion limit on arbitrary polygons, allowing this method to be applied to radiative transfer problems with these types of meshes. We also present numerical results for multiple problems on quadrilateral grids and compare these results to the well-known bi-linear discontinuous finite element method.
An IPOT meshless method using DC PSE approximation for fluid flow equations in 2D and 3D geometries
NASA Astrophysics Data System (ADS)
Bourantas, G. C.; Loukopoulos, V. C.; Skouras, E. D.; Burganos, V. N.; Nikiforidis, G. C.
2016-06-01
Navier-Stokes (N-S) equations, in their primitive variable (u-v-p) formulation, are numerically solved using the Implicit Potential (IPOT) numerical scheme in the context of strong form Meshless Point Collocation (MPC) method. The unknown field functions are computed using the Discretization Correction Particle Strength Exchange (DC PSE) approximation method. The latter makes use of discrete moment conditions to derive the operator kernels, which leads to low condition number for the moment matrix compared to other meshless interpolation methods and increased stability for the numerical solution. The proposed meshless scheme is applied on 2D and 3D spatial domains, using uniform or irregular set of nodes to represent the domain. The numerical results obtained are compared against those obtained using well-established methods.
Fevotte, F.; Lathuiliere, B.
2013-07-01
The large increase in computing power over the past few years now makes it possible to consider developing 3D full-core heterogeneous deterministic neutron transport solvers for reference calculations. Among all approaches presented in the literature, the method first introduced in [1] seems very promising. It consists in iterating over resolutions of 2D and ID MOC problems by taking advantage of prismatic geometries without introducing approximations of a low order operator such as diffusion. However, before developing a solver with all industrial options at EDF, several points needed to be clarified. In this work, we first prove the convergence of this iterative process, under some assumptions. We then present our high-performance, parallel implementation of this algorithm in the MICADO solver. Benchmarking the solver against the Takeda case shows that the 2D-1D coupling algorithm does not seem to affect the spatial convergence order of the MOC solver. As for performance issues, our study shows that even though the data distribution is suited to the 2D solver part, the efficiency of the ID part is sufficient to ensure a good parallel efficiency of the global algorithm. After this study, the main remaining difficulty implementation-wise is about the memory requirement of a vector used for initialization. An efficient acceleration operator will also need to be developed. (authors)
Blodwell, J.F.
1987-10-01
It is argued that the point structure of space and time must be constructed from the primitive extensional character of space and time. A procedure for doing this is laid down and applied to one-dimensional and two-dimensional systems of abstract extensions. Topological and metrical properties of the constructed point systems, which differ nontrivially from the usual R and R/sup 2/ models, are examined. Briefly, constructed points are associated with directions and the Cartesian point is split. In one-dimension each point splits into a point pair compatible with the linear ordering. An application to one-dimensional particle motion is given, with the result that natural topological assumptions force the number of left point, right point transitions to remain locally finite in a continuous motion. In general, Cartesian points are seen to correspond to certain filters on a suitable Boolean algebra. Constructed points correspond to ultrafilters. Thus, point construction gives a natural refinement of the Cartesian systems.
Electronic Absolute Cartesian Autocollimator
NASA Technical Reports Server (NTRS)
Leviton, Douglas B.
2006-01-01
An electronic absolute Cartesian autocollimator performs the same basic optical function as does a conventional all-optical or a conventional electronic autocollimator but differs in the nature of its optical target and the manner in which the position of the image of the target is measured. The term absolute in the name of this apparatus reflects the nature of the position measurement, which, unlike in a conventional electronic autocollimator, is based absolutely on the position of the image rather than on an assumed proportionality between the position and the levels of processed analog electronic signals. The term Cartesian in the name of this apparatus reflects the nature of its optical target. Figure 1 depicts the electronic functional blocks of an electronic absolute Cartesian autocollimator along with its basic optical layout, which is the same as that of a conventional autocollimator. Referring first to the optical layout and functions only, this or any autocollimator is used to measure the compound angular deviation of a flat datum mirror with respect to the optical axis of the autocollimator itself. The optical components include an illuminated target, a beam splitter, an objective or collimating lens, and a viewer or detector (described in more detail below) at a viewing plane. The target and the viewing planes are focal planes of the lens. Target light reflected by the datum mirror is imaged on the viewing plane at unit magnification by the collimating lens. If the normal to the datum mirror is parallel to the optical axis of the autocollimator, then the target image is centered on the viewing plane. Any angular deviation of the normal from the optical axis manifests itself as a lateral displacement of the target image from the center. The magnitude of the displacement is proportional to the focal length and to the magnitude (assumed to be small) of the angular deviation. The direction of the displacement is perpendicular to the axis about which the
NASA Astrophysics Data System (ADS)
Shea, Thomas; Krimer, Daniel; Costa, Fidel; Hammer, Julia
2014-05-01
One of the achievements in recent years in volcanology is the determination of time-scales of magmatic processes via diffusion in minerals and its addition to the petrologists' and volcanologists' toolbox. The method typically requires one-dimensional modeling of randomly cut crystals from two-dimensional thin sections. Here we address the question whether using 1D (traverse) or 2D (surface) datasets exploited from randomly cut 3D crystals introduces a bias or dispersion in the time-scales estimated, and how this error can be improved or eliminated. Computational simulations were performed using a concentration-dependent, finite-difference solution to the diffusion equation in 3D. The starting numerical models involved simple geometries (spheres, parallelepipeds), Mg/Fe zoning patterns (either normal or reverse), and isotropic diffusion coefficients. Subsequent models progressively incorporated more complexity, 3D olivines possessing representative polyhedral morphologies, diffusion anisotropy along the different crystallographic axes, and more intricate core-rim zoning patterns. Sections and profiles used to compare 1, 2 and 3D diffusion models were selected to be (1) parallel to the crystal axes, (2) randomly oriented but passing through the olivine center, or (3) randomly oriented and sectioned. Results show that time-scales estimated on randomly cut traverses (1D) or surfaces (2D) can be widely distributed around the actual durations of 3D diffusion (~0.2 to 10 times the true diffusion time). The magnitude over- or underestimations of duration are a complex combination of the geometry of the crystal, the zoning pattern, the orientation of the cuts with respect to the crystallographic axes, and the degree of diffusion anisotropy. Errors on estimated time-scales retrieved from such models may thus be significant. Drastic reductions in the uncertainty of calculated diffusion times can be obtained by following some simple guidelines during the course of data
A Cartesian grid approach with hierarchical refinement for compressible flows
NASA Technical Reports Server (NTRS)
Quirk, James J.
1994-01-01
Many numerical studies of flows that involve complex geometries are limited by the difficulties in generating suitable grids. We present a Cartesian boundary scheme for two-dimensional, compressible flows that is unfettered by the need to generate a computational grid and so it may be used, routinely, even for the most awkward of geometries. In essence, an arbitrary-shaped body is allowed to blank out some region of a background Cartesian mesh and the resultant cut-cells are singled out for special treatment. This is done within a finite-volume framework and so, in principle, any explicit flux-based integration scheme can take advantage of this method for enforcing solid boundary conditions. For best effect, the present Cartesian boundary scheme has been combined with a sophisticated, local mesh refinement scheme, and a number of examples are shown in order to demonstrate the efficacy of the combined algorithm for simulations of shock interaction phenomena.
Development and Applications of 3D Cartesian CFD Technology
NASA Technical Reports Server (NTRS)
Melton, John E.; Berger, Marsha J.; VanDalsem, William (Technical Monitor)
1994-01-01
The urgent need for dramatic reductions in aircraft design cycle time is focusing scrutiny upon all aspects of computational fluid dynamics (CFD). These reductions will most likely come not from increased reliance upon user-interactive (and therefore time-expensive) methods, but instead from methods that can be fully automated and incorporated into 'black box' solutions. In comparison with tetrahedral methods, three-dimensional Cartesian grid approaches are in relative infancy, but initial experiences with automated Cartesian techniques are quite promising. Our research is targeted at furthering the development of Cartesian methods so that they can become key elements of a completely automatic grid generation/flow solution procedure applicable to the Euler analysis of complex aircraft geometries.
GSRP/David Marshall: Fully Automated Cartesian Grid CFD Application for MDO in High Speed Flows
NASA Technical Reports Server (NTRS)
2003-01-01
With the renewed interest in Cartesian gridding methodologies for the ease and speed of gridding complex geometries in addition to the simplicity of the control volumes used in the computations, it has become important to investigate ways of extending the existing Cartesian grid solver functionalities. This includes developing methods of modeling the viscous effects in order to utilize Cartesian grids solvers for accurate drag predictions and addressing the issues related to the distributed memory parallelization of Cartesian solvers. This research presents advances in two areas of interest in Cartesian grid solvers, viscous effects modeling and MPI parallelization. The development of viscous effects modeling using solely Cartesian grids has been hampered by the widely varying control volume sizes associated with the mesh refinement and the cut cells associated with the solid surface. This problem is being addressed by using physically based modeling techniques to update the state vectors of the cut cells and removing them from the finite volume integration scheme. This work is performed on a new Cartesian grid solver, NASCART-GT, with modifications to its cut cell functionality. The development of MPI parallelization addresses issues associated with utilizing Cartesian solvers on distributed memory parallel environments. This work is performed on an existing Cartesian grid solver, CART3D, with modifications to its parallelization methodology.
ERIC Educational Resources Information Center
Planinsic, G.; Kos, M.; Jerman, R.
2004-01-01
It is quite easy to make a version of the well known Cartesian diver experiment that uses two immiscible liquids. This allows students to test their knowledge of density and pressure in explaining the diver's behaviour. Construction details are presented here together with a mathematical model to explain the observations.
Cartesian control of redundant robots
NASA Technical Reports Server (NTRS)
Colbaugh, R.; Glass, K.
1989-01-01
A Cartesian-space position/force controller is presented for redundant robots. The proposed control structure partitions the control problem into a nonredundant position/force trajectory tracking problem and a redundant mapping problem between Cartesian control input F is a set member of the set R(sup m) and robot actuator torque T is a set member of the set R(sup n) (for redundant robots, m is less than n). The underdetermined nature of the F yields T map is exploited so that the robot redundancy is utilized to improve the dynamic response of the robot. This dynamically optimal F yields T map is implemented locally (in time) so that it is computationally efficient for on-line control; however, it is shown that the map possesses globally optimal characteristics. Additionally, it is demonstrated that the dynamically optimal F yields T map can be modified so that the robot redundancy is used to simultaneously improve the dynamic response and realize any specified kinematic performance objective (e.g., manipulability maximization or obstacle avoidance). Computer simulation results are given for a four degree of freedom planar redundant robot under Cartesian control, and demonstrate that position/force trajectory tracking and effective redundancy utilization can be achieved simultaneously with the proposed controller.
A one-loop test for construction of 4D N = 4 SYM from 2D SYM via fuzzy-sphere geometry
NASA Astrophysics Data System (ADS)
Matsuura, So; Sugino, Fumihiko
2016-04-01
As a perturbative check of the construction of 4D N=4 supersymmetric Yang-Mills theory (SYM) from mass-deformed N=(8,8) SYM on the 2D lattice, the one-loop effective action for scalar kinetic terms is computed in N=4 U(k) SYM on R^2 × (fuzzy S^2), which is obtained by expanding 2D N=(8,8) U(N) SYM with mass deformation around its fuzzy-sphere classical solution. The radius of the fuzzy sphere is proportional to the inverse of the mass. We consider two successive limits: (1) decompactify the fuzzy sphere to a noncommutative (Moyal) plane and (2) turn off the noncommutativity of the Moyal plane. It is straightforward at the classical level to obtain the ordinary N=4 SYM on R^4 in the limits, while it is nontrivial at the quantum level. The one-loop effective action for the SU(k) sector of the gauge group U(k) coincides with that of the ordinary 4D N=4 SYM in the above limits. Although a "noncommutative anomaly" appears in the overall U(1) sector of the U(k) gauge group, this can be expected to be a gauge artifact not affecting gauge-invariant observables.
A Cartesian embedded boundary method for hyperbolic conservation laws
Sjogreen, B; Petersson, N A
2006-12-04
The authors develop an embedded boundary finite difference technique for solving the compressible two- or three-dimensional Euler equations in complex geometries on a Cartesian grid. The method is second order accurate with an explicit time step determined by the grid size away from the boundary. Slope limiters are used on the embedded boundary to avoid non-physical oscillations near shock waves. They show computed examples of supersonic flow past a cylinder and compare with results computed on a body fitted grid. Furthermore, they discuss the implementation of the method for thin geometries, and show computed examples of transonic flow past an airfoil.
Quality-based generation of weather radar Cartesian products
NASA Astrophysics Data System (ADS)
Osrodka, K.; Szturc, J.
2015-05-01
Weather radar data volumes are commonly processed to obtain various 2-D Cartesian products based on the transfer from polar to Cartesian representations through a certain interpolation method. In this research an algorithm of the spatial interpolation of polar reflectivity data employing quality index data is applied to find the Cartesian reflectivity as plan position indicator products. On this basis, quality-based versions of standard algorithms for the generation of the following products have been developed: ETOP (echo top), MAX (maximum of reflectivity), and VIL (vertically integrated liquid water). Moreover, as an example of a higher-level product, a CONVECTION (detection of convection) has been defined as a specific combination of the above-listed standard products. A corresponding quality field is determined for each generated product, taking into account the quality of the pixels from which a given product was determined and how large a fraction of the investigated heights was scanned. Examples of such quality-based products are presented in the paper.
NASA Astrophysics Data System (ADS)
Ouadah, S.; Stayman, J. W.; Gang, G.; Uneri, A.; Ehtiati, T.; Siewerdsen, J. H.
2015-03-01
Purpose: Robotic C-arm systems are capable of general noncircular orbits whose trajectories can be driven by the particular imaging task. However obtaining accurate calibrations for reconstruction in such geometries can be a challenging problem. This work proposes a method to perform a unique geometric calibration of an arbitrary C-arm orbit by registering 2D projections to a previously acquired 3D image to determine the transformation parameters representing the system geometry. Methods: Experiments involved a cone-beam CT (CBCT) bench system, a robotic C-arm, and three phantoms. A robust 3D-2D registration process was used to compute the 9 degree of freedom (DOF) transformation between each projection and an existing 3D image by maximizing normalized gradient information with a digitally reconstructed radiograph (DRR) of the 3D volume. The quality of the resulting "self-calibration" was evaluated in terms of the agreement with an established calibration method using a BB phantom as well as image quality in the resulting CBCT reconstruction. Results: The self-calibration yielded CBCT images without significant difference in spatial resolution from the standard ("true") calibration methods (p-value >0.05 for all three phantoms), and the differences between CBCT images reconstructed using the "self" and "true" calibration methods were on the order of 10-3 mm-1. Maximum error in magnification was 3.2%, and back-projection ray placement was within 0.5 mm. Conclusion: The proposed geometric "self" calibration provides a means for 3D imaging on general noncircular orbits in CBCT systems for which a geometric calibration is either not available or not reproducible. The method forms the basis of advanced "task-based" 3D imaging methods now in development for robotic C-arms.
NASA Astrophysics Data System (ADS)
Ge, Wenjun; Modest, Michael F.; Roy, Somesh P.
2016-03-01
The high-order spherical harmonics (PN) method for 2-D Cartesian domains is extracted from the 3-D formulation. The number of equations and intensity coefficients reduces to (N + 1)2 / 4 in the 2-D Cartesian formulation compared with N(N + 1) / 2 for the general 3-D PN formulation. The Marshak boundary conditions are extended to solve problems with nonblack and mixed diffuse-specular surfaces. Additional boundary conditions for specified radiative wall flux, for symmetry/specular reflection boundaries have also been developed. The mathematical details of the formulations and their implementation in the OpenFOAM finite volume based CFD software platform are presented. The accuracy and computational cost of the 2-D Cartesian PN are compared with that of the 3-D PN solver and a Photon Monte Carlo solver for a square enclosure, as well as a 45° wedge geometry with variable radiative properties. The new boundary conditions have been applied for both test cases, and the boundary condition for mixed diffuse-specular surfaces is further illustrated by numerical examples of a rectangular geometry enclosed by walls with different surface characteristics.
Development and application of a 3D Cartesian grid Euler method
NASA Technical Reports Server (NTRS)
Melton, John E.; Aftosmis, Michael J.; Berger, Marsha J.; Wong, Michael D.
1995-01-01
This report describes recent progress in the development and application of 3D Cartesian grid generation and Euler flow solution techniques. Improvements to flow field grid generation algorithms, geometry representations, and geometry refinement criteria are presented, including details of a procedure for correctly identifying and resolving extremely thin surface features. An initial implementation of automatic flow field refinement is also presented. Results for several 3D multi-component configurations are provided and discussed.
Turing instabilities on Cartesian product networks
Asllani, Malbor; Busiello, Daniel M.; Carletti, Timoteo; Fanelli, Duccio; Planchon, Gwendoline
2015-01-01
The problem of Turing instabilities for a reaction-diffusion system defined on a complex Cartesian product network is considered. To this end we operate in the linear regime and expand the time dependent perturbation on a basis formed by the tensor product of the eigenvectors of the discrete Laplacian operators, associated to each of the individual networks that build the Cartesian product. The dispersion relation which controls the onset of the instability depends on a set of discrete wavelengths, the eigenvalues of the aforementioned Laplacians. Patterns can develop on the Cartesian network, if they are supported on at least one of its constitutive sub-graphs. Multiplex networks are also obtained under specific prescriptions. In this case, the criteria for the instability reduce to compact explicit formulae. Numerical simulations carried out for the Mimura-Murray reaction kinetics confirm the adequacy of the proposed theory. PMID:26245138
Turing instabilities on Cartesian product networks
NASA Astrophysics Data System (ADS)
Asllani, Malbor; Busiello, Daniel M.; Carletti, Timoteo; Fanelli, Duccio; Planchon, Gwendoline
2015-08-01
The problem of Turing instabilities for a reaction-diffusion system defined on a complex Cartesian product network is considered. To this end we operate in the linear regime and expand the time dependent perturbation on a basis formed by the tensor product of the eigenvectors of the discrete Laplacian operators, associated to each of the individual networks that build the Cartesian product. The dispersion relation which controls the onset of the instability depends on a set of discrete wavelengths, the eigenvalues of the aforementioned Laplacians. Patterns can develop on the Cartesian network, if they are supported on at least one of its constitutive sub-graphs. Multiplex networks are also obtained under specific prescriptions. In this case, the criteria for the instability reduce to compact explicit formulae. Numerical simulations carried out for the Mimura-Murray reaction kinetics confirm the adequacy of the proposed theory.
Turing instabilities on Cartesian product networks.
Asllani, Malbor; Busiello, Daniel M; Carletti, Timoteo; Fanelli, Duccio; Planchon, Gwendoline
2015-01-01
The problem of Turing instabilities for a reaction-diffusion system defined on a complex Cartesian product network is considered. To this end we operate in the linear regime and expand the time dependent perturbation on a basis formed by the tensor product of the eigenvectors of the discrete Laplacian operators, associated to each of the individual networks that build the Cartesian product. The dispersion relation which controls the onset of the instability depends on a set of discrete wavelengths, the eigenvalues of the aforementioned Laplacians. Patterns can develop on the Cartesian network, if they are supported on at least one of its constitutive sub-graphs. Multiplex networks are also obtained under specific prescriptions. In this case, the criteria for the instability reduce to compact explicit formulae. Numerical simulations carried out for the Mimura-Murray reaction kinetics confirm the adequacy of the proposed theory. PMID:26245138
High Energy Boundary Conditions for a Cartesian Mesh Euler Solver
NASA Technical Reports Server (NTRS)
Pandya, Shishir; Murman, Scott; Aftosmis, Michael
2003-01-01
Inlets and exhaust nozzles are common place in the world of flight. Yet, many aerodynamic simulation packages do not provide a method of modelling such high energy boundaries in the flow field. For the purposes of aerodynamic simulation, inlets and exhausts are often fared over and it is assumed that the flow differences resulting from this assumption are minimal. While this is an adequate assumption for the prediction of lift, the lack of a plume behind the aircraft creates an evacuated base region thus effecting both drag and pitching moment values. In addition, the flow in the base region is often mis-predicted resulting in incorrect base drag. In order to accurately predict these quantities, a method for specifying inlet and exhaust conditions needs to be available in aerodynamic simulation packages. A method for a first approximation of a plume without accounting for chemical reactions is added to the Cartesian mesh based aerodynamic simulation package CART3D. The method consists of 3 steps. In the first step, a components approach where each triangle is assigned a component number is used. Here, a method for marking the inlet or exhaust plane triangles as separate components is discussed. In step two, the flow solver is modified to accept a reference state for the components marked inlet or exhaust. In the third step, the flow solver uses these separated components and the reference state to compute the correct flow condition at that triangle. The present method is implemented in the CART3D package which consists of a set of tools for generating a Cartesian volume mesh from a set of component triangulations. The Euler equations are solved on the resulting unstructured Cartesian mesh. The present methods is implemented in this package and its usefulness is demonstrated with two validation cases. A generic missile body is also presented to show the usefulness of the method on a real world geometry.
The 3D Euler solutions using automated Cartesian grid generation
NASA Technical Reports Server (NTRS)
Melton, John E.; Enomoto, Francis Y.; Berger, Marsha J.
1993-01-01
Viewgraphs on 3-dimensional Euler solutions using automated Cartesian grid generation are presented. Topics covered include: computational fluid dynamics (CFD) and the design cycle; Cartesian grid strategy; structured body fit; grid generation; prolate spheroid; and ONERA M6 wing.
Static Aeroelastic Analysis with an Inviscid Cartesian Method
NASA Technical Reports Server (NTRS)
Rodriguez, David L.; Aftosmis, Michael J.; Nemec, Marian; Smith, Stephen C.
2014-01-01
An embedded-boundary Cartesian-mesh flow solver is coupled with a three degree-offreedom structural model to perform static, aeroelastic analysis of complex aircraft geometries. The approach solves the complete system of aero-structural equations using a modular, loosely-coupled strategy which allows the lower-fidelity structural model to deform the highfidelity CFD model. The approach uses an open-source, 3-D discrete-geometry engine to deform a triangulated surface geometry according to the shape predicted by the structural model under the computed aerodynamic loads. The deformation scheme is capable of modeling large deflections and is applicable to the design of modern, very-flexible transport wings. The interface is modular so that aerodynamic or structural analysis methods can be easily swapped or enhanced. This extended abstract includes a brief description of the architecture, along with some preliminary validation of underlying assumptions and early results on a generic 3D transport model. The final paper will present more concrete cases and validation of the approach. Preliminary results demonstrate convergence of the complete aero-structural system and investigate the accuracy of the approximations used in the formulation of the structural model.
Static Aeroelastic Analysis with an Inviscid Cartesian Method
NASA Technical Reports Server (NTRS)
Rodriguez, David L.; Aftosmis, Michael J.; Nemec, Marian; Smith, Stephen C.
2014-01-01
An embedded-boundary, Cartesian-mesh flow solver is coupled with a three degree-of-freedom structural model to perform static, aeroelastic analysis of complex aircraft geometries. The approach solves a nonlinear, aerostructural system of equations using a loosely-coupled strategy. An open-source, 3-D discrete-geometry engine is utilized to deform a triangulated surface geometry according to the shape predicted by the structural model under the computed aerodynamic loads. The deformation scheme is capable of modeling large deflections and is applicable to the design of modern, very-flexible transport wings. The coupling interface is modular so that aerodynamic or structural analysis methods can be easily swapped or enhanced. After verifying the structural model with comparisons to Euler beam theory, two applications of the analysis method are presented as validation. The first is a relatively stiff, transport wing model which was a subject of a recent workshop on aeroelasticity. The second is a very flexible model recently tested in a low speed wind tunnel. Both cases show that the aeroelastic analysis method produces results in excellent agreement with experimental data.
On the Use of Parmetric-CAD Systems and Cartesian Methods for Aerodynamic Design
NASA Technical Reports Server (NTRS)
Nemec, Marian; Aftosmis, Michael J.; Pulliam, Thomas H.
2004-01-01
Automated, high-fidelity tools for aerodynamic design face critical issues in attempting to optimize real-life geometry arid in permitting radical design changes. Success in these areas promises not only significantly shorter design- cycle times, but also superior and unconventional designs. To address these issues, we investigate the use of a parmetric-CAD system in conjunction with an embedded-boundary Cartesian method. Our goal is to combine the modeling capabilities of feature-based CAD with the robustness and flexibility of component-based Cartesian volume-mesh generation for complex geometry problems. We present the development of an automated optimization frame-work with a focus on the deployment of such a CAD-based design approach in a heterogeneous parallel computing environment.
Transonic airfoil design using Cartesian coordinates
NASA Technical Reports Server (NTRS)
Carlson, L. A.
1976-01-01
A numerical technique for designing transonic airfoils having a prescribed pressure distribution (the inverse problem) is presented. The method employs the basic features of Jameson's iterative solution for the full potential equation, except that inverse boundary conditions and Cartesian coordinates are used. The method is a direct-inverse approach that controls trailing-edge closure. Examples show the application of the method to design aft-cambered and other airfoils specifically for transonic flight.
Equivalence of the Path Integral for Fermions in Cartesian and Spherical Coordinates
NASA Astrophysics Data System (ADS)
Briggs, Andrew; Camblong, Horacio E.; Ordóñez, Carlos R.
2013-06-01
The path integral calculation for the free energy of a spin-1/2 Dirac-fermion gas is performed in spherical polar coordinates for a flat space-time geometry. Its equivalence with the Cartesian-coordinate representation is explicitly established. This evaluation involves a relevant limiting case of the fermionic path integral in a Schwarzschild background, whose near-horizon limit has been shown to be related to black hole thermodynamics.
Adjoint Formulation for an Embedded-Boundary Cartesian Method
NASA Technical Reports Server (NTRS)
Nemec, Marian; Aftosmis, Michael J.; Murman, Scott M.; Pulliam, Thomas H.
2004-01-01
Many problems in aerodynamic design can be characterized by smooth and convex objective functions. This motivates the use of gradient-based algorithms, particularly for problems with a large number of design variables, to efficiently determine optimal shapes and configurations that maximize aerodynamic performance. Accurate and efficient computation of the gradient, however, remains a challenging task. In optimization problems where the number of design variables dominates the number of objectives and flow- dependent constraints, the cost of gradient computations can be significantly reduced by the use of the adjoint method. The problem of aerodynamic optimization using the adjoint method has been analyzed and validated for both structured and unstructured grids. The method has been applied to design problems governed by the potential, Euler, and Navier-Stokes equations and can be subdivided into the continuous and discrete formulations. Giles and Pierce provide a detailed review of both approaches. Most implementations rely on grid-perturbation or mapping procedures during the gradient computation that explicitly couple changes in the surface shape to the volume grid. The solution of the adjoint equation is usually accomplished using the same scheme that solves the governing flow equations. Examples of such code reuse include multistage Runge-Kutta schemes coupled with multigrid, approximate-factorization, line-implicit Gauss-Seidel, and also preconditioned GMRES. The development of the adjoint method for aerodynamic optimization problems on Cartesian grids has been limited. In contrast to implementations on structured and unstructured grids, Cartesian grid methods decouple the surface discretization from the volume grid. This feature makes Cartesian methods well suited for the automated analysis of complex geometry problems, and consequently a promising approach to aerodynamic optimization. Melvin e t al. developed an adjoint formulation for the TRANAIR code
Schwartz, Peter; Barad, Michael; Colella, Phillip; Ligocki, Terry
2004-11-02
We present an algorithm for solving Poisson's equation and the heat equation on irregular domains in three dimensions. Our work uses the Cartesian grid embedded boundary algorithm for 2D problems of Johansen and Colella (1998, J. Comput. Phys. 147(2):60-85) and extends work of McCorquodale, Colella and Johansen (2001, J. Comput. Phys. 173(2):60-85). Our method is based on a finite-volume discretization of the operator, on the control volumes formed by intersecting the Cartesian grid cells with the domain, combined with a second-order accurate discretization of the fluxes. The resulting method provides uniformly second-order accurate solutions and gradients and is amenable to geometric multigrid solvers.
A Cartesian Adaptive Level Set Method for Two-Phase Flows
NASA Technical Reports Server (NTRS)
Ham, F.; Young, Y.-N.
2003-01-01
In the present contribution we develop a level set method based on local anisotropic Cartesian adaptation as described in Ham et al. (2002). Such an approach should allow for the smallest possible Cartesian grid capable of resolving a given flow. The remainder of the paper is organized as follows. In section 2 the level set formulation for free surface calculations is presented and its strengths and weaknesses relative to the other free surface methods reviewed. In section 3 the collocated numerical method is described. In section 4 the method is validated by solving the 2D and 3D drop oscilation problem. In section 5 we present some results from more complex cases including the 3D drop breakup in an impulsively accelerated free stream, and the 3D immiscible Rayleigh-Taylor instability. Conclusions are given in section 6.
Influence of the Geometry on Mantle Convection Models
NASA Astrophysics Data System (ADS)
Noack, L.; Tosi, N.
2012-04-01
Modelling of geodynamic processes like mantle or core convection has strongly improved over the last two decades thanks to the steady development of numerical codes that tend to incorporate a more and more realistic physics. High-performance parallel computations allow the simulation of complex problems, such as the self-consistent generation of tectonic plates or the formation of planetary magnetic fields. However, the need to perform broad explorations of the parameter space and the large computational demands imposed by the non-linear, multi-scale nature of convection require several simplifications, in the domain geometry as well as in the physical complexity of the problem. A straightforward approach to limit the computational complexity of the simulations is to decrease the total number of degrees of freedom of the problem by reducing either the number of dimensions or the size of the model domain. On the one hand, for a given resolution, a 3D spherical shell clearly needs a much larger number of grid points than a 2D cylindrical shell or a 2D Cartesian box. At the resolutions typically employed to solve mantle convection problems, this difference amounts to at least a factor of a few hundreds. On the other hand, for certain problems, only a relatively small part of the mantle may be of interest, as in the case of the modelling of subduction [1], mid-ocean ridges or transform faults [2]. We adapted the code GAIA [3] to solve the Stokes problem in several different geometries (Cartesian box, cylindrical, spherical and regional-spherical) and dimensions (2D and 3D) and started a benchmark along the lines of [4] to assess the loss of accuracy when using reduced domains instead of a 3D spherical shell [5]. In general, upwellings in Cartesian geometry are rather flat, whereas the spherical geometry changes their shape to more mushroom-like structures. Furthermore, the number of plumes, which is representative of the characteristic wavelength of convection, varies
CAD-Based Aerodynamic Design of Complex Configurations using a Cartesian Method
NASA Technical Reports Server (NTRS)
Nemec, Marian; Aftosmis, Michael J.; Pulliam, Thomas H.
2003-01-01
A modular framework for aerodynamic optimization of complex geometries is developed. By working directly with a parametric CAD system, complex-geometry models are modified nnd tessellated in an automatic fashion. The use of a component-based Cartesian method significantly reduces the demands on the CAD system, and also provides for robust and efficient flowfield analysis. The optimization is controlled using either a genetic or quasi-Newton algorithm. Parallel efficiency of the framework is maintained even when subject to limited CAD resources by dynamically re-allocating the processors of the flow solver. Overall, the resulting framework can explore designs incorporating large shape modifications and changes in topology.
Explicitly computing geodetic coordinates from Cartesian coordinates
NASA Astrophysics Data System (ADS)
Zeng, Huaien
2013-04-01
This paper presents a new form of quartic equation based on Lagrange's extremum law and a Groebner basis under the constraint that the geodetic height is the shortest distance between a given point and the reference ellipsoid. A very explicit and concise formulae of the quartic equation by Ferrari's line is found, which avoids the need of a good starting guess for iterative methods. A new explicit algorithm is then proposed to compute geodetic coordinates from Cartesian coordinates. The convergence region of the algorithm is investigated and the corresponding correct solution is given. Lastly, the algorithm is validated with numerical experiments.
NASA Astrophysics Data System (ADS)
Leboeuf, Jean-Noel; Decyk, Viktor; Newman, David; Sanchez, Raul
2012-03-01
The massively parallel, nonlinear, 3D, toroidal, electrostatic, gyrokinetic, PIC, Cartesian geometry UCAN code, with particle ions and adiabatic electrons, has been successfully exercised to identify non-diffusive transport characteristics in DIII-D-like tokamak discharges. The limitation in applying UCAN to larger scale discharges is the 1D domain decomposition in the toroidal (or z-) direction for massively parallel implementation using MPI which has restricted the calculations to a few hundred ion Larmor radii per minor radius. To exceed these sizes, we have implemented 2D domain decomposition in UCAN with the addition of the y-direction to the processor mix. This has been facilitated by use of relevant components in the 2D domain decomposed PLIB2 library of field and particle management routines developed for UCLA's UPIC framework of conventional PIC codes. The gyro-averaging in gyrokinetic codes has necessitated the use of replicated arrays for efficient charge accumulation and particle push. The 2D domain-decomposed UCAN2 code reproduces the original 1D domain results within roundoff. Production calculations at large system sizes have been performed with UCAN2 on 131072 processors of the Cray XE6 at NERSC.
Software for Automated Generation of Cartesian Meshes
NASA Technical Reports Server (NTRS)
Aftosmis, Michael J.; Melton, John E.; Berger, Marshal J.
2006-01-01
Cart3D is a collection of computer programs for generating Cartesian meshes [for computational fluid dynamics (CFD) and other applications] in volumes bounded by solid objects. Aspects of Cart3D at earlier stages of development were reported in "Robust and Efficient Generation of Cartesian Meshes for CFD" (ARC-14275), NASA Tech Briefs, Vol. 23, No. 8 (August 1999), page 30. The geometric input to Cart3D comprises surface triangulations like those commonly generated by computer-aided-design programs. Complexly shaped objects can be represented as assemblies of simpler ones. Cart3D deletes all portions of such an assembled object that are not on the exterior surface. Intersections between components are preserved in the resulting triangulation. A tie-breaking routine unambiguously resolves geometric degeneracies. Then taking the intersected surface triangulation as input, the volume mesh is generated through division of cells of an initially coarse hexahedral grid. Cells are subdivided to refine the grid in regions of increased surface curvature and/or increased flow gradients. Cells that become split into multiple unconnected regions by thin pieces of surface are identified.
Adjoint Algorithm for CAD-Based Shape Optimization Using a Cartesian Method
NASA Technical Reports Server (NTRS)
Nemec, Marian; Aftosmis, Michael J.
2004-01-01
Adjoint solutions of the governing flow equations are becoming increasingly important for the development of efficient analysis and optimization algorithms. A well-known use of the adjoint method is gradient-based shape optimization. Given an objective function that defines some measure of performance, such as the lift and drag functionals, its gradient is computed at a cost that is essentially independent of the number of design variables (geometric parameters that control the shape). More recently, emerging adjoint applications focus on the analysis problem, where the adjoint solution is used to drive mesh adaptation, as well as to provide estimates of functional error bounds and corrections. The attractive feature of this approach is that the mesh-adaptation procedure targets a specific functional, thereby localizing the mesh refinement and reducing computational cost. Our focus is on the development of adjoint-based optimization techniques for a Cartesian method with embedded boundaries.12 In contrast t o implementations on structured and unstructured grids, Cartesian methods decouple the surface discretization from the volume mesh. This feature makes Cartesian methods well suited for the automated analysis of complex geometry problems, and consequently a promising approach to aerodynamic optimization. Melvin et developed an adjoint formulation for the TRANAIR code, which is based on the full-potential equation with viscous corrections. More recently, Dadone and Grossman presented an adjoint formulation for the Euler equations. In both approaches, a boundary condition is introduced to approximate the effects of the evolving surface shape that results in accurate gradient computation. Central to automated shape optimization algorithms is the issue of geometry modeling and control. The need to optimize complex, "real-life" geometry provides a strong incentive for the use of parametric-CAD systems within the optimization procedure. In previous work, we presented
A Two-dimensional Cartesian and Axisymmetric Study of Combustion-acoustic Interaction
NASA Technical Reports Server (NTRS)
Hood, Caroline; Frendi, Abdelkader
2006-01-01
This paper describes a study of a lean premixed (LP) methane-air combustion wave in a two-dimensional Cartesian and axisymmetric coordinate system. Lean premixed combustors provide low emission and high efficiency; however, they are susceptible to combustion instabilities. The present study focuses on the behavior of the flame as it interacts with an external acoustic disturbance. It was found that the flame oscillations increase as the disturbance amplitude is increased. Furthermore, when the frequency of the disturbance is at resonance with a chamber frequency, the instabilities increase. For the axisymmetric geometry, the flame is found to be more unstable compared to the Cartesian case. In some cases, these instabilities were severe and led to flame extinction. In the axisymmetric case, several passive control devices were tested to assess their effectiveness. It is found that an acoustic cavity is better able at controlling the pressure fluctuations in the chamber.
OTAHAL,THOMAS J.; GALLIS,MICHAIL A.; BARTEL,TIMOTHY J.
2000-06-27
This paper presents an investigation of a technique for using two-dimensional bodies composed of simple polygons with a body decoupled uniform Cmtesian grid in the Direct Simulation Monte Carlo method (DSMC). The method employs an automated grid pre-processing scheme beginning form a CAD geometry definition file, and is based on polygon triangulation using a trapezoid algorithm. A particle-body intersection time comparison is presented between the Icarus DSMC code using a body-fitted structured grid and using a structured body-decoupled Cartesian grid with both linear and logarithmic search techniques. A comparison of neutral flow over a cylinder is presented using the structured body fitted grid and the Cartesian body de-coupled grid.
Elking, Dennis M
2016-08-15
New equations for torque and atomic force are derived for use in flexible molecule force fields with atomic multipoles. The expressions are based on Cartesian tensors with arbitrary multipole rank. The standard method for rotating Cartesian tensor multipoles and calculating torque is to first represent the tensor with n indexes and 3(n) redundant components. In this work, new expressions for directly rotating the unique (n + 1)(n + 2)/2 Cartesian tensor multipole components Θpqr are given by introducing Cartesian tensor rotation matrix elements X(R). A polynomial expression and a recursion relation for X(R) are derived. For comparison, the analogous rotation matrix for spherical tensor multipoles are the Wigner functions D(R). The expressions for X(R) are used to derive simple equations for torque and atomic force. The torque and atomic force equations are applied to the geometry optimization of small molecule crystal unit cells. In addition, a discussion of computational efficiency as a function of increasing multipole rank is given for Cartesian tensors. © 2016 Wiley Periodicals, Inc. PMID:27349179
ARC2D - EFFICIENT SOLUTION METHODS FOR THE NAVIER-STOKES EQUATIONS (DEC RISC ULTRIX VERSION)
NASA Technical Reports Server (NTRS)
Biyabani, S. R.
1994-01-01
ARC2D is a computational fluid dynamics program developed at the NASA Ames Research Center specifically for airfoil computations. The program uses implicit finite-difference techniques to solve two-dimensional Euler equations and thin layer Navier-Stokes equations. It is based on the Beam and Warming implicit approximate factorization algorithm in generalized coordinates. The methods are either time accurate or accelerated non-time accurate steady state schemes. The evolution of the solution through time is physically realistic; good solution accuracy is dependent on mesh spacing and boundary conditions. The mathematical development of ARC2D begins with the strong conservation law form of the two-dimensional Navier-Stokes equations in Cartesian coordinates, which admits shock capturing. The Navier-Stokes equations can be transformed from Cartesian coordinates to generalized curvilinear coordinates in a manner that permits one computational code to serve a wide variety of physical geometries and grid systems. ARC2D includes an algebraic mixing length model to approximate the effect of turbulence. In cases of high Reynolds number viscous flows, thin layer approximation can be applied. ARC2D allows for a variety of solutions to stability boundaries, such as those encountered in flows with shocks. The user has considerable flexibility in assigning geometry and developing grid patterns, as well as in assigning boundary conditions. However, the ARC2D model is most appropriate for attached and mildly separated boundary layers; no attempt is made to model wake regions and widely separated flows. The techniques have been successfully used for a variety of inviscid and viscous flowfield calculations. The Cray version of ARC2D is written in FORTRAN 77 for use on Cray series computers and requires approximately 5Mb memory. The program is fully vectorized. The tape includes variations for the COS and UNICOS operating systems. Also included is a sample routine for CONVEX
ARC2D - EFFICIENT SOLUTION METHODS FOR THE NAVIER-STOKES EQUATIONS (CRAY VERSION)
NASA Technical Reports Server (NTRS)
Pulliam, T. H.
1994-01-01
ARC2D is a computational fluid dynamics program developed at the NASA Ames Research Center specifically for airfoil computations. The program uses implicit finite-difference techniques to solve two-dimensional Euler equations and thin layer Navier-Stokes equations. It is based on the Beam and Warming implicit approximate factorization algorithm in generalized coordinates. The methods are either time accurate or accelerated non-time accurate steady state schemes. The evolution of the solution through time is physically realistic; good solution accuracy is dependent on mesh spacing and boundary conditions. The mathematical development of ARC2D begins with the strong conservation law form of the two-dimensional Navier-Stokes equations in Cartesian coordinates, which admits shock capturing. The Navier-Stokes equations can be transformed from Cartesian coordinates to generalized curvilinear coordinates in a manner that permits one computational code to serve a wide variety of physical geometries and grid systems. ARC2D includes an algebraic mixing length model to approximate the effect of turbulence. In cases of high Reynolds number viscous flows, thin layer approximation can be applied. ARC2D allows for a variety of solutions to stability boundaries, such as those encountered in flows with shocks. The user has considerable flexibility in assigning geometry and developing grid patterns, as well as in assigning boundary conditions. However, the ARC2D model is most appropriate for attached and mildly separated boundary layers; no attempt is made to model wake regions and widely separated flows. The techniques have been successfully used for a variety of inviscid and viscous flowfield calculations. The Cray version of ARC2D is written in FORTRAN 77 for use on Cray series computers and requires approximately 5Mb memory. The program is fully vectorized. The tape includes variations for the COS and UNICOS operating systems. Also included is a sample routine for CONVEX
Mackie, Cameron J; Candian, Alessandra; Huang, Xinchuan; Lee, Timothy J; Tielens, Alexander G G M
2015-06-28
A full derivation of the analytic transformation of the quadratic, cubic, and quartic force constants from normal coordinates to Cartesian coordinates is given. Previous attempts at this transformation have resulted in non-linear transformations; however, for the first time, a simple linear transformation is presented here. Two different approaches have been formulated and implemented, one of which does not require prior knowledge of the translation-rotation eigenvectors from diagonalization of the Hessian matrix. The validity of this method is tested using two molecules H2O and c-C3H2D(+). PMID:26133410
Electrostatic PIC with adaptive Cartesian mesh
NASA Astrophysics Data System (ADS)
Kolobov, Vladimir; Arslanbekov, Robert
2016-05-01
We describe an initial implementation of an electrostatic Particle-in-Cell (ES-PIC) module with adaptive Cartesian mesh in our Unified Flow Solver framework. Challenges of PIC method with cell-based adaptive mesh refinement (AMR) are related to a decrease of the particle-per-cell number in the refined cells with a corresponding increase of the numerical noise. The developed ES-PIC solver is validated for capacitively coupled plasma, its AMR capabilities are demonstrated for simulations of streamer development during high-pressure gas breakdown. It is shown that cell-based AMR provides a convenient particle management algorithm for exponential multiplications of electrons and ions in the ionization events.
An adaptive Cartesian control scheme for manipulators
NASA Technical Reports Server (NTRS)
Seraji, H.
1987-01-01
A adaptive control scheme for direct control of manipulator end-effectors to achieve trajectory tracking in Cartesian space is developed. The control structure is obtained from linear multivariable theory and is composed of simple feedforward and feedback controllers and an auxiliary input. The direct adaptation laws are derived from model reference adaptive control theory and are not based on parameter estimation of the robot model. The utilization of feedforward control and the inclusion of auxiliary input are novel features of the present scheme and result in improved dynamic performance over existing adaptive control schemes. The adaptive controller does not require the complex mathematical model of the robot dynamics or any knowledge of the robot parameters or the payload, and is computationally fast for online implementation with high sampling rates.
Energy Science and Technology Software Center (ESTSC)
2005-07-01
Aniso2d is a two-dimensional seismic forward modeling code. The earth is parameterized by an X-Z plane in which the seismic properties Can have monoclinic with x-z plane symmetry. The program uses a user define time-domain wavelet to produce synthetic seismograms anrwhere within the two-dimensional media.
NASA Astrophysics Data System (ADS)
Martínez, F.; Maksymowicz, A.; Ochoa, H.; Díaz, D.
2015-08-01
This paper discusses an integrated approach that provides new ideas about the structural geometry of the NNE-striking, Cretaceous Chañarcillo Basin located along the eastern Coastal Cordillera in the western Central Andes of northern Chile (27-28° S). The results obtained from the integration of two transverse (E-W) gravity profiles with previous geological information, show that the architecture of this basin is defined by a large NNE-SSE-trending and east-vergent anticline ("Tierra Amarilla Anticlinorium"), which is related to the positive reactivation of a former Cretaceous normal fault (Elisa de Bordos Master Fault). Moreover, intercalations of high and low gravity anomalies and steep gravity gradients reveal a set of buried, west-tilted half-grabens associated with a synthetic normal fault pattern. These results, together with the uplift and folding style of the Cretaceous syn-rift recognized within the basin, suggest that their complete structural geometry could be explained by an inverted fault system linked to the shortening of pre-existing Cretaceous normal fault systems. Ages of the synorogenic deposits exposed unconformably over the frontal limb of the Tierra Amarilla Anticlinorium confirm a Late Cretaceous age for the Andean deformation and tectonic inversion of the basin.
Greg Flach, Frank Smith
2011-12-31
Mesh2d is a Fortran90 program designed to generate two-dimensional structured grids of the form [x(i),y(i,j)] where [x,y] are grid coordinates identified by indices (i,j). The x(i) coordinates alone can be used to specify a one-dimensional grid. Because the x-coordinates vary only with the i index, a two-dimensional grid is composed in part of straight vertical lines. However, the nominally horizontal y(i,j0) coordinates along index i are permitted to undulate or otherwise vary. Mesh2d also assigns an integer material type to each grid cell, mtyp(i,j), in a user-specified manner. The complete grid is specified through three separate input files defining the x(i), y(i,j), and mtyp(i,j) variations.
Energy Science and Technology Software Center (ESTSC)
2011-12-31
Mesh2d is a Fortran90 program designed to generate two-dimensional structured grids of the form [x(i),y(i,j)] where [x,y] are grid coordinates identified by indices (i,j). The x(i) coordinates alone can be used to specify a one-dimensional grid. Because the x-coordinates vary only with the i index, a two-dimensional grid is composed in part of straight vertical lines. However, the nominally horizontal y(i,j0) coordinates along index i are permitted to undulate or otherwise vary. Mesh2d also assignsmore » an integer material type to each grid cell, mtyp(i,j), in a user-specified manner. The complete grid is specified through three separate input files defining the x(i), y(i,j), and mtyp(i,j) variations.« less
NASA Astrophysics Data System (ADS)
Lotsch, Bettina V.
2015-07-01
Graphene's legacy has become an integral part of today's condensed matter science and has equipped a whole generation of scientists with an armory of concepts and techniques that open up new perspectives for the postgraphene area. In particular, the judicious combination of 2D building blocks into vertical heterostructures has recently been identified as a promising route to rationally engineer complex multilayer systems and artificial solids with intriguing properties. The present review highlights recent developments in the rapidly emerging field of 2D nanoarchitectonics from a materials chemistry perspective, with a focus on the types of heterostructures available, their assembly strategies, and their emerging properties. This overview is intended to bridge the gap between two major—yet largely disjunct—developments in 2D heterostructures, which are firmly rooted in solid-state chemistry or physics. Although the underlying types of heterostructures differ with respect to their dimensions, layer alignment, and interfacial quality, there is common ground, and future synergies between the various assembly strategies are to be expected.
Efficient and Robust Cartesian Mesh Generation for Building-Cube Method
NASA Astrophysics Data System (ADS)
Ishida, Takashi; Takahashi, Shun; Nakahashi, Kazuhiro
In this study, an efficient and robust Cartesian mesh generation method for Building-Cube Method (BCM) is proposed. It can handle “dirty” geometry data whose surface has cracks, overlaps, and reverse of triangle. BCM mesh generation is implemented by two procedures; cube generation and cell generation in each cube. The cell generation procedure in this study is managed in each cube individually, and parallelized by OpenMP. Efficiency of the parallelized BCM mesh generation is demonstrated for several three-dimensional test cases using a multi-core PC.
The adaptive, cut-cell Cartesian approach (warts and all)
NASA Technical Reports Server (NTRS)
Powell, Kenneth G.
1995-01-01
Solution-adaptive methods based on cutting bodies out of Cartesian grids are gaining popularity now that the ways of circumventing the accuracy problems associated with small cut cells have been developed. Researchers are applying Cartesian-based schemes to a broad class of problems now, and, although there is still development work to be done, it is becoming clearer which problems are best suited to the approach (and which are not). The purpose of this paper is to give a candid assessment, based on applying Cartesian schemes to a variety of problems, of the strengths and weaknesses of the approach as it is currently implemented.
NASA Astrophysics Data System (ADS)
Martínez, F.; Maksymowicz, A.; Ochoa, H.; Díaz, D.
2015-12-01
This paper discusses an integrated approach that provides new ideas about the structural geometry of the NNE-striking, Cretaceous Chañarcillo Basin located along the eastern Coastal Cordillera in the western Central Andes of northern Chile (27-28° S). The results obtained from the integration of two transverse (E-W) gravity profiles with previous geological information show that the architecture of this basin is defined by a large NNE-SSE-trending and east-vergent anticline ("Tierra Amarilla Anticlinorium"), which is related to the positive reactivation of a former Cretaceous normal fault (Elisa de Bordos Master Fault). Moreover, intercalations of high and low gravity anomalies and steep gravity gradients reveal a set of buried, west-tilted half-grabens associated with a synthetic normal fault pattern. These results, together with the uplift and folding style of the Cretaceous synextensional deposits recognized within the basin, suggest that its structure could be explained by an inverted fault system linked to the shortening of pre-existing Cretaceous normal fault systems. Ages of the synorogenic deposits exposed unconformably over the frontal limb of the Tierra Amarilla Anticlinorium confirm a Late Cretaceous age for the Andean deformation and tectonic inversion of the basin.
Dominant 2D magnetic turbulence in the solar wind
Bieber, John W.; Wanner, Wolfgang; Matthaeus, William H.
1996-07-20
There have been recent suggestions that solar wind magnetic turbulence may be a composite of slab geometry (wavevectors aligned with the mean magnetic field) and 2D geometry (wavevectors perpendicular to the mean field). We report results of two new tests of this hypothesis using Helios measurements of mid-inertial range magnetic spectra in the solar wind. The first test is based upon a characteristic difference between reduced magnetic power spectra in the two different directions perpendicular to the mean field. Such a difference is expected for 2D geometry but not for slab geometry. The second test examines the dependence of power spectrum density upon the magnetic field angle (i.e., the angle between the mean magnetic field and the radial direction), a relationship which is expected to be in opposite directions for the slab and 2D components. Both tests support the presence of a dominant ({approx}85% by energy) 2D component in solar wind magnetic turbulence.
2D constant-loss taper for mode conversion
NASA Astrophysics Data System (ADS)
Horth, Alexandre; Kashyap, Raman; Quitoriano, Nathaniel J.
2015-03-01
Proposed in this manuscript is a novel taper geometry, the constant-loss taper (CLT). This geometry is derived with 1D slabs of silicon embedded in silicon dioxide using coupled-mode theory (CMT). The efficiency of the CLT is compared to both linear and parabolic tapers using CMT and 2D finite-difference time-domain simulations. It is shown that over a short 2D, 4.45 μm long taper the CLT's mode conversion efficiency is ~90% which is 10% and 18% more efficient than a 2D parabolic or linear taper, respectively.
Solwnd: A 3D Compressible MHD Code for Solar Wind Studies. Version 1.0: Cartesian Coordinates
NASA Technical Reports Server (NTRS)
Deane, Anil E.
1996-01-01
Solwnd 1.0 is a three-dimensional compressible MHD code written in Fortran for studying the solar wind. Time-dependent boundary conditions are available. The computational algorithm is based on Flux Corrected Transport and the code is based on the existing code of Zalesak and Spicer. The flow considered is that of shear flow with incoming flow that perturbs this base flow. Several test cases corresponding to pressure balanced magnetic structures with velocity shear flow and various inflows including Alfven waves are presented. Version 1.0 of solwnd considers a rectangular Cartesian geometry. Future versions of solwnd will consider a spherical geometry. Some discussions of this issue is presented.
Frequency-Offset Cartesian Feedback Based on Polyphase Difference Amplifiers
Zanchi, Marta G.; Pauly, John M.; Scott, Greig C.
2010-01-01
A modified Cartesian feedback method called “frequency-offset Cartesian feedback” and based on polyphase difference amplifiers is described that significantly reduces the problems associated with quadrature errors and DC-offsets in classic Cartesian feedback power amplifier control systems. In this method, the reference input and feedback signals are down-converted and compared at a low intermediate frequency (IF) instead of at DC. The polyphase difference amplifiers create a complex control bandwidth centered at this low IF, which is typically offset from DC by 200–1500 kHz. Consequently, the loop gain peak does not overlap DC where voltage offsets, drift, and local oscillator leakage create errors. Moreover, quadrature mismatch errors are significantly attenuated in the control bandwidth. Since the polyphase amplifiers selectively amplify the complex signals characterized by a +90° phase relationship representing positive frequency signals, the control system operates somewhat like single sideband (SSB) modulation. However, the approach still allows the same modulation bandwidth control as classic Cartesian feedback. In this paper, the behavior of the polyphase difference amplifier is described through both the results of simulations, based on a theoretical analysis of their architecture, and experiments. We then describe our first printed circuit board prototype of a frequency-offset Cartesian feedback transmitter and its performance in open and closed loop configuration. This approach should be especially useful in magnetic resonance imaging transmit array systems. PMID:20814450
Frequency-Offset Cartesian Feedback Based on Polyphase Difference Amplifiers.
Zanchi, Marta G; Pauly, John M; Scott, Greig C
2010-05-01
A modified Cartesian feedback method called "frequency-offset Cartesian feedback" and based on polyphase difference amplifiers is described that significantly reduces the problems associated with quadrature errors and DC-offsets in classic Cartesian feedback power amplifier control systems.In this method, the reference input and feedback signals are down-converted and compared at a low intermediate frequency (IF) instead of at DC. The polyphase difference amplifiers create a complex control bandwidth centered at this low IF, which is typically offset from DC by 200-1500 kHz. Consequently, the loop gain peak does not overlap DC where voltage offsets, drift, and local oscillator leakage create errors. Moreover, quadrature mismatch errors are significantly attenuated in the control bandwidth. Since the polyphase amplifiers selectively amplify the complex signals characterized by a +90° phase relationship representing positive frequency signals, the control system operates somewhat like single sideband (SSB) modulation. However, the approach still allows the same modulation bandwidth control as classic Cartesian feedback.In this paper, the behavior of the polyphase difference amplifier is described through both the results of simulations, based on a theoretical analysis of their architecture, and experiments. We then describe our first printed circuit board prototype of a frequency-offset Cartesian feedback transmitter and its performance in open and closed loop configuration. This approach should be especially useful in magnetic resonance imaging transmit array systems. PMID:20814450
A Geometric Boolean Library for 2D Objects
Energy Science and Technology Software Center (ESTSC)
2006-01-05
The 2D Boolean Library is a collection of C++ classes -- which primarily represent 2D geometric data and relationships, and routines -- which contain algorithms for 2D geometric Boolean operations and utility functions. Classes are provided for 2D points, lines, arcs, edgeuses, loops, surfaces and mask sets. Routines are provided that incorporate the Boolean operations Union(OR), XOR, Intersection and Difference. Various analytical geometry routines and routines for importing and exporting the data in various filemore » formats, are also provided in the library.« less
A Geometric Boolean Library for 2D Objects
McBride, Corey L.; Yarberry, Victor; Jorgensen, Craig
2006-01-05
The 2D Boolean Library is a collection of C++ classes -- which primarily represent 2D geometric data and relationships, and routines -- which contain algorithms for 2D geometric Boolean operations and utility functions. Classes are provided for 2D points, lines, arcs, edgeuses, loops, surfaces and mask sets. Routines are provided that incorporate the Boolean operations Union(OR), XOR, Intersection and Difference. Various analytical geometry routines and routines for importing and exporting the data in various file formats, are also provided in the library.
Efficient Fluid Dynamic Design Optimization Using Cartesian Grids
NASA Technical Reports Server (NTRS)
Dadone, A.; Grossman, B.; Sellers, Bill (Technical Monitor)
2004-01-01
This report is subdivided in three parts. The first one reviews a new approach to the computation of inviscid flows using Cartesian grid methods. The crux of the method is the curvature-corrected symmetry technique (CCST) developed by the present authors for body-fitted grids. The method introduces ghost cells near the boundaries whose values are developed from an assumed flow-field model in vicinity of the wall consisting of a vortex flow, which satisfies the normal momentum equation and the non-penetration condition. The CCST boundary condition was shown to be substantially more accurate than traditional boundary condition approaches. This improved boundary condition is adapted to a Cartesian mesh formulation, which we call the Ghost Body-Cell Method (GBCM). In this approach, all cell centers exterior to the body are computed with fluxes at the four surrounding cell edges. There is no need for special treatment corresponding to cut cells which complicate other Cartesian mesh methods.
A Cartesian parametrization for the numerical analysis of material instability
Mota, Alejandro; Chen, Qiushi; Foulk, III, James W.; Ostien, Jakob T.; Lai, Zhengshou
2016-02-25
We examine four parametrizations of the unit sphere in the context of material stability analysis by means of the singularity of the acoustic tensor. We then propose a Cartesian parametrization for vectors that lie a cube of side length two and use these vectors in lieu of unit normals to test for the loss of the ellipticity condition. This parametrization is then used to construct a tensor akin to the acoustic tensor. It is shown that both of these tensors become singular at the same time and in the same planes in the presence of a material instability. Furthermore, themore » performance of the Cartesian parametrization is compared against the other parametrizations, with the results of these comparisons showing that in general, the Cartesian parametrization is more robust and more numerically efficient than the others.« less
On differential transformations between Cartesian and curvilinear (geodetic) coordinates
NASA Technical Reports Server (NTRS)
Soler, T.
1976-01-01
Differential transformations are developed between Cartesian and curvilinear orthogonal coordinates. Only matrix algebra is used for the presentation of the basic concepts. After defining the reference systems used the rotation (R), metric (H), and Jacobian (J) matrices of the transformations between cartesian and curvilinear coordinate systems are introduced. A value of R as a function of H and J is presented. Likewise an analytical expression for J(-1) as a function of H(-2) and R is obtained. Emphasis is placed on showing that differential equations are equivalent to conventional similarity transformations. Scaling methods are discussed along with ellipsoidal coordinates. Differential transformations between elipsoidal and geodetic coordinates are established.
NASA Astrophysics Data System (ADS)
Wang, Jin; Ma, Jianyong; Zhou, Changhe
2014-11-01
A 3×3 high divergent 2D-grating with period of 3.842μm at wavelength of 850nm under normal incidence is designed and fabricated in this paper. This high divergent 2D-grating is designed by the vector theory. The Rigorous Coupled Wave Analysis (RCWA) in association with the simulated annealing (SA) is adopted to calculate and optimize this 2D-grating.The properties of this grating are also investigated by the RCWA. The diffraction angles are more than 10 degrees in the whole wavelength band, which are bigger than the traditional 2D-grating. In addition, the small period of grating increases the difficulties of fabrication. So we fabricate the 2D-gratings by direct laser writing (DLW) instead of traditional manufacturing method. Then the method of ICP etching is used to obtain the high divergent 2D-grating.
Polar versus Cartesian velocity models for maneuvering target tracking with IMM
NASA Astrophysics Data System (ADS)
Laneuville, Dann
This paper compares various model sets in different IMM filters for the maneuvering target tracking problem. The aim is to see whether we can improve the tracking performance of what is certainly the most widely used model set in the literature for the maneuvering target tracking problem: a Nearly Constant Velocity model and a Nearly Coordinated Turn model. Our new challenger set consists of a mixed Cartesian position and polar velocity state vector to describe the uniform motion segments and is augmented with the turn rate to obtain the second model for the maneuvering segments. This paper also gives a general procedure to discretize up to second order any non-linear continuous time model with linear diffusion. Comparative simulations on an air defence scenario with a 2D radar, show that this new approach improves significantly the tracking performance in this case.
NASA Technical Reports Server (NTRS)
Finley, Dennis B.
1995-01-01
This report documents results from the Euler Technology Assessment program. The objective was to evaluate the efficacy of Euler computational fluid dynamics (CFD) codes for use in preliminary aircraft design. Both the accuracy of the predictions and the rapidity of calculations were to be assessed. This portion of the study was conducted by Lockheed Fort Worth Company, using a recently developed in-house Cartesian-grid code called SPLITFLOW. The Cartesian grid technique offers several advantages for this study, including ease of volume grid generation and reduced number of cells compared to other grid schemes. SPLITFLOW also includes grid adaptation of the volume grid during the solution convergence to resolve high-gradient flow regions. This proved beneficial in resolving the large vortical structures in the flow for several configurations examined in the present study. The SPLITFLOW code predictions of the configuration forces and moments are shown to be adequate for preliminary design analysis, including predictions of sideslip effects and the effects of geometry variations at low and high angles of attack. The time required to generate the results from initial surface definition is on the order of several hours, including grid generation, which is compatible with the needs of the design environment.
Applications of Space-Filling-Curves to Cartesian Methods for CFD
NASA Technical Reports Server (NTRS)
Aftosmis, Michael J.; Berger, Marsha J.; Murman, Scott M.
2003-01-01
The proposed paper presents a variety novel uses of Space-Filling-Curves (SFCs) for Cartesian mesh methods in 0. While these techniques will be demonstrated using non-body-fitted Cartesian meshes, most are applicable on general body-fitted meshes -both structured and unstructured. We demonstrate the use of single O(N log N) SFC-based reordering to produce single-pass (O(N)) algorithms for mesh partitioning, multigrid coarsening, and inter-mesh interpolation. The intermesh interpolation operator has many practical applications including warm starts on modified geometry, or as an inter-grid transfer operator on remeshed regions in moving-body simulations. Exploiting the compact construction of these operators, we further show that these algorithms are highly amenable to parallelization. Examples using the SFC-based mesh partitioner show nearly linear speedup to 512 CPUs even when using multigrid as a smoother. Partition statistics are presented showing that the SFC partitions are, on-average, within 10% of ideal even with only around 50,000 cells in each subdomain. The inter-mesh interpolation operator also has linear asymptotic complexity and can be used to map a solution with N unknowns to another mesh with M unknowns with O(max(M,N)) operations. This capability is demonstrated both on moving-body simulations and in mapping solutions to perturbed meshes for finite-difference-based gradient design methods.
Applications of Space-Filling-Curves to Cartesian Methods for CFD
NASA Technical Reports Server (NTRS)
Aftosmis, M. J.; Murman, S. M.; Berger, M. J.
2003-01-01
This paper presents a variety of novel uses of space-filling-curves (SFCs) for Cartesian mesh methods in CFD. While these techniques will be demonstrated using non-body-fitted Cartesian meshes, many are applicable on general body-fitted meshes-both structured and unstructured. We demonstrate the use of single theta(N log N) SFC-based reordering to produce single-pass (theta(N)) algorithms for mesh partitioning, multigrid coarsening, and inter-mesh interpolation. The intermesh interpolation operator has many practical applications including warm starts on modified geometry, or as an inter-grid transfer operator on remeshed regions in moving-body simulations Exploiting the compact construction of these operators, we further show that these algorithms are highly amenable to parallelization. Examples using the SFC-based mesh partitioner show nearly linear speedup to 640 CPUs even when using multigrid as a smoother. Partition statistics are presented showing that the SFC partitions are, on-average, within 15% of ideal even with only around 50,000 cells in each sub-domain. The inter-mesh interpolation operator also has linear asymptotic complexity and can be used to map a solution with N unknowns to another mesh with M unknowns with theta(M + N) operations. This capability is demonstrated both on moving-body simulations and in mapping solutions to perturbed meshes for control surface deflection or finite-difference-based gradient design methods.
2D stepping drive for hyperspectral systems
NASA Astrophysics Data System (ADS)
Endrödy, Csaba; Mehner, Hannes; Grewe, Adrian; Sinzinger, Stefan; Hoffmann, Martin
2015-07-01
We present the design, fabrication and characterization of a compact 2D stepping microdrive for pinhole array positioning. The miniaturized solution enables a highly integrated compact hyperspectral imaging system. Based on the geometry of the pinhole array, an inch-worm drive with electrostatic actuators was designed resulting in a compact (1 cm2) positioning system featuring a step size of about 15 µm in a 170 µm displacement range. The high payload (20 mg) as required for the pinhole array and the compact system design exceed the known electrostatic inch-worm-based microdrives.
The Cartesian Diver as an Aid for Teaching Respiratory Physiology
ERIC Educational Resources Information Center
Fitch, Greg K.
2004-01-01
The mechanism by which air enters the mammalian lung is difficult for many students of physiology. In particular, some students have trouble seeing how pressure can be transmitted through a fluid such as the intrapleural fluid and how the magnitude of that pressure can change. A Cartesian diver, an old-time child's toy, may be used as a visual aid…
The Cartesian Diver, Surface Tension and the Cheerios Effect
ERIC Educational Resources Information Center
Chen, Chi-Tung; Lee, Wen-Tang; Kao, Sung-Kai
2014-01-01
A Cartesian diver can be used to measure the surface tension of a liquid to a certain extent. The surface tension measurement is related to the two critical pressures at which the diver is about to sink and about to emerge. After sinking because of increasing pressure, the diver is repulsed to the centre of the vessel. After the pressure is…
A Lot of Good Physics in the Cartesian Diver
ERIC Educational Resources Information Center
De Luca, Roberto; Ganci, Salvatore
2011-01-01
The Cartesian diver experiment certainly occupies a place of honour in old physics textbooks as a vivid demonstration of Archimedes' buoyancy. The original experiment, as described in old textbooks, shows Archimedes buoyancy qualitatively: when the increased weight of the diver is not counterbalanced by Archimedes' buoyancy, the diver sinks. When…
Dominant 2D magnetic turbulence in the solar wind
NASA Technical Reports Server (NTRS)
Bieber, John W.; Wanner, Wolfgang; Matthaeus, William H.
1995-01-01
There have been recent suggestions that solar wind magnetic turbulence may be a composite of slab geometry (wavevector aligned with the mean magnetic field) and 2D geometry (wavevectors perpendicular to the mean field). We report results of two new tests of this hypothesis using Helios measurements of inertial ranged magnetic spectra in the solar wind. The first test is based upon a characteristic difference between perpendicular and parallel reduced power spectra which is expected for the 2D component but not for the slab component. The second test examines the dependence of power spectrum density upon the magnetic field angle (i.e., the angle between the mean magnetic field and the radial direction), a relationship which is expected to be in opposite directions for the slab and 2D components. Both tests support the presence of a dominant (approximately 85 percent by energy) 2D component in solar wind magnetic turbulence.
Validation of Inlet and Exhaust Boundary Conditions for a Cartesian Method
NASA Technical Reports Server (NTRS)
Pandya, Shishir A.; Murman, Scott M.; Aftosmis, Michael J.
2004-01-01
Inlets and exhaust nozzles are often omitted in aerodynamic simulations of aircraft due to the complexities involved in the modeling of engine details and flow physics. However, the omission is often improper since inlet or plume flows may have a substantial effect on vehicle aerodynamics. A method for modeling the effect of inlets and exhaust plumes using boundary conditions within an inviscid Cartesian flow solver is presented. This approach couples with both CAD systems and legacy geometry to provide an automated tool suitable for parameter studies. The method is validated using two and three-dimensional test problems which are compared with both theoretical and experimental results. The numerical results demonstrate excellent agreement with theory and available data, even for extremely strong jets and very sensitive inlets.
2D Regimes of Non-Fourier Convection
NASA Astrophysics Data System (ADS)
Papanicolaou, N. C.
2010-11-01
In this work, we investigate the 2D flow in a rectangular cavity subject to both vertical and horizontal temperature gradients. The linearized model is studied and the effect of thermal relaxation, as described by the Maxwell-Cattaneo law of heat conduction is examined. To this end, a spectral numerical model is created based on a Galerkin expansion. The basis is the Cartesian product of systems of beam functions and trigonometric functions. The natural modes of the system are derived for both the Fourier and non-Fourier models. The results are compared to earlier works for the plain Fourier law. Our computations show that for the same set of parameters, the Maxwell-Cattaneo law yields modes which are quantitatively different from the Fourier. It is found that the real parts of the eigenvalues increase with the Straughan number Sg, which quantifies the non-Fourier effects. This confirms the destabilizing effect of the MC-law on the convective flow.
Generation and Radiation of Acoustic Waves from a 2-D Shear Layer
NASA Technical Reports Server (NTRS)
Agarwal, Anurag; Morris, Philip J.
2000-01-01
A parallel numerical simulation of the radiation of sound from an acoustic source inside a 2-D jet is presented in this paper. This basic benchmark problem is used as a test case for scattering problems that are presently being solved by using the Impedance Mismatch Method (IMM). In this technique, a solid body in the domain is represented by setting the acoustic impedance of each medium, encountered by a wave, to a different value. This impedance discrepancy results in reflected and scattered waves with appropriate amplitudes. The great advantage of the use of this method is that no modifications to a simple Cartesian grid need to be made for complicated geometry bodies. Thus, high order finite difference schemes may be applied simply to all parts of the domain. In the IMM, the total perturbation field is split into incident and scattered fields. The incident pressure is assumed to be known and the equivalent sources for the scattered field are associated with the presence of the scattering body (through the impedance mismatch) and the propagation of the incident field through a non-uniform flow. An earlier version of the technique could only handle uniform flow in the vicinity of the source and at the outflow boundary. Scattering problems in non-uniform mean flow are of great practical importance (for example, scattering from a high lift device in a non-uniform mean flow or the effects of a fuselage boundary layer). The solution to this benchmark problem, which has an acoustic wave propagating through a non-uniform mean flow, serves as a test case for the extensions of the IMM technique.
NASA Astrophysics Data System (ADS)
Yin, Bo; Luo, Haoxiang
2010-11-01
Two-phase flows in various industrial applications often occur in complex geometries. To simulate this type of flows, we have combined two different immersed-boundary methods to handle the fluid-fluid and the fluid-solid interfaces separately. For the fluid-fluid interface, a diffuse-interface method is employed where the discontinuities of the material properties and the traction jump are all regularized using an approximate Dirac's Delta function. For the fluid-solid interface, ghost nodes and a local flow reconstruction are employed to complement the finite-difference discretization and to incorporate the boundary conditions. A single-block Cartesian mesh is used to discretize the entire domain. Both 2D and 3D codes have been implemented. Validation and code applications will be demonstrated at the conference. *Supported by the ACS Petroleum Research Fund.
Energy Science and Technology Software Center (ESTSC)
2004-08-01
AnisWave2D is a 2D finite-difference code for a simulating seismic wave propagation in fully anisotropic materials. The code is implemented to run in parallel over multiple processors and is fully portable. A mesh refinement algorithm has been utilized to allow the grid-spacing to be tailored to the velocity model, avoiding the over-sampling of high-velocity materials that usually occurs in fixed-grid schemes.
On the Use of CAD and Cartesian Methods for Aerodynamic Optimization
NASA Technical Reports Server (NTRS)
Nemec, M.; Aftosmis, M. J.; Pulliam, T. H.
2004-01-01
The objective for this paper is to present the development of an optimization capability for Curt3D, a Cartesian inviscid-flow analysis package. We present the construction of a new optimization framework and we focus on the following issues: 1) Component-based geometry parameterization approach using parametric-CAD models and CAPRI. A novel geometry server is introduced that addresses the issue of parallel efficiency while only sparingly consuming CAD resources; 2) The use of genetic and gradient-based algorithms for three-dimensional aerodynamic design problems. The influence of noise on the optimization methods is studied. Our goal is to create a responsive and automated framework that efficiently identifies design modifications that result in substantial performance improvements. In addition, we examine the architectural issues associated with the deployment of a CAD-based approach in a heterogeneous parallel computing environment that contains both CAD workstations and dedicated compute engines. We demonstrate the effectiveness of the framework for a design problem that features topology changes and complex geometry.
Simulations of 6-DOF Motion with a Cartesian Method
NASA Technical Reports Server (NTRS)
Murman, Scott M.; Aftosmis, Michael J.; Berger, Marsha J.; Kwak, Dochan (Technical Monitor)
2003-01-01
Coupled 6-DOF/CFD trajectory predictions using an automated Cartesian method are demonstrated by simulating a GBU-32/JDAM store separating from an F-18C aircraft. Numerical simulations are performed at two Mach numbers near the sonic speed, and compared with flight-test telemetry and photographic-derived data. Simulation results obtained with a sequential-static series of flow solutions are contrasted with results using a time-dependent flow solver. Both numerical methods show good agreement with the flight-test data through the first half of the simulations. The sequential-static and time-dependent methods diverge over the last half of the trajectory prediction. after the store produces peak angular rates. A cost comparison for the Cartesian method is included, in terms of absolute cost and relative to computing uncoupled 6-DOF trajectories. A detailed description of the 6-DOF method, as well as a verification of its accuracy, is provided in an appendix.
Cartesian to geodetic coordinates conversion on a triaxial ellipsoid
NASA Astrophysics Data System (ADS)
Ligas, Marcin
2012-04-01
A new method of transforming Cartesian to geodetic (or planetographic) coordinates on a triaxial ellipsoid is presented. The method is based on simple reasoning coming from essentials of vector calculus. The reasoning results in solving a nonlinear system of equations for coordinates of the point being the projection of a point located outside or inside a triaxial ellipsoid along the normal to the ellipsoid. The presented method has been compared to a vector method of Feltens (J Geod 83:129-137, 2009) who claims that no other methods are available in the literature. Generally, our method turns out to be more accurate, faster and applicable to celestial bodies characterized by different geometric parameters. The presented method also fits to the classical problem of converting Cartesian to geodetic coordinates on the ellipsoid of revolution.
Irreducible Cartesian tensors of highest weight, for arbitrary order
NASA Astrophysics Data System (ADS)
Mane, S. R.
2016-03-01
A closed form expression is presented for the irreducible Cartesian tensor of highest weight, for arbitrary order. Two proofs are offered, one employing bookkeeping of indices and, after establishing the connection with the so-called natural tensors and their projection operators, the other one employing purely coordinate-free tensor manipulations. Some theorems and formulas in the published literature are generalized from SO(3) to SO(n), for dimensions n ≥ 3.
Errors of Remapping of Radar Estimates onto Cartesian Coordinates
NASA Astrophysics Data System (ADS)
Sharif, H. O.; Ogden, F. L.
2014-12-01
Recent upgrades to operational radar rainfall products in terms of quality and resolution call for re-examination of the factors that contribute to the uncertainty of radar rainfall estimation. Remapping or gridding of radar polar observations onto Cartesian coordinates is implemented using various methods, and is often applied when radar estimates are compared against rain gauge observations, in hydrologic applications, or for merging data from different radars. However, assuming perfect radar observations, many of the widely used remapping methodologies do not conserve mass for the rainfall rate field. Research has suggested that optimal remapping should select all polar bins falling within or intersecting a Cartesian grid and assign them weights based on the proportion of each individual bin's area falling within the grid. However, to reduce computational demand practitioners use a variety of approximate remapping approaches. The most popular approximate approaches used are those based on extracting information from radar bins whose centers fall within a certain distance from the center of the Cartesian grid. This paper introduces a mass-conserving method for remapping, which we call "precise remapping", and evaluates it by comparing against two other commonly used remapping methods based on areal weighting and distance. Results show that the choice of the remapping method can lead to large errors in grid-averaged rainfall accumulations.
Frequency-offset Cartesian feedback for MRI power amplifier linearization.
Zanchi, Marta G; Stang, Pascal; Kerr, Adam; Pauly, John M; Scott, Greig C
2011-02-01
High-quality magnetic resonance imaging (MRI) requires precise control of the transmit radio-frequency (RF) field. In parallel excitation applications such as transmit SENSE, high RF power linearity is essential to cancel aliased excitations. In widely-employed class AB power amplifiers, gain compression, cross-over distortion, memory effects, and thermal drift all distort the RF field modulation and can degrade image quality. Cartesian feedback (CF) linearization can mitigate these effects in MRI, if the quadrature mismatch and dc offset imperfections inherent in the architecture can be minimized. In this paper, we present a modified Cartesian feedback technique called "frequency-offset Cartesian feedback" (FOCF) that significantly reduces these problems. In the FOCF architecture, the feedback control is performed at a low intermediate frequency rather than dc, so that quadrature ghosts and dc errors are shifted outside the control bandwidth. FOCF linearization is demonstrated with a variety of typical MRI pulses. Simulation of the magnetization obtained with the Bloch equation demonstrates that high-fidelity RF reproduction can be obtained even with inexpensive class AB amplifiers. Finally, the enhanced RF fidelity of FOCF over CF is demonstrated with actual images obtained in a 1.5 T MRI system. PMID:20959264
Topology preserving advection of implicit interfaces on Cartesian grids
NASA Astrophysics Data System (ADS)
Qin, Zhipeng; Delaney, Keegan; Riaz, Amir; Balaras, Elias
2015-06-01
Accurate representation of implicit interface topology is important for the numerical computation of two phase flow on Cartesian grids. A new method is proposed for the construction of signed distance function by geometrically projecting interface topology onto the Cartesian grid using a multi-level projection framework. The method involves a stepwise improvement in the approximation to the signed distance function based on pointwise, piecewise and locally smooth reconstructions of the interface. We show that this approach provides accurate representation of the projected interface and its topology on the Cartesian grid, including the distance from the interface and the interface normal and curvature. The projected interface can be in the form of either a connected set of marker particles that evolve with Lagrangian advection, or a discrete set of points associated with an implicit interface that evolves with the advection of a scalar function. The signed distance function obtained with geometric projection is independent of the details of the scaler field, in contrast to the conventional approach where advection and reinitialization cannot be decoupled. As a result, errors introduced by reinitialization do not amplify advection errors, which leads to substantial improvement in both volume conservation and topology representation.
Frequency-Offset Cartesian Feedback for MRI Power Amplifier Linearization
Zanchi, Marta Gaia; Stang, Pascal; Kerr, Adam; Pauly, John Mark; Scott, Greig Cameron
2011-01-01
High-quality magnetic resonance imaging (MRI) requires precise control of the transmit radio-frequency field. In parallel excitation applications such as transmit SENSE, high RF power linearity is essential to cancel aliased excitations. In widely-employed class AB power amplifiers, gain compression, cross-over distortion, memory effects, and thermal drift all distort the RF field modulation and can degrade image quality. Cartesian feedback (CF) linearization can mitigate these effects in MRI, if the quadrature mismatch and DC offset imperfections inherent in the architecture can be minimized. In this paper, we present a modified Cartesian feedback technique called “frequency-offset Cartesian feedback” (FOCF) that significantly reduces these problems. In the FOCF architecture, the feedback control is performed at a low intermediate frequency rather than DC, so that quadrature ghosts and DC errors are shifted outside the control bandwidth. FOCF linearization is demonstrated with a variety of typical MRI pulses. Simulation of the magnetization obtained with the Bloch equation demonstrates that high-fidelity RF reproduction can be obtained even with inexpensive class AB amplifiers. Finally, the enhanced RF fidelity of FOCF over CF is demonstrated with actual images obtained in a 1.5 T MRI system. PMID:20959264
NASA Astrophysics Data System (ADS)
Bakhsh, Abeer; Samtaney, Ravi
2015-11-01
Numerical simulations and analysis in Cartesian slab geometry for nonlinear ideal magnetohydrodynamics (MHD) indicate that the Richtmyer-Meshkov instability (RMI) is suppressed in the presence of a magnetic field. An analytical solution of incompressible 2-D MHD RMI of an impulsively accelerated interface was investigated by Wheatley et al. (Phys. Rev. Lett. 2005; J. Fluid Mech. 2005) who found that, for a finite magnetic field, although the initial growth rate of the interface is unaffected by the presence of magnetic field, the late-time amplitude of the interface asymptotes to a constant value. In the framework of incompressible MHD, we examine analytically the behavior of an impulsively accelerated interface separating conducting fluids of different densities in cylindrical geometry. We investigate the stability properties of such a system and study the influence of the magnetic field on the growth rate of the interface. In converging cylindrical geometry, the RMI is followed by a Rayleigh-Taylor (RT) phase. Our analysis does not account for the RT phase of the instability but is valid for the duration of the RMI phase. We compare results of the incompressible analysis with linear compressible MHD simulations. Supported by the KAUST Office of Competitive Research Funds under Award No. URF/1/2162-01.
A transient, quadratic nodal method for triangular-Z geometry
DeLorey, T.F.
1993-06-01
Many systematically-derived nodal methods have been developed for Cartesian geometry due to the extensive interest in Light Water Reactors. These methods typically model the transverse-integrated flux as either an analytic or low order polynomial function of position within the node. Recently, quadratic nodal methods have been developed for R-Z and hexagonal geometry. A static and transient quadratic nodal method is developed for triangular-Z geometry. This development is particularly challenging because the quadratic expansion in each node must be performed between the node faces and the triangular points. As a consequence, in the 2-D plane, the flux and current at the points of the triangles must be treated. Quadratic nodal equations are solved using a non-linear iteration scheme, which utilizes the corrected, mesh-centered finite difference equations, and forces these equations to match the quadratic equations by computing discontinuity factors during the solution. Transient nodal equations are solved using the improved quasi-static method, which has been shown to be a very efficient solution method for transient problems. Several static problems are used to compare the quadratic nodal method to the Coarse Mesh Finite Difference (CMFD) method. The quadratic method is shown to give more accurate node-averaged fluxes. However, it appears that the method has difficulty predicting node leakages near reactor boundaries and severe material interfaces. The consequence is that the eigenvalue may be poorly predicted for certain reactor configurations. The transient methods are tested using a simple analytic test problem, a heterogeneous heavy water reactor benchmark problem, and three thermal hydraulic test problems. Results indicate that the transient methods have been implemented correctly.
On the phase diagram of 2d Lorentzian Quantum Gravity
NASA Astrophysics Data System (ADS)
Ambjørn, Jan; Anagnostopoulos, K. N.; Loll, R.
The phase diagram of 2d Lorentzian quantum gravity (LQG) coupled to conformal matter is studied. A phase transition is observed at c = c crit ( {1}/{2} < c crit < 4) which can be thought of as the analogue of the c = 1 barrier of Euclidean quantum gravity (EQG). The non-trivial properties of the quantum geometry are discussed.
NASA Astrophysics Data System (ADS)
Mayor, Louise
2016-05-01
Graphene might be the most famous example, but there are other 2D materials and compounds too. Louise Mayor explains how these atomically thin sheets can be layered together to create flexible “van der Waals heterostructures”, which could lead to a range of novel applications.
Double-diffusive magnetic buoyancy instability in a quasi-two-dimensional Cartesian geometry
NASA Astrophysics Data System (ADS)
Skinner, D. M.; Silvers, L. J.
2013-11-01
Magnetic buoyancy, believed to occur in the solar tachocline, is both an important part of large-scale solar dynamo models and the picture of how sunspots are formed. Given that in the tachocline region the ratio of magnetic diffusivity to thermal diffusivity is small it is important, for both the dynamo and sunspot formation pictures, to understand magnetic buoyancy in this regime. Furthermore, the tachocline is a region of strong shear and such investigations must involve structures that become buoyant in the double-diffusive regime which are generated entirely from a shear flow. In a previous study, we have illustrated that shear-generated double-diffusive magnetic buoyancy instability is possible in the tachocline. However, this study was severely limited due to the computational requirements of running three-dimensional magnetohydrodynamic simulations over diffusive time-scales. A more comprehensive investigation is required to fully understand the double-diffusive magnetic buoyancy instability and its dependency on a number of key parameters; such an investigation requires the consideration of a reduced model. Here we consider a quasi-two-dimensional model where all gradients in the x direction are set to zero. We show how the instability is sensitive to changes in the thermal diffusivity and also show how different initial configurations of the forced shear flow affect the behaviour of the instability. Finally, we conclude that if the tachocline is thinner than currently stated then the double-diffusive magnetic buoyancy instability can more easily occur.
SU-E-I-41: Non-Cartesian MR Image Reconstruction with Integrated Gradient Non-Linearity Correction
Tao, S; Trzasko, JD; Polley, TW; Shu, Y; Bernstein, MA
2014-06-01
Purpose: Nonlinearities in the spatial encoding gradients of MRI systems cause geometric distortion in images. Typically, this is retrospectively corrected via image-domain interpolation (a.k.a., “gradwarp”) albeit with a loss of spatial resolution. For non-Cartesian MRI, the latter problem is exaggerated by noise and undersampling artifact. In this study, we describe a novel correction strategy that accounts for gradient nonlinearities during — rather than after — non-Cartesian MRI reconstruction, and demonstrate that this approach mitigates the resolution loss that can occur with standard methods. Methods: To test the proposed method, the American College of Radiology (ACR) quality control phantom was scanned on at 1.5 T (General Electric, v16.0, “zoom” gradient) using a 1.6x undersampled 3D non- Cartesian Shells trajectory (GRE, FOV=24 cm3, 120 shells, 16552 shots, 512 readout, matrix=2403). Image reconstruction was first performed via standard k-space density-compensated gridding and retrospectively corrected via cubic spline interpolation. Image reconstruction was then separately performed using a k-space and image-domain densitycompensated type-3 non-uniform fast Fourier transform (NUFFT), which provides a direct mapping between non-Cartesian k-space samples and warped image space voxel locations. Thus, no separate distortion correction procedure is needed for the proposed approach. The gradient distortion field was determined using vendor provided calibration data. Results: Phantom scan results show that both processing approaches successfully correct geometric distortion. However, visual inspection of the ACR phantom spatial resolution inserts shows that the proposed strategy preserves the resolution of the nominal (uncorrected) reconstruction while “gradwarp” imparts marked spatial blurring (especially for the 1.0 and 1.1 mm inserts) and thus resolution loss. Conclusion: We've presented a novel reconstruction strategy for non-Cartesian MRI
Resistivity inversion in 2-D anisotropic media: numerical experiments
NASA Astrophysics Data System (ADS)
Wiese, Timothy; Greenhalgh, Stewart; Zhou, Bing; Greenhalgh, Mark; Marescot, Laurent
2015-04-01
Many rocks and layered/fractured sequences have a clearly expressed electrical anisotropy although it is rare in practice to incorporate anisotropy into resistivity inversion. In this contribution, we present a series of 2.5-D synthetic inversion experiments for various electrode configurations and 2-D anisotropic models. We examine and compare the image reconstructions obtained using the correct anisotropic inversion code with those obtained using the false but widely used isotropic assumption. Superior reconstruction in terms of reduced data misfit, true anomaly shape and position, and anisotropic background parameters were obtained when the correct anisotropic assumption was employed for medium to high coefficients of anisotropy. However, for low coefficient values the isotropic assumption produced better-quality results. When an erroneous isotropic inversion is performed on medium to high level anisotropic data, the images are dominated by patterns of banded artefacts and high data misfits. Various pole-pole, pole-dipole and dipole-dipole data sets were investigated and evaluated for the accuracy of the inversion result. The eigenvalue spectra of the pseudo-Hessian matrix and the formal resolution matrix were also computed to determine the information content and goodness of the results. We also present a data selection strategy based on high sensitivity measurements which drastically reduces the number of data to be inverted but still produces comparable results to that of the comprehensive data set. Inversion was carried out using transversely isotropic model parameters described in two different co-ordinate frames for the conductivity tensor, namely Cartesian versus natural or eigenframe. The Cartesian frame provided a more stable inversion product. This can be simply explained from inspection of the eigenspectra of the pseudo-Hessian matrix for the two model descriptions.
Nested Cartesian grid method in incompressible viscous fluid flow
NASA Astrophysics Data System (ADS)
Peng, Yih-Ferng; Mittal, Rajat; Sau, Amalendu; Hwang, Robert R.
2010-09-01
In this work, the local grid refinement procedure is focused by using a nested Cartesian grid formulation. The method is developed for simulating unsteady viscous incompressible flows with complex immersed boundaries. A finite-volume formulation based on globally second-order accurate central-difference schemes is adopted here in conjunction with a two-step fractional-step procedure. The key aspects that needed to be considered in developing such a nested grid solver are proper imposition of interface conditions on the nested-block boundaries, and accurate discretization of the governing equations in cells that are with block-interface as a control-surface. The interpolation procedure adopted in the study allows systematic development of a discretization scheme that preserves global second-order spatial accuracy of the underlying solver, and as a result high efficiency/accuracy nested grid discretization method is developed. Herein the proposed nested grid method has been widely tested through effective simulation of four different classes of unsteady incompressible viscous flows, thereby demonstrating its performance in the solution of various complex flow-structure interactions. The numerical examples include a lid-driven cavity flow and Pearson vortex problems, flow past a circular cylinder symmetrically installed in a channel, flow past an elliptic cylinder at an angle of attack, and flow past two tandem circular cylinders of unequal diameters. For the numerical simulations of flows past bluff bodies an immersed boundary (IB) method has been implemented in which the solid object is represented by a distributed body force in the Navier-Stokes equations. The main advantages of the implemented immersed boundary method are that the simulations could be performed on a regular Cartesian grid and applied to multiple nested-block (Cartesian) structured grids without any difficulty. Through the numerical experiments the strength of the solver in effectively
Claes Hellerström and Cartesian diver microrespirometry
Welsh, Michael
2016-01-01
Cartesian diver microrespirometry was introduced by Claes Hellerström at the Department of Histology/Medical Cell Biology at Uppsala University, Sweden, to determine rates of oxygen consumption in islets of Langerhans. The theory behind this method is touched upon and the main findings described. Glucose-stimulated beta cell respiration significantly contributes to increased ATP generation, which is a prerequisite for stimulated insulin secretion and synthesis. This has had major implications for understanding the beta cell stimulus–secretion coupling. PMID:27181825
Ion Transport in 2-D Graphene Nanochannels
NASA Astrophysics Data System (ADS)
Xie, Quan; Foo, Elbert; Duan, Chuanhua
2015-11-01
Graphene membranes have recently attracted wide attention due to its great potential in water desalination and selective molecular sieving. Further developments of these membranes, including enhancing their mass transport rate and/or molecular selectivity, rely on the understanding of fundamental transport mechanisms through graphene membranes, which has not been studied experimentally before due to fabrication and measurement difficulties. Herein we report the fabrication of the basic constituent of graphene membranes, i.e. 2-D single graphene nanochannels (GNCs) and the study of ion transport in these channels. A modified bonding technique was developed to form GNCs with well-defined geometry and uniform channel height. Ion transport in such GNCs was studied using DC conductance measurement. Our preliminary results showed that the ion transport in GNCs is still governed by surface charge at low concentrations (10-6M to 10-4M). However, GNCs exhibits much higher ionic conductances than silica nanochannels with the same geometries in the surface-charge-governed regime. This conductance enhancement can be attributed to the pre-accumulation of charges on graphene surfaces. The work is supported by the Faculty Startup Fund (Boston University, USA).
Energy Science and Technology Software Center (ESTSC)
2001-01-31
This software reduces the data from two-dimensional kSA MOS program, k-Space Associates, Ann Arbor, MI. Initial MOS data is recorded without headers in 38 columns, with one row of data per acquisition per lase beam tracked. The final MOSS 2d data file is reduced, graphed, and saved in a tab-delimited column format with headers that can be plotted in any graphing software.
Maintain rigid structures in Verlet based Cartesian molecular dynamics simulations
NASA Astrophysics Data System (ADS)
Tao, Peng; Wu, Xiongwu; Brooks, Bernard R.
2012-10-01
An algorithm is presented to maintain rigid structures in Verlet based Cartesian molecular dynamics (MD) simulations. After each unconstrained MD step, the coordinates of selected particles are corrected to maintain rigid structures through an iterative procedure of rotation matrix computation. This algorithm, named as SHAPE and implemented in CHARMM program suite, avoids the calculations of Lagrange multipliers, so that the complexity of computation does not increase with the number of particles in a rigid structure. The implementation of this algorithm does not require significant modification of propagation integrator, and can be plugged into any Cartesian based MD integration scheme. A unique feature of the SHAPE method is that it is interchangeable with SHAKE for any object that can be constrained as a rigid structure using multiple SHAKE constraints. Unlike SHAKE, the SHAPE method can be applied to large linear (with three or more centers) and planar (with four or more centers) rigid bodies. Numerical tests with four model systems including two proteins demonstrate that the accuracy and reliability of the SHAPE method are comparable to the SHAKE method, but with much more applicability and efficiency.
[Cartesian misunderstanding as a cause of therapeutic failure].
Isler, H
1986-01-01
Headache patients disassociate themselves from their own automatic responses, relying on the traditional separation of body and mind. On the other hand, patients who obtain voluntary control of automatic functions by biofeedback training modify not only vegetative but also voluntary behaviour patterns, losing "neurotic" traits. The basic misconception of the separation of body and mind, Cartesian dualism, is now ingrained in our culture. In the 17th century Descartes asserted that concepts applied to the soul must be entirely different from those used for the body in order to improve comprehension of the immortality of the soul. This dualism also led to "enlightenment" and to many later social and philosophical developments. But his basic neurophysiology was obsolete when he wrote it down. Other models from mainstream natural philosophy were better compatible with observation and experiments. Gassendi assumed a "body soul" consisting of energy as the functional principle of the nervous system, and Willis accommodated a series of anticipations of 19th century discoveries within this model. No comparable progress resulted from Descartes' own medieval model. Cartesian dualism has become untenable in view of recent neuropsychology but it still obstructs our management of functional patients. Instead of reinforcing the delusion of separation of psyche and soma, we ought to encourage patients to understand that their malfunctioning organs are on-line with their emotions, and with their mind. PMID:2420000
Nanoimprint lithography: 2D or not 2D? A review
NASA Astrophysics Data System (ADS)
Schift, Helmut
2015-11-01
Nanoimprint lithography (NIL) is more than a planar high-end technology for the patterning of wafer-like substrates. It is essentially a 3D process, because it replicates various stamp topographies by 3D displacement of material and takes advantage of the bending of stamps while the mold cavities are filled. But at the same time, it keeps all assets of a 2D technique being able to pattern thin masking layers like in photon- and electron-based traditional lithography. This review reports about 20 years of development of replication techniques at Paul Scherrer Institut, with a focus on 3D aspects of molding, which enable NIL to stay 2D, but at the same time enable 3D applications which are "more than Moore." As an example, the manufacturing of a demonstrator for backlighting applications based on thermally activated selective topography equilibration will be presented. This technique allows generating almost arbitrary sloped, convex and concave profiles in the same polymer film with dimensions in micro- and nanometer scale.
A general time element for orbit integration in Cartesian coordinates
NASA Technical Reports Server (NTRS)
Janin, G.; Bond, V. R.
1981-01-01
Two techniques are discussed for increasing the accuracy of the numerical integration of eccentric orbits in Cartesian coordinates. One involves the use of an independent variable different from time; this increases the efficiency of the numerical integration. The other uses a time element, which reduces the in-track error. A general expression is given of a time element valid for an arbitrary independent variable. It is pointed out that this time element makes it possible to switch the independent variable merely by applying a scaling factor; there is no need to change the differential equations of the motion. Eccentric, true, and elliptic anomalies are used as independent variables in the case of a transfer orbit for a geosynchronous orbit. The elliptic anomaly is shown to perform much better than the other classical anomalies.
Direct adaptive control of manipulators in Cartesian space
NASA Technical Reports Server (NTRS)
Seraji, H.
1987-01-01
A new adaptive-control scheme for direct control of manipulator end effector to achieve trajectory tracking in Cartesian space is developed in this article. The control structure is obtained from linear multivariable theory and is composed of simple feedforward and feedback controllers and an auxiliary input. The direct adaptation laws are derived from model reference adaptive control theory and are not based on parameter estimation of the robot model. The utilization of adaptive feedforward control and the inclusion of auxiliary input are novel features of the present scheme and result in improved dynamic performance over existing adaptive control schemes. The adaptive controller does not require the complex mathematical model of the robot dynamics or any knowledge of the robot parameters or the payload, and is computationally fast for on-line implementation with high sampling rates. The control scheme is applied to a two-link manipulator for illustration.
NASA Technical Reports Server (NTRS)
Finley, Dennis B.; Karman, Steve L., Jr.
1996-01-01
The objective of the second phase of the Euler Technology Assessment program was to evaluate the ability of Euler computational fluid dynamics codes to predict compressible flow effects over a generic fighter wind tunnel model. This portion of the study was conducted by Lockheed Martin Tactical Aircraft Systems, using an in-house Cartesian-grid code called SPLITFLOW. The Cartesian grid technique offers several advantages, including ease of volume grid generation and reduced number of cells compared to other grid schemes. SPLITFLOW also includes grid adaption of the volume grid during the solution to resolve high-gradient regions. The SPLITFLOW code predictions of configuration forces and moments are shown to be adequate for preliminary design, including predictions of sideslip effects and the effects of geometry variations at low and high angles-of-attack. The transonic pressure prediction capabilities of SPLITFLOW are shown to be improved over subsonic comparisons. The time required to generate the results from initial surface data is on the order of several hours, including grid generation, which is compatible with the needs of the design environment.
NASA Astrophysics Data System (ADS)
Miao, Sha; Hendrickson, Kelli; Liu, Yuming; Subramani, Hariprasad
2015-11-01
This work presents a novel and efficient Cartesian-grid based simulation capability for the study of an incompressible, turbulent gas layer over a liquid flow with disparate Reynolds numbers in two phases. This capability couples a turbulent gas-flow solver and a liquid-layer based on a second-order accurate Boundary Data Immersion Method (BDIM) at the deformable interface. The turbulent gas flow solver solves the incompressible Navier-Stokes equations via direct numerical simulation or through turbulence closure (unsteady Reynolds-Averaged Navier-Stokes Models) for Reynolds numbers O(106). In this application, a laminar liquid layer solution is obtained from depth-integrated Navier-Stokes equations utilizing shallow water wave assumptions. The immersed boundary method (BDIM) enforces the coupling at the deformable interface, the boundary conditions to turbulence closure equations and defines the domain geometry on the Cartesian grid. Validations are made for the turbulent gas channel flow over high-viscosity liquid. This simulation capability can be applied to problems in the oil and industrial sector such as channel and pipe flows with heavy oils as well as wind wave generation in shallow waters. Sponsored by the Chevron Energy Technology Company.
Multi-fault Tolerance for Cartesian Data Distributions
Ali, Nawab; Krishnamoorthy, Sriram; Halappanavar, Mahantesh; Daily, Jeffrey A.
2013-06-01
Faults are expected to play an increasingly important role in how algorithms and applications are designed to run on future extreme-scale sys- tems. Algorithm-based fault tolerance (ABFT) is a promising approach that involves modications to the algorithm to recover from faults with lower over- heads than replicated storage and a signicant reduction in lost work compared to checkpoint-restart techniques. Fault-tolerant linear algebra (FTLA) algo- rithms employ additional processors that store parities along the dimensions of a matrix to tolerate multiple, simultaneous faults. Existing approaches as- sume regular data distributions (blocked or block-cyclic) with the failures of each data block being independent. To match the characteristics of failures on parallel computers, we extend these approaches to mapping parity blocks in several important ways. First, we handle parity computation for generalized Cartesian data distributions with each processor holding arbitrary subsets of blocks in a Cartesian-distributed array. Second, techniques to handle corre- lated failures, i.e., multiple processors that can be expected to fail together, are presented. Third, we handle the colocation of parity blocks with the data blocks and do not require them to be on additional processors. Several al- ternative approaches, based on graph matching, are presented that attempt to balance the memory overhead on processors while guaranteeing the same fault tolerance properties as existing approaches that assume independent fail- ures on regular blocked data distributions. The evaluation of these algorithms demonstrates that the additional desirable properties are provided by the pro- posed approach with minimal overhead.
NASA Technical Reports Server (NTRS)
Aftosmis, M. J.; Berger, M. J.; Adomavicius, G.
2000-01-01
Preliminary verification and validation of an efficient Euler solver for adaptively refined Cartesian meshes with embedded boundaries is presented. The parallel, multilevel method makes use of a new on-the-fly parallel domain decomposition strategy based upon the use of space-filling curves, and automatically generates a sequence of coarse meshes for processing by the multigrid smoother. The coarse mesh generation algorithm produces grids which completely cover the computational domain at every level in the mesh hierarchy. A series of examples on realistically complex three-dimensional configurations demonstrate that this new coarsening algorithm reliably achieves mesh coarsening ratios in excess of 7 on adaptively refined meshes. Numerical investigations of the scheme's local truncation error demonstrate an achieved order of accuracy between 1.82 and 1.88. Convergence results for the multigrid scheme are presented for both subsonic and transonic test cases and demonstrate W-cycle multigrid convergence rates between 0.84 and 0.94. Preliminary parallel scalability tests on both simple wing and complex complete aircraft geometries shows a computational speedup of 52 on 64 processors using the run-time mesh partitioner.
NASA Technical Reports Server (NTRS)
Aftosmis, M. J.; Berger, M. J.; Adomavicius, G.; Nixon, David (Technical Monitor)
1998-01-01
The work presents a new method for on-the-fly domain decomposition technique for mapping grids and solution algorithms to parallel machines, and is applicable to both shared-memory and message-passing architectures. It will be demonstrated on the Cray T3E, HP Exemplar, and SGI Origin 2000. Computing time has been secured on all these platforms. The decomposition technique is an outgrowth of techniques used in computational physics for simulations of N-body problems and the event horizons of black holes, and has not been previously used by the CFD community. Since the technique offers on-the-fly partitioning, it offers a substantial increase in flexibility for computing in heterogeneous environments, where the number of available processors may not be known at the time of job submission. In addition, since it is dynamic it permits the job to be repartitioned without global communication in cases where additional processors become available after the simulation has begun, or in cases where dynamic mesh adaptation changes the mesh size during the course of a simulation. The platform for this partitioning strategy is a completely new Cartesian Euler solver tarcreted at parallel machines which may be used in conjunction with Ames' "Cart3D" arbitrary geometry simulation package.
Viability of Bioprinted Cellular Constructs Using a Three Dispenser Cartesian Printer
Dennis, SG.; Trusk, T.; Richards, D.; Jia, J.; Tan, Y.; Mei, Y.; Fann, S.; Markwald, R.; Yost, M.
2016-01-01
Tissue engineering has centralized its focus on the construction of replacements for non-functional or damaged tissue. The utilization of three-dimensional bioprinting in tissue engineering has generated new methods for the printing of cells and matrix to fabricate biomimetic tissue constructs. The solid freeform fabrication (SFF) method developed for three-dimensional bioprinting uses an additive manufacturing approach by depositing droplets of cells and hydrogels in a layer-by-layer fashion. Bioprinting fabrication is dependent on the specific placement of biological materials into three-dimensional architectures, and the printed constructs should closely mimic the complex organization of cells and extracellular matrices in native tissue. This paper highlights the use of the Palmetto Printer, a Cartesian bioprinter, as well as the process of producing spatially organized, viable constructs while simultaneously allowing control of environmental factors. This methodology utilizes computer-aided design and computer-aided manufacturing to produce these specific and complex geometries. Finally, this approach allows for the reproducible production of fabricated constructs optimized by controllable printing parameters. PMID:26436877
Viability of Bioprinted Cellular Constructs Using a Three Dispenser Cartesian Printer.
Dennis, Sarah Grace; Trusk, Thomas; Richards, Dylan; Jia, Jia; Tan, Yu; Mei, Ying; Fann, Stephen; Markwald, Roger; Yost, Michael
2015-01-01
Tissue engineering has centralized its focus on the construction of replacements for non-functional or damaged tissue. The utilization of three-dimensional bioprinting in tissue engineering has generated new methods for the printing of cells and matrix to fabricate biomimetic tissue constructs. The solid freeform fabrication (SFF) method developed for three-dimensional bioprinting uses an additive manufacturing approach by depositing droplets of cells and hydrogels in a layer-by-layer fashion. Bioprinting fabrication is dependent on the specific placement of biological materials into three-dimensional architectures, and the printed constructs should closely mimic the complex organization of cells and extracellular matrices in native tissue. This paper highlights the use of the Palmetto Printer, a Cartesian bioprinter, as well as the process of producing spatially organized, viable constructs while simultaneously allowing control of environmental factors. This methodology utilizes computer-aided design and computer-aided manufacturing to produce these specific and complex geometries. Finally, this approach allows for the reproducible production of fabricated constructs optimized by controllable printing parameters. PMID:26436877
Extending a CAD-Based Cartesian Mesh Generator for the Lattice Boltzmann Method
Cantrell, J Nathan; Inclan, Eric J; Joshi, Abhijit S; Popov, Emilian L; Jain, Prashant K
2012-01-01
This paper describes the development of a custom preprocessor for the PaRAllel Thermal Hydraulics simulations using Advanced Mesoscopic methods (PRATHAM) code based on an open-source mesh generator, CartGen [1]. PRATHAM is a three-dimensional (3D) lattice Boltzmann method (LBM) based parallel flow simulation software currently under development at the Oak Ridge National Laboratory. The LBM algorithm in PRATHAM requires a uniform, coordinate system-aligned, non-body-fitted structured mesh for its computational domain. CartGen [1], which is a GNU-licensed open source code, already comes with some of the above needed functionalities. However, it needs to be further extended to fully support the LBM specific preprocessing requirements. Therefore, CartGen is being modified to (i) be compiler independent while converting a neutral-format STL (Stereolithography) CAD geometry to a uniform structured Cartesian mesh, (ii) provide a mechanism for PRATHAM to import the mesh and identify the fluid/solid domains, and (iii) provide a mechanism to visually identify and tag the domain boundaries on which to apply different boundary conditions.
Competing coexisting phases in 2D water
Zanotti, Jean-Marc; Judeinstein, Patrick; Dalla-Bernardina, Simona; Creff, Gaëlle; Brubach, Jean-Blaise; Roy, Pascale; Bonetti, Marco; Ollivier, Jacques; Sakellariou, Dimitrios; Bellissent-Funel, Marie-Claire
2016-01-01
The properties of bulk water come from a delicate balance of interactions on length scales encompassing several orders of magnitudes: i) the Hydrogen Bond (HBond) at the molecular scale and ii) the extension of this HBond network up to the macroscopic level. Here, we address the physics of water when the three dimensional extension of the HBond network is frustrated, so that the water molecules are forced to organize in only two dimensions. We account for the large scale fluctuating HBond network by an analytical mean-field percolation model. This approach provides a coherent interpretation of the different events experimentally (calorimetry, neutron, NMR, near and far infra-red spectroscopies) detected in interfacial water at 160, 220 and 250 K. Starting from an amorphous state of water at low temperature, these transitions are respectively interpreted as the onset of creation of transient low density patches of 4-HBonded molecules at 160 K, the percolation of these domains at 220 K and finally the total invasion of the surface by them at 250 K. The source of this surprising behaviour in 2D is the frustration of the natural bulk tetrahedral local geometry and the underlying very significant increase in entropy of the interfacial water molecules. PMID:27185018
Parallel map analysis on 2-D grids
Berry, M.; Comiskey, J.; Minser, K.
1993-12-31
In landscape ecology, computer modeling is used to assess habitat fragmentation and its ecological iMPLications. Specifically, maps (2-D grids) of habitat clusters must be analyzed to determine number, sizes and geometry of clusters. Models prior to this study relied upon sequential Fortran-77 programs which limited the sizes of maps and densities of clusters which could be analyzed. In this paper, we present more efficient computer models which can exploit recursion or parallelism. Significant improvements over the original Fortran-77 programs have been achieved using both recursive and nonrecursive C implementations on a variety of workstations such as the Sun Sparc 2, IBM RS/6000-350, and HP 9000-750. Parallel implementations on a 4096-processor MasPar MP-1 and a 32-processor CM-5 are also studied. Preliminary experiments suggest that speed improvements for the parallel model on the MasPar MP-1 (written in MPL) and on the CM-5 (written in C using CMMD) can be as much as 39 and 34 times faster, respectively, than the most efficient sequential C program on a Sun Sparc 2 for a 512 map. An important goal in this research effort is to produce a scalable map analysis algorithm for the identification and characterization of clusters for relatively large maps on massively-parallel computers.
Competing coexisting phases in 2D water
NASA Astrophysics Data System (ADS)
Zanotti, Jean-Marc; Judeinstein, Patrick; Dalla-Bernardina, Simona; Creff, Gaëlle; Brubach, Jean-Blaise; Roy, Pascale; Bonetti, Marco; Ollivier, Jacques; Sakellariou, Dimitrios; Bellissent-Funel, Marie-Claire
2016-05-01
The properties of bulk water come from a delicate balance of interactions on length scales encompassing several orders of magnitudes: i) the Hydrogen Bond (HBond) at the molecular scale and ii) the extension of this HBond network up to the macroscopic level. Here, we address the physics of water when the three dimensional extension of the HBond network is frustrated, so that the water molecules are forced to organize in only two dimensions. We account for the large scale fluctuating HBond network by an analytical mean-field percolation model. This approach provides a coherent interpretation of the different events experimentally (calorimetry, neutron, NMR, near and far infra-red spectroscopies) detected in interfacial water at 160, 220 and 250 K. Starting from an amorphous state of water at low temperature, these transitions are respectively interpreted as the onset of creation of transient low density patches of 4-HBonded molecules at 160 K, the percolation of these domains at 220 K and finally the total invasion of the surface by them at 250 K. The source of this surprising behaviour in 2D is the frustration of the natural bulk tetrahedral local geometry and the underlying very significant increase in entropy of the interfacial water molecules.
Competing coexisting phases in 2D water.
Zanotti, Jean-Marc; Judeinstein, Patrick; Dalla-Bernardina, Simona; Creff, Gaëlle; Brubach, Jean-Blaise; Roy, Pascale; Bonetti, Marco; Ollivier, Jacques; Sakellariou, Dimitrios; Bellissent-Funel, Marie-Claire
2016-01-01
The properties of bulk water come from a delicate balance of interactions on length scales encompassing several orders of magnitudes: i) the Hydrogen Bond (HBond) at the molecular scale and ii) the extension of this HBond network up to the macroscopic level. Here, we address the physics of water when the three dimensional extension of the HBond network is frustrated, so that the water molecules are forced to organize in only two dimensions. We account for the large scale fluctuating HBond network by an analytical mean-field percolation model. This approach provides a coherent interpretation of the different events experimentally (calorimetry, neutron, NMR, near and far infra-red spectroscopies) detected in interfacial water at 160, 220 and 250 K. Starting from an amorphous state of water at low temperature, these transitions are respectively interpreted as the onset of creation of transient low density patches of 4-HBonded molecules at 160 K, the percolation of these domains at 220 K and finally the total invasion of the surface by them at 250 K. The source of this surprising behaviour in 2D is the frustration of the natural bulk tetrahedral local geometry and the underlying very significant increase in entropy of the interfacial water molecules. PMID:27185018
Crockett, Robert; Graves, Daniel; Colella, Phillip
2009-10-23
We present a method for solving Poisson and heat equations with discon- tinuous coefficients in two- and three-dimensions. It uses a Cartesian cut-cell/embedded boundary method to represent the interface between materi- als, as described in Johansen& Colella (1998). Matching conditions across the interface are enforced using an approximation to fluxes at the boundary. Overall second order accuracy is achieved, as indicated by an array of tests using non-trivial interface geometries. Both the elliptic and heat solvers are shown to remain stable and efficient for material coefficient contrasts up to 106, thanks in part to the use of geometric multigrid. A test of accuracy when adaptive mesh refinement capabilities are utilized is also performed. An example problem relevant to nuclear reactor core simulation is presented, demonstrating the ability of the method to solve problems with realistic physical parameters.
Conversion of Cartesian coordinates from and to Generalized Balanced Ternary addresses
van Roessel, Jan W.
1988-01-01
Hexagonal grids have several advantages over square grids, such as a greater angular resolution and unambiguous connectivity. The Generalized Balanced Ternary (GBT) system is a spatial addressing method for hexagonal grids in which the hexagons are arranged in hierarchical aggregates, and which accommodates vector operations in GBT space. Efficient algorithms for converting Cartesian coordinates from and to GBT addresses are based on the dual representation of the hexagonal tessellation. The GBT-to-Cartesian algorithm is an order of magnitude faster than the Cartesian-to-GBT algorithm, the latter requiring interpolation and GBT addition for each digit of the generated GBT address.
On NUFFT-based gridding for non-Cartesian MRI
Fessler, Jeffrey A.
2007-01-01
For MRI with non-Cartesian sampling, the conventional approach to reconstructing images is to use the gridding method with a Kaiser-Bessel (KB) interpolation kernel. Recently, Sha et al. [1] proposed an alternative method based on a nonuniform FFT (NUFFT) with least-squares (LS) design of the interpolation coefficients. They described this LS_NUFFT method as shift variant and reported that it yielded smaller reconstruction approximation errors than the conventional shift-invariant KB approach. This paper analyzes the LS_NUFFT approach in detail. We show that when one accounts for a certain linear phase factor, the core of the LS_NUFFT interpolator is in fact real and shift invariant. Furthermore, we find that the KB approach yields smaller errors than the original LS_NUFFT approach. We show that optimizing certain scaling factors can lead to a somewhat improved LS_NUFFT approach, but the high computation cost seems to outweigh the modest reduction in reconstruction error. We conclude that the standard KB approach, with appropriate parameters as described in the literature, remains the practical method of choice for gridding reconstruction in MRI. PMID:17689121
Shared Memory Parallelism for 3D Cartesian Discrete Ordinates Solver
NASA Astrophysics Data System (ADS)
Moustafa, Salli; Dutka-Malen, Ivan; Plagne, Laurent; Ponçot, Angélique; Ramet, Pierre
2014-06-01
This paper describes the design and the performance of DOMINO, a 3D Cartesian SN solver that implements two nested levels of parallelism (multicore+SIMD) on shared memory computation nodes. DOMINO is written in C++, a multi-paradigm programming language that enables the use of powerful and generic parallel programming tools such as Intel TBB and Eigen. These two libraries allow us to combine multi-thread parallelism with vector operations in an efficient and yet portable way. As a result, DOMINO can exploit the full power of modern multi-core processors and is able to tackle very large simulations, that usually require large HPC clusters, using a single computing node. For example, DOMINO solves a 3D full core PWR eigenvalue problem involving 26 energy groups, 288 angular directions (S16), 46 × 106 spatial cells and 1 × 1012 DoFs within 11 hours on a single 32-core SMP node. This represents a sustained performance of 235 GFlops and 40:74% of the SMP node peak performance for the DOMINO sweep implementation. The very high Flops/Watt ratio of DOMINO makes it a very interesting building block for a future many-nodes nuclear simulation tool.
A multilevel Cartesian non-uniform grid time domain algorithm
Meng Jun; Boag, Amir; Lomakin, Vitaliy; Michielssen, Eric
2010-11-01
A multilevel Cartesian non-uniform grid time domain algorithm (CNGTDA) is introduced to rapidly compute transient wave fields radiated by time dependent three-dimensional source constellations. CNGTDA leverages the observation that transient wave fields generated by temporally bandlimited and spatially confined source constellations can be recovered via interpolation from appropriately delay- and amplitude-compensated field samples. This property is used in conjunction with a multilevel scheme, in which the computational domain is hierarchically decomposed into subdomains with sparse non-uniform grids used to obtain the fields. For both surface and volumetric source distributions, the computational cost of CNGTDA to compute the transient field at N{sub s} observation locations from N{sub s} collocated sources for N{sub t} discrete time instances scales as O(N{sub t}N{sub s}logN{sub s}) and O(N{sub t}N{sub s}log{sup 2}N{sub s}) in the low- and high-frequency regimes, respectively. Coupled with marching-on-in-time (MOT) time domain integral equations, CNGTDA can facilitate efficient analysis of large scale time domain electromagnetic and acoustic problems.
On NUFFT-based gridding for non-Cartesian MRI
NASA Astrophysics Data System (ADS)
Fessler, Jeffrey A.
2007-10-01
For MRI with non-Cartesian sampling, the conventional approach to reconstructing images is to use the gridding method with a Kaiser-Bessel (KB) interpolation kernel. Recently, Sha et al. [L. Sha, H. Guo, A.W. Song, An improved gridding method for spiral MRI using nonuniform fast Fourier transform, J. Magn. Reson. 162(2) (2003) 250-258] proposed an alternative method based on a nonuniform FFT (NUFFT) with least-squares (LS) design of the interpolation coefficients. They described this LS_NUFFT method as shift variant and reported that it yielded smaller reconstruction approximation errors than the conventional shift-invariant KB approach. This paper analyzes the LS_NUFFT approach in detail. We show that when one accounts for a certain linear phase factor, the core of the LS_NUFFT interpolator is in fact real and shift invariant. Furthermore, we find that the KB approach yields smaller errors than the original LS_NUFFT approach. We show that optimizing certain scaling factors can lead to a somewhat improved LS_NUFFT approach, but the high computation cost seems to outweigh the modest reduction in reconstruction error. We conclude that the standard KB approach, with appropriate parameters as described in the literature, remains the practical method of choice for gridding reconstruction in MRI.
On NUFFT-based gridding for non-Cartesian MRI.
Fessler, Jeffrey A
2007-10-01
For MRI with non-Cartesian sampling, the conventional approach to reconstructing images is to use the gridding method with a Kaiser-Bessel (KB) interpolation kernel. Recently, Sha et al. [L. Sha, H. Guo, A.W. Song, An improved gridding method for spiral MRI using nonuniform fast Fourier transform, J. Magn. Reson. 162(2) (2003) 250-258] proposed an alternative method based on a nonuniform FFT (NUFFT) with least-squares (LS) design of the interpolation coefficients. They described this LS_NUFFT method as shift variant and reported that it yielded smaller reconstruction approximation errors than the conventional shift-invariant KB approach. This paper analyzes the LS_NUFFT approach in detail. We show that when one accounts for a certain linear phase factor, the core of the LS_NUFFT interpolator is in fact real and shift invariant. Furthermore, we find that the KB approach yields smaller errors than the original LS_NUFFT approach. We show that optimizing certain scaling factors can lead to a somewhat improved LS_NUFFT approach, but the high computation cost seems to outweigh the modest reduction in reconstruction error. We conclude that the standard KB approach, with appropriate parameters as described in the literature, remains the practical method of choice for gridding reconstruction in MRI. PMID:17689121
NKG2D ligands as therapeutic targets
Spear, Paul; Wu, Ming-Ru; Sentman, Marie-Louise; Sentman, Charles L.
2013-01-01
The Natural Killer Group 2D (NKG2D) receptor plays an important role in protecting the host from infections and cancer. By recognizing ligands induced on infected or tumor cells, NKG2D modulates lymphocyte activation and promotes immunity to eliminate ligand-expressing cells. Because these ligands are not widely expressed on healthy adult tissue, NKG2D ligands may present a useful target for immunotherapeutic approaches in cancer. Novel therapies targeting NKG2D ligands for the treatment of cancer have shown preclinical success and are poised to enter into clinical trials. In this review, the NKG2D receptor and its ligands are discussed in the context of cancer, infection, and autoimmunity. In addition, therapies targeting NKG2D ligands in cancer are also reviewed. PMID:23833565
Geometry of Thin Nematic Elastomer Sheets
NASA Astrophysics Data System (ADS)
Aharoni, Hillel; Sharon, Eran; Kupferman, Raz
A thin sheet of nematic elastomer attains 3D configurations depending on the nematic director field upon heating. In this talk we describe the intrinsic geometry of such a sheet, and derive an expression for the metric induced by general smooth nematic director fields. Furthermore, we investigate the reverse problem of constructing a director field that induces a specified 2D geometry. We provide an explicit analytical recipe for constructing any surface of revolution using this method. We demonstrate how the design of an arbitrary 2D geometry is accessible using approximate numerical methods.
3D Euler flow solutions using unstructured Cartesian and prismatic grids
NASA Technical Reports Server (NTRS)
Melton, John E.; Pandya, Shishir A.; Steger, Joseph L.
1993-01-01
A hyperbolic prismatic grid generation technique is combined with a background Cartesian grid for the study of inviscid three-dimensional flows. The mathematics of the hyperbolic prismatic grid generation algorithm are described, and some simple inviscid demonstration cases are presented. By combining the simplicity of the Cartesian background grid with the geometric flexibility and computational efficiencies inherent to prismatic grids, this approach shows promise for improving computational aerodynamic simulations.
Choi, Cheol Ho
2004-02-22
A new way of generating the multipole moments of Cartesian Gaussian functions in spherical polar coordinates has been established, bypassing the intermediary of Cartesian moment tensors. A new set of recurrence relations have also been derived for the resulting analytic integral values. The new method furnishes a conceptually simple and numerically efficient evaluation procedure for the multipole moments. The advantages over existing methods are documented. The results are relevant for the linear scaling quantum theories based on the fast multipole method. PMID:15268515
Rheological Properties of Quasi-2D Fluids in Microgravity
NASA Technical Reports Server (NTRS)
Stannarius, Ralf; Trittel, Torsten; Eremin, Alexey; Harth, Kirsten; Clark, Noel; Maclennan, Joseph; Glaser, Matthew; Park, Cheol; Hall, Nancy; Tin, Padetha
2015-01-01
In recent years, research on complex fluids and fluids in restricted geometries has attracted much attention in the scientific community. This can be attributed not only to the development of novel materials based on complex fluids but also to a variety of important physical phenomena which have barely been explored. One example is the behavior of membranes and thin fluid films, which can be described by two-dimensional (2D) rheology behavior that is quite different from 3D fluids. In this study, we have investigated the rheological properties of freely suspended films of a thermotropic liquid crystal in microgravity experiments. This model system mimics isotropic and anisotropic quasi 2D fluids [46]. We use inkjet printing technology to dispense small droplets (inclusions) onto the film surface. The motion of these inclusions provides information on the rheological properties of the films and allows the study of a variety of flow instabilities. Flat films have been investigated on a sub-orbital rocket flight and curved films (bubbles) have been studied in the ISS project OASIS. Microgravity is essential when the films are curved in order to avoid sedimentation. The experiments yield the mobility of the droplets in the films as well as the mutual mobility of pairs of particles. Experimental results will be presented for 2D-isotropic (smectic-A) and 2D-nematic (smectic-C) phases.
Multiscale geometric modeling of macromolecules I: Cartesian representation
NASA Astrophysics Data System (ADS)
Xia, Kelin; Feng, Xin; Chen, Zhan; Tong, Yiying; Wei, Guo-Wei
2014-01-01
This paper focuses on the geometric modeling and computational algorithm development of biomolecular structures from two data sources: Protein Data Bank (PDB) and Electron Microscopy Data Bank (EMDB) in the Eulerian (or Cartesian) representation. Molecular surface (MS) contains non-smooth geometric singularities, such as cusps, tips and self-intersecting facets, which often lead to computational instabilities in molecular simulations, and violate the physical principle of surface free energy minimization. Variational multiscale surface definitions are proposed based on geometric flows and solvation analysis of biomolecular systems. Our approach leads to geometric and potential driven Laplace-Beltrami flows for biomolecular surface evolution and formation. The resulting surfaces are free of geometric singularities and minimize the total free energy of the biomolecular system. High order partial differential equation (PDE)-based nonlinear filters are employed for EMDB data processing. We show the efficacy of this approach in feature-preserving noise reduction. After the construction of protein multiresolution surfaces, we explore the analysis and characterization of surface morphology by using a variety of curvature definitions. Apart from the classical Gaussian curvature and mean curvature, maximum curvature, minimum curvature, shape index, and curvedness are also applied to macromolecular surface analysis for the first time. Our curvature analysis is uniquely coupled to the analysis of electrostatic surface potential, which is a by-product of our variational multiscale solvation models. As an expository investigation, we particularly emphasize the numerical algorithms and computational protocols for practical applications of the above multiscale geometric models. Such information may otherwise be scattered over the vast literature on this topic. Based on the curvature and electrostatic analysis from our multiresolution surfaces, we introduce a new concept, the
Multiscale geometric modeling of macromolecules I: Cartesian representation
Xia, Kelin; Feng, Xin; Chen, Zhan; Tong, Yiying; Wei, Guo-Wei
2014-01-15
This paper focuses on the geometric modeling and computational algorithm development of biomolecular structures from two data sources: Protein Data Bank (PDB) and Electron Microscopy Data Bank (EMDB) in the Eulerian (or Cartesian) representation. Molecular surface (MS) contains non-smooth geometric singularities, such as cusps, tips and self-intersecting facets, which often lead to computational instabilities in molecular simulations, and violate the physical principle of surface free energy minimization. Variational multiscale surface definitions are proposed based on geometric flows and solvation analysis of biomolecular systems. Our approach leads to geometric and potential driven Laplace–Beltrami flows for biomolecular surface evolution and formation. The resulting surfaces are free of geometric singularities and minimize the total free energy of the biomolecular system. High order partial differential equation (PDE)-based nonlinear filters are employed for EMDB data processing. We show the efficacy of this approach in feature-preserving noise reduction. After the construction of protein multiresolution surfaces, we explore the analysis and characterization of surface morphology by using a variety of curvature definitions. Apart from the classical Gaussian curvature and mean curvature, maximum curvature, minimum curvature, shape index, and curvedness are also applied to macromolecular surface analysis for the first time. Our curvature analysis is uniquely coupled to the analysis of electrostatic surface potential, which is a by-product of our variational multiscale solvation models. As an expository investigation, we particularly emphasize the numerical algorithms and computational protocols for practical applications of the above multiscale geometric models. Such information may otherwise be scattered over the vast literature on this topic. Based on the curvature and electrostatic analysis from our multiresolution surfaces, we introduce a new concept, the
Multiscale geometric modeling of macromolecules I: Cartesian representation
Xia, Kelin; Feng, Xin; Chen, Zhan; Tong, Yiying; Wei, Guo Wei
2013-01-01
This paper focuses on the geometric modeling and computational algorithm development of biomolecular structures from two data sources: Protein Data Bank (PDB) and Electron Microscopy Data Bank (EMDB) in the Eulerian (or Cartesian) representation. Molecular surface (MS) contains non-smooth geometric singularities, such as cusps, tips and self-intersecting facets, which often lead to computational instabilities in molecular simulations, and violate the physical principle of surface free energy minimization. Variational multiscale surface definitions are proposed based on geometric flows and solvation analysis of biomolecular systems. Our approach leads to geometric and potential driven Laplace-Beltrami flows for biomolecular surface evolution and formation. The resulting surfaces are free of geometric singularities and minimize the total free energy of the biomolecular system. High order partial differential equation (PDE)-based nonlinear filters are employed for EMDB data processing. We show the efficacy of this approach in feature-preserving noise reduction. After the construction of protein multiresolution surfaces, we explore the analysis and characterization of surface morphology by using a variety of curvature definitions. Apart from the classical Gaussian curvature and mean curvature, maximum curvature, minimum curvature, shape index, and curvedness are also applied to macromolecular surface analysis for the first time. Our curvature analysis is uniquely coupled to the analysis of electrostatic surface potential, which is a by-product of our variational multiscale solvation models. As an expository investigation, we particularly emphasize the numerical algorithms and computational protocols for practical applications of the above multiscale geometric models. Such information may otherwise be scattered over the vast literature on this topic. Based on the curvature and electrostatic analysis from our multiresolution surfaces, we introduce a new concept, the
Multiscale geometric modeling of macromolecules I: Cartesian representation.
Xia, Kelin; Feng, Xin; Chen, Zhan; Tong, Yiying; Wei, Guo Wei
2014-01-01
This paper focuses on the geometric modeling and computational algorithm development of biomolecular structures from two data sources: Protein Data Bank (PDB) and Electron Microscopy Data Bank (EMDB) in the Eulerian (or Cartesian) representation. Molecular surface (MS) contains non-smooth geometric singularities, such as cusps, tips and self-intersecting facets, which often lead to computational instabilities in molecular simulations, and violate the physical principle of surface free energy minimization. Variational multiscale surface definitions are proposed based on geometric flows and solvation analysis of biomolecular systems. Our approach leads to geometric and potential driven Laplace-Beltrami flows for biomolecular surface evolution and formation. The resulting surfaces are free of geometric singularities and minimize the total free energy of the biomolecular system. High order partial differential equation (PDE)-based nonlinear filters are employed for EMDB data processing. We show the efficacy of this approach in feature-preserving noise reduction. After the construction of protein multiresolution surfaces, we explore the analysis and characterization of surface morphology by using a variety of curvature definitions. Apart from the classical Gaussian curvature and mean curvature, maximum curvature, minimum curvature, shape index, and curvedness are also applied to macromolecular surface analysis for the first time. Our curvature analysis is uniquely coupled to the analysis of electrostatic surface potential, which is a by-product of our variational multiscale solvation models. As an expository investigation, we particularly emphasize the numerical algorithms and computational protocols for practical applications of the above multiscale geometric models. Such information may otherwise be scattered over the vast literature on this topic. Based on the curvature and electrostatic analysis from our multiresolution surfaces, we introduce a new concept, the
Extension and application of the Preissmann slot model to 2D transient mixed flows
NASA Astrophysics Data System (ADS)
Maranzoni, Andrea; Dazzi, Susanna; Aureli, Francesca; Mignosa, Paolo
2015-08-01
This paper presents an extension of the Preissmann slot concept for the modeling of highly transient two-dimensional (2D) mixed flows. The classic conservative formulation of the 2D shallow water equations for free surface flows is adapted by assuming that two fictitious vertical slots, aligned along the two Cartesian plane directions and normally intersecting, are added on the ceiling of each integration element. Accordingly, transitions between free surface and pressurized flow can be handled in a natural and straightforward way by using the same set of governing equations. The opportunity of coupling free surface and pressurized flows is actually useful not only in one-dimensional (1D) problems concerning sewer systems but also for modeling 2D flooding phenomena in which the pressurization of bridges, culverts, or other crossing hydraulic structures can be expected. Numerical simulations are performed by using a shock-capturing MUSCL-Hancock finite volume scheme combined with the FORCE (First-Order Centred) solver for the evaluation of the numerical fluxes. The validation of the mathematical model is accomplished on the basis of both exact solutions of 1D discontinuous initial value problems and reference radial solutions of idealized test cases with cylindrical symmetry. Furthermore, the capability of the model to deal with practical field-scale applications is assessed by simulating the transit of a bore under an arch bridge. Numerical results show that the proposed model is suitable for the prediction of highly transient 2D mixed flows.
NASA Astrophysics Data System (ADS)
Civicioglu, Pinar
2012-09-01
In order to solve numerous practical navigational, geodetic and astro-geodetic problems, it is necessary to transform geocentric cartesian coordinates into geodetic coordinates or vice versa. It is very easy to solve the problem of transforming geodetic coordinates into geocentric cartesian coordinates. On the other hand, it is rather difficult to solve the problem of transforming geocentric cartesian coordinates into geodetic coordinates as it is very hard to define a mathematical relationship between the geodetic latitude (φ) and the geocentric cartesian coordinates (X, Y, Z). In this paper, a new algorithm, the Differential Search Algorithm (DS), is presented to solve the problem of transforming the geocentric cartesian coordinates into geodetic coordinates and its performance is compared with the performances of the classical methods (i.e., Borkowski, 1989; Bowring, 1976; Fukushima, 2006; Heikkinen, 1982; Jones, 2002; Zhang, 2005; Borkowski, 1987; Shu, 2010 and Lin, 1995) and Computational-Intelligence algorithms (i.e., ABC, JDE, JADE, SADE, EPSDE, GSA, PSO2011, and CMA-ES). The statistical tests realized for the comparison of performances indicate that the problem-solving success of DS algorithm in transforming the geocentric cartesian coordinates into geodetic coordinates is higher than those of all classical methods and Computational-Intelligence algorithms used in this paper.
Differential geometry measures of nonlinearity for the bearing-only tracking problem
NASA Astrophysics Data System (ADS)
Mallick, Mahendra; La Scala, Barbara F.; Arulampalam, M. S.
2005-05-01
The bearing-only tracking problem arises in many radar and sonar tracking applications. Since the bearing measurement model is a nonlinear function of the target state, the filtering problem is nonlinear in nature. A great deal of attention has been focused on this problem due to the difficulty posed by the so-called high degree of nonlinearity (DoN) in the problem. However, a quantitative measure of the DoN is not calculated in previous works. It has been observed that the extended Kalman filter (EKF) in which the state vector consists of the Cartesian components of position and velocity is unstable and diverges in some cases. The range parametrized EKF (RPEKF) and particle filter (PF) have been shown to produce improved estimates for the bearing-only tracking problem. In this paper, we calculate two measures of nonlinearity, (1) the parameter-effects curvature and (2) intrinsic curvature for the bearing-only tracking problem using the differential geometry measures of nonlinearity. We present numerical results using simulated data for the constant velocity motion of a target in 2D with bearing-only measurements where the sensor platform uses a higher order motion than the target to achieve observability. We analyze the DoN by varying the distance between the target and sensor.
Perspectives for spintronics in 2D materials
NASA Astrophysics Data System (ADS)
Han, Wei
2016-03-01
The past decade has been especially creative for spintronics since the (re)discovery of various two dimensional (2D) materials. Due to the unusual physical characteristics, 2D materials have provided new platforms to probe the spin interaction with other degrees of freedom for electrons, as well as to be used for novel spintronics applications. This review briefly presents the most important recent and ongoing research for spintronics in 2D materials.
2D FEM Heat Transfer & E&M Field Code
Energy Science and Technology Software Center (ESTSC)
1992-04-02
TOPAZ and TOPAZ2D are two-dimensional implicit finite element computer codes for heat transfer analysis. TOPAZ2D can also be used to solve electrostatic and magnetostatic problems. The programs solve for the steady-state or transient temperature or electrostatic and magnetostatic potential field on two-dimensional planar or axisymmetric geometries. Material properties may be temperature or potential-dependent and either isotropic or orthotropic. A variety of time and temperature-dependent boundary conditions can be specified including temperature, flux, convection, and radiation.more » By implementing the user subroutine feature, users can model chemical reaction kinetics and allow for any type of functional representation of boundary conditions and internal heat generation. The programs can solve problems of diffuse and specular band radiation in an enclosure coupled with conduction in the material surrounding the enclosure. Additional features include thermal contact resistance across an interface, bulk fluids, phase change, and energy balances.« less
2-D and 3-D computations of curved accelerator magnets
Turner, L.R.
1991-01-01
In order to save computer memory, a long accelerator magnet may be computed by treating the long central region and the end regions separately. The dipole magnets for the injector synchrotron of the Advanced Photon Source (APS), now under construction at Argonne National Laboratory (ANL), employ magnet iron consisting of parallel laminations, stacked with a uniform radius of curvature of 33.379 m. Laplace's equation for the magnetic scalar potential has a different form for a straight magnet (x-y coordinates), a magnet with surfaces curved about a common center (r-{theta} coordinates), and a magnet with parallel laminations like the APS injector dipole. Yet pseudo 2-D computations for the three geometries give basically identical results, even for a much more strongly curved magnet. Hence 2-D (x-y) computations of the central region and 3-D computations of the end regions can be combined to determine the overall magnetic behavior of the magnets. 1 ref., 6 figs.
2D FEM Heat Transfer & E&M Field Code
1992-04-02
TOPAZ and TOPAZ2D are two-dimensional implicit finite element computer codes for heat transfer analysis. TOPAZ2D can also be used to solve electrostatic and magnetostatic problems. The programs solve for the steady-state or transient temperature or electrostatic and magnetostatic potential field on two-dimensional planar or axisymmetric geometries. Material properties may be temperature or potential-dependent and either isotropic or orthotropic. A variety of time and temperature-dependent boundary conditions can be specified including temperature, flux, convection, and radiation. By implementing the user subroutine feature, users can model chemical reaction kinetics and allow for any type of functional representation of boundary conditions and internal heat generation. The programs can solve problems of diffuse and specular band radiation in an enclosure coupled with conduction in the material surrounding the enclosure. Additional features include thermal contact resistance across an interface, bulk fluids, phase change, and energy balances.
2D resistivity inversion using conjugate gradients for a finite element discretization
NASA Astrophysics Data System (ADS)
Bortolozo, C. A.; Santos, F. M.; Porsani, J. L.
2014-12-01
In this work we present a DC 2D inversion algorithm using conjugate gradients relaxation to solve the maximum likelihood inverse equations. We apply, according to Zhang (1995), the maximum likelihood inverse theory developed by Tarantola and Valette (1982) to our 2D resistivity inversion. This algorithm was chosen to this research because it doesn't need to calculate the field's derivatives. Since conjugate gradient techniques only need the results of the sensitivity matrix Ã or its transpose ÃT multiplying a vector, the actual computation of the sensitivity matrix are not performed, according to the methodology described in Zhang (1995). In Zhang (1995), the terms Ãx and ÃTy, are dependent of the stiffness matrix K and its partial derivative ∂K⁄∂ρ. The inversion methodology described in Zhang (1995) is for the case of 3D electrical resistivity by finite differences discretization. So it was necessary to make a series of adjustments to obtain a satisfactory result for 2D electrical inversion using finite element method. The difference between the modeling of 3D resistivity with finite difference and the 2D finite element method are in the integration variable, used in the 2D case. In the 2D case the electrical potential are initially calculated in the transformed domain, including the stiffness matrix, and only in the end is transformed in Cartesian domain. In the case of 3D, described by Zhang (1995) this is done differently, the calculation is done directly in the Cartesian domain. In the literature was not found any work describing how to deal with this problem. Because the calculations of Ãx and ÃTy must be done without having the real stiffness matrix, the adaptation consist in calculate the stiffness matrix and its partial derivative using a set of integration variables. We transform those matrix in the same form has in the potential case, but with different sets of variables. The results will be presented and are very promising.
ERIC Educational Resources Information Center
Desseyn, H. O.; And Others
1985-01-01
Compares linear-nonlinear and planar-nonplanar geometry through the valence-shell electron pairs repulsion (V.S.E.P.R.), Mulliken-Walsh, and electrostatic force theories. Indicates that although the V.S.E.P.R. theory has more advantages for elementary courses, an explanation of the best features of the different theories offers students a better…
Finite Element Analysis of 2-D Elastic Contacts Involving FGMs
NASA Astrophysics Data System (ADS)
Abhilash, M. N.; Murthy, H.
2014-05-01
The response of elastic indenters in contact with Functionally Graded Material (FGM) coated homogeneous elastic half space has been presented in the current paper. Finite element analysis has been used due to its ability to handle complex geometry, material, and boundary conditions. Indenters of different typical surface profiles have been considered and the problem has been idealized as a two-dimensional (2D) plane strain problem considering only normal loads. Initially, indenters were considered to be rigid and the results were validated with the solutions presented in the literature. The analysis has then been extended to the case of elastic indenters on FGM-coated half spaces and the results are discussed.
VECTUM. Irregular 2D Velocity Vector Field Plotting Package
McClurg, F.R.; Mousseau, V.A.
1992-05-04
VECTUM is a NCAR Graphics based package, for generating a plot of an irregular 2D velocity vector field. The program reads an ASCII database of x, y, u, v, data pairs and produces a plot in Computer Graphics Metafile (CGM) format. The program also uses an ASCII parameter file for controlling annotation details such as the plot title, arrowhead style, scale of vectors, windowing, etc. Simple geometry (i.e. lines, arcs, splines) can be defined to be included with the velocity vectors. NCAR Graphics drivers can be used to display the CGM file into PostScript, HPGL, HDF, etc, output.
Annotated Bibliography of EDGE2D Use
J.D. Strachan and G. Corrigan
2005-06-24
This annotated bibliography is intended to help EDGE2D users, and particularly new users, find existing published literature that has used EDGE2D. Our idea is that a person can find existing studies which may relate to his intended use, as well as gain ideas about other possible applications by scanning the attached tables.
Staring 2-D hadamard transform spectral imager
Gentry, Stephen M.; Wehlburg, Christine M.; Wehlburg, Joseph C.; Smith, Mark W.; Smith, Jody L.
2006-02-07
A staring imaging system inputs a 2D spatial image containing multi-frequency spectral information. This image is encoded in one dimension of the image with a cyclic Hadamarid S-matrix. The resulting image is detecting with a spatial 2D detector; and a computer applies a Hadamard transform to recover the encoded image.
Analysis of Crustal Magnetisation in Cartesian Vector Harmonics
NASA Astrophysics Data System (ADS)
Gubbins, D.; Ivers, D.; Williams, S.
2015-12-01
We present a new set of functions, Vector Cartesian Harmonics (VCH), analogous to the Vector Spherical Harmonics that we have applied recently to global models of crustal and lithospheric magnetisation. Like their spherical counterpart, the VCH form a complete, orthogonal set: planar models of magnetisation can be expanded in them. There are 3 distinct types of VCH, one representing that part of the magnetisation which generates the potential magnetic field above the surface, another the potential magnetic field below the surface, and a toroidal function that generates only a non-potential field. One function therefore describes the magnetisation detected by observations of the magnetic anomaly while the other two describe the null space of an inversion of magnetic observations for magnetisation. The formalism is therefore ideal for analysing the results of inversions for magnetic structures in plane layers such as local or regional surveys where Earth's curvature can be ignored. The null space is in general very large, being an arbitrary combination of a doubly-infinite set of vector functions. However, in the absence of remanence and when the inducing field is uniform the null space reduces to only 2 types of structure, uniform susceptibility (Runcorn's Theorem) and a pattern of susceptibility induced by a uniform field, the null space is restricted to uniform magnetisation and 1D patterns of susceptibility aligned with a horizontal inducing field. Both these cases are already well known, but this analysis shows them to be the ONLY members of the null space. We also give results for familiar text-book structures to show the nature of the null space in each case. Curiously, inversion of the magnetic field from a buried dipole returns exactly half the correct magnitude plus a spurious distributed magnetisation. A more complex application is the topographic structure based on the Bishop formation in California (Fairhead and Williams, SEG exp. abstr. 25, 845, 2006
NASA Astrophysics Data System (ADS)
Cembranos, J. A. R.; Dobado, A.; Maroto, A. L.
Extra-dimensional theories contain additional degrees of freedom related to the geometry of the extra space which can be interpreted as new particles. Such theories allow to reformulate most of the fundamental problems of physics from a completely different point of view. In this essay, we concentrate on the brane fluctuations which are present in brane-worlds, and how such oscillations of the own space-time geometry along curved extra dimensions can help to resolve the Universe missing mass problem. The energy scales involved in these models are low compared to the Planck scale, and this means that some of the brane fluctuations distinctive signals could be detected in future colliders and in direct or indirect dark matter searches.
Light field morphing using 2D features.
Wang, Lifeng; Lin, Stephen; Lee, Seungyong; Guo, Baining; Shum, Heung-Yeung
2005-01-01
We present a 2D feature-based technique for morphing 3D objects represented by light fields. Existing light field morphing methods require the user to specify corresponding 3D feature elements to guide morph computation. Since slight errors in 3D specification can lead to significant morphing artifacts, we propose a scheme based on 2D feature elements that is less sensitive to imprecise marking of features. First, 2D features are specified by the user in a number of key views in the source and target light fields. Then the two light fields are warped view by view as guided by the corresponding 2D features. Finally, the two warped light fields are blended together to yield the desired light field morph. Two key issues in light field morphing are feature specification and warping of light field rays. For feature specification, we introduce a user interface for delineating 2D features in key views of a light field, which are automatically interpolated to other views. For ray warping, we describe a 2D technique that accounts for visibility changes and present a comparison to the ideal morphing of light fields. Light field morphing based on 2D features makes it simple to incorporate previous image morphing techniques such as nonuniform blending, as well as to morph between an image and a light field. PMID:15631126
2D materials for nanophotonic devices
NASA Astrophysics Data System (ADS)
Xu, Renjing; Yang, Jiong; Zhang, Shuang; Pei, Jiajie; Lu, Yuerui
2015-12-01
Two-dimensional (2D) materials have become very important building blocks for electronic, photonic, and phononic devices. The 2D material family has four key members, including the metallic graphene, transition metal dichalcogenide (TMD) layered semiconductors, semiconducting black phosphorous, and the insulating h-BN. Owing to the strong quantum confinements and defect-free surfaces, these atomically thin layers have offered us perfect platforms to investigate the interactions among photons, electrons and phonons. The unique interactions in these 2D materials are very important for both scientific research and application engineering. In this talk, I would like to briefly summarize and highlight the key findings, opportunities and challenges in this field. Next, I will introduce/highlight our recent achievements. We demonstrated atomically thin micro-lens and gratings using 2D MoS2, which is the thinnest optical component around the world. These devices are based on our discovery that the elastic light-matter interactions in highindex 2D materials is very strong. Also, I would like to introduce a new two-dimensional material phosphorene. Phosphorene has strongly anisotropic optical response, which creates 1D excitons in a 2D system. The strong confinement in phosphorene also enables the ultra-high trion (charged exciton) binding energies, which have been successfully measured in our experiments. Finally, I will briefly talk about the potential applications of 2D materials in energy harvesting.
Inertial solvation in femtosecond 2D spectra
NASA Astrophysics Data System (ADS)
Hybl, John; Albrecht Ferro, Allison; Farrow, Darcie; Jonas, David
2001-03-01
We have used 2D Fourier transform spectroscopy to investigate polar solvation. 2D spectroscopy can reveal molecular lineshapes beneath ensemble averaged spectra and freeze molecular motions to give an undistorted picture of the microscopic dynamics of polar solvation. The transition from "inhomogeneous" to "homogeneous" 2D spectra is governed by both vibrational relaxation and solvent motion. Therefore, the time dependence of the 2D spectrum directly reflects the total response of the solvent-solute system. IR144, a cyanine dye with a dipole moment change upon electronic excitation, was used to probe inertial solvation in methanol and propylene carbonate. Since the static Stokes' shift of IR144 in each of these solvents is similar, differences in the 2D spectra result from solvation dynamics. Initial results indicate that the larger propylene carbonate responds more slowly than methanol, but appear to be inconsistent with rotational estimates of the inertial response. To disentangle intra-molecular vibrations from solvent motion, the 2D spectra of IR144 will be compared to the time-dependent 2D spectra of the structurally related nonpolar cyanine dye HDITCP.
Internal Photoemission Spectroscopy of 2-D Materials
NASA Astrophysics Data System (ADS)
Nguyen, Nhan; Li, Mingda; Vishwanath, Suresh; Yan, Rusen; Xiao, Shudong; Xing, Huili; Cheng, Guangjun; Hight Walker, Angela; Zhang, Qin
Recent research has shown the great benefits of using 2-D materials in the tunnel field-effect transistor (TFET), which is considered a promising candidate for the beyond-CMOS technology. The on-state current of TFET can be enhanced by engineering the band alignment of different 2D-2D or 2D-3D heterostructures. Here we present the internal photoemission spectroscopy (IPE) approach to determine the band alignments of various 2-D materials, in particular SnSe2 and WSe2, which have been proposed for new TFET designs. The metal-oxide-2-D semiconductor test structures are fabricated and characterized by IPE, where the band offsets from the 2-D semiconductor to the oxide conduction band minimum are determined by the threshold of the cube root of IPE yields as a function of photon energy. In particular, we find that SnSe2 has a larger electron affinity than most semiconductors and can be combined with other semiconductors to form near broken-gap heterojunctions with low barrier heights which can produce a higher on-state current. The details of data analysis of IPE and the results from Raman spectroscopy and spectroscopic ellipsometry measurements will also be presented and discussed.
Geometry Dependence of Stellarator Turbulence
H.E. Mynick, P. Xanthopoulos and A.H. Boozer
2009-08-10
Using the nonlinear gyrokinetic code package GENE/GIST, we study the turbulent transport in a broad family of stellarator designs, to understand the geometry-dependence of the microturbulence. By using a set of flux tubes on a given flux surface, we construct a picture of the 2D structure of the microturbulence over that surface, and relate this to relevant geometric quantities, such as the curvature, local shear, and effective potential in the Schrodinger-like equation governing linear drift modes.
Rabow, A. A.; Scheraga, H. A.
1996-01-01
We have devised a Cartesian combination operator and coding scheme for improving the performance of genetic algorithms applied to the protein folding problem. The genetic coding consists of the C alpha Cartesian coordinates of the protein chain. The recombination of the genes of the parents is accomplished by: (1) a rigid superposition of one parent chain on the other, to make the relation of Cartesian coordinates meaningful, then, (2) the chains of the children are formed through a linear combination of the coordinates of their parents. The children produced with this Cartesian combination operator scheme have similar topology and retain the long-range contacts of their parents. The new scheme is significantly more efficient than the standard genetic algorithm methods for locating low-energy conformations of proteins. The considerable superiority of genetic algorithms over Monte Carlo optimization methods is also demonstrated. We have also devised a new dynamic programming lattice fitting procedure for use with the Cartesian combination operator method. The procedure finds excellent fits of real-space chains to the lattice while satisfying bond-length, bond-angle, and overlap constraints. PMID:8880904
Freitas, Andreia C; Wylezinska, Marzena; Birch, Malcolm J; Petersen, Steffen E; Miquel, Marc E
2016-01-01
Dynamic imaging of the vocal tract using real-time MRI has been an active and growing area of research, having demonstrated great potential to become routinely performed in the clinical evaluation of speech and swallowing disorders. Although many technical advances have been made in regards to acquisition and reconstruction methodologies, there is still no consensus in best practice protocols. This study aims to compare Cartesian and non-Cartesian real-time MRI sequences, regarding image quality and temporal resolution trade-off, for dynamic speech imaging. Five subjects were imaged at 1.5T, while performing normal phonation, in order to assess velar motion and velopharyngeal closure. Data was acquired using both Cartesian and non-Cartesian (spiral and radial) real-time sequences at five different spatial-temporal resolution sets, between 10 fps (1.7×1.7×10 mm3) and 25 fps (1.5×1.5×10 mm3). Only standard scanning resources provided by the MRI scanner manufacturer were used to ensure easy applicability to clinical evaluation and reproducibility. Data sets were evaluated by comparing measurements of the velar structure, dynamic contrast-to-noise ratio and image quality visual scoring. Results showed that for all proposed sequences, FLASH spiral acquisitions provided higher contrast-to-noise ratio, up to a 170.34% increase at 20 fps, than equivalent bSSFP Cartesian acquisitions for the same spatial-temporal resolution. At higher frame rates (22 and 25 fps), spiral protocols were optimal and provided higher CNR and visual scoring than equivalent radial protocols. Comparison of dynamic imaging at 10 and 22 fps for radial and spiral acquisitions revealed no significant difference in CNR performance, thus indicating that temporal resolution can be doubled without compromising spatial resolution (1.9×1.9 mm2) or CNR. In summary, this study suggests that the use of FLASH spiral protocols should be preferred over bSSFP Cartesian for the dynamic imaging of velopharyngeal
Freitas, Andreia C.; Wylezinska, Marzena; Birch, Malcolm J.; Petersen, Steffen E.; Miquel, Marc E.
2016-01-01
Dynamic imaging of the vocal tract using real-time MRI has been an active and growing area of research, having demonstrated great potential to become routinely performed in the clinical evaluation of speech and swallowing disorders. Although many technical advances have been made in regards to acquisition and reconstruction methodologies, there is still no consensus in best practice protocols. This study aims to compare Cartesian and non-Cartesian real-time MRI sequences, regarding image quality and temporal resolution trade-off, for dynamic speech imaging. Five subjects were imaged at 1.5T, while performing normal phonation, in order to assess velar motion and velopharyngeal closure. Data was acquired using both Cartesian and non-Cartesian (spiral and radial) real-time sequences at five different spatial-temporal resolution sets, between 10 fps (1.7×1.7×10 mm3) and 25 fps (1.5×1.5×10 mm3). Only standard scanning resources provided by the MRI scanner manufacturer were used to ensure easy applicability to clinical evaluation and reproducibility. Data sets were evaluated by comparing measurements of the velar structure, dynamic contrast-to-noise ratio and image quality visual scoring. Results showed that for all proposed sequences, FLASH spiral acquisitions provided higher contrast-to-noise ratio, up to a 170.34% increase at 20 fps, than equivalent bSSFP Cartesian acquisitions for the same spatial-temporal resolution. At higher frame rates (22 and 25 fps), spiral protocols were optimal and provided higher CNR and visual scoring than equivalent radial protocols. Comparison of dynamic imaging at 10 and 22 fps for radial and spiral acquisitions revealed no significant difference in CNR performance, thus indicating that temporal resolution can be doubled without compromising spatial resolution (1.9×1.9 mm2) or CNR. In summary, this study suggests that the use of FLASH spiral protocols should be preferred over bSSFP Cartesian for the dynamic imaging of velopharyngeal
Design Application Translates 2-D Graphics to 3-D Surfaces
NASA Technical Reports Server (NTRS)
2007-01-01
Fabric Images Inc., specializing in the printing and manufacturing of fabric tension architecture for the retail, museum, and exhibit/tradeshow communities, designed software to translate 2-D graphics for 3-D surfaces prior to print production. Fabric Images' fabric-flattening design process models a 3-D surface based on computer-aided design (CAD) specifications. The surface geometry of the model is used to form a 2-D template, similar to a flattening process developed by NASA's Glenn Research Center. This template or pattern is then applied in the development of a 2-D graphic layout. Benefits of this process include 11.5 percent time savings per project, less material wasted, and the ability to improve upon graphic techniques and offer new design services. Partners include Exhibitgroup/Giltspur (end-user client: TAC Air, a division of Truman Arnold Companies Inc.), Jack Morton Worldwide (end-user client: Nickelodeon), as well as 3D Exhibits Inc., and MG Design Associates Corp.
Brittle damage models in DYNA2D
Faux, D.R.
1997-09-01
DYNA2D is an explicit Lagrangian finite element code used to model dynamic events where stress wave interactions influence the overall response of the system. DYNA2D is often used to model penetration problems involving ductile-to-ductile impacts; however, with the advent of the use of ceramics in the armor-anti-armor community and the need to model damage to laser optics components, good brittle damage models are now needed in DYNA2D. This report will detail the implementation of four brittle damage models in DYNA2D, three scalar damage models and one tensor damage model. These new brittle damage models are then used to predict experimental results from three distinctly different glass damage problems.
Ginsparg, P.
1991-01-01
These are introductory lectures for a general audience that give an overview of the subject of matrix models and their application to random surfaces, 2d gravity, and string theory. They are intentionally 1.5 years out of date.
Ginsparg, P.
1991-12-31
These are introductory lectures for a general audience that give an overview of the subject of matrix models and their application to random surfaces, 2d gravity, and string theory. They are intentionally 1.5 years out of date.
2D electronic materials for army applications
NASA Astrophysics Data System (ADS)
O'Regan, Terrance; Perconti, Philip
2015-05-01
The record electronic properties achieved in monolayer graphene and related 2D materials such as molybdenum disulfide and hexagonal boron nitride show promise for revolutionary high-speed and low-power electronic devices. Heterogeneous 2D-stacked materials may create enabling technology for future communication and computation applications to meet soldier requirements. For instance, transparent, flexible and even wearable systems may become feasible. With soldier and squad level electronic power demands increasing, the Army is committed to developing and harnessing graphene-like 2D materials for compact low size-weight-and-power-cost (SWAP-C) systems. This paper will review developments in 2D electronic materials at the Army Research Laboratory over the last five years and discuss directions for future army applications.
2-d Finite Element Code Postprocessor
Energy Science and Technology Software Center (ESTSC)
1996-07-15
ORION is an interactive program that serves as a postprocessor for the analysis programs NIKE2D, DYNA2D, TOPAZ2D, and CHEMICAL TOPAZ2D. ORION reads binary plot files generated by the two-dimensional finite element codes currently used by the Methods Development Group at LLNL. Contour and color fringe plots of a large number of quantities may be displayed on meshes consisting of triangular and quadrilateral elements. ORION can compute strain measures, interface pressures along slide lines, reaction forcesmore » along constrained boundaries, and momentum. ORION has been applied to study the response of two-dimensional solids and structures undergoing finite deformations under a wide variety of large deformation transient dynamic and static problems and heat transfer analyses.« less
Bottom-up design of 2D organic photocatalysts for visible-light driven hydrogen evolution.
Wang, Peng; Jiang, Xue; Zhao, Jijun
2016-01-27
To design two-dimensional (2D) organocatalysts, three series of covalent organic frameworks (COFs) are constructed using bottom-up strategies, i.e. molecular selection, tunable linkage, and functionalization. First-principles calculations are performed to confirm their photocatalytic activity under visible light. Two of our constructed 2D COF models (B1 and C3) are identified as a sufficiently efficient organocatalyst for visible light water splitting. The controllable construction of such COFs from suitable organic subunit, linkage, and functional groups paves the way for correlating band edge alignments and geometry parameters of 2D organic materials. Our theoretical prediction not only provides essential insights into designing 2D-COF photocatalysts for water splitting, but also sparks other technological applications for 2D organic materials. PMID:26704386
Bottom-up design of 2D organic photocatalysts for visible-light driven hydrogen evolution
NASA Astrophysics Data System (ADS)
Wang, Peng; Jiang, Xue; Zhao, Jijun
2016-01-01
To design two-dimensional (2D) organocatalysts, three series of covalent organic frameworks (COFs) are constructed using bottom-up strategies, i.e. molecular selection, tunable linkage, and functionalization. First-principles calculations are performed to confirm their photocatalytic activity under visible light. Two of our constructed 2D COF models (B1 and C3) are identified as a sufficiently efficient organocatalyst for visible light water splitting. The controllable construction of such COFs from suitable organic subunit, linkage, and functional groups paves the way for correlating band edge alignments and geometry parameters of 2D organic materials. Our theoretical prediction not only provides essential insights into designing 2D-COF photocatalysts for water splitting, but also sparks other technological applications for 2D organic materials.
Chemical Approaches to 2D Materials.
Samorì, Paolo; Palermo, Vincenzo; Feng, Xinliang
2016-08-01
Chemistry plays an ever-increasing role in the production, functionalization, processing and applications of graphene and other 2D materials. This special issue highlights a selection of enlightening chemical approaches to 2D materials, which nicely reflect the breadth of the field and convey the excitement of the individuals involved in it, who are trying to translate graphene and related materials from the laboratory into a real, high-impact technology. PMID:27478083
TOPAZ2D heat transfer code users manual and thermal property data base
NASA Astrophysics Data System (ADS)
Shapiro, A. B.; Edwards, A. L.
1990-05-01
TOPAZ2D is a two dimensional implicit finite element computer code for heat transfer analysis. This user's manual provides information on the structure of a TOPAZ2D input file. Also included is a material thermal property data base. This manual is supplemented with The TOPAZ2D Theoretical Manual and the TOPAZ2D Verification Manual. TOPAZ2D has been implemented on the CRAY, SUN, and VAX computers. TOPAZ2D can be used to solve for the steady state or transient temperature field on two dimensional planar or axisymmetric geometries. Material properties may be temperature dependent and either isotropic or orthotropic. A variety of time and temperature dependent boundary conditions can be specified including temperature, flux, convection, and radiation. Time or temperature dependent internal heat generation can be defined locally be element or globally by material. TOPAZ2D can solve problems of diffuse and specular band radiation in an enclosure coupled with conduction in material surrounding the enclosure. Additional features include thermally controlled reactive chemical mixtures, thermal contact resistance across an interface, bulk fluid flow, phase change, and energy balances. Thermal stresses can be calculated using the solid mechanics code NIKE2D which reads the temperature state data calculated by TOPAZ2D. A three dimensional version of the code, TOPAZ3D is available.
Extended 2D generalized dilaton gravity theories
NASA Astrophysics Data System (ADS)
de Mello, R. O.
2008-09-01
We show that an anomaly-free description of matter in (1+1) dimensions requires a deformation of the 2D relativity principle, which introduces a non-trivial centre in the 2D Poincaré algebra. Then we work out the reduced phase space of the anomaly-free 2D relativistic particle, in order to show that it lives in a noncommutative 2D Minkowski space. Moreover, we build a Gaussian wave packet to show that a Planck length is well defined in two dimensions. In order to provide a gravitational interpretation for this noncommutativity, we propose to extend the usual 2D generalized dilaton gravity models by a specific Maxwell component, which guages the extra symmetry associated with the centre of the 2D Poincaré algebra. In addition, we show that this extension is a high energy correction to the unextended dilaton theories that can affect the topology of spacetime. Further, we couple a test particle to the general extended dilaton models with the purpose of showing that they predict a noncommutativity in curved spacetime, which is locally described by a Moyal star product in the low energy limit. We also conjecture a probable generalization of this result, which provides strong evidence that the noncommutativity is described by a certain star product which is not of the Moyal type at high energies. Finally, we prove that the extended dilaton theories can be formulated as Poisson Sigma models based on a nonlinear deformation of the extended Poincaré algebra.
High Energy Boundary Conditions for a Cartesian Mesh Euler Solver
NASA Technical Reports Server (NTRS)
Pandya, Shishir A.; Murman, Scott M.; Aftosmis, Michael J.
2004-01-01
Inlets and exhaust nozzles are often omitted or fared over in aerodynamic simulations of aircraft due to the complexities involving in the modeling of engine details such as complex geometry and flow physics. However, the assumption is often improper as inlet or plume flows have a substantial effect on vehicle aerodynamics. A tool for specifying inlet and exhaust plume conditions through the use of high-energy boundary conditions in an established inviscid flow solver is presented. The effects of the plume on the flow fields near the inlet and plume are discussed.
NASA Astrophysics Data System (ADS)
Rousseau, Guy; Gay, David; Piché, Michel
2006-09-01
A recent analysis [G. Rousseau, D. Gay and M. Piché, One-dimensional description of cylindrically symmetric laser beams: application to Bessel-type nondiffracting beams, J. Opt. Soc. Am. A, 22 (2005) 1274] has shown that any cylindrically symmetric laser beam can be synthesized from a single wave called a constituent wave. This representation allows the introduction of one-dimensional Cartesian root-mean-square (rms) parameters to describe the conical structure of cylindrically symmetric laser beams. In this paper, we compare the rms characterization of Bessel-Gauss beams in polar coordinates with that of their respective constituent waves in Cartesian coordinates. Numerical results reveal an asymptotic correspondence between polar and Cartesian rms parameters of Bessel-Gauss beams propagating in a nondiffracting regime. Such a correspondence eliminates misleading interpretations about the propagation factor and the Rayleigh range of nondiffracting Bessel-type beams characterized in terms of polar rms parameters.
An accuracy assessment of Cartesian-mesh approaches for the Euler equations
NASA Technical Reports Server (NTRS)
Coirier, William J.; Powell, Kenneth G.
1995-01-01
A critical assessment of the accuracy of Cartesian-mesh approaches for steady, transonic solutions of the Euler equations of gas dynamics is made. An exact solution of the Euler equations (Ringleb's flow) is used not only to infer the order of the truncation error of the Cartesian-mesh approaches, but also to compare the magnitude of the discrete error directly to that obtained with a structured mesh approach. Uniformly and adaptively refined solutions using a Cartesian-mesh approach are obtained and compared to each other and to uniformly refined structured mesh results. The effect of cell merging is investigated as well as the use of two different K-exact reconstruction procedures. The solution methodology of the schemes is explained and tabulated results are presented to compare the solution accuracies.
Smith, C U
2001-08-01
J. C. Eccles (1903-1997) had a highly distinguished career in neurophysiology, being awarded the Nobel Prize for Medicine or Physiology in 1963. This paper sets him within the Cartesian tradition of British neurophysiology initiated by Thomas Henry Huxley in the mid-19th century. It shows how the mind-brain problematique of the Cartesian tradition troubled him throughout his career, leading him finally to a solution in terms of quantum microphysics and microphysiology. This position, which has subsequently become fashionable, is discussed and shown (at least in the form Eccles espoused) to provide no solution to the problem posed by Descartes in the early 17th century. PMID:11487286
Random subspaces for encryption based on a private shared Cartesian frame
Bartlett, Stephen D.; Hayden, Patrick; Spekkens, Robert W.
2005-11-15
A private shared Cartesian frame is a novel form of private shared correlation that allows for both private classical and quantum communication. Cryptography using a private shared Cartesian frame has the remarkable property that asymptotically, if perfect privacy is demanded, the private classical capacity is three times the private quantum capacity. We demonstrate that if the requirement for perfect privacy is relaxed, then it is possible to use the properties of random subspaces to nearly triple the private quantum capacity, almost closing the gap between the private classical and quantum capacities.
A Cartesian cut cell method for rarefied flow simulations around moving obstacles
NASA Astrophysics Data System (ADS)
Dechristé, G.; Mieussens, L.
2016-06-01
For accurate simulations of rarefied gas flows around moving obstacles, we propose a cut cell method on Cartesian grids: it allows exact conservation and accurate treatment of boundary conditions. Our approach is designed to treat Cartesian cells and various kinds of cut cells by the same algorithm, with no need to identify the specific shape of each cut cell. This makes the implementation quite simple, and allows a direct extension to 3D problems. Such simulations are also made possible by using an adaptive mesh refinement technique and a hybrid parallel implementation. This is illustrated by several test cases, including a 3D unsteady simulation of the Crookes radiometer.
NASA Astrophysics Data System (ADS)
Rhee, Young Min
2000-10-01
A modified method to construct an accurate potential energy surface by interpolation is presented. The modification is based on the use of Cartesian coordinates in the weighting function. The translational and rotational invariance of the potential is incorporated by a proper definition of the distance between two Cartesian configurations. A numerical algorithm to find the distance is developed. It is shown that the present method is more exact in describing a planar system compared to the previous methods with weightings in internal coordinates. The applicability of the method to reactive systems is also demonstrated by performing classical trajectory simulations on the surface.
Parameter Studies, time-dependent simulations and design with automated Cartesian methods
NASA Technical Reports Server (NTRS)
Aftosmis, Michael
2005-01-01
Over the past decade, NASA has made a substantial investment in developing adaptive Cartesian grid methods for aerodynamic simulation. Cartesian-based methods played a key role in both the Space Shuttle Accident Investigation and in NASA's return to flight activities. The talk will provide an overview of recent technological developments focusing on the generation of large-scale aerodynamic databases, automated CAD-based design, and time-dependent simulations with of bodies in relative motion. Automation, scalability and robustness underly all of these applications and research in each of these topics will be presented.
A Fast and Robust Poisson-Boltzmann Solver Based on Adaptive Cartesian Grids.
Boschitsch, Alexander H; Fenley, Marcia O
2011-05-10
An adaptive Cartesian grid (ACG) concept is presented for the fast and robust numerical solution of the 3D Poisson-Boltzmann Equation (PBE) governing the electrostatic interactions of large-scale biomolecules and highly charged multi-biomolecular assemblies such as ribosomes and viruses. The ACG offers numerous advantages over competing grid topologies such as regular 3D lattices and unstructured grids. For very large biological molecules and multi-biomolecule assemblies, the total number of grid-points is several orders of magnitude less than that required in a conventional lattice grid used in the current PBE solvers thus allowing the end user to obtain accurate and stable nonlinear PBE solutions on a desktop computer. Compared to tetrahedral-based unstructured grids, ACG offers a simpler hierarchical grid structure, which is naturally suited to multigrid, relieves indirect addressing requirements and uses fewer neighboring nodes in the finite difference stencils. Construction of the ACG and determination of the dielectric/ionic maps are straightforward, fast and require minimal user intervention. Charge singularities are eliminated by reformulating the problem to produce the reaction field potential in the molecular interior and the total electrostatic potential in the exterior ionic solvent region. This approach minimizes grid-dependency and alleviates the need for fine grid spacing near atomic charge sites. The technical portion of this paper contains three parts. First, the ACG and its construction for general biomolecular geometries are described. Next, a discrete approximation to the PBE upon this mesh is derived. Finally, the overall solution procedure and multigrid implementation are summarized. Results obtained with the ACG-based PBE solver are presented for: (i) a low dielectric spherical cavity, containing interior point charges, embedded in a high dielectric ionic solvent - analytical solutions are available for this case, thus allowing rigorous
Optical modulators with 2D layered materials
NASA Astrophysics Data System (ADS)
Sun, Zhipei; Martinez, Amos; Wang, Feng
2016-04-01
Light modulation is an essential operation in photonics and optoelectronics. With existing and emerging technologies increasingly demanding compact, efficient, fast and broadband optical modulators, high-performance light modulation solutions are becoming indispensable. The recent realization that 2D layered materials could modulate light with superior performance has prompted intense research and significant advances, paving the way for realistic applications. In this Review, we cover the state of the art of optical modulators based on 2D materials, including graphene, transition metal dichalcogenides and black phosphorus. We discuss recent advances employing hybrid structures, such as 2D heterostructures, plasmonic structures, and silicon and fibre integrated structures. We also take a look at the future perspectives and discuss the potential of yet relatively unexplored mechanisms, such as magneto-optic and acousto-optic modulation.
Large Area Synthesis of 2D Materials
NASA Astrophysics Data System (ADS)
Vogel, Eric
Transition metal dichalcogenides (TMDs) have generated significant interest for numerous applications including sensors, flexible electronics, heterostructures and optoelectronics due to their interesting, thickness-dependent properties. Despite recent progress, the synthesis of high-quality and highly uniform TMDs on a large scale is still a challenge. In this talk, synthesis routes for WSe2 and MoS2 that achieve monolayer thickness uniformity across large area substrates with electrical properties equivalent to geological crystals will be described. Controlled doping of 2D semiconductors is also critically required. However, methods established for conventional semiconductors, such as ion implantation, are not easily applicable to 2D materials because of their atomically thin structure. Redox-active molecular dopants will be demonstrated which provide large changes in carrier density and workfunction through the choice of dopant, treatment time, and the solution concentration. Finally, several applications of these large-area, uniform 2D materials will be described including heterostructures, biosensors and strain sensors.
2D microwave imaging reflectometer electronics
Spear, A. G.; Domier, C. W. Hu, X.; Muscatello, C. M.; Ren, X.; Luhmann, N. C.; Tobias, B. J.
2014-11-15
A 2D microwave imaging reflectometer system has been developed to visualize electron density fluctuations on the DIII-D tokamak. Simultaneously illuminated at four probe frequencies, large aperture optics image reflections from four density-dependent cutoff surfaces in the plasma over an extended region of the DIII-D plasma. Localized density fluctuations in the vicinity of the plasma cutoff surfaces modulate the plasma reflections, yielding a 2D image of electron density fluctuations. Details are presented of the receiver down conversion electronics that generate the in-phase (I) and quadrature (Q) reflectometer signals from which 2D density fluctuation data are obtained. Also presented are details on the control system and backplane used to manage the electronics as well as an introduction to the computer based control program.
2D microwave imaging reflectometer electronics
NASA Astrophysics Data System (ADS)
Spear, A. G.; Domier, C. W.; Hu, X.; Muscatello, C. M.; Ren, X.; Tobias, B. J.; Luhmann, N. C.
2014-11-01
A 2D microwave imaging reflectometer system has been developed to visualize electron density fluctuations on the DIII-D tokamak. Simultaneously illuminated at four probe frequencies, large aperture optics image reflections from four density-dependent cutoff surfaces in the plasma over an extended region of the DIII-D plasma. Localized density fluctuations in the vicinity of the plasma cutoff surfaces modulate the plasma reflections, yielding a 2D image of electron density fluctuations. Details are presented of the receiver down conversion electronics that generate the in-phase (I) and quadrature (Q) reflectometer signals from which 2D density fluctuation data are obtained. Also presented are details on the control system and backplane used to manage the electronics as well as an introduction to the computer based control program.
2D microwave imaging reflectometer electronics.
Spear, A G; Domier, C W; Hu, X; Muscatello, C M; Ren, X; Tobias, B J; Luhmann, N C
2014-11-01
A 2D microwave imaging reflectometer system has been developed to visualize electron density fluctuations on the DIII-D tokamak. Simultaneously illuminated at four probe frequencies, large aperture optics image reflections from four density-dependent cutoff surfaces in the plasma over an extended region of the DIII-D plasma. Localized density fluctuations in the vicinity of the plasma cutoff surfaces modulate the plasma reflections, yielding a 2D image of electron density fluctuations. Details are presented of the receiver down conversion electronics that generate the in-phase (I) and quadrature (Q) reflectometer signals from which 2D density fluctuation data are obtained. Also presented are details on the control system and backplane used to manage the electronics as well as an introduction to the computer based control program. PMID:25430247
2D-Crystal-Based Functional Inks.
Bonaccorso, Francesco; Bartolotta, Antonino; Coleman, Jonathan N; Backes, Claudia
2016-08-01
The possibility to produce and process graphene, related 2D crystals, and heterostructures in the liquid phase makes them promising materials for an ever-growing class of applications as composite materials, sensors, in flexible optoelectronics, and energy storage and conversion. In particular, the ability to formulate functional inks with on-demand rheological and morphological properties, i.e., lateral size and thickness of the dispersed 2D crystals, is a step forward toward the development of industrial-scale, reliable, inexpensive printing/coating processes, a boost for the full exploitation of such nanomaterials. Here, the exfoliation strategies of graphite and other layered crystals are reviewed, along with the advances in the sorting of lateral size and thickness of the exfoliated sheets together with the formulation of functional inks and the current development of printing/coating processes of interest for the realization of 2D-crystal-based devices. PMID:27273554
Toward a 2-D magneto-optical trap for polar molecules
NASA Astrophysics Data System (ADS)
Hummon, Matthew; Stuhl, Benjamin; Yeo, Mark; Collopy, Alejandra; Ye, Jun
2012-06-01
The additional structure that arises from the rotational degree of freedom in diatomic molecules makes difficult the adaptation of a traditional atomic magneto-optical trap (MOT) for use with molecules. We describe progress toward development of a 2-D MOT for laser cooled yttrium monoxide molecules based on a resonant LC baseball coil geometry.
Parallel algorithms for 2-D cylindrical transport equations of Eigenvalue problem
Wei, J.; Yang, S.
2013-07-01
In this paper, aimed at the neutron transport equations of eigenvalue problem under 2-D cylindrical geometry on unstructured grid, the discrete scheme of Sn discrete ordinate and discontinuous finite is built, and the parallel computation for the scheme is realized on MPI systems. Numerical experiments indicate that the designed parallel algorithm can reach perfect speedup, it has good practicality and scalability. (authors)
2D and 3D Method of Characteristic Tools for Complex Nozzle Development
NASA Technical Reports Server (NTRS)
Rice, Tharen
2003-01-01
This report details the development of a 2D and 3D Method of Characteristic (MOC) tool for the design of complex nozzle geometries. These tools are GUI driven and can be run on most Windows-based platforms. The report provides a user's manual for these tools as well as explains the mathematical algorithms used in the MOC solutions.
The Use of Geometric Properties of 2D Arrays across Development
ERIC Educational Resources Information Center
Gibson, Brett M.; Leichtman, Michelle D.; Costa, Rachel; Bemis, Rhyannon
2009-01-01
Four- to 10-year-old children (n = 50) participated in a 2D search task that included geometry (with- and without lines) and feature conditions. During each of 27 trials, participants watched as a cartoon character hid behind one of three landmarks arranged in a triangle on a computer screen. During feature condition trials, participants could use…
The 2D lingual appliance system.
Cacciafesta, Vittorio
2013-09-01
The two-dimensional (2D) lingual bracket system represents a valuable treatment option for adult patients seeking a completely invisible orthodontic appliance. The ease of direct or simplified indirect bonding of 2D lingual brackets in combination with low friction mechanics makes it possible to achieve a good functional and aesthetic occlusion, even in the presence of a severe malocclusion. The use of a self-ligating bracket significantly reduces chair-side time for the orthodontist, and the low-profile bracket design greatly improves patient comfort. PMID:24005953
Inkjet printing of 2D layered materials.
Li, Jiantong; Lemme, Max C; Östling, Mikael
2014-11-10
Inkjet printing of 2D layered materials, such as graphene and MoS2, has attracted great interests for emerging electronics. However, incompatible rheology, low concentration, severe aggregation and toxicity of solvents constitute critical challenges which hamper the manufacturing efficiency and product quality. Here, we introduce a simple and general technology concept (distillation-assisted solvent exchange) to efficiently overcome these challenges. By implementing the concept, we have demonstrated excellent jetting performance, ideal printing patterns and a variety of promising applications for inkjet printing of 2D layered materials. PMID:25169938
Measurement of 2D birefringence distribution
NASA Astrophysics Data System (ADS)
Noguchi, Masato; Ishikawa, Tsuyoshi; Ohno, Masahiro; Tachihara, Satoru
1992-10-01
A new measuring method of 2-D birefringence distribution has been developed. It has not been an easy job to get a birefringence distribution in an optical element with conventional ellipsometry because of its lack of scanning means. Finding an analogy between the rotating analyzer method in ellipsometry and the phase-shifting method in recently developed digital interferometry, we have applied the phase-shifting algorithm to ellipsometry, and have developed a new method that makes the measurement of 2-D birefringence distribution easy and possible. The system contains few moving parts, assuring reliability, and measures a large area of a sample at one time, making the measuring time very short.
Topological Invariants and CW Complexes of Cartesian Product and Hexagonal Tiling Paces
NASA Astrophysics Data System (ADS)
Escudero, Juan García
2011-09-01
The cohomology of a class of cartesian product tiling spaces in N dimensions when the inflation factor is a Pisot-Vijayaraghavan unit is analyzed. A CW complex for an hexagonal tiling space is defined in terms of collared tiles for the study of its topological invariants.
ERIC Educational Resources Information Center
Earnest, Darrell Steven
2012-01-01
This dissertation explores fifth and eighth grade students' interpretations of three kinds of mathematical representations: number lines, the Cartesian plane, and graphs of linear functions. Two studies were conducted. In Study 1, I administered the paper-and-pencil Linear Representations Assessment (LRA) to examine students'…
Real-time cartesian force feedback control of a teleoperated robot
NASA Technical Reports Server (NTRS)
Campbell, Perry
1989-01-01
Active cartesian force control of a teleoperated robot is investigated. An economical microcomputer based control method was tested. Limitations are discussed and methods of performance improvement suggested. To demonstrate the performance of this technique, a preliminary test was performed with success. A general purpose bilateral force reflecting hand controller is currently being constructed based on this control method.
"Mens Sana in Corpore Sano": Cartesian Dualism and the Marginalisation of Sex Education
ERIC Educational Resources Information Center
Paechter, Carrie
2004-01-01
Cartesian dualism has left a heavy legacy in terms of how we think about ourselves, so that we treat humans as minds within bodies rather than mind/body unities. This has far-reaching effects on our conceptualisation of the sex/gender distinction and on the relationship between bodies and identities. Related to this is a dualism that is embedded…
NASA Astrophysics Data System (ADS)
Le Chenadec, Vincent; Bay, Yong Yi
2015-11-01
The treatment of complex geometries in Computational Fluid Dynamics applications is a challenging endeavor, which immersed boundary and cut-cell techniques can significantly simplify by alleviating the meshing process required by body-fitted meshes. These methods also introduce new challenges, in that the formulation of accurate and well-posed discrete operators is not trivial. A cut-cell method for the solution of the incompressible Navier-Stokes equation is proposed for staggered Cartesian grids. In both scalar and vector cases, the emphasis is set on the structure of the discrete operators, designed to mimic the properties of the continuous ones while retaining a nearest-neighbor stencil. For convective transport, different forms are proposed (divergence, advective and skew-symmetric), and shown to be equivalent when the discrete continuity equation is satisfied. This ensures mass, momentum and kinetic energy conservation. For diffusive transport, conservative and symmetric operators are proposed for both Dirichlet and Neumann boundary conditions. Symmetry ensures the existence of a sink term (viscous dissipation) in the discrete kinetic energy budget, which is beneficial for stability. The accuracy of method is finally assessed in standard test cases.
NASA Astrophysics Data System (ADS)
Trost, Nico; Jiménez, Javier; Imke, Uwe; Sanchez, Victor
2014-06-01
TWOPORFLOW is a thermo-hydraulic code based on a porous media approach to simulate single- and two-phase flow including boiling. It is under development at the Institute for Neutron Physics and Reactor Technology (INR) at KIT. The code features a 3D transient solution of the mass, momentum and energy conservation equations for two inter-penetrating fluids with a semi-implicit continuous Eulerian type solver. The application domain of TWOPORFLOW includes the flow in standard porous media and in structured porous media such as micro-channels and cores of nuclear power plants. In the latter case, the fluid domain is coupled to a fuel rod model, describing the heat flow inside the solid structure. In this work, detailed profiling tools have been utilized to determine the optimization potential of TWOPORFLOW. As a result, bottle-necks were identified and reduced in the most feasible way, leading for instance to an optimization of the water-steam property computation. Furthermore, an OpenMP implementation addressing the routines in charge of inter-phase momentum-, energy- and mass-coupling delivered good performance together with a high scalability on shared memory architectures. In contrast to that, the approach for distributed memory systems was to solve sub-problems resulting by the decomposition of the initial Cartesian geometry. Thread communication for the sub-problem boundary updates was accomplished by the Message Passing Interface (MPI) standard.
NASA Technical Reports Server (NTRS)
Aftosmis, M. J.; Berger, M. J.; Murman, S. M.; Kwak, Dochan (Technical Monitor)
2002-01-01
The proposed paper will present recent extensions in the development of an efficient Euler solver for adaptively-refined Cartesian meshes with embedded boundaries. The paper will focus on extensions of the basic method to include solution adaptation, time-dependent flow simulation, and arbitrary rigid domain motion. The parallel multilevel method makes use of on-the-fly parallel domain decomposition to achieve extremely good scalability on large numbers of processors, and is coupled with an automatic coarse mesh generation algorithm for efficient processing by a multigrid smoother. Numerical results are presented demonstrating parallel speed-ups of up to 435 on 512 processors. Solution-based adaptation may be keyed off truncation error estimates using tau-extrapolation or a variety of feature detection based refinement parameters. The multigrid method is extended to for time-dependent flows through the use of a dual-time approach. The extension to rigid domain motion uses an Arbitrary Lagrangian-Eulerlarian (ALE) formulation, and results will be presented for a variety of two- and three-dimensional example problems with both simple and complex geometry.
Parallel stitching of 2D materials
Ling, Xi; Wu, Lijun; Lin, Yuxuan; Ma, Qiong; Wang, Ziqiang; Song, Yi; Yu, Lili; Huang, Shengxi; Fang, Wenjing; Zhang, Xu; et al
2016-01-27
Diverse parallel stitched 2D heterostructures, including metal–semiconductor, semiconductor–semiconductor, and insulator–semiconductor, are synthesized directly through selective “sowing” of aromatic molecules as the seeds in the chemical vapor deposition (CVD) method. Lastly, the methodology enables the large-scale fabrication of lateral heterostructures, which offers tremendous potential for its application in integrated circuits.
Parallel Stitching of 2D Materials.
Ling, Xi; Lin, Yuxuan; Ma, Qiong; Wang, Ziqiang; Song, Yi; Yu, Lili; Huang, Shengxi; Fang, Wenjing; Zhang, Xu; Hsu, Allen L; Bie, Yaqing; Lee, Yi-Hsien; Zhu, Yimei; Wu, Lijun; Li, Ju; Jarillo-Herrero, Pablo; Dresselhaus, Mildred; Palacios, Tomás; Kong, Jing
2016-03-01
Diverse parallel stitched 2D heterostructures, including metal-semiconductor, semiconductor-semiconductor, and insulator-semiconductor, are synthesized directly through selective "sowing" of aromatic molecules as the seeds in the chemical vapor deposition (CVD) method. The methodology enables the large-scale fabrication of lateral heterostructures, which offers tremendous potential for its application in integrated circuits. PMID:26813882
Baby universes in 2d quantum gravity
NASA Astrophysics Data System (ADS)
Ambjørn, Jan; Jain, Sanjay; Thorleifsson, Gudmar
1993-06-01
We investigate the fractal structure of 2d quantum gravity, both for pure gravity and for gravity coupled to multiple gaussian fields and for gravity coupled to Ising spins. The roughness of the surfaces is described in terms of baby universes and using numerical simulations we measure their distribution which is related to the string susceptibility exponent γstring.
Steady 2D Detonations and the DSD Theory
NASA Astrophysics Data System (ADS)
Lubyatinsky, S. N.; Loboiko, B. G.; Filin, V. P.; Kostitsin, O. V.; Smirnov, E. B.
2006-07-01
In the framework of the simplest DSD theory we obtained ODEs describing steady 2D detonation front shapes for slab, cylinder, and rib geometries. It was found that one solution (a steady detonation front shape) corresponds to several combinations of the confinement material and the defining charge dimension (slab thickness, cylinder radius, or inner rib radius). Comparing experimental data for these combinations and analyzing the shape difference at the edge provide information on the D(κ) relation at low D. The analysis of the data on IHE ribs detonation indicates that as D decreases total curvature κ tends to a limit of about 0.1 mm-1, i.e., double the reciprocal critical diameter. This correction makes the DSD theory consistent with the experimental critical diameter.
2D numerical modelling of meandering channel formation
NASA Astrophysics Data System (ADS)
XIAO, Y.; ZHOU, G.; YANG, F. S.
2016-03-01
A 2D depth-averaged model for hydrodynamic sediment transport and river morphological adjustment was established. The sediment transport submodel takes into account the influence of non-uniform sediment with bed surface armoring and considers the impact of secondary flow in the direction of bed-load transport and transverse slope of the river bed. The bank erosion submodel incorporates a simple simulation method for updating bank geometry during either degradational or aggradational bed evolution. Comparison of the results obtained by the extended model with experimental and field data, and numerical predictions validate that the proposed model can simulate grain sorting in river bends and duplicate the characteristics of meandering river and its development. The results illustrate that by using its control factors, the improved numerical model can be applied to simulate channel evolution under different scenarios and improve understanding of patterning processes.
Novel antenna coupled 2D plasmonic terahertz detection.
Allen, Jim; Dyer, Greg; Reno, John Louis; Shaner, Eric Arthur
2010-03-01
Resonant plasmonic detectors are potentially important for terahertz (THz) spectroscopic imaging. We have fabricated and characterized antenna coupled detectors that integrate a broad-band antenna, which improves coupling of THz radiation. The vertex of the antenna contains the tuning gates and the bolometric barrier gate. Incident THz radiation may excite 2D plasmons with wave-vectors defined by either a periodic grating gate or a plasmonic cavity determined by ohmic contacts and gate terminals. The latter approach of exciting plasmons in a cavity defined by a short micron-scale channel appears most promising. With this short-channel geometry, we have observed multiple harmonics of THz plasmons. At 20 K with detector bias optimized we report responsivity on resonance of 2.5 kV/W and an NEP of 5 x 10{sup -10} W/Hz{sup 1/2}.
Interpretive 2-D treatment of scrape-off-layer plasmas
Umansky, M.; Allen, A.; Daughton, W.
1996-12-31
The width of the scrape-off-layer in a tokamak is determined by cross field transport. In Alcator C-mod the plasma parameters in the scrape-off-layer are measured at upstream and divertor plate locations. We solve a 2-D scrape-off-layer heat conduction equation in the flux geometry (as determined by EFIT) of the C-mod experiment. Bolometric measurements are utilized for the radiative loss term. We use the end wall probe measurements of electron temperature as a boundary condition and the fast scanning probe measurements of upstream temperature are treated as constraints to determine the cross field transport and thermal conductivity. Results are compared with 1-D onion-skin-model predictions.
2D foam coarsening in a microfluidic system
NASA Astrophysics Data System (ADS)
Marchalot, J.; Lambert, J.; Cantat, I.; Tabeling, P.; Jullien, M.-C.
2008-09-01
We report an experimental study of 2D microfoam coarsening confined in a micrometer scale geometry, the typical bubbles diameter being of the order of 50-100 μm. These experiments raise both fundamental and applicative issues. For applicative issues: what is the typical time of foam ageing (for a polydisperse foam) in microsystems in scope of gas pocket storage in lab-on-a-chips? Experimental results show that a typical time of 2-3 mn is found, leading to the possibility of short-time storing, depending on the application. For fundamental interests, 2D foam ageing is generally described by von Neumann's law (von Neumann J., Metal Interfaces (American Society of Metals, Cleveland) 1952, p. 108) which is based on the hypothesis that bubbles are separated by thin films. Does this hypothesis still hold for foams confined in a 40 μm height geometry? This problematic is analyzed and it is shown that von Neumann's law still holds but that the diffusion coefficient involved in this law is modified by the confinement which imposes a curvature radius at Plateau borders. More precisely, it is shown that the liquid fraction is high on a film cross-section, in contrast with macrometric experiments where drainage occurs. An analytical description of the diffusion is developped taking into account the fact that soap film height is only a fraction of the cell height. While most of microfoams are flowing, the experimental set-up we describe leads to the achievement of a motionless confined microfoam.
Modelling RF sources using 2-D PIC codes
Eppley, K.R.
1993-03-01
In recent years, many types of RF sources have been successfully modelled using 2-D PIC codes. Both cross field devices (magnetrons, cross field amplifiers, etc.) and pencil beam devices (klystrons, gyrotrons, TWT`S, lasertrons, etc.) have been simulated. All these devices involve the interaction of an electron beam with an RF circuit. For many applications, the RF structure may be approximated by an equivalent circuit, which appears in the simulation as a boundary condition on the electric field (``port approximation``). The drive term for the circuit is calculated from the energy transfer between beam and field in the drift space. For some applications it may be necessary to model the actual geometry of the structure, although this is more expensive. One problem not entirely solved is how to accurately model in 2-D the coupling to an external waveguide. Frequently this is approximated by a radial transmission line, but this sometimes yields incorrect results. We also discuss issues in modelling the cathode and injecting the beam into the PIC simulation.
Modelling RF sources using 2-D PIC codes
Eppley, K.R.
1993-03-01
In recent years, many types of RF sources have been successfully modelled using 2-D PIC codes. Both cross field devices (magnetrons, cross field amplifiers, etc.) and pencil beam devices (klystrons, gyrotrons, TWT'S, lasertrons, etc.) have been simulated. All these devices involve the interaction of an electron beam with an RF circuit. For many applications, the RF structure may be approximated by an equivalent circuit, which appears in the simulation as a boundary condition on the electric field ( port approximation''). The drive term for the circuit is calculated from the energy transfer between beam and field in the drift space. For some applications it may be necessary to model the actual geometry of the structure, although this is more expensive. One problem not entirely solved is how to accurately model in 2-D the coupling to an external waveguide. Frequently this is approximated by a radial transmission line, but this sometimes yields incorrect results. We also discuss issues in modelling the cathode and injecting the beam into the PIC simulation.
3D track initiation in clutter using 2D measurements
NASA Astrophysics Data System (ADS)
Lin, Lin; Kirubarajan, Thiagalingam; Bar-Shalom, Yaakov
2001-11-01
In this paper we present an algorithm for initiating 3-D tracks using range and azimuth (bearing) measurements from a 2-D radar on a moving platform. The work is motivated by the need to track possibly low-flying targets, e.g., cruise missiles, using reports from an aircraft-based surveillance radar. Previous work on this problem considered simple linear motion in a flat earth coordinate frame. Our research extends this to a more realistic scenario where the earth"s curvature is also considered. The target is assumed to be moving along a great circle at a constant altitude. After the necessary coordinate transformations, the measurements are nonlinear functions of the target state and the observability of target altitude is severely limited. The observability, quantified by the Cramer-Rao Lower Bound (CRLB), is very sensitive to the sensor-to-target geometry. The paper presents a Maximum Likelihood (ML) estimator for estimating the target motion parameters in the Earth Centered Earth Fixed coordinate frame from 2-D range and angle measurements. In order to handle the possibility of false measurements and missed detections, which was not considered in, we use the Probabilistic Data Association (PDA) algorithm to weight the detections in a frame. The PDA-based modified global likelihood is optimized using a numerical search. The accuracies obtained by the resulting ML-PDA estimator are quantified using the CRLB for different sensor-target configurations. It is shown that the proposed estimator is efficient, that is, it meets the CRLB. Of particular interest is the achievable accuracy for estimating the target altitude, which is not observed directly by the 2-D radar, but can be only inferred from the range and bearing observations.