Science.gov

Sample records for quadratic tetrahedral elements

  1. Finite element simulation of articular contact mechanics with quadratic tetrahedral elements.

    PubMed

    Maas, Steve A; Ellis, Benjamin J; Rawlins, David S; Weiss, Jeffrey A

    2016-03-21

    Although it is easier to generate finite element discretizations with tetrahedral elements, trilinear hexahedral (HEX8) elements are more often used in simulations of articular contact mechanics. This is due to numerical shortcomings of linear tetrahedral (TET4) elements, limited availability of quadratic tetrahedron elements in combination with effective contact algorithms, and the perceived increased computational expense of quadratic finite elements. In this study we implemented both ten-node (TET10) and fifteen-node (TET15) quadratic tetrahedral elements in FEBio (www.febio.org) and compared their accuracy, robustness in terms of convergence behavior and computational cost for simulations relevant to articular contact mechanics. Suitable volume integration and surface integration rules were determined by comparing the results of several benchmark contact problems. The results demonstrated that the surface integration rule used to evaluate the contact integrals for quadratic elements affected both convergence behavior and accuracy of predicted stresses. The computational expense and robustness of both quadratic tetrahedral formulations compared favorably to the HEX8 models. Of note, the TET15 element demonstrated superior convergence behavior and lower computational cost than both the TET10 and HEX8 elements for meshes with similar numbers of degrees of freedom in the contact problems that we examined. Finally, the excellent accuracy and relative efficiency of these quadratic tetrahedral elements was illustrated by comparing their predictions with those for a HEX8 mesh for simulation of articular contact in a fully validated model of the hip. These results demonstrate that TET10 and TET15 elements provide viable alternatives to HEX8 elements for simulation of articular contact mechanics. PMID:26900037

  2. An Improved Linear Tetrahedral Element for Plasticity

    SciTech Connect

    Puso, M

    2005-04-25

    A stabilized, nodally integrated linear tetrahedral is formulated and analyzed. It is well known that linear tetrahedral elements perform poorly in problems with plasticity, nearly incompressible materials, and acute bending. For a variety of reasons, linear tetrahedral elements are preferable to quadratic tetrahedral elements in most nonlinear problems. Whereas, mixed methods work well for linear hexahedral elements, they don't for linear tetrahedrals. On the other hand, automatic mesh generation is typically not feasible for building many 3D hexahedral meshes. A stabilized, nodally integrated linear tetrahedral is developed and shown to perform very well in problems with plasticity, nearly incompressible materials and acute bending. Furthermore, the formulation is analytically and numerically shown to be stable and optimally convergent. The element is demonstrated to perform well in several standard linear and nonlinear benchmarks.

  3. An 8-node tetrahedral finite element suitable for explicit transient dynamic simulations

    SciTech Connect

    Key, S.W.; Heinstein, M.W.; Stone, C.M.

    1997-12-31

    Considerable effort has been expended in perfecting the algorithmic properties of 8-node hexahedral finite elements. Today the element is well understood and performs exceptionally well when used in modeling three-dimensional explicit transient dynamic events. However, the automatic generation of all-hexahedral meshes remains an elusive achievement. The alternative of automatic generation for all-tetrahedral finite element is a notoriously poor performer, and the 10-node quadratic tetrahedral finite element while a better performer numerically is computationally expensive. To use the all-tetrahedral mesh generation extant today, the authors have explored the creation of a quality 8-node tetrahedral finite element (a four-node tetrahedral finite element enriched with four midface nodal points). The derivation of the element`s gradient operator, studies in obtaining a suitable mass lumping and the element`s performance in applications are presented. In particular, they examine the 80node tetrahedral finite element`s behavior in longitudinal plane wave propagation, in transverse cylindrical wave propagation, and in simulating Taylor bar impacts. The element only samples constant strain states and, therefore, has 12 hourglass modes. In this regard, it bears similarities to the 8-node, mean-quadrature hexahedral finite element. Given automatic all-tetrahedral meshing, the 8-node, constant-strain tetrahedral finite element is a suitable replacement for the 8-node hexahedral finite element and handbuilt meshes.

  4. A Family of Uniform Strain Tetrahedral Elements and a Method for Connecting Dissimilar Finite Element Meshes

    SciTech Connect

    Dohrmann, C.R.; Heinstein, M.W.; Jung, J.; Key, S.W.

    1999-01-01

    This report documents a collection of papers on a family of uniform strain tetrahedral finite elements and their connection to different element types. Also included in the report are two papers which address the general problem of connecting dissimilar meshes in two and three dimensions. Much of the work presented here was motivated by the development of the tetrahedral element described in the report "A Suitable Low-Order, Eight-Node Tetrahedral Finite Element For Solids," by S. W. Key {ital et al.}, SAND98-0756, March 1998. Two basic issues addressed by the papers are: (1) the performance of alternative tetrahedral elements with uniform strain and enhanced uniform strain formulations, and (2) the proper connection of tetrahedral and other element types when two meshes are "tied" together to represent a single continuous domain.

  5. A computational study of nodal-based tetrahedral element behavior.

    SciTech Connect

    Gullerud, Arne S.

    2010-09-01

    This report explores the behavior of nodal-based tetrahedral elements on six sample problems, and compares their solution to that of a corresponding hexahedral mesh. The problems demonstrate that while certain aspects of the solution field for the nodal-based tetrahedrons provide good quality results, the pressure field tends to be of poor quality. Results appear to be strongly affected by the connectivity of the tetrahedral elements. Simulations that rely on the pressure field, such as those which use material models that are dependent on the pressure (e.g. equation-of-state models), can generate erroneous results. Remeshing can also be strongly affected by these issues. The nodal-based test elements as they currently stand need to be used with caution to ensure that their numerical deficiencies do not adversely affect critical values of interest.

  6. Tetrahedral mesh improvement via optimization of the element condition number

    SciTech Connect

    FREITAG,LORI A.; KNUPP,PATRICK

    2000-05-22

    The authors present a new shape measure for tetrahedral elements that is optimal in that it gives the distance of a tetrahedron from the set of inverted elements. This measure is constructed from the condition number of the linear transformation between a unit equilateral tetrahedron and any tetrahedron with positive volume. Using this shape measure, they formulate two optimization objective functions that are differentiated by their goal: the first seeks to improve the average quality of the tetrahedral mesh; the second aims to improve the worst-quality element in the mesh. They review the optimization techniques used with each objective function and presents experimental results that demonstrate the effectiveness of the mesh improvement methods. They show that a combined optimization approach that uses both objective functions obtains the best-quality meshes for several complex geometries.

  7. A finite element boundary integral formulation for radiation and scattering by cavity antennas using tetrahedral elements

    NASA Technical Reports Server (NTRS)

    Gong, J.; Volakis, J. L.; Chatterjee, A.; Jin, J. M.

    1992-01-01

    A hybrid finite element boundary integral formulation is developed using tetrahedral and/or triangular elements for discretizing the cavity and/or aperture of microstrip antenna arrays. The tetrahedral elements with edge based linear expansion functions are chosen for modeling the volume region and triangular elements are used for discretizing the aperture. The edge based expansion functions are divergenceless thus removing the requirement to introduce a penalty term and the tetrahedral elements permit greater geometrical adaptability than the rectangular bricks. The underlying theory and resulting expressions are discussed in detail together with some numerical scattering examples for comparison and demonstration.

  8. A mesh generator for tetrahedral elements using Delaunay triangulation

    SciTech Connect

    Yuan, J.S.; Fitzsimons, C.J. )

    1993-03-01

    A tetrahedral mesh generator has been developed. The generator is based on the Delaunay triangulation which is implemented by employing the insertion polyhedron algorithm. In this paper some new methods to deal with the problems associated with the three-dimensional Delaunay triangulation and the insertion polyhedron algorithm are presented: degeneracy, the crossing situation, identification of the internal elements and internal point generation. The generator works both for convex and non-convex domains, including those with high aspect-ratio subdomains. Some examples are given in this paper to illustrate the capability of the generator.

  9. A suitable low-order, eight-node tetrahedral finite element for solids

    SciTech Connect

    Key, S.W.; Heinstein, M.S.; Stone, C.M.; Mello, F.J.; Blanford, M.L.; Budge, K.G.

    1998-03-01

    To use the all-tetrahedral mesh generation existing today, the authors have explored the creation of a computationally efficient eight-node tetrahedral finite element (a four-node tetrahedral finite element enriched with four mid-face nodal points). The derivation of the element`s gradient operator, studies in obtaining a suitable mass lumping, and the element`s performance in applications are presented. In particular they examine the eight-node tetrahedral finite element`s behavior in longitudinal plane wave propagation, in transverse cylindrical wave propagation, and in simulating Taylor bar impacts. The element samples only constant strain states and, therefore, has 12 hour-glass modes. In this regard it bears similarities to the eight-node, mean-quadrature hexahedral finite element. Comparisons with the results obtained from the mean-quadrature eight-node hexahedral finite element and the four-node tetrahedral finite element are included. Given automatic all-tetrahedral meshing, the eight-node, constant-strain tetrahedral finite element is a suitable replacement for the eight-node hexahedral finite element in those cases where mesh generation requires an inordinate amount of user intervention and direction to obtain acceptable mesh properties.

  10. A 3D Frictional Segment-to-Segment Contact Method for Large Deformations and Quadratic Elements

    SciTech Connect

    Puso, M; Laursen, T; Solberg, J

    2004-04-01

    Node-on-segment contact is the most common form of contact used today but has many deficiencies ranging from potential locking to non-smooth behavior with large sliding. Furthermore, node-on-segment approaches are not at all applicable to higher order discretizations (e.g. quadratic elements). In a previous work, [3, 4] we developed a segment-to-segment contact approach for eight node hexahedral elements based on the mortar method that was applicable to large deformation mechanics. The approach proved extremely robust since it eliminated the over-constraint that caused 'locking' and provided smooth force variations in large sliding. Here, we extend this previous approach to treat frictional contact problems. In addition, the method is extended to 3D quadratic tetrahedrals and hexahedrals. The proposed approach is then applied to several challenging frictional contact problems that demonstrate its effectiveness.

  11. Stabilized tetrahedral elements for crystal plasticity finite element analysis overcoming volumetric locking

    NASA Astrophysics Data System (ADS)

    Cheng, Jiahao; Shahba, Ahmad; Ghosh, Somnath

    2016-05-01

    Image-based CPFE modeling involves computer generation of virtual polycrystalline microstructures from experimental data, followed by discretization into finite element meshes. Discretization is commonly accomplished using three-dimensional four-node tetrahedral or TET4 elements, which conform to the complex geometries. It has been commonly observed that TET4 elements suffer from severe volumetric locking when simulating deformation of incompressible or nearly incompressible materials. This paper develops and examines three locking-free stabilized finite element formulations in the context of crystal plasticity finite element analysis. They include a node-based uniform strain (NUS) element, a locally integrated B-bar (LIB) based element and a F-bar patch (FP) based element. All three formulations are based on the partitioning of TET4 element meshes and integrating over patches to obtain favorable incompressibility constraint ratios without adding large degrees of freedom. The results show that NUS formulation introduces unstable spurious energy modes, while the LIB and FP elements stabilize the solutions and are preferred for reliable CPFE analysis. The FP element is found to be computationally efficient over the LIB element.

  12. Uniform Strain Elements for Three-Node Triangular and Four-Node Tetrahedral Meshes

    SciTech Connect

    Dohrmann, C.R.; Heinstein, M.W.; Jung, J.; Key, S.W.; Witkowski, W.R.

    1999-03-02

    A family of uniform strain elements is presented for three-node triangular and four-node tetrahedral meshes. The elements use the linear interpolation functions of the original mesh, but each element is associated with a single node. As a result, a favorable constraint ratio for the volumetric response is obtained for problems in solid mechanics. The uniform strain elements do not require the introduction of additional degrees of freedom and their performance is shown to be significantly better than that of three-node triangular or four-node tetrahedral elements. In addition, nodes inside the boundary of the mesh are observed to exhibit superconvergent behavior for a set of example problems.

  13. Tetrahedral element shape optimization via the Jacobian determinant and condition number.

    SciTech Connect

    Freitag, L. A.; Knupp, P. M.

    1999-07-30

    We present a new shape measure for tetrahedral elements that is optimal in the sense that it gives the distance of a tetrahedron from the set of inverted elements. This measure is constructed from the condition number of the linear transformation between a unit equilateral tetrahedron and any tetrahedron with positive volume. We use this shape measure to formulate two optimization objective functions that are differentiated by their goal: the first seeks to improve the average quality of the tetrahedral mesh; the second aims to improve the worst-quality element in the mesh. Because the element condition number is not defined for tetrahedral with negative volume, these objective functions can be used only when the initial mesh is valid. Therefore, we formulate a third objective function using the determinant of the element Jacobian that is suitable for mesh untangling. We review the optimization techniques used with each objective function and present experimental results that demonstrate the effectiveness of the mesh improvement and untangling methods. We show that a combined optimization approach that uses both condition number objective functions obtains the best-quality meshes.

  14. QUADRATIC SERENDIPITY FINITE ELEMENTS ON POLYGONS USING GENERALIZED BARYCENTRIC COORDINATES.

    PubMed

    Rand, Alexander; Gillette, Andrew; Bajaj, Chandrajit

    2014-01-01

    We introduce a finite element construction for use on the class of convex, planar polygons and show it obtains a quadratic error convergence estimate. On a convex n-gon, our construction produces 2n basis functions, associated in a Lagrange-like fashion to each vertex and each edge midpoint, by transforming and combining a set of n(n + 1)/2 basis functions known to obtain quadratic convergence. The technique broadens the scope of the so-called 'serendipity' elements, previously studied only for quadrilateral and regular hexahedral meshes, by employing the theory of generalized barycentric coordinates. Uniform a priori error estimates are established over the class of convex quadrilaterals with bounded aspect ratio as well as over the class of convex planar polygons satisfying additional shape regularity conditions to exclude large interior angles and short edges. Numerical evidence is provided on a trapezoidal quadrilateral mesh, previously not amenable to serendipity constructions, and applications to adaptive meshing are discussed. PMID:25301974

  15. QUADRATIC SERENDIPITY FINITE ELEMENTS ON POLYGONS USING GENERALIZED BARYCENTRIC COORDINATES

    PubMed Central

    RAND, ALEXANDER; GILLETTE, ANDREW; BAJAJ, CHANDRAJIT

    2013-01-01

    We introduce a finite element construction for use on the class of convex, planar polygons and show it obtains a quadratic error convergence estimate. On a convex n-gon, our construction produces 2n basis functions, associated in a Lagrange-like fashion to each vertex and each edge midpoint, by transforming and combining a set of n(n + 1)/2 basis functions known to obtain quadratic convergence. The technique broadens the scope of the so-called ‘serendipity’ elements, previously studied only for quadrilateral and regular hexahedral meshes, by employing the theory of generalized barycentric coordinates. Uniform a priori error estimates are established over the class of convex quadrilaterals with bounded aspect ratio as well as over the class of convex planar polygons satisfying additional shape regularity conditions to exclude large interior angles and short edges. Numerical evidence is provided on a trapezoidal quadrilateral mesh, previously not amenable to serendipity constructions, and applications to adaptive meshing are discussed. PMID:25301974

  16. Quadratic Finite Element Method for 1D Deterministic Transport

    SciTech Connect

    Tolar, Jr., D R; Ferguson, J M

    2004-01-06

    In the discrete ordinates, or SN, numerical solution of the transport equation, both the spatial ({und r}) and angular ({und {Omega}}) dependences on the angular flux {psi}{und r},{und {Omega}}are modeled discretely. While significant effort has been devoted toward improving the spatial discretization of the angular flux, we focus on improving the angular discretization of {psi}{und r},{und {Omega}}. Specifically, we employ a Petrov-Galerkin quadratic finite element approximation for the differencing of the angular variable ({mu}) in developing the one-dimensional (1D) spherical geometry S{sub N} equations. We develop an algorithm that shows faster convergence with angular resolution than conventional S{sub N} algorithms.

  17. Workflows for generating tetrahedral meshes for finite element simulations on complex geological structures

    NASA Astrophysics Data System (ADS)

    Zehner, Björn; Börner, Jana H.; Görz, Ines; Spitzer, Klaus

    2015-06-01

    Subsurface processing numerical simulations require accurate discretization of the modeling domain such that the geological units are represented correctly. Unstructured tetrahedral grids are particularly flexible in adapting to the shape of geo-bodies and are used in many finite element codes. In order to generate a tetrahedral mesh on a 3D geological model, the tetrahedrons have to belong completely to one geological unit and have to describe geological boundaries by connected facets of tetrahedrons. This is especially complicated at the contact points between several units and for irregular sharp-shaped bodies, especially in case of faulted zones. This study develops, tests and validates three workflows to generate a good tetrahedral mesh from a geological basis model. The tessellation of the model needs (i) to be of good quality to guarantee a stable calculation, (ii) to include certain nodes to apply boundary conditions for the numerical solution, and (iii) support local mesh refinement. As a test case we use the simulation of a transient electromagnetic measurement above a salt diapir. We can show that the suggested workflows lead to a tessellation of the structure on which the simulation can be run robustly. All workflows show advantages and disadvantages with respect to the workload, the control the user has over the resulting mesh and the skills in software handling that are required.

  18. 3-D inversion of magnetotelluric data using unstructured tetrahedral elements: applicability to data affected by topography

    NASA Astrophysics Data System (ADS)

    Usui, Yoshiya

    2015-08-01

    A 3-D magnetotelluric (MT) inversion code using unstructured tetrahedral elements has been developed in order to correct the topographic effect by directly incorporating it into computational grids. The electromagnetic field and response functions get distorted at the observation sites of MT surveys because of the undulating surface topography, and without correcting this distortion, the subsurface structure can be misinterpreted. Of the two methods proposed to correct the topographic effect, the method incorporating topography explicitly in the inversion is applicable to a wider range of surveys. For forward problems, it has been shown that the finite element method using unstructured tetrahedral elements is useful for the incorporation of topography. Therefore, this paper shows the applicability of unstructured tetrahedral elements in MT inversion using the newly developed code. The inversion code is capable of using the impedance tensor, the vertical magnetic transfer function (VMTF), and the phase tensor as observational data, and it estimates the subsurface resistivity values and the distortion tensor of each observation site. The forward part of the code was verified using two test models, one incorporating topographic effect and one without, and the verifications showed that the results were almost the same as those of previous works. The developed inversion code was then applied to synthetic data from a MT survey, and was verified as being able to recover the resistivity structure as well as other inversion codes. Finally, to confirm its applicability to the data affected by topography, inversion was performed using the synthetic data of the model that included two overlapping mountains. In each of the cases using the impedance tensor, the VMTF and the phase tensor, by including the topography in the mesh, the subsurface resistivity was determined more proficiently than in the case using the flat-surface mesh. Although the locations of the anomalies were

  19. A three-dimensional FE analysis of large deformations for impact loadings using tetrahedral elements

    NASA Astrophysics Data System (ADS)

    Yoo, Y. H.; Lee, M.

    A three-dimensional dynamic program for the anaysis of large deformations in contact-penetration problems is developed using the finite element Lagrangian method with explicit time integration. By incorporating a tetrahedral element, which allows a single-point integration without a special hourglass control scheme, this program can be more effective to the present problem. The position code algorithm is used to search contact surface. Eroding surfaces are also considered. The defense node algorithm was slightly modified for the calculation of contact forces. A study of obliquity effects on metallic plate perforation and ricochet processes in thin plates impacted by a sphere was conducted. It is well simulated that on separation of two parts of the sphere, the portion still within the crater tends to perforate, while the portion in contact with the plate surface ricochets. This deformation pattern is observed in experiments, especially at high obliquities. A long rod that impacts an oblique steel plate at high impact velocity was also simulated in order to study the dynamics of the rod caused by the three dimensional asymmetric contact. The agreement between simulated and experimental results is quite good. Fracture phenomena occuring at high obliquity deserves further investigations.

  20. Fast calculation of the sensitivity matrix in magnetic induction tomography by tetrahedral edge finite elements and the reciprocity theorem.

    PubMed

    Hollaus, K; Magele, C; Merwa, R; Scharfetter, H

    2004-02-01

    Magnetic induction tomography of biological tissue is used to reconstruct the changes in the complex conductivity distribution by measuring the perturbation of an alternating primary magnetic field. To facilitate the sensitivity analysis and the solution of the inverse problem a fast calculation of the sensitivity matrix, i.e. the Jacobian matrix, which maps the changes of the conductivity distribution onto the changes of the voltage induced in a receiver coil, is needed. The use of finite differences to determine the entries of the sensitivity matrix does not represent a feasible solution because of the high computational costs of the basic eddy current problem. Therefore, the reciprocity theorem was exploited. The basic eddy current problem was simulated by the finite element method using symmetric tetrahedral edge elements of second order. To test the method various simulations were carried out and discussed. PMID:15005313

  1. Finite-element analysis of earing using non-quadratic yield surfaces

    SciTech Connect

    Logan, R.W.

    1995-06-18

    During deep draw cupping, the phenomenon known as earing may occur as the cup wall is formed, resulting in a periodic variation of cup wall height around the perimeter of the finished cup. This is generally due to planar anisotropy of flow in rolled sheet product. It is generally observed that the anisotropy parameter R will vary in the plane of the sheet when ears are observed in cupping, with a parameter {Delta}R describing the variation of R in the plane of the sheet. For many common textures in face-centered and body-centered materials, the ears form relative to the sheet rolling direction at 0{degrees} and 90{degrees} around the perimeter if {Delta}R>0, and at -45{degrees} and +45{degrees} if {Delta}R<0. There is extensive experimental evidence that ear height shows a linear correlation with {Delta}R/R, but attempts to duplicate this using the finite-element method are highly dependent on both the methodology and yield surface used. It was shown previously that using a coarse mesh and the quadratic Hill yield surface tends to greatly under predict earing. In this study, we have used two different finite-element codes developed at LLNL to examine the predicted earing using both quadratic Hill and alternative non-quadratic yield surfaces. These results are compared to experimental data and conclusions drawn about the most desirable closed-form yield surfaces to duplicate the observed earing phenomena.

  2. Solving the transport equation with quadratic finite elements: Theory and applications

    SciTech Connect

    Ferguson, J.M.

    1997-12-31

    At the 4th Joint Conference on Computational Mathematics, the author presented a paper introducing a new quadratic finite element scheme (QFEM) for solving the transport equation. In the ensuing year the author has obtained considerable experience in the application of this method, including solution of eigenvalue problems, transmission problems, and solution of the adjoint form of the equation as well as the usual forward solution. He will present detailed results, and will also discuss other refinements of his transport codes, particularly for 3-dimensional problems on rectilinear and non-rectilinear grids.

  3. Parallel Anisotropic Tetrahedral Adaptation

    NASA Technical Reports Server (NTRS)

    Park, Michael A.; Darmofal, David L.

    2008-01-01

    An adaptive method that robustly produces high aspect ratio tetrahedra to a general 3D metric specification without introducing hybrid semi-structured regions is presented. The elemental operators and higher-level logic is described with their respective domain-decomposed parallelizations. An anisotropic tetrahedral grid adaptation scheme is demonstrated for 1000-1 stretching for a simple cube geometry. This form of adaptation is applicable to more complex domain boundaries via a cut-cell approach as demonstrated by a parallel 3D supersonic simulation of a complex fighter aircraft. To avoid the assumptions and approximations required to form a metric to specify adaptation, an approach is introduced that directly evaluates interpolation error. The grid is adapted to reduce and equidistribute this interpolation error calculation without the use of an intervening anisotropic metric. Direct interpolation error adaptation is illustrated for 1D and 3D domains.

  4. Performance of continuous mass-lumped tetrahedral elements for elastic wave propagation with and without global assembly

    NASA Astrophysics Data System (ADS)

    Mulder, W. A.; Shamasundar, R.

    2016-07-01

    We consider isotropic elastic wave propagation with continuous mass-lumped finite elements on tetrahedra with explicit time stepping. These elements require higher-order polynomials in their interior to preserve accuracy after mass lumping and are only known up to degree 3. Global assembly of the symmetric stiffness matrix is a natural approach but requires large memory. Local assembly on the fly, in the form of matrix-vector products per element at each time step, has a much smaller memory footprint. With dedicated expressions for local assembly, our code ran about 1.3 times faster for degree 2 and 1.9 times for degree 3 on a simple homogeneous test problem, using 24 cores. This is similar to the acoustic case. For a more realistic problem, the gain in efficiency was a factor 2.5 for degree 2 and 3 for degree 3. For the lowest degree, the linear element, the expressions for both the global and local assembly can be further simplified. In that case, global assembly is more efficient than local assembly. Among the three degrees, the element of degree 3 is the most efficient in terms of accuracy at a given cost.

  5. Quadratic Damping

    ERIC Educational Resources Information Center

    Fay, Temple H.

    2012-01-01

    Quadratic friction involves a discontinuous damping term in equations of motion in order that the frictional force always opposes the direction of the motion. Perhaps for this reason this topic is usually omitted from beginning texts in differential equations and physics. However, quadratic damping is more realistic than viscous damping in many…

  6. Quadratic elongation: A quantitative measure of distortion in coordination polyhedra

    USGS Publications Warehouse

    Robinson, Kelly F.; Gibbs, G.V.; Ribbe, P.H.

    1971-01-01

    Quadratic elongation and the variance of bond angles are linearly correlated for distorted octahedral and tetrahedral coordination complexes, both of which show variations in bond length and bond angle. The quadratic elonga tion is dimensionless, giving a quantitative measure of polyhedral distortion which is independent of the effective size of the polyhedron.

  7. Optimization of tetrahedral meshes

    SciTech Connect

    Briere De L`Isle, E.; George, P.L.

    1995-12-31

    Finite element computations are all the more exact if we start from {open_quotes}good{close_quotes} elements. We are interested in meshes where the elements are tetrahedra and we shall develop utilities allowing us to improve the quality of these meshes.

  8. Parallel tetrahedral mesh refinement with MOAB.

    SciTech Connect

    Thompson, David C.; Pebay, Philippe Pierre

    2008-12-01

    In this report, we present the novel functionality of parallel tetrahedral mesh refinement which we have implemented in MOAB. This report details work done to implement parallel, edge-based, tetrahedral refinement into MOAB. The theoretical basis for this work is contained in [PT04, PT05, TP06] while information on design, performance, and operation specific to MOAB are contained herein. As MOAB is intended mainly for use in pre-processing and simulation (as opposed to the post-processing bent of previous papers), the primary use case is different: rather than refining elements with non-linear basis functions, the goal is to increase the number of degrees of freedom in some region in order to more accurately represent the solution to some system of equations that cannot be solved analytically. Also, MOAB has a unique mesh representation which impacts the algorithm. This introduction contains a brief review of streaming edge-based tetrahedral refinement. The remainder of the report is broken into three sections: design and implementation, performance, and conclusions. Appendix A contains instructions for end users (simulation authors) on how to employ the refiner.

  9. Lattice Cleaving: A Multimaterial Tetrahedral Meshing Algorithm with Guarantees

    PubMed Central

    Bronson, Jonathan; Levine, Joshua A.; Whitaker, Ross

    2014-01-01

    We introduce a new algorithm for generating tetrahedral meshes that conform to physical boundaries in volumetric domains consisting of multiple materials. The proposed method allows for an arbitrary number of materials, produces high-quality tetrahedral meshes with upper and lower bounds on dihedral angles, and guarantees geometric fidelity. Moreover, the method is combinatoric so its implementation enables rapid mesh construction. These meshes are structured in a way that also allows grading, to reduce element counts in regions of homogeneity. Additionally, we provide proofs showing that both element quality and geometric fidelity are bounded using this approach. PMID:24356365

  10. Self-Replicating Quadratics

    ERIC Educational Resources Information Center

    Withers, Christopher S.; Nadarajah, Saralees

    2012-01-01

    We show that there are exactly four quadratic polynomials, Q(x) = x [superscript 2] + ax + b, such that (x[superscript 2] + ax + b) (x[superscript 2] - ax + b) = (x[superscript 4] + ax[superscript 2] + b). For n = 1, 2, ..., these quadratic polynomials can be written as the product of N = 2[superscript n] quadratic polynomials in x[superscript…

  11. Feature-Sensitive Tetrahedral Mesh Generation with Guaranteed Quality

    PubMed Central

    Wang, Jun; Yu, Zeyun

    2012-01-01

    Tetrahedral meshes are being extensively used in finite element methods (FEM). This paper proposes an algorithm to generate feature-sensitive and high-quality tetrahedral meshes from an arbitrary surface mesh model. A top-down octree subdivision is conducted on the surface mesh and a set of tetrahedra are constructed using adaptive body-centered cubic (BCC) lattices. Special treatments are given to the tetrahedra near the surface such that the quality of the resulting tetrahedral mesh is provably guaranteed: the smallest dihedral angle is always greater than 5.71°. The meshes generated by our method are not only adaptive from the interior to the boundary, but also feature-sensitive on the surface with denser elements in high-curvature regions where geometric feature most likely reside. A variety of experimental results are presented to demonstrate the effectiveness and robustness of this algorithm. PMID:22328787

  12. Au20: A Tetrahedral Cluster

    SciTech Connect

    Li, Jun; Li, Xi; Zhai, Hua Jin; Wang, Lai S.

    2003-02-07

    Photoelectron spectroscopy revealed that a 20 atom gold cluster has an extremely large energy gap, which is even greater than that of C60, and an electron affinity comparable with that of C60. This observation suggests that the Au20 cluster must be extremely stable and chemically inert. Using relativistic density functional calculations, we found that Au20 possesses a remarkable tetrahedral structure, which is a fragment of the bulk face-centered cubic lattice of gold with a small structural relaxation. Au20 is thus a true cluster molecule, while at the same time it is exactly part of the bulk, but with very different properties. The tetrahedral Au20 may possess interesting catalytic properties and may be synthesized in bulk quantity or assembled on non-interacting surfaces.

  13. A tetrahedral entropy for water.

    PubMed

    Kumar, Pradeep; Buldyrev, Sergey V; Stanley, H Eugene

    2009-12-29

    We introduce the space-dependent correlation function C (Q)(r) and time-dependent autocorrelation function C (Q)(t) of the local tetrahedral order parameter Q identical with Q(r,t). By using computer simulations of 512 waterlike particles interacting through the transferable interaction potential with five points (TIP5 potential), we investigate C (Q)(r) in a broad region of the phase diagram. We find that at low temperatures C (Q)(t) exhibits a two-step time-dependent decay similar to the self-intermediate scattering function and that the corresponding correlation time tau(Q) displays a dynamic cross-over from non-Arrhenius behavior for T > T (W) to Arrhenius behavior for T < T (W), where T (W) denotes the Widom temperature where the correlation length has a maximum as T is decreased along a constant-pressure path. We define a tetrahedral entropy S (Q) associated with the local tetrahedral order of water molecules and find that it produces a major contribution to the specific heat maximum at the Widom line. Finally, we show that tau(Q) can be extracted from S (Q) by using an analog of the Adam-Gibbs relation. PMID:20018692

  14. Quadratic eigenvalue problems.

    SciTech Connect

    Walsh, Timothy Francis; Day, David Minot

    2007-04-01

    In this report we will describe some nonlinear eigenvalue problems that arise in the areas of solid mechanics, acoustics, and coupled structural acoustics. We will focus mostly on quadratic eigenvalue problems, which are a special case of nonlinear eigenvalue problems. Algorithms for solving the quadratic eigenvalue problem will be presented, along with some example calculations.

  15. Water's non-tetrahedral side.

    PubMed

    Henchman, Richard H; Cockram, Stuart J

    2013-01-01

    The case for liquid water having non-tetrahedral as well as tetrahedral coordination is put forward. Given the dependence of structure on the hydrogen bond definition, a recent conceptual breakthrough has been the topological hydrogen bond definition which overcomes the shortcomings of traditional cut-off-based hydrogen bond definitions. It identifies the hydrogen bonds in water's first coordination shell using assumed transition states as boundaries instead of fixed cut-offs. Here, the topological definition is applied to liquid water to characterise the distances, angles and energies of the hydrogen bonds for the different types of coordinations found. These coordinations include bent, trigonal, tetrahedral, trigonal bipyramidal, and octahedral structures, as well as bifurcated hydrogens, bifurcated oxygens and cyclic dimers, and larger polygons. All species are shown to have properties consistent with their classification, justifying their assignments, and supporting the structure of water as a continuous, single phase mixture. However, a detailed analysis to assess the existence of the assumed transition states reveals the remarkable finding that hydrogen bond switching via a bifurcated hydrogen under certain circumstances is a barrierless process. The likelihood of a switch depends on both the acceptor numbers and on the proximity of a donor to its acceptor. Specifically, a donor in an acceptor's outermost subshell switches uphill to an acceptor of the same or higher coordination to the starting acceptor, downhill to an acceptor of lower coordination by two or more, or sits bifurcated between two acceptors if the new acceptor has a coordination lower by only one. Which it is depends intimately on the donor molecule's oscillations and on other hydrogen bond switches that control the nearby acceptors' coordinations. Finally, a search is conducted for long-range structure in water in terms of asymmetry in the distribution of the donor-acceptor bias but none is

  16. Self-replicating quadratics

    NASA Astrophysics Data System (ADS)

    Withers, Christopher S.; Nadarajah, Saralees

    2012-06-01

    We show that there are exactly four quadratic polynomials, Q(x) = x 2 + ax + b, such that For n = 1, 2, … , these quadratic polynomials can be written as the product of N = 2 n quadratic polynomials in x 1/N , namely, ? , where w N is the Nth root of 1.

  17. Resolvability and the Tetrahedral Configuration of Carbon.

    ERIC Educational Resources Information Center

    Kauffman, George B.

    1983-01-01

    Discusses evidence for the tetrahedral configuration of the carbon atom, indicating that three symmetrical configurations are theoretically possible for coordination number four. Includes table indicating that resolvability of compounds of type CR'R"R"'R"" is a necessary but not sufficient condition for proving tetrahedral configuration. (JN)

  18. Building Tetrahedral Kites. Grades 6-8.

    ERIC Educational Resources Information Center

    Rushton, Erik; Ryan, Emily; Swift, Charles

    Working in teams of four, students build a tetrahedral kite following a specific set of directions and using specific provided materials. Students use basic processes of manufacturing systems-- cutting, shaping, forming, conditioning, assembling, joining, finishing, and quality control--to manufacture a complete tetrahedral kite within a given…

  19. Sequentially deployable maneuverable tetrahedral beam

    NASA Technical Reports Server (NTRS)

    Mikulas, M. M., Jr.; Crawford, R. F. (Inventor)

    1985-01-01

    A tetrahedral beam that can be compactly stowed, sequentially deployed, and widely manipulated to provide a structurally sound yet highly maneuverable truss structure is comprised of a number of repeating units of tandem tetralhedral sharing common sides. Fixed length battens are jointed into equilateral triangles called batten frames. Apexes of adjacent triangles are interconnected by longerons having a mid-point folding hinge. Joints, comprised of gussets pivotabley connected by links, permit two independent degrees of rotational freedom between joined adjacent batten frames, and provide a stable structure from packaged configuration to complete deployment. The longerons and joints can be actuated in any sequence, independently of one another. The beam is suited to remote actuation. Longerons may be provided with powered mid-point hinges enabling beam erection and packaging under remote control. Providing one or more longerons with powered telescoping segments permits the shape of the beam central axis to be remotely manipulated so that the beam may function as a remote manipulator arm.

  20. Sequentially deployable maneuverable tetrahedral beam

    NASA Astrophysics Data System (ADS)

    Mikulas, M. M., Jr.; Crawford, R. F.

    1985-12-01

    A tetrahedral beam that can be compactly stowed, sequentially deployed, and widely manipulated to provide a structurally sound yet highly maneuverable truss structure is comprised of a number of repeating units of tandem tetralhedral sharing common sides. Fixed length battens are jointed into equilateral triangles called batten frames. Apexes of adjacent triangles are interconnected by longerons having a mid-point folding hinge. Joints, comprised of gussets pivotabley connected by links, permit two independent degrees of rotational freedom between joined adjacent batten frames, and provide a stable structure from packaged configuration to complete deployment. The longerons and joints can be actuated in any sequence, independently of one another. The beam is suited to remote actuation. Longerons may be provided with powered mid-point hinges enabling beam erection and packaging under remote control. Providing one or more longerons with powered telescoping segments permits the shape of the beam central axis to be remotely manipulated so that the beam may function as a remote manipulator arm.

  1. Au40: A Large Tetrahedral Magic Cluster

    SciTech Connect

    Jiang, Deen; Walter, Michael

    2011-01-01

    40 is a magic number for tetrahedral symmetry predicted in both nuclear physics and the electronic jellium model. We show that Au{sub 40} could be such a magic cluster from density functional theory-based basin hopping for global minimization. The putative global minimum found for Au{sub 40} has a twisted pyramid structure, reminiscent of the famous tetrahedral Au{sub 20}, and a sizable HOMO-LUMO gap of 0.69 eV, indicating its molecular nature. Analysis of the electronic states reveals that the gap is related to shell closings of the metallic electrons in a tetrahedrally distorted effective potential.

  2. Au40: A large tetrahedral magic cluster

    NASA Astrophysics Data System (ADS)

    Jiang, De-En; Walter, Michael

    2011-11-01

    40 is a magic number for tetrahedral symmetry predicted in both nuclear physics and the electronic jellium model. We show that Au40 could be such a a magic cluster from density functional theory-based basin hopping for global minimization. The putative global minimum found for Au40 has a twisted pyramid structure, reminiscent of the famous tetrahedral Au20, and a sizable HOMO-LUMO gap of 0.69 eV, indicating its molecular nature. Analysis of the electronic states reveals that the gap is related to shell closings of the metallic electrons in a tetrahedrally distorted effective potential.

  3. A velocity-pressure integrated, mixed interpolation, Galerkin finite element method for high Reynolds number laminar flows

    NASA Technical Reports Server (NTRS)

    Kim, Sang-Wook

    1988-01-01

    A velocity-pressure integrated, mixed interpolation, Galerkin finite element method for the Navier-Stokes equations is presented. In the method, the velocity variables were interpolated using complete quadratic shape functions and the pressure was interpolated using linear shape functions. For the two dimensional case, the pressure is defined on a triangular element which is contained inside the complete biquadratic element for velocity variables; and for the three dimensional case, the pressure is defined on a tetrahedral element which is again contained inside the complete tri-quadratic element. Thus the pressure is discontinuous across the element boundaries. Example problems considered include: a cavity flow for Reynolds number of 400 through 10,000; a laminar backward facing step flow; and a laminar flow in a square duct of strong curvature. The computational results compared favorable with those of the finite difference methods as well as experimental data available. A finite elememt computer program for incompressible, laminar flows is presented.

  4. The Mystical "Quadratic Formula."

    ERIC Educational Resources Information Center

    March, Robert H.

    1993-01-01

    Uses projectile motion to explain the two roots found when using the quadratic formula. An example is provided for finding the time of flight for a projectile which has a negative root implying a negative time of flight. This negative time of flight also has a useful physical meaning. (MVL)

  5. A Quadratic Spring Equation

    ERIC Educational Resources Information Center

    Fay, Temple H.

    2010-01-01

    Through numerical investigations, we study examples of the forced quadratic spring equation [image omitted]. By performing trial-and-error numerical experiments, we demonstrate the existence of stability boundaries in the phase plane indicating initial conditions yielding bounded solutions, investigate the resonance boundary in the [omega]…

  6. Details of tetrahedral anisotropic mesh adaptation

    NASA Astrophysics Data System (ADS)

    Jensen, Kristian Ejlebjerg; Gorman, Gerard

    2016-04-01

    We have implemented tetrahedral anisotropic mesh adaptation using the local operations of coarsening, swapping, refinement and smoothing in MATLAB without the use of any for- N loops, i.e. the script is fully vectorised. In the process of doing so, we have made three observations related to details of the implementation: 1. restricting refinement to a single edge split per element not only simplifies the code, it also improves mesh quality, 2. face to edge swapping is unnecessary, and 3. optimising for the Vassilevski functional tends to give a little higher value for the mean condition number functional than optimising for the condition number functional directly. These observations have been made for a uniform and a radial shock metric field, both starting from a structured mesh in a cube. Finally, we compare two coarsening techniques and demonstrate the importance of applying smoothing in the mesh adaptation loop. The results pertain to a unit cube geometry, but we also show the effect of corners and edges by applying the implementation in a spherical geometry.

  7. MMS Spacecraft Transition to Tetrahedral Flying Formation

    NASA Video Gallery

    In the latter half of July 2015, the four satellites of the Magnetosphere Multi-scale (MMS) mission move into their tetrahedral formation flying configuration as part of the checkout for the scienc...

  8. Tetrahedrality and hydrogen bonds in water

    NASA Astrophysics Data System (ADS)

    Székely, Eszter; Varga, Imre K.; Baranyai, András

    2016-06-01

    We carried out extensive calculations of liquid water at different temperatures and pressures using the BK3 model suggested recently [P. T. Kiss and A. Baranyai, J. Chem. Phys. 138, 204507 (2013)]. In particular, we were interested in undercooled regions to observe the propensity of water to form tetrahedral coordination of closest neighbors around a central molecule. We compared the found tetrahedral order with the number of hydrogen bonds and with the partial pair correlation functions unfolded as distributions of the closest, the second closest, etc. neighbors. We found that contrary to the number of hydrogen bonds, tetrahedrality changes substantially with state variables. Not only the number of tetrahedral arrangements increases with lowering the pressure, the density, and the temperature but the domain size of connecting tetrahedral structures as well. The difference in tetrahedrality is very pronounced between the two sides of the Widom line and even more so between the low density amorphous (LDA) and high density amorphous (HDA) phases. We observed that in liquid water and in HDA, the 5th water molecule, contrary to ice and LDA, is positioned between the first and the second coordination shell. We found no convincing evidence of structural heterogeneity or regions referring to structural transition.

  9. Motion compensation for PET image reconstruction using deformable tetrahedral meshes

    NASA Astrophysics Data System (ADS)

    Manescu, P.; Ladjal, H.; Azencot, J.; Beuve, M.; Shariat, B.

    2015-12-01

    Respiratory-induced organ motion is a technical challenge to PET imaging. This motion induces displacements and deformation of the organs tissues, which need to be taken into account when reconstructing the spatial radiation activity. Classical image-based methods that describe motion using deformable image registration (DIR) algorithms cannot fully take into account the non-reproducibility of the respiratory internal organ motion nor the tissue volume variations that occur during breathing. In order to overcome these limitations, various biomechanical models of the respiratory system have been developed in the past decade as an alternative to DIR approaches. In this paper, we describe a new method of correcting motion artefacts in PET image reconstruction adapted to motion estimation models such as those based on the finite element method. In contrast with the DIR-based approaches, the radiation activity was reconstructed on deforming tetrahedral meshes. For this, we have re-formulated the tomographic reconstruction problem by introducing a time-dependent system matrix based calculated using tetrahedral meshes instead of voxelized images. The MLEM algorithm was chosen as the reconstruction method. The simulations performed in this study show that the motion compensated reconstruction based on tetrahedral deformable meshes has the capability to correct motion artefacts. Results demonstrate that, in the case of complex deformations, when large volume variations occur, the developed tetrahedral based method is more appropriate than the classical DIR-based one. This method can be used, together with biomechanical models controlled by external surrogates, to correct motion artefacts in PET images and thus reducing the need for additional internal imaging during the acquisition.

  10. A finite element computational method for high Reynolds number laminar flows

    NASA Technical Reports Server (NTRS)

    Kim, Sang-Wook

    1987-01-01

    A velocity-pressure integrated, mixed interpolation, Galerkin finite element method for the Navier-Stokes equations is presented. In the method, the velocity variables are interpolated using complete quadratic shape functions, and the pressure is interpolated using linear shape functions which are defined on a triangular element for the two-dimensional case and on a tetrahedral element for the three-dimensional case. The triangular element and the tetrahedral element are contained inside the complete bi- and tri-quadratic elements for velocity variables for two and three dimensional cases, respectively, so that the pressure is discontinuous across the element boundaries. Example problems considered include: a cavity flow of Reynolds numbers 400 through 10,000; a laminar backward facing step flow; and a laminar flow in a square duct of strong curvature. The computational results compared favorably with the finite difference computational results and/or experimental data available. It was found that the present method can capture the delicate pressure driven recirculation zones, that the method did not yield any spurious pressure modes, and that the method requires fewer grid points than the finite difference methods to obtain comparable computational results.

  11. Nonzero Quadrupole Moments of Candidate Tetrahedral Bands

    SciTech Connect

    Bark, R. A.; Lawrie, E. A.; Lawrie, J. J.; Mullins, S. M.; Murray, S. H. T.; Ncapayi, N. J.; Smit, F. D.; Sharpey-Schafer, J. F.; Lindsay, R.

    2010-01-15

    Negative-parity bands in the vicinity of {sup 156}Gd and {sup 160}Yb have been suggested as candidates for the rotation of tetrahedral nuclei. We report the observation of the odd and even-spin members of the lowest energy negative-parity bands in {sup 160}Yb and {sup 154}Gd. The properties of these bands are similar to the proposed tetrahedral band of {sup 156}Gd and its even-spin partner. Band-mixing calculations are performed and absolute and relative quadrupole moments deduced for {sup 160}Yb and {sup 154}Gd. The values are inconsistent with zero, as required for tetrahedral shape, and the bands are interpreted as octupole vibrational bands. The failure to observe the in-band E2 transitions of the bands at low spins can be understood using the measured B(E1) and B(E2) values.

  12. Solitons in quadratic media

    NASA Astrophysics Data System (ADS)

    Colin, M.; Di Menza, L.; Saut, J. C.

    2016-03-01

    In this paper, we investigate the properties of solitonic structures arising in quadratic media. First, we recall the derivation of systems governing the interaction process for waves propagating in such media and we check the local and global well-posedness of the corresponding Cauchy problem. Then, we look for stationary states in the context of normal or anomalous dispersion regimes, that lead us to either elliptic or non-elliptic systems and we address the problem of orbital stability. Finally, some numerical experiments are carried out in order to compute localized states for several regimes and to study dynamic stability as well as long-time asymptotics.

  13. Evidence for tetrahedral symmetry in (16)O.

    PubMed

    Bijker, R; Iachello, F

    2014-04-18

    We derive the rotation-vibration spectrum of a 4α configuration with tetrahedral symmetry Td and show evidence for the occurrence of this symmetry in the low-lying spectrum of (16)O. All vibrational states with A, E, and F symmetry appear to have been observed as well as the rotational bands with LP=0+, 3-, 4+, 6+ on the A states and part of the rotational bands built on the E, F states. We derive analytic expressions for the form factors and B(EL) values of the ground-state rotational band and show that the measured values support the tetrahedral symmetry of this band. PMID:24785032

  14. The NASA Tetrahedral Unstructured Software System (TETRUSS)

    NASA Technical Reports Server (NTRS)

    Frink, Neal T.; Pirzadeh, Shahyar Z.; Parikh, Paresh C.; Pandya, Mohagna J.; Bhat, M. K.

    2000-01-01

    The NASA Tetrahedral Unstructured Software System (TetrUSS) was developed during the 1990's to provide a rapid aerodynamic analysis and design capability to applied aerodynamicists. The system is comprised of loosely integrated, user-friendly software that enables the application of advanced Euler and Navier-Stokes tetrahedral finite volume technology to complex aerodynamic problems. TetrUSS has matured well because of the generous feedback from many willing users representing a broad cross-section of background and skill levels. This paper presents an overview of the current capabilities of the TetrUSS system along with some representative results from selected applications.

  15. Quadratic spatial soliton interactions

    NASA Astrophysics Data System (ADS)

    Jankovic, Ladislav

    Quadratic spatial soliton interactions were investigated in this Dissertation. The first part deals with characterizing the principal features of multi-soliton generation and soliton self-reflection. The second deals with two beam processes leading to soliton interactions and collisions. These subjects were investigated both theoretically and experimentally. The experiments were performed by using potassium niobate (KNBO 3) and periodically poled potassium titanyl phosphate (KTP) crystals. These particular crystals were desirable for these experiments because of their large nonlinear coefficients and, more importantly, because the experiments could be performed under non-critical-phase-matching (NCPM) conditions. The single soliton generation measurements, performed on KNBO3 by launching the fundamental component only, showed a broad angular acceptance bandwidth which was important for the soliton collisions performed later. Furthermore, at high input intensities multi-soliton generation was observed for the first time. The influence on the multi-soliton patterns generated of the input intensity and beam symmetry was investigated. The combined experimental and theoretical efforts indicated that spatial and temporal noise on the input laser beam induced multi-soliton patterns. Another research direction pursued was intensity dependent soliton routing by using of a specially engineered quadratically nonlinear interface within a periodically poled KTP sample. This was the first time demonstration of the self-reflection phenomenon in a system with a quadratic nonlinearity. The feature investigated is believed to have a great potential for soliton routing and manipulation by engineered structures. A detailed investigation was conducted on two soliton interaction and collision processes. Birth of an additional soliton resulting from a two soliton collision was observed and characterized for the special case of a non-planar geometry. A small amount of spiraling, up to 30

  16. Tetrahedral boron in naturally occurring tourmaline

    SciTech Connect

    Tagg, S.L.; Cho, H.; Dyar, M.D.; Grew, E.S.

    1999-09-01

    Evidence for boron in both trigonal and tetrahedral coordination has been found in {sup 11}B magic-angle-spinning (MAS) nuclear magnetic resonance (NMR) spectra of natural, inclusion-free specimens of aluminum-rich lithian tourmaline from granitic pregmatites.

  17. Quadratic soliton self-reflection at a quadratically nonlinear interface

    NASA Astrophysics Data System (ADS)

    Jankovic, Ladislav; Kim, Hongki; Stegeman, George; Carrasco, Silvia; Torner, Lluis; Katz, Mordechai

    2003-11-01

    The reflection of bulk quadratic solutions incident onto a quadratically nonlinear interface in periodically poled potassium titanyl phosphate was observed. The interface consisted of the boundary between two quasi-phase-matched regions displaced from each other by a half-period. At high intensities and small angles of incidence the soliton is reflected.

  18. On a combined adaptive tetrahedral tracing and edge diffraction model

    NASA Astrophysics Data System (ADS)

    Hart, Carl R.

    A major challenge in architectural acoustics is the unification of diffraction models and geometric acoustics. For example, geometric acoustics is insufficient to quantify the scattering characteristics of acoustic diffusors. Typically the time-independent boundary element method (BEM) is the method of choice. In contrast, time-domain computations are of interest for characterizing both the spatial and temporal scattering characteristics of acoustic diffusors. Hence, a method is sought that predicts acoustic scattering in the time-domain. A prediction method, which combines an advanced image source method and an edge diffraction model, is investigated for the prediction of time-domain scattering. Adaptive tetrahedral tracing is an advanced image source method that generates image sources through an adaptive process. Propagating tetrahedral beams adapt to ensonified geometry mapping the geometric sound field in space and along boundaries. The edge diffraction model interfaces with the adaptive tetrahedral tracing process by the transfer of edge geometry and visibility information. Scattering is quantified as the contribution of secondary sources along a single or multiple interacting edges. Accounting for a finite number of diffraction permutations approximates the scattered sound field. Superposition of the geometric and scattered sound fields results in a synthesized impulse response between a source and a receiver. Evaluation of the prediction technique involves numerical verification and numerical validation. Numerical verification is based upon a comparison with analytic and numerical (BEM) solutions for scattering geometries. Good agreement is shown for the selected scattering geometries. Numerical validation is based upon experimentally determined scattered impulse responses of acoustic diffusors. Experimental data suggests that the predictive model is appropriate for high-frequency predictions. For the experimental determination of the scattered impulse

  19. Students' Understanding of Quadratic Equations

    ERIC Educational Resources Information Center

    López, Jonathan; Robles, Izraim; Martínez-Planell, Rafael

    2016-01-01

    Action-Process-Object-Schema theory (APOS) was applied to study student understanding of quadratic equations in one variable. This required proposing a detailed conjecture (called a genetic decomposition) of mental constructions students may do to understand quadratic equations. The genetic decomposition which was proposed can contribute to help…

  20. Search for Tetrahedral Symmetry in 70Ge

    NASA Astrophysics Data System (ADS)

    Le, Khanh; Haring-Kaye, R. A.; Elder, R. M.; Jones, K. D.; Morrow, S. I.; Tabor, S. L.; Tripathi, V.; Bender, P. C.; Allegro, P. R. P.; Medina, N. H.; Oliveira, J. R. B.; Doring, J.

    2014-09-01

    The even-even Ge isotopes have recently become an active testing ground for a variety of exotic structural characteristics, including the existence of tetrahedral symmetry (pyramid-like shapes). Although theoretical shape calculations predict the onset of tetrahedral symmetry near 72Ge, the experimental signatures (including vanishing quadrupole moments within high-spin bands) remain elusive. This study searched for possible experimental evidence of tetrahedral symmetry in 70Ge. Excited states in 70Ge were populated at Florida State University using the 55Mn(18O,p2n) fusion-evaporation reaction at 50 MeV. Prompt γ- γ coincidences were measured with a Compton-suppressed Ge array consisting of three Clover detectors and seven single-crystal detectors. The existing level scheme was enhanced through the addition of 20 new transitions and the rearrangement of five others based on the measured coincidence relations and relative intensities. Lifetimes of 24 states were measured using the Doppler-shift attenuation method, from which transition quadrupole moments were inferred. These results will be compared with those obtained from cranked Woods-Saxon calculations. The even-even Ge isotopes have recently become an active testing ground for a variety of exotic structural characteristics, including the existence of tetrahedral symmetry (pyramid-like shapes). Although theoretical shape calculations predict the onset of tetrahedral symmetry near 72Ge, the experimental signatures (including vanishing quadrupole moments within high-spin bands) remain elusive. This study searched for possible experimental evidence of tetrahedral symmetry in 70Ge. Excited states in 70Ge were populated at Florida State University using the 55Mn(18O,p2n) fusion-evaporation reaction at 50 MeV. Prompt γ- γ coincidences were measured with a Compton-suppressed Ge array consisting of three Clover detectors and seven single-crystal detectors. The existing level scheme was enhanced through the addition

  1. Influence of the Rigid Connection Between Discs in the Tetrahedral Prisms on Equivalent Stresses When Cutting Work Faces

    NASA Astrophysics Data System (ADS)

    Khoreshok, A. A.; Mametyev, L. E.; Borisov, A. Yu; Vorobyev, A. V.

    2016-04-01

    This paper presents the results of modeling of the stressed state of structural elements of the paired fastening points of the two disc tools to the tetrahedral prisms of the working bodies of the roadheaders of selective action when cutting work faces of heterogeneous structure. The advantages of cooperative mode rotation to separate two disc tools on each of the tetrahedral prisms placed between the axial cutting crowns.

  2. Lattice Cleaving: Conforming Tetrahedral Meshes of Multimaterial Domains with Bounded Quality

    PubMed Central

    Bronson, Jonathan R.; Levine, Joshua A.; Whitaker, Ross T.

    2013-01-01

    Summary We introduce a new algorithm for generating tetrahedral meshes that conform to physical boundaries in volumetric domains consisting of multiple materials. The proposed method allows for an arbitrary number of materials, produces high-quality tetrahedral meshes with upper and lower bounds on dihedral angles, and guarantees geometric fidelity. Moreover, the method is combinatoric so its implementation enables rapid mesh construction. These meshes are structured in a way that also allows grading, in order to reduce element counts in regions of homogeneity. PMID:25309969

  3. Kinetically Trapped Tetrahedral Cages via Alkyne Metathesis.

    PubMed

    Lee, Semin; Yang, Anna; Moneypenny, Timothy P; Moore, Jeffrey S

    2016-02-24

    In dynamic covalent synthesis, kinetic traps are perceived as disadvantageous, hindering the system from reaching its thermodynamic equilibrium. Here we present the near-quantitative preparation of tetrahedral cages from simple tritopic precursors using alkyne metathesis. While the cages are the presumed thermodynamic sink, we experimentally demonstrate that the products no longer exchange their vertices once they have formed. The example reported here illustrates that kinetically trapped products may facilitate high yields of complex products from dynamic covalent synthesis. PMID:26854552

  4. Tetrahedrally coordinated carbonates in Earth's lower mantle.

    PubMed

    Boulard, Eglantine; Pan, Ding; Galli, Giulia; Liu, Zhenxian; Mao, Wendy L

    2015-01-01

    Carbonates are the main species that bring carbon deep into our planet through subduction. They are an important rock-forming mineral group, fundamentally distinct from silicates in the Earth's crust in that carbon binds to three oxygen atoms, while silicon is bonded to four oxygens. Here we present experimental evidence that under the sufficiently high pressures and high temperatures existing in the lower mantle, ferromagnesian carbonates transform to a phase with tetrahedrally coordinated carbons. Above 80 GPa, in situ synchrotron infrared experiments show the unequivocal spectroscopic signature of the high-pressure phase of (Mg,Fe)CO3. Using ab-initio calculations, we assign the new infrared signature to C-O bands associated with tetrahedrally coordinated carbon with asymmetric C-O bonds. Tetrahedrally coordinated carbonates are expected to exhibit substantially different reactivity than low-pressure threefold coordinated carbonates, as well as different chemical properties in the liquid state. Hence, this may have significant implications for carbon reservoirs and fluxes, and the global geodynamic carbon cycle. PMID:25692448

  5. Enhanced temperature uniformity by tetrahedral laser heating

    SciTech Connect

    Schroers, Jan; Bossuyt, Sven; Rhim, Won-Kyu; Li Jianzhong; Zhou Zhenhua; Johnson, William L.

    2004-11-01

    Temperature profile on a spherical sample that is heated by laser beams in various geometries while processed in vacuum is analyzed. Sample heating by one or four laser beams was considered. An analytical expression was derived for directional sample heating cases. It suggests an enhanced temperature uniformity over the samples when heated with four diffuse laser beams arranged in a tetrahedral geometry. This was experimentally verified by heating a spherical stainless steel sample by laser beams. Both the calculated and experimentally determined temperature variations over the sample suggest that use of diffuse four beams arranged in tetrahedral geometry would be effective in reducing temperature variation to within 1 K. The enhancement in the temperature uniformity for four diffuse beams arranged in a tetrahedral geometry by a factor of 50 over a single focused beam is promising to accurately measure of thermophysical properties. This drastic improvement in temperature uniformity might even enable atomic diffusion measurements in the undercooled liquid states of the bulk glass forming alloys since Marangoni and gravity driven convection will be substantially reduced.

  6. Methods for prismatic/tetrahedral grid generation and adaptation

    NASA Astrophysics Data System (ADS)

    Kallinderis, Y.

    1995-10-01

    The present work involves generation of hybrid prismatic/tetrahedral grids for complex 3-D geometries including multi-body domains. The prisms cover the region close to each body's surface, while tetrahedra are created elsewhere. Two developments are presented for hybrid grid generation around complex 3-D geometries. The first is a new octree/advancing front type of method for generation of the tetrahedra of the hybrid mesh. The main feature of the present advancing front tetrahedra generator that is different from previous such methods is that it does not require the creation of a background mesh by the user for the determination of the grid-spacing and stretching parameters. These are determined via an automatically generated octree. The second development is a method for treating the narrow gaps in between different bodies in a multiply-connected domain. This method is applied to a two-element wing case. A High Speed Civil Transport (HSCT) type of aircraft geometry is considered. The generated hybrid grid required only 170 K tetrahedra instead of an estimated two million had a tetrahedral mesh been used in the prisms region as well. A solution adaptive scheme for viscous computations on hybrid grids is also presented. A hybrid grid adaptation scheme that employs both h-refinement and redistribution strategies is developed to provide optimum meshes for viscous flow computations. Grid refinement is a dual adaptation scheme that couples 3-D, isotropic division of tetrahedra and 2-D, directional division of prisms.

  7. Boundary Recovery For Delaunay Tetrahedral Meshes Using Local Topological Transformations

    PubMed Central

    Ghadyani, Hamid; Sullivan, John; Wu, Ziji

    2009-01-01

    Numerous high-quality, volume mesh-generation systems exist. However, no strategy can address all geometry situations without some element qualities being compromised. Many 3D mesh generation algorithms are based on Delaunay tetrahedralization which frequently fails to preserve the input boundary surface topology. For biomedical applications, this surface preservation can be critical as they usually contain multiple material regions of interest coherently connected. In this paper we present an algorithm as a post-processing method that optimizes local regions of compromised element quality and recovers the original boundary surface facets (triangles) regardless of the original mesh generation strategy. The algorithm carves out a small sub-volume in the vicinity of the missing boundary facet or compromised element, creating a cavity. If the task is to recover a surface boundary facet, a natural exit hole in the cavity will be present. This hole is patched with the missing boundary surface face first followed by other patches to seal the cavity. If the task was to improve a compromised region, then the cavity is already sealed. Every triangular facet of the cavity shell is classified as an active face and can be connected to another shell node creating a tetrahedron. In the process the base of the tetrahedron is removed from the active face list and potentially 3 new active faces are created. This methodology is the underpinnings of our last resort method. Each active face can be viewed as the trunk of a tree. An exhaustive breath and depth search will identify all possible tetrahedral combinations to uniquely fill the cavity. We have streamlined this recursive process reducing the time complexity by orders of magnitude. The original surfaces boundaries (internal and external) are fully restored and the quality of compromised regions improved. PMID:20305743

  8. Small Power Technology for Tetrahedral Rovers

    NASA Astrophysics Data System (ADS)

    Clark, P. E.; Floyd, S. R.; Butler, C. D.; Flom, Y.

    2006-01-01

    The Small Power Technology (SPOT) being studied at GSFC has the potential to be an efficient and compact radioisotope based electrical power system. Such a system would provide power for innovative tetrahedral robotic arms and walkers to support the lunar exploration initiative within the next decade. Presently, NASA has designated two flight qualified Radioisotope Power Supplies (RPS): the Multi-Mission RTG (MMRTG) which uses thermocouple technology and the more efficient but more massive Stirling RTG (SRTG) which uses a mechanical heat (Stirling) engine technology. With SPOT, thermal output from a radioisotope source is converted to electrical power using a combination of shape memory material and piezoelectric crystals. The SPOT combined energy conversion technologies are potentially more efficient than thermocouples and do not require moving parts, thus keeping efficiency high with an excellent mass to power ratio. Applications of particular interest are highly modular, addressable, reconfigurable arrays of tetrahedral structural components designed to be arms or rovers with high mobility in rough terrain. Such prototypes are currently being built at GSFC. Missions requiring long-lived operation in unilluminated environments preclude the use of solar cells as the main power source and must rely on the use of RPS technology. The design concept calls for a small motor and battery assembly for each strut, and thus a distributed power system. We estimate, based on performance of our current tetrahedral prototypes and power scaling for small motors, that such devices require tens of watts of power output per kilogram of power supply. For these reasons, SPOT is a good candidate for the ART (addressable Reconfigurable Technology) baseline power system.

  9. Extreme Mobility: Next Generation Tetrahedral Rovers

    NASA Astrophysics Data System (ADS)

    Clark, P. E.; Curtis, S. A.; Rilee, M. L.; Cheung, C. Y.; Wesenberg, R.; Brown, G.; Cooperrider, C.

    2007-01-01

    This paper describes the development and testing of a patented rover concept called Tetrahedral Explorer Technologies (TET), designed to provide extreme mobility and plug-and-play utility through reconfigurable addressable architecture. Here, we present the results of preliminary lab and field tests of Prototype III. Reconfigurable architecture is essential in exploration because reaching features of the great potential interest will require crossing a wide range of terrains largely inaccessible to permanently appendaged vehicles. One surface might be relatively flat and navigable, while another could be rough, variably sloping, broken, or dominated by unconsolidated debris. To be totally functional, structures must form pseudo-appendages varying in size, rate, and manner of deployment (gait) and moving at a speed approaching that of a human in rugged terrain. TET architecture is based on the tetrahedron, the basic space-filling shape, as building block. Tetrahedra are interconnected, their apices acting as nodes from which struts reversibly deploy. The tetrahedral framework acts as a simple skeletal muscular structure. Two simple robotic walker prototypes have already been developed from a single reconfigurable tetrahedron capable of tumbling. This paper presents the results of our attempts to simulate motions, improve the hardware, and develop gaits for a more evolved 12Tetrahedral Walker (Prototype 3) which high degrees of freedom locomotion commandable through a user friendly interface. Our rover is an early level mission concept, realized as an electromechanical system at present, which would allow autonomous in situ exploration of lunar sites when we return to the Moon. Such a rover could carry into inaccessible terrain an in situ analysis payload designed to provide not only details of composition of traversed terrain, but the identification of sites with resources useful for permanent bases, including water and high Ti glass.

  10. Dark Matter from Binary Tetrahedral Flavor Symmetry

    NASA Astrophysics Data System (ADS)

    Eby, David; Frampton, Paul

    2012-03-01

    Binary Tetrahedral Flavor Symmetry, originally developed as a quark family symmetry and later adapted to leptons, has proved both resilient and versatile over the past decade. In 2008 a minimal T' model was developed to accommodate quark and lepton masses and mixings using a family symmetry of (T'xZ2). We examine an expansion of this earlier model using an additional Z2 group that facilitates predictions of WIMP dark matter, the Cabibbo angle, and deviations from Tribimaximal Mixing, while giving hints at the nature of leptogenesis.

  11. A quadratic analog-to-digital converter

    NASA Technical Reports Server (NTRS)

    Harrison, D. C.; Staples, M. H.

    1980-01-01

    An analog-to-digital converter with a square root transfer function has been developed for use with a pair of CCD imaging detectors in the White Light Coronagraph/X-ray XUV Telescope experiment to be flown as part of the Internal Solar Polar Mission. It is shown that in background-noise-limited instrumentation systems a quadratic analog-to-digital converter will allow a maximum dynamic range with a fixed number of data bits. Low power dissipation, moderately fast conversion time, and reliability are achieved in the proposed design using standard components and avoiding nonlinear elements.

  12. Students' understanding of quadratic equations

    NASA Astrophysics Data System (ADS)

    López, Jonathan; Robles, Izraim; Martínez-Planell, Rafael

    2016-05-01

    Action-Process-Object-Schema theory (APOS) was applied to study student understanding of quadratic equations in one variable. This required proposing a detailed conjecture (called a genetic decomposition) of mental constructions students may do to understand quadratic equations. The genetic decomposition which was proposed can contribute to help students achieve an understanding of quadratic equations with improved interrelation of ideas and more flexible application of solution methods. Semi-structured interviews with eight beginning undergraduate students explored which of the mental constructions conjectured in the genetic decomposition students could do, and which they had difficulty doing. Two of the mental constructions that form part of the genetic decomposition are highlighted and corresponding further data were obtained from the written work of 121 undergraduate science and engineering students taking a multivariable calculus course. The results suggest the importance of explicitly considering these two highlighted mental constructions.

  13. A bicontinuous tetrahedral structure in a liquid-crystalline lipid

    NASA Astrophysics Data System (ADS)

    Longley, William; McIntosh, Thomas J.

    1983-06-01

    The structure of most lipid-water phases can be visualized as an ordered distribution of two liquid media, water and hydrocarbons, separated by a continuous surface covered by the polar groups of the lipid molecules1. In the cubic phases in particular, rod-like elements are linked into three-dimensional networks1,2. Two of these phases (space groups Ia3d and Pn3m) contain two such three-dimensional networks mutually inter-woven and unconnected. Under the constraints of energy minimization3, the interface between the components in certain of these `porous fluids' may well resemble one of the periodic minimal surface structures of the type described mathematically by Schwarz4,5. A structure of this sort has been proposed for the viscous isotropic (cubic) form of glycerol monooleate (GMO) by Larsson et al.6 who suggested that the X-ray diagrams of Lindblom et al.7 indicated a body-centred crystal structure in which lipid bilayers might be arranged as in Schwarz's octahedral surface4. We have now found that at high water contents, a primitive cubic lattice better fits the X-ray evidence with the material in the crystal arranged in a tetrahedral way. The lipid appears to form a single bilayer, continuous in three dimensions, separating two continuous interlinked networks of water. Each of the water networks has the symmetry of the diamond crystal structure and the bilayer lies in the space between them following a surface resembling Schwarz's tetrahedral surface4.

  14. Phase diagram of a truncated tetrahedral model.

    PubMed

    Krcmar, Roman; Gendiar, Andrej; Nishino, Tomotoshi

    2016-08-01

    Phase diagram of a discrete counterpart of the classical Heisenberg model, the truncated tetrahedral model, is analyzed on the square lattice, when the interaction is ferromagnetic. Each spin is represented by a unit vector that can point to one of the 12 vertices of the truncated tetrahedron, which is a continuous interpolation between the tetrahedron and the octahedron. Phase diagram of the model is determined by means of the statistical analog of the entanglement entropy, which is numerically calculated by the corner transfer matrix renormalization group method. The obtained phase diagram consists of four different phases, which are separated by five transition lines. In the parameter region, where the octahedral anisotropy is dominant, a weak first-order phase transition is observed. PMID:27627273

  15. Artificial intelligence approach to planning the robotic assembly of large tetrahedral truss structures

    NASA Technical Reports Server (NTRS)

    Homemdemello, Luiz S.

    1992-01-01

    An assembly planner for tetrahedral truss structures is presented. To overcome the difficulties due to the large number of parts, the planner exploits the simplicity and uniformity of the shapes of the parts and the regularity of their interconnection. The planning automation is based on the computational formalism known as production system. The global data base consists of a hexagonal grid representation of the truss structure. This representation captures the regularity of tetrahedral truss structures and their multiple hierarchies. It maps into quadratic grids and can be implemented in a computer by using a two-dimensional array data structure. By maintaining the multiple hierarchies explicitly in the model, the choice of a particular hierarchy is only made when needed, thus allowing a more informed decision. Furthermore, testing the preconditions of the production rules is simple because the patterned way in which the struts are interconnected is incorporated into the topology of the hexagonal grid. A directed graph representation of assembly sequences allows the use of both graph search and backtracking control strategies.

  16. Updates to Multi-Dimensional Flux Reconstruction for Hypersonic Simulations on Tetrahedral Grids

    NASA Technical Reports Server (NTRS)

    Gnoffo, Peter A.

    2010-01-01

    The quality of simulated hypersonic stagnation region heating with tetrahedral meshes is investigated by using an updated three-dimensional, upwind reconstruction algorithm for the inviscid flux vector. An earlier implementation of this algorithm provided improved symmetry characteristics on tetrahedral grids compared to conventional reconstruction methods. The original formulation however displayed quantitative differences in heating and shear that were as large as 25% compared to a benchmark, structured-grid solution. The primary cause of this discrepancy is found to be an inherent inconsistency in the formulation of the flux limiter. The inconsistency is removed by employing a Green-Gauss formulation of primitive gradients at nodes to replace the previous Gram-Schmidt algorithm. Current results are now in good agreement with benchmark solutions for two challenge problems: (1) hypersonic flow over a three-dimensional cylindrical section with special attention to the uniformity of the solution in the spanwise direction and (2) hypersonic flow over a three-dimensional sphere. The tetrahedral cells used in the simulation are derived from a structured grid where cell faces are bisected across the diagonal resulting in a consistent pattern of diagonals running in a biased direction across the otherwise symmetric domain. This grid is known to accentuate problems in both shock capturing and stagnation region heating encountered with conventional, quasi-one-dimensional inviscid flux reconstruction algorithms. Therefore the test problems provide a sensitive indicator for algorithmic effects on heating. Additional simulations on a sharp, double cone and the shuttle orbiter are then presented to demonstrate the capabilities of the new algorithm on more geometrically complex flows with tetrahedral grids. These results provide the first indication that pure tetrahedral elements utilizing the updated, three-dimensional, upwind reconstruction algorithm may be used for the

  17. Linear state feedback, quadratic weights, and closed loop eigenstructures. M.S. Thesis

    NASA Technical Reports Server (NTRS)

    Thompson, P. M.

    1979-01-01

    Results are given on the relationships between closed loop eigenstructures, state feedback gain matrices of the linear state feedback problem, and quadratic weights of the linear quadratic regulator. Equations are derived for the angles of general multivariable root loci and linear quadratic optimal root loci, including angles of departure and approach. The generalized eigenvalue problem is used for the first time to compute angles of approach. Equations are also derived to find the sensitivity of closed loop eigenvalues and the directional derivatives of closed loop eigenvectors (with respect to a scalar multiplying the feedback gain matrix or the quadratic control weight). An equivalence class of quadratic weights that produce the same asymptotic eigenstructure is defined, sufficient conditions to be in it are given, a canonical element is defined, and an algorithm to find it is given. The behavior of the optimal root locus in the nonasymptotic region is shown to be different for quadratic weights with the same asymptotic properties.

  18. More About the Tetrahedral Unstructured Software System

    NASA Technical Reports Server (NTRS)

    Abdol-Hamid, Khaled S.; Frink, Neal T.; Hunter, Craig A.; Parikh, Paresh C.; Pizadeh, Shalyar Z.; Samareh, Jamshid A.; Bhat, Maharaj K.; Pandya, Mohagna J.; Grismer, Matthew J.

    2006-01-01

    TetrUSS is a comprehensive suite of computational fluid dynamics (CFD) programs that won the Software of the Year award in 1996 and has found increasing use in government, academia, and industry for solving realistic flow problems (especially in aerodynamics and aeroelastics of aircraft having complex shapes). TetrUSS includes not only programs for solving basic equations of flow but also programs that afford capabilities for efficient generation and utilization of computational grids and for graphical representation of computed flows (see figure). The 2004 version of the Tetrahedral Unstructured Software System (TetrUSS), which is one of two software systems reported in "NASA s 2004 Software of the Year," NASA Tech Briefs, Vol. 28, No. 10 (October 2004), page 18, has been improved greatly since 1996. These improvements include (1) capabilities to simulate viscous flow by solving the Navier-Stokes equations on unstructured grids, (2) portability to personal computers from diverse manufacturers, (3) advanced models of turbulence, (4) a parallel-processing version of one of the unstructured-grid Navier-Stokes-equation-solving programs, and (5) advanced programs for generating unstructured grids.

  19. Streaming Compression of Tetrahedral Volume Meshes

    SciTech Connect

    Isenburg, M; Lindstrom, P; Gumhold, S; Shewchuk, J

    2005-11-21

    Geometry processing algorithms have traditionally assumed that the input data is entirely in main memory and available for random access. This assumption does not scale to large data sets, as exhausting the physical memory typically leads to IO-inefficient thrashing. Recent works advocate processing geometry in a 'streaming' manner, where computation and output begin as soon as possible. Streaming is suitable for tasks that require only local neighbor information and batch process an entire data set. We describe a streaming compression scheme for tetrahedral volume meshes that encodes vertices and tetrahedra in the order they are written. To keep the memory footprint low, the compressor is informed when vertices are referenced for the last time (i.e. are finalized). The compression achieved depends on how coherent the input order is and how many tetrahedra are buffered for local reordering. For reasonably coherent orderings and a buffer of 10,000 tetrahedra, we achieve compression rates that are only 25 to 40 percent above the state-of-the-art, while requiring drastically less memory resources and less than half the processing time.

  20. Quadratic expressions by means of `summing all the matchsticks'

    NASA Astrophysics Data System (ADS)

    Faaiz Gierdien, M.

    2012-09-01

    This note presents demonstrations of quadratic expressions that come about when particular problems are posed with respect to matchsticks that form regular triangles, squares, pentagons and so on. Usually when such 'matchstick' problems are used as ways to foster algebraic thinking, the expressions for the number of matchstick quantities are linear and not quadratic. It will be shown that a pedagogy of 'summing all the matchsticks' is central to the emergence of quadratic expressions. This pedagogy involves generational and transformational activities which are considered as some of the main activities of algebra. Key elements to these activities are processes such as recognizing and extending patterns, and specializing and generalizing particular functional relationships. Implications of these processes in terms of algebraic thinking are considered.

  1. Orthogonality preserving infinite dimensional quadratic stochastic operators

    SciTech Connect

    Akın, Hasan; Mukhamedov, Farrukh

    2015-09-18

    In the present paper, we consider a notion of orthogonal preserving nonlinear operators. We introduce π-Volterra quadratic operators finite and infinite dimensional settings. It is proved that any orthogonal preserving quadratic operator on finite dimensional simplex is π-Volterra quadratic operator. In infinite dimensional setting, we describe all π-Volterra operators in terms orthogonal preserving operators.

  2. Preliminary design of a large tetrahedral truss/hexagonal heatshield panel aerobrake

    NASA Technical Reports Server (NTRS)

    Dorsey, John T.; Mikulas, Martin M., Jr.

    1989-01-01

    An aerobrake structural concept is introduced which consists of two primary components: (1) a lightweight erectable tetrahedral support truss; and (2) sandwich hexagonal heatshield panels which, when attached to the truss, form a continuous impermeable aerobraking surface. Generic finite element models and a general analysis procedure to design tetrahedral truss/hexagonal heatshield panel aerobrakes is developed, and values of the aerobrake design parameters which minimize mass and packaging volume for a 120-foot-diameter aerobrake are determined. Sensitivity of the aerobrake design to variations in design parameters is also assessed. The results show that a 120-foot-diameter aerobrake is viable using the concept presented (i.e., the aerobrake mass is less than or equal to 15 percent of the payload spacecraft mass). Minimizing the aerobrake mass (by increasing the number of rings in the support truss) however, leads to aerobrakes with the highest part count.

  3. Adaptive hybrid prismatic-tetrahedral grids for viscous flows

    NASA Technical Reports Server (NTRS)

    Kallinderis, Yannis; Khawaja, Aly; Mcmorris, Harlan

    1995-01-01

    The paper presents generation of adaptive hybrid prismatic/tetrahedral grids for complex 3-D geometries including multi-body domains. The prisms cover the region close to each body's surface, while tetrahedra are created elsewhere. Two developments are presented for hybrid grid generation around complex 3-D geometries. The first is a new octree/advancing front type of method for generation of the tetrahedra of the hybrid mesh. The main feature of the present advancing front tetrahedra generator that is different from previous such methods is that it does not require the creation of a background mesh by the user for the determination of the grid-spacing and stretching parameters. These are determined via an automatically generated octree. The second development is an Automatic Receding Method (ARM) for treating the narrow gaps in between different bodies in a multiply-connected domain. This method is applied to a two-element wing case. A hybrid grid adaptation scheme that employs both h-refinement and redistribution strategies is developed to provide optimum meshes for viscous flow computations. Grid refinement is a dual adaptation scheme that couples division of tetrahedra, as well as 2-D directional division of prisms.

  4. Adaptive hybrid prismatic-tetrahedral grids for viscous flows

    NASA Astrophysics Data System (ADS)

    Kallinderis, Yannis; Khawaja, Aly; McMorris, Harlan

    1995-03-01

    The paper presents generation of adaptive hybrid prismatic/tetrahedral grids for complex 3-D geometries including multi-body domains. The prisms cover the region close to each body's surface, while tetrahedra are created elsewhere. Two developments are presented for hybrid grid generation around complex 3-D geometries. The first is a new octree/advancing front type of method for generation of the tetrahedra of the hybrid mesh. The main feature of the present advancing front tetrahedra generator that is different from previous such methods is that it does not require the creation of a background mesh by the user for the determination of the grid-spacing and stretching parameters. These are determined via an automatically generated octree. The second development is an Automatic Receding Method (ARM) for treating the narrow gaps in between different bodies in a multiply-connected domain. This method is applied to a two-element wing case. A hybrid grid adaptation scheme that employs both h-refinement and redistribution strategies is developed to provide optimum meshes for viscous flow computations. Grid refinement is a dual adaptation scheme that couples division of tetrahedra, as well as 2-D directional division of prisms.

  5. Dynamic Modulation of DNA Hybridization Using Allosteric DNA Tetrahedral Nanostructures.

    PubMed

    Song, Ping; Li, Min; Shen, Juwen; Pei, Hao; Chao, Jie; Su, Shao; Aldalbahi, Ali; Wang, Lihua; Shi, Jiye; Song, Shiping; Wang, Lianhui; Fan, Chunhai; Zuo, Xiaolei

    2016-08-16

    The fixed dynamic range of traditional biosensors limits their utility in several real applications. For example, viral load monitoring requires the dynamic range spans several orders of magnitude; whereas, monitoring of drugs requires extremely narrow dynamic range. To overcome this limitation, here, we devised tunable biosensing interface using allosteric DNA tetrahedral bioprobes to tune the dynamic range of DNA biosensors. Our strategy takes the advantage of the readily and flexible structure design and predictable geometric reconfiguration of DNA nanotechnology. We reconfigured the DNA tetrahedral bioprobes by inserting the effector sequence into the DNA tetrahedron, through which, the binding affinity of DNA tetrahedral bioprobes can be tuned. As a result, the detection limit of DNA biosensors can be programmably regulated. The dynamic range of DNA biosensors can be tuned (narrowed or extended) for up to 100-fold. Using the regulation of binding affinity, we realized the capture and release of biomolecules by tuning the binding behavior of DNA tetrahedral bioprobes. PMID:27435955

  6. A boundary recovery algorithm for Delaunay tetrahedral meshing

    SciTech Connect

    Sharov, D.; Nakahashi, K.

    1996-12-31

    A method for automatic generation of unstructured grids comprised of tetrahedra is discussed. Delaunay approach for tetrahedral grid generation is used. Particular attention is given to the boundary constraining problem. A simple and robust algorithm for the boundary constraining by successive use of boundary edge swapping, tetrahedral edge swapping and direct subdivision of tetrahedra is used. Small modifications allow to apply the method for viscous grid generation as well. Grid examples demonstrate efficiency of the method.

  7. Coherent states for quadratic Hamiltonians

    NASA Astrophysics Data System (ADS)

    Contreras-Astorga, Alonso; Fernández C, David J.; Velázquez, Mercedes

    2011-01-01

    The coherent states for a set of quadratic Hamiltonians in the trap regime are constructed. A matrix technique which allows us to directly identify the creation and annihilation operators will be presented. Then, the coherent states as simultaneous eigenstates of the annihilation operators will be derived, and will be compared with those attained through the displacement operator method. The corresponding wavefunction will be found, and a general procedure for obtaining several mean values involving the canonical operators in these states will be described. The results will be illustrated through the asymmetric Penning trap.

  8. Geometrical and Graphical Solutions of Quadratic Equations.

    ERIC Educational Resources Information Center

    Hornsby, E. John, Jr.

    1990-01-01

    Presented are several geometrical and graphical methods of solving quadratic equations. Discussed are Greek origins, Carlyle's method, von Staudt's method, fixed graph methods and imaginary solutions. (CW)

  9. Single-photon quadratic optomechanics

    PubMed Central

    Liao, Jie-Qiao; Nori, Franco

    2014-01-01

    We present exact analytical solutions to study the coherent interaction between a single photon and the mechanical motion of a membrane in quadratic optomechanics. We consider single-photon emission and scattering when the photon is initially inside the cavity and in the fields outside the cavity, respectively. Using our solutions, we calculate the single-photon emission and scattering spectra, and find relations between the spectral features and the system's inherent parameters, such as: the optomechanical coupling strength, the mechanical frequency, and the cavity-field decay rate. In particular, we clarify the conditions for the phonon sidebands to be visible. We also study the photon-phonon entanglement for the long-time emission and scattering states. The linear entropy is employed to characterize this entanglement by treating it as a bipartite one between a single mode of phonons and a single photon. PMID:25200128

  10. THE GENERATION OF TETRAHEDRAL MESH MODELS FOR NEUROANATOMICAL MRI

    PubMed Central

    Lederman, Carl; Joshi, Anand; Dinov, Ivo; Vese, Luminita; Toga, Arthur; Van Horn, John Darrell

    2010-01-01

    In this article, we describe a detailed method for automatically generating tetrahedral meshes from 3D images having multiple region labels. An adaptively sized tetrahedral mesh modeling approach is described that is capable of producing meshes conforming precisely to the voxelized regions in the image. Efficient tetrahedral construction is performed minimizing an energy function containing three terms: a smoothing term to remove the voxelization, a fidelity term to maintain continuity with the image data, and a novel elasticity term to prevent the tetrahedra from becoming flattened or inverted as the mesh deforms while allowing the voxelization to be removed entirely. The meshing algorithm is applied to structural MR image data that has been automatically segmented into 56 neuroanatomical sub-divisions as well as on two other examples. The resulting tetrahedral representation has several desirable properties such as tetrahedra with dihedral angles away from 0 and 180 degrees, smoothness, and a high resolution. Tetrahedral modeling via the approach described here has applications in modeling brain structure in normal as well as diseased brain in human and non-human data and facilitates examination of 3D object deformations resulting from neurological illness (e.g. Alzheimer’s Disease), development, and/or aging. PMID:21073968

  11. A comparison of tetrahedral mesh improvement techniques

    SciTech Connect

    Freitag, L.A.; Ollivier-Gooch, C.

    1996-12-01

    Automatic mesh generation and adaptive refinement methods for complex three-dimensional domains have proven to be very successful tools for the efficient solution of complex applications problems. These methods can, however, produce poorly shaped elements that cause the numerical solution to be less accurate and more difficult to compute. Fortunately, the shape of the elements can be improved through several mechanisms, including face-swapping techniques that change local connectivity and optimization-based mesh smoothing methods that adjust grid point location. The authors consider several criteria for each of these two methods and compare the quality of several meshes obtained by using different combinations of swapping and smoothing. Computational experiments show that swapping is critical to the improvement of general mesh quality and that optimization-based smoothing is highly effective in eliminating very small and very large angles. The highest quality meshes are obtained by using a combination of swapping and smoothing techniques.

  12. Binary Quadratic Forms: A Historical View

    ERIC Educational Resources Information Center

    Khosravani, Azar N.; Beintema, Mark B.

    2006-01-01

    We present an expository account of the development of the theory of binary quadratic forms. Beginning with the formulation and proof of the Two-Square Theorem, we show how the study of forms of the type x[squared] + ny[squared] led to the discovery of the Quadratic Reciprocity Law, and how this theorem, along with the concept of reduction relates…

  13. An Unexpected Influence on a Quadratic

    ERIC Educational Resources Information Center

    Davis, Jon D.

    2013-01-01

    Using technology to explore the coefficients of a quadratic equation can lead to an unexpected result. This article describes an investigation that involves sliders and dynamically linked representations. It guides students to notice the effect that the parameter "a" has on the graphical representation of a quadratic function in the form…

  14. Factorising a Quadratic Expression with Geometric Insights

    ERIC Educational Resources Information Center

    Joarder, Anwar H.

    2015-01-01

    An algorithm is presented for factorising a quadratic expression to facilitate instruction and learning. It appeals to elementary geometry which may provide better insights to some students or teachers. There have been many methods for factorising a quadratic expression described in school text books. However, students often seem to struggle with…

  15. Molecular origin of auxetic behavior in tetrahedral framework silicates.

    PubMed

    Alderson, Andrew; Evans, Kenneth E

    2002-11-25

    Recent analytical models for the Poisson's ratios (nu(ij)) of tetrahedral frameworks are applied to alpha-cristobalite and alpha-quartz for the first time. Rotation and dilation of the SiO4 tetrahedral subunits are considered. Each mechanism leads to negative nu(31) values, whereas negative and positive values are possible when they act concurrently. The concurrent model is in excellent agreement with experiment and explains the dichotomy between negative and positive nu(31) values in alpha-cristobalite and alpha-quartz, respectively. The predicted strain-dependent trends confirm those from molecular modeling. PMID:12485081

  16. The response of cranial biomechanical finite element models to variations in mesh density.

    PubMed

    Bright, Jen A; Rayfield, Emily J

    2011-04-01

    Finite element (FE) models provide discrete solutions to continuous problems. Therefore, to arrive at the correct solution, it is vital to ensure that FE models contain a sufficient number of elements to fully resolve all the detail encountered in a continuum structure. Mesh convergence testing is the process of comparing successively finer meshes to identify the point of diminishing returns; where increasing resolution has marginal effects on results and further detail would become costly and unnecessary. Historically, convergence has not been considered in most CT-based biomechanical reconstructions involving complex geometries like the skull, as generating such models has been prohibitively time-consuming. To assess how mesh convergence influences results, 18 increasingly refined CT-based models of a domestic pig skull were compared to identify the point of convergence for strain and displacement, using both linear and quadratic tetrahedral elements. Not all regions of the skull converged at the same rate, and unexpectedly, areas of high strain converged faster than low-strain regions. Linear models were slightly stiffer than their quadratic counterparts, but did not converge less rapidly. As expected, insufficiently dense models underestimated strain and displacement, and failed to resolve strain "hot-spots" notable in contour plots. In addition to quantitative differences, visual assessments of such plots often inform conclusions drawn in many comparative studies, highlighting that mesh convergence should be performed on all finite element models before further analysis takes place. PMID:21370496

  17. L -functions and class numbers of imaginary quadratic fields and of quadratic extensions of an imaginary quadratic field

    NASA Astrophysics Data System (ADS)

    Louboutin, Stephane

    1992-07-01

    Starting from the analytic class number formula involving its L-function, we first give an expression for the class number of an imaginary quadratic field which, in the case of large discriminants, provides us with a much more powerful numerical technique than that of counting the number of reduced definite positive binary quadratic forms, as has been used by Buell in order to compute his class number tables. Then, using class field theory, we will construct a periodic character &chi , defined on the ring of integers of a field K that is a quadratic extension of a principal imaginary quadratic field k, such that the zeta function of K is the product of the zeta function of k and of the L-function L(s,χ) . We will then determine an integral representation of this L-function that enables us to calculate the class number of K numerically, as soon as its regulator is known. It will also provide us with an upper bound for these class numbers, showing that Hua's bound for the class numbers of imaginary and real quadratic fields is not the best that one could expect. We give statistical results concerning the class numbers of the first 50000 quadratic extensions of {Q}(i) with prime relative discriminant (and with K/Q a non-Galois quartic extension). Our analytic calculation improves the algebraic calculation used by Lakein in the same way as the analytic calculation of the class numbers of real quadratic fields made by Williams and Broere improved the algebraic calculation consisting in counting the number of cycles of reduced ideals. Finally, we give upper bounds for class numbers of K that is a quadratic extension of an imaginary quadratic field k which is no longer assumed to be of class number one.

  18. Tetrahedral Models of Learning: Application to College Reading.

    ERIC Educational Resources Information Center

    Nist, Sherrie L.

    J. D. Bransford's tetrahedral model of learning considers four variables: (1) learning activities, (2) characteristics of the learner, (3) criterial tasks, and (4) the nature of the materials. Bransford's model provides a research-based theoretical framework that can be used to teach, model, and have students apply a variety of study strategies to…

  19. Selective refinement queries for volume visualization of unstructured tetrahedral meshes.

    PubMed

    Cignoni, Paolo; De Floriani, Leila; Magillo, Paola; Puppo, Enrico; Scopigno, Roberto

    2004-01-01

    In this paper, we address the problem of the efficient visualization of large irregular volume data sets by exploiting a multiresolution model based on tetrahedral meshes. Multiresolution models, also called Level-Of-Detail (LOD) models, allow encoding the whole data set at a virtually continuous range of different resolutions. We have identified a set of queries for extracting meshes at variable resolution from a multiresolution model, based on field values, domain location, or opacity of the transfer function. Such queries allow trading off between resolution and speed in visualization. We define a new compact data structure for encoding a multiresolution tetrahedral mesh built through edge collapses to support selective refinement efficiently and show that such a structure has a storage cost from 3 to 5.5 times lower than standard data structures used for tetrahedral meshes. The data structures and variable resolution queries have been implemented together with state-of-the art visualization techniques in a system for the interactive visualization of three-dimensional scalar fields defined on tetrahedral meshes. Experimental results show that selective refinement queries can support interactive visualization of large data sets. PMID:15382696

  20. Hinge specification for a square-faceted tetrahedral truss

    NASA Technical Reports Server (NTRS)

    Adams, L. R.

    1984-01-01

    A square-faceted tetrahedral truss is geometrically analyzed. Expressions are developed for single degree of freedom hinges which allow packaging of the structure into a configuration in which all members are parallel and closely packed in a square pattern. Deployment is sequential, thus providing control over the structure during deployment.

  1. Quadratic boundedness of uncertain nonlinear dynamic systems

    NASA Astrophysics Data System (ADS)

    Brockman, Mark Lawrence

    Physical systems are often perturbed by unknown external disturbances or contain important system parameters which are difficult to model exactly. However, engineers are expected to design systems which perform well even in the presence of uncertainties. For example, an airplane designer can never know the precise direction or magnitude of wind gusts, or the exact mass distribution inside the aircraft, but passengers expect to arrive on time after a smooth ride. This thesis will first present the concept of quadratic boundedness of an uncertain nonlinear dynamic system, and then develop analysis techniques and control design methods for systems containing unknown disturbances and parameters. For a class of nonlinear systems, conditions for quadratic boundedness are given, and the relationship between quadratic boundedness and quadratic stability is explored. An important consequence of quadratic boundedness is the ability to calculate an upper bound on the system gain of an uncertain nonlinear system. For nominally linear systems, necessary and sufficient conditions for quadratic boundedness are given. The innovative use of linear matrix inequalities in an iterative algorithm provides a means to analyze the quadratic boundedness properties of systems containing parameter uncertainties. The analysis results establish a framework for the development of design methods which integrate performance specifications into the control design process for all the types of systems considered. Numerous examples illustrate the major results of the thesis.

  2. Maintaining tetrahedral mesh quality in response to time-dependent topological and geometrical deformation

    SciTech Connect

    Kuprat, A.; George, D.

    1998-12-01

    When modeling deformation of geometrically complex regions, unstructured tetrahedral meshes provide the flexibility necessary to track interfaces as they change geometrically and topologically. In the class of time-dependent simulations considered in this paper, multimaterial interfaces are represented by sets of triangular facets, and motion of the interfaces is controlled by physical considerations. The motion of interior points in the conforming tetrahedral mesh (i.e., points not on interfaces) is arbitrary and may be chosen to produce good element shapes. In the context of specified boundary motion driven by physical considerations, they have found that a rather large glossary of mesh changes is required to allow the simulation to survive all the transitions of interface geometry and topology that occur during time evolution. This paper will describe mesh changes required to maintain good element quality as the geometry evolves, as well as mesh changes required to capture changes i n topology that occur when material regions collapse or pinch off. This paper will present a detailed description of mesh changes necessary for capturing the aforementioned geometrical and topological changes, as implemented in the code GRAIN3D, and will provide examples from a metallic grain growth simulation in which the normal velocity of the grain boundary is proportional to mean curvature.

  3. Seven Wonders of the Ancient and Modern Quadratic World.

    ERIC Educational Resources Information Center

    Taylor, Sharon E.; Mittag, Kathleen Cage

    2001-01-01

    Presents four methods for solving a quadratic equation using graphing calculator technology: (1) graphing with the CALC feature; (2) quadratic formula program; (3) table; and (4) solver. Includes a worksheet for a lab activity on factoring quadratic equations. (KHR)

  4. Tetrahedral-Mesh Simulation of Turbulent Flows with the Space-Time Conservative Schemes

    NASA Technical Reports Server (NTRS)

    Chang, Chau-Lyan; Venkatachari, Balaji; Cheng, Gary C.

    2015-01-01

    Direct numerical simulations of turbulent flows are predominantly carried out using structured, hexahedral meshes despite decades of development in unstructured mesh methods. Tetrahedral meshes offer ease of mesh generation around complex geometries and the potential of an orientation free grid that would provide un-biased small-scale dissipation and more accurate intermediate scale solutions. However, due to the lack of consistent multi-dimensional numerical formulations in conventional schemes for triangular and tetrahedral meshes at the cell interfaces, numerical issues exist when flow discontinuities or stagnation regions are present. The space-time conservative conservation element solution element (CESE) method - due to its Riemann-solver-free shock capturing capabilities, non-dissipative baseline schemes, and flux conservation in time as well as space - has the potential to more accurately simulate turbulent flows using unstructured tetrahedral meshes. To pave the way towards accurate simulation of shock/turbulent boundary-layer interaction, a series of wave and shock interaction benchmark problems that increase in complexity, are computed in this paper with triangular/tetrahedral meshes. Preliminary computations for the normal shock/turbulence interactions are carried out with a relatively coarse mesh, by direct numerical simulations standards, in order to assess other effects such as boundary conditions and the necessity of a buffer domain. The results indicate that qualitative agreement with previous studies can be obtained for flows where, strong shocks co-exist along with unsteady waves that display a broad range of scales, with a relatively compact computational domain and less stringent requirements for grid clustering near the shock. With the space-time conservation properties, stable solutions without any spurious wave reflections can be obtained without a need for buffer domains near the outflow/farfield boundaries. Computational results for the

  5. Quadratic Stochastic Operators with Countable State Space

    NASA Astrophysics Data System (ADS)

    Ganikhodjaev, Nasir

    2016-03-01

    In this paper, we provide the classes of Poisson and Geometric quadratic stochastic operators with countable state space, study the dynamics of these operators and discuss their application to economics.

  6. Schur Stability Regions for Complex Quadratic Polynomials

    ERIC Educational Resources Information Center

    Cheng, Sui Sun; Huang, Shao Yuan

    2010-01-01

    Given a quadratic polynomial with complex coefficients, necessary and sufficient conditions are found in terms of the coefficients such that all its roots have absolute values less than 1. (Contains 3 figures.)

  7. Quantum integrability of quadratic Killing tensors

    SciTech Connect

    Duval, C.; Valent, G.

    2005-05-01

    Quantum integrability of classical integrable systems given by quadratic Killing tensors on curved configuration spaces is investigated. It is proven that, using a 'minimal' quantization scheme, quantum integrability is ensured for a large class of classic examples.

  8. Weight of quadratic forms and graph states

    NASA Astrophysics Data System (ADS)

    Cosentino, Alessandro; Severini, Simone

    2009-11-01

    We prove a connection between Schmidt rank and weight of quadratic forms. This provides a new tool for the classification of graph states based on entanglement. Our main tool arises from a reformulation of previously known results concerning the weight of quadratic forms in terms of graph states properties. As a byproduct, we obtain a straightforward characterization of the weight of functions associated with pivot-minor of bipartite graphs.

  9. Quadratic mutual information for dimensionality reduction and classification

    NASA Astrophysics Data System (ADS)

    Gray, David M.; Principe, José C.

    2010-04-01

    A research area based on the application of information theory to machine learning has attracted considerable interest in the last few years. This research area has been coined information-theoretic learning within the community. In this paper we apply elements of information-theoretic learning to the problem of automatic target recognition (ATR). A number of researchers have previously shown the benefits of designing classifiers based on maximizing the mutual information between the class data and the class labels. Following prior research in information-theoretic learning, in the current results we show that quadratic mutual information, derived using a special case of the more general Renyi's entropy, can be used for classifier design. In this implementation, a simple subspace projection classifier is formulated to find the optimal projection weights such that the quadratic mutual information between the class data and the class labels is maximized. This subspace projection accomplishes a dimensionality reduction of the raw data set wherein information about the class membership is retained while irrelevant information is discarded. A subspace projection based on this criterion preserves as much class discriminability as possible within the subspace. For this paper, laser radar images are used to demonstrate the results. Classification performance against this data set is compared for a gradient descent MLP classifier and a quadratic mutual information MLP classifier.

  10. Structural stiffness, strength and dynamic characteristics of large tetrahedral space truss structures

    NASA Technical Reports Server (NTRS)

    Mikulas, M. M., Jr.; Bush, H. G.; Card, M. F.

    1977-01-01

    Physical characteristics of large skeletal frameworks for space applications are investigated by analyzing one concept: the tetrahedral truss, which is idealized as a sandwich plate with isotropic faces. Appropriate analytical relations are presented in terms of the truss column element properties which for calculations were taken as slender graphite/epoxy tubes. Column loads, resulting from gravity gradient control and orbital transfer, are found to be small for the class structure investigated. Fundamental frequencies of large truss structures are shown to be an order of magnitude lower than large earth based structures. Permissible loads are shown to result in small lateral deflections of the truss due to low-strain at Euler buckling of the slender graphite/epoxy truss column elements. Lateral thermal deflections are found to be a fraction of the truss depth using graphite/epoxy columns.

  11. On ideal structure in quadratic DDS in R{sup 2}

    SciTech Connect

    Kutnjak, Milan

    2008-11-13

    We consider the dynamics in a special case of two-dimensional quadratic homogeneous discrete dynamical systems. It is well known (c.f. [1, 2]) that homogeneous quadratic maps are in one to one correspondence with two-dimensional commutative (nonassociative) algebras. Algebraic concepts (such as the structure of algebra and existence of special elements like idempotents and nilpotents) help us to study the dynamics of the corresponding discrete homogeneous quadratic maps. It is well-known that such systems can exhibit chaotic behavior [3], In this article we consider the influence of the existence of an algebra ideal on the dynamics of the corresponding discrete homogeneous quadratic system. We also present some examples in the plane.

  12. Quality Tetrahedral Mesh Smoothing via Boundary-Optimized Delaunay Triangulation.

    PubMed

    Gao, Zhanheng; Yu, Zeyun; Holst, Michael

    2012-12-01

    Despite its great success in improving the quality of a tetrahedral mesh, the original optimal Delaunay triangulation (ODT) is designed to move only inner vertices and thus cannot handle input meshes containing "bad" triangles on boundaries. In the current work, we present an integrated approach called boundary-optimized Delaunay triangulation (B-ODT) to smooth (improve) a tetrahedral mesh. In our method, both inner and boundary vertices are repositioned by analytically minimizing the error between a paraboloid function and its piecewise linear interpolation over the neighborhood of each vertex. In addition to the guaranteed volume-preserving property, the proposed algorithm can be readily adapted to preserve sharp features in the original mesh. A number of experiments are included to demonstrate the performance of our method. PMID:23144522

  13. Transition Strength Ratios in the Tetrahedral Candidate ^156Dy

    NASA Astrophysics Data System (ADS)

    Hartley, D. J.; Riedinger, L. L.; Curien, D.; Dudek, J.; Gall, B.; Allmond, J. M.; Beausang, C. W.; Carpenter, M. P.; Chiara, C. J.; Janssens, R. V. F.; Kondev, F. G.; Lauritsen, T.; McCutchan, E. A.; Stefanescu, I.; Zhu, S.; Garrett, P. E.; Kulp, W. D.; Wood, J. L.; Mazurek, K.; Riley, M. A.; Wang, X.; Schunck, N.; Yu, C.-H.; Sharpey-Schafer, J.; Simpson, J.

    2009-10-01

    A new symmetry has been recently proposed where nuclei may stabilize in a tetrahedral (pyramid) shape. One of the consequences of this symmetry is that the transition strength, B(E2), of the inband transitions should approach zero in the ideal case. Thus, one signal of this exotic shape would be a rotational band where the inband E2 transitions are extremely weak or nonexistent. Such bands exist in many of the lowest negative-parity bands in the N 90 nuclei, which is also a predicted ``magic" region for tetrahedral symmetry. A Gammasphere experiment was performed to measure the B(E2)/B(E1) ratios of such a negative-parity band in ^156Dy. The results (which are consistent with the theory) will be presented, as well as a discussion of the proposed follow-up experiment to directly measure the B(E2) rates.

  14. Quality Tetrahedral Mesh Smoothing via Boundary-Optimized Delaunay Triangulation

    PubMed Central

    Gao, Zhanheng; Yu, Zeyun; Holst, Michael

    2012-01-01

    Despite its great success in improving the quality of a tetrahedral mesh, the original optimal Delaunay triangulation (ODT) is designed to move only inner vertices and thus cannot handle input meshes containing “bad” triangles on boundaries. In the current work, we present an integrated approach called boundary-optimized Delaunay triangulation (B-ODT) to smooth (improve) a tetrahedral mesh. In our method, both inner and boundary vertices are repositioned by analytically minimizing the error between a paraboloid function and its piecewise linear interpolation over the neighborhood of each vertex. In addition to the guaranteed volume-preserving property, the proposed algorithm can be readily adapted to preserve sharp features in the original mesh. A number of experiments are included to demonstrate the performance of our method. PMID:23144522

  15. Coarse-grained theory of a realistic tetrahedral liquid model

    NASA Astrophysics Data System (ADS)

    Procaccia, I.; Regev, I.

    2012-02-01

    Tetrahedral liquids such as water and silica-melt show unusual thermodynamic behavior such as a density maximum and an increase in specific heat when cooled to low temperatures. Previous work had shown that Monte Carlo and mean-field solutions of a lattice model can exhibit these anomalous properties with or without a phase transition, depending on the values of the different terms in the Hamiltonian. Here we use a somewhat different approach, where we start from a very popular empirical model of tetrahedral liquids —the Stillinger-Weber model— and construct a coarse-grained theory which directly quantifies the local structure of the liquid as a function of volume and temperature. We compare the theory to molecular-dynamics simulations and show that the theory can rationalize the simulation results and the anomalous behavior.

  16. Optimal channels for channelized quadratic estimators.

    PubMed

    Kupinski, Meredith K; Clarkson, Eric

    2016-06-01

    We present a new method for computing optimized channels for estimation tasks that is feasible for high-dimensional image data. Maximum-likelihood (ML) parameter estimates are challenging to compute from high-dimensional likelihoods. The dimensionality reduction from M measurements to L channels is a critical advantage of channelized quadratic estimators (CQEs), since estimating likelihood moments from channelized data requires smaller sample sizes and inverting a smaller covariance matrix is easier. The channelized likelihood is then used to form ML estimates of the parameter(s). In this work we choose an imaging example in which the second-order statistics of the image data depend upon the parameter of interest: the correlation length. Correlation lengths are used to approximate background textures in many imaging applications, and in these cases an estimate of the correlation length is useful for pre-whitening. In a simulation study we compare the estimation performance, as measured by the root-mean-squared error (RMSE), of correlation length estimates from CQE and power spectral density (PSD) distribution fitting. To abide by the assumptions of the PSD method we simulate an ergodic, isotropic, stationary, and zero-mean random process. These assumptions are not part of the CQE formalism. The CQE method assumes a Gaussian channelized likelihood that can be a valid for non-Gaussian image data, since the channel outputs are formed from weighted sums of the image elements. We have shown that, for three or more channels, the RMSE of CQE estimates of correlation length is lower than conventional PSD estimates. We also show that computing CQE by using a standard nonlinear optimization method produces channels that yield RMSE within 2% of the analytic optimum. CQE estimates of anisotropic correlation length estimation are reported to demonstrate this technique on a two-parameter estimation problem. PMID:27409452

  17. Mesh quality control for multiply-refined tetrahedral grids

    NASA Technical Reports Server (NTRS)

    Biswas, Rupak; Strawn, Roger

    1994-01-01

    A new algorithm for controlling the quality of multiply-refined tetrahedral meshes is presented in this paper. The basic dynamic mesh adaption procedure allows localized grid refinement and coarsening to efficiently capture aerodynamic flow features in computational fluid dynamics problems; however, repeated application of the procedure may significantly deteriorate the quality of the mesh. Results presented show the effectiveness of this mesh quality algorithm and its potential in the area of helicopter aerodynamics and acoustics.

  18. ATTILA: A three-dimensional, unstructured tetrahedral mesh discrete ordinates transport code

    SciTech Connect

    Wareing, T.A.; McGhee, J.M.; Morel, J.E.

    1996-12-31

    Many applications of radiation transport require the accurate modeling of complex three-dimensional geometries. Historically, Monte Carlo codes have been used for such applications. Existing deterministic transport codes were not applied to such problems because of the difficulties of modeling complex three-dimensional geometries with rectangular meshes. The authors have developed a three-dimensional discrete ordinates (S{sub n}) code, ATTILA, which uses linear-discontinuous finite element spatial differencing in conjunction with diffusion-synthetic acceleration (DSA) on an unstructured tetrahedral mesh. This tetrahedral mesh capability enables the authors to efficiently model complex three-dimensional geometries. One interesting and challenging application of neutron and/or gamma-ray transport is nuclear well-logging applications. Nuclear well-logging problems usually involve a complex geometry with fixed sources and one or more detectors. Detector responses must generally be accurate to within {approx}1%. The combination of complex three-dimensional geometries and high accuracy requirements makes it difficult to perform logging problems with traditional S{sub n} differencing schemes and rectangular meshes. Hence, it is not surprising that deterministic S{sub n} codes have seen limited use in nuclear well-logging applications. The geometric modeling capabilities and the advanced spatial differencing of ATTILA give it a significant advantage, relative to traditional S{sub n} codes, for performing nuclear well-logging calculations.

  19. Schwarz and multilevel methods for quadratic spline collocation

    SciTech Connect

    Christara, C.C.; Smith, B.

    1994-12-31

    Smooth spline collocation methods offer an alternative to Galerkin finite element methods, as well as to Hermite spline collocation methods, for the solution of linear elliptic Partial Differential Equations (PDEs). Recently, optimal order of convergence spline collocation methods have been developed for certain degree splines. Convergence proofs for smooth spline collocation methods are generally more difficult than for Galerkin finite elements or Hermite spline collocation, and they require stronger assumptions and more restrictions. However, numerical tests indicate that spline collocation methods are applicable to a wider class of problems, than the analysis requires, and are very competitive to finite element methods, with respect to efficiency. The authors will discuss Schwarz and multilevel methods for the solution of elliptic PDEs using quadratic spline collocation, and compare these with domain decomposition methods using substructuring. Numerical tests on a variety of parallel machines will also be presented. In addition, preliminary convergence analysis using Schwarz and/or maximum principle techniques will be presented.

  20. The Factorability of Quadratics: Motivation for More Techniques

    ERIC Educational Resources Information Center

    Bosse, Michael J.; Nandakumar, N. R.

    2005-01-01

    Typically, secondary and college algebra students attempt to utilize either completing the square or the quadratic formula as techniques to solve a quadratic equation only after frustration with factoring has arisen. While both completing the square and the quadratic formula are techniques which can determine solutions for all quadratic equations,…

  1. Fast approximate quadratic programming for graph matching.

    PubMed

    Vogelstein, Joshua T; Conroy, John M; Lyzinski, Vince; Podrazik, Louis J; Kratzer, Steven G; Harley, Eric T; Fishkind, Donniell E; Vogelstein, R Jacob; Priebe, Carey E

    2015-01-01

    Quadratic assignment problems arise in a wide variety of domains, spanning operations research, graph theory, computer vision, and neuroscience, to name a few. The graph matching problem is a special case of the quadratic assignment problem, and graph matching is increasingly important as graph-valued data is becoming more prominent. With the aim of efficiently and accurately matching the large graphs common in big data, we present our graph matching algorithm, the Fast Approximate Quadratic assignment algorithm. We empirically demonstrate that our algorithm is faster and achieves a lower objective value on over 80% of the QAPLIB benchmark library, compared with the previous state-of-the-art. Applying our algorithm to our motivating example, matching C. elegans connectomes (brain-graphs), we find that it efficiently achieves performance. PMID:25886624

  2. Fast Approximate Quadratic Programming for Graph Matching

    PubMed Central

    Vogelstein, Joshua T.; Conroy, John M.; Lyzinski, Vince; Podrazik, Louis J.; Kratzer, Steven G.; Harley, Eric T.; Fishkind, Donniell E.; Vogelstein, R. Jacob; Priebe, Carey E.

    2015-01-01

    Quadratic assignment problems arise in a wide variety of domains, spanning operations research, graph theory, computer vision, and neuroscience, to name a few. The graph matching problem is a special case of the quadratic assignment problem, and graph matching is increasingly important as graph-valued data is becoming more prominent. With the aim of efficiently and accurately matching the large graphs common in big data, we present our graph matching algorithm, the Fast Approximate Quadratic assignment algorithm. We empirically demonstrate that our algorithm is faster and achieves a lower objective value on over 80% of the QAPLIB benchmark library, compared with the previous state-of-the-art. Applying our algorithm to our motivating example, matching C. elegans connectomes (brain-graphs), we find that it efficiently achieves performance. PMID:25886624

  3. Limit cycles near hyperbolas in quadratic systems

    NASA Astrophysics Data System (ADS)

    Artés, Joan C.; Dumortier, Freddy; Llibre, Jaume

    In this paper we introduce the notion of infinity strip and strip of hyperbolas as organizing centers of limit cycles in polynomial differential systems on the plane. We study a strip of hyperbolas occurring in some quadratic systems. We deal with the cyclicity of the degenerate graphics DI2a from the programme, set up in [F. Dumortier, R. Roussarie, C. Rousseau, Hilbert's 16th problem for quadratic vector fields, J. Differential Equations 110 (1994) 86-133], to solve the finiteness part of Hilbert's 16th problem for quadratic systems. Techniques from geometric singular perturbation theory are combined with the use of the Bautin ideal. We also rely on the theory of Darboux integrability.

  4. Quadratic forms of projective spaces over rings

    NASA Astrophysics Data System (ADS)

    Levchuk, V. M.; Starikova, O. A.

    2006-06-01

    In the passage from fields to rings of coefficients quadratic forms with invertible matrices lose their decisive role. It turns out that if all quadratic forms over a ring are diagonalizable, then in effect this is always a local principal ideal ring R with 2\\in R^*. The problem of the construction of a `normal' diagonal form of a quadratic form over a ring R faces obstacles in the case of indices \\vert R^*:R^{*2}\\vert greater than 1. In the case of index 2 this problem has a solution given in Theorem 2.1 for 1+R^{*2}\\subseteq R^{*2} (an extension of the law of inertia for real quadratic forms) and in Theorem 2.2 for 1+R^2 containing an invertible non-square. Under the same conditions on a ring R with nilpotent maximal ideal the number of classes of projectively congruent quadratic forms of the projective space associated with a free R-module of rank n is explicitly calculated (Proposition 3.2). Up to projectivities, the list of forms is presented for the projective plane over R and also (Theorem 3.3) over the local ring F\\lbrack\\lbrack x,y\\rbrack\\rbrack/\\langle x^{2},xy,y^{2}\\rangle with non-principal maximal ideal, where F=2F is a field with an invertible non-square in 1+F^{2} and \\vert F^{*}:F^{*2}\\vert=2. In the latter case the number of classes of non-diagonalizable quadratic forms of rank 0 depends on one's choice of the field F and is not even always finite; all the other forms make up 21 classes.

  5. Quintessence with quadratic coupling to dark matter

    SciTech Connect

    Boehmer, Christian G.; Chan, Nyein; Caldera-Cabral, Gabriela; Lazkoz, Ruth; Maartens, Roy

    2010-04-15

    We introduce a new form of coupling between dark energy and dark matter that is quadratic in their energy densities. Then we investigate the background dynamics when dark energy is in the form of exponential quintessence. The three types of quadratic coupling all admit late-time accelerating critical points, but these are not scaling solutions. We also show that two types of coupling allow for a suitable matter era at early times and acceleration at late times, while the third type of coupling does not admit a suitable matter era.

  6. Heredity in one-dimensional quadratic maps

    NASA Astrophysics Data System (ADS)

    Romera, M.; Pastor, G.; Alvarez, G.; Montoya, F.

    1998-12-01

    In an iterative process, as is the case of a one-dimensional quadratic map, heredity has never been mentioned. In this paper we show that the pattern of a superstable orbit of a one-dimensional quadratic map can be expressed as the sum of the gene of the chaotic band where the pattern is to be found, and the ancestral path that joins all its ancestors. The ancestral path holds all the needed genetic information to calculate the descendants of the pattern. The ancestral path and successive descendant generations of the pattern constitute the family tree of the pattern, which is important to study and understand the orbit's ordering.

  7. On orthogonality preserving quadratic stochastic operators

    SciTech Connect

    Mukhamedov, Farrukh; Taha, Muhammad Hafizuddin Mohd

    2015-05-15

    A quadratic stochastic operator (in short QSO) is usually used to present the time evolution of differing species in biology. Some quadratic stochastic operators have been studied by Lotka and Volterra. In the present paper, we first give a simple characterization of Volterra QSO in terms of absolutely continuity of discrete measures. Further, we introduce a notion of orthogonal preserving QSO, and describe such kind of operators defined on two dimensional simplex. It turns out that orthogonal preserving QSOs are permutations of Volterra QSO. The associativity of genetic algebras generated by orthogonal preserving QSO is studied too.

  8. Quadratic-Like Dynamics of Cubic Polynomials

    NASA Astrophysics Data System (ADS)

    Blokh, Alexander; Oversteegen, Lex; Ptacek, Ross; Timorin, Vladlen

    2016-02-01

    A small perturbation of a quadratic polynomial f with a non-repelling fixed point gives a polynomial g with an attracting fixed point and a Jordan curve Julia set, on which g acts like angle doubling. However, there are cubic polynomials with a non-repelling fixed point, for which no perturbation results into a polynomial with Jordan curve Julia set. Motivated by the study of the closure of the Cubic Principal Hyperbolic Domain, we describe such polynomials in terms of their quadratic-like restrictions.

  9. Guises and disguises of quadratic divergences

    SciTech Connect

    Cherchiglia, A.L.; Vieira, A.R.; Hiller, Brigitte; Baêta Scarpelli, A.P.; Sampaio, Marcos

    2014-12-15

    In this contribution, we present a new perspective on the control of quadratic divergences in quantum field theory, in general, and in the Higgs naturalness problem, in particular. Our discussion is essentially based on an approach where UV divergences are parameterized, after being reduced to basic divergent integrals (BDI) in one internal momentum, as functions of a cutoff and a renormalization group scale λ. We illustrate our proposal with well-known examples, such as the gluon vacuum self energy of QCD and the Higgs decay in two photons within this approach. We also discuss frameworks in effective low-energy QCD models, where quadratic divergences are indeed fundamental.

  10. An assessment of the adaptive unstructured tetrahedral grid, Euler Flow Solver Code FELISA

    NASA Technical Reports Server (NTRS)

    Djomehri, M. Jahed; Erickson, Larry L.

    1994-01-01

    A three-dimensional solution-adaptive Euler flow solver for unstructured tetrahedral meshes is assessed, and the accuracy and efficiency of the method for predicting sonic boom pressure signatures about simple generic models are demonstrated. Comparison of computational and wind tunnel data and enhancement of numerical solutions by means of grid adaptivity are discussed. The mesh generation is based on the advancing front technique. The FELISA code consists of two solvers, the Taylor-Galerkin and the Runge-Kutta-Galerkin schemes, both of which are spacially discretized by the usual Galerkin weighted residual finite-element methods but with different explicit time-marching schemes to steady state. The solution-adaptive grid procedure is based on either remeshing or mesh refinement techniques. An alternative geometry adaptive procedure is also incorporated.

  11. Analytic relations for reconstructing piecewise linear interfaces in triangular and tetrahedral grids

    NASA Astrophysics Data System (ADS)

    Yang, Xiaofeng; James, Ashley J.

    2006-05-01

    In volume of fluid methods for interfacial flow simulations, one essential process is the so-called interface reconstruction, in which an approximate interface is reconstructed from a given discrete volume fraction field. In [J. Comput. Phys. 164 (2000) 228-237], Scardovelli and Zaleski presented analytical relations connecting linear interfaces and volume fractions in rectangular grids. Here, we present analytical relations connecting linear interfaces and volume fractions in triangular and tetrahedral grids. For computing the volume of fluid in an arbitrary polygonal or polyhedral fluid element, we also cite some of the most efficient formulas for polygon area and polyhedron volume computations. Simple test cases show that this analytic method of interface reconstruction is about 18 times faster than an iterative method in two dimensions, and four to six times faster in three dimensions. The results can be in general applied to other fields as well.

  12. Platelet adhesion on phosphorus-incorporated tetrahedral amorphous carbon films

    NASA Astrophysics Data System (ADS)

    Liu, Aiping; Zhu, Jiaqi; Liu, Meng; Dai, Zhifei; Han, Xiao; Han, Jiecai

    2008-11-01

    The haemocompatibility of phosphorus-incorporated tetrahedral amorphous carbon (ta-C:P) films, synthesized by filtered cathodic vacuum arc technique with PH 3 as the dopant source, was assessed by in vitro platelet adhesion tests. Results based on scanning electron microscopy and contact angle measurements reveal that phosphorus incorporation improves the wettability and blood compatibility of ta-C film. Our studies may provide a novel approach for the design and synthesis of doped ta-C films to repel platelet adhesion and reduce thrombosis risk.

  13. Novel biomedical tetrahedral mesh methods: algorithms and applications

    NASA Astrophysics Data System (ADS)

    Yu, Xiao; Jin, Yanfeng; Chen, Weitao; Huang, Pengfei; Gu, Lixu

    2007-12-01

    Tetrahedral mesh generation algorithm, as a prerequisite of many soft tissue simulation methods, becomes very important in the virtual surgery programs because of the real-time requirement. Aiming to speed up the computation in the simulation, we propose a revised Delaunay algorithm which makes a good balance of quality of tetrahedra, boundary preservation and time complexity, with many improved methods. Another mesh algorithm named Space-Disassembling is also presented in this paper, and a comparison of Space-Disassembling, traditional Delaunay algorithm and the revised Delaunay algorithm is processed based on clinical soft-tissue simulation projects, including craniofacial plastic surgery and breast reconstruction plastic surgery.

  14. Evaluation of a Kinematically-Driven Finite Element Footstrike Model.

    PubMed

    Hannah, Iain; Harland, Andy; Price, Dan; Schlarb, Heiko; Lucas, Tim

    2016-06-01

    A dynamic finite element model of a shod running footstrike was developed and driven with 6 degree of freedom foot segment kinematics determined from a motion capture running trial. Quadratic tetrahedral elements were used to mesh the footwear components with material models determined from appropriate mechanical tests. Model outputs were compared with experimental high-speed video (HSV) footage, vertical ground reaction force (GRF), and center of pressure (COP) excursion to determine whether such an approach is appropriate for the development of athletic footwear. Although unquantified, good visual agreement to the HSV footage was observed but significant discrepancies were found between the model and experimental GRF and COP readings (9% and 61% of model readings outside of the mean experimental reading ± 2 standard deviations, respectively). Model output was also found to be highly sensitive to input kinematics with a 120% increase in maximum GRF observed when translating the force platform 2 mm vertically. While representing an alternative approach to existing dynamic finite element footstrike models, loading highly representative of an experimental trial was not found to be achievable when employing exclusively kinematic boundary conditions. This significantly limits the usefulness of employing such an approach in the footwear development process. PMID:26671721

  15. Quadratic minima and modular forms II

    NASA Astrophysics Data System (ADS)

    Brent, Barry

    We give upper bounds on the size of the gap between a non-zero constant term and the next non-zero Fourier coefficient of an entire level two modular form. We give upper bounds for the minimum positive integer represented by a level two even positive-definite quadratic form. These bounds extend partial results in part I.

  16. Curious Consequences of a Miscopied Quadratic

    ERIC Educational Resources Information Center

    Poet, Jeffrey L.; Vestal, Donald L., Jr.

    2005-01-01

    The starting point of this article is a search for pairs of quadratic polynomials x[superscript 2] + bx plus or minus c with the property that they both factor over the integers. The search leads quickly to some number theory in the form of primitive Pythagorean triples, and this paper develops the connection between these two topics.

  17. Integration of the Quadratic Function and Generalization

    ERIC Educational Resources Information Center

    Mitsuma, Kunio

    2011-01-01

    We will first recall useful formulas in integration that simplify the calculation of certain definite integrals with the quadratic function. A main formula relies only on the coefficients of the function. We will then explore a geometric proof of one of these formulas. Finally, we will extend the formulas to more general cases. (Contains 3…

  18. Fourier analysis of quadratic phase interferograms

    NASA Astrophysics Data System (ADS)

    Muñoz-Maciel, Jesús; Mora-González, Miguel; Casillas-Rodríguez, Francisco J.; Peña-Lecona, Francisco G.

    2015-06-01

    A phase demodulation method from a single interferogram with a quadratic phase term is developed. The fringe pattern being analysed may contain circular, elliptic or astigmatic fringes. The Fourier transform of such interferograms is seen to be also a sine or a cosine of a second order polynomial in both the real and imaginary parts. In this work we take a discrete Fourier transform of the fringe patterns and then we take separate inverse discrete transforms of the real and imaginary parts of the frequency spectrum. This results in two new interferograms corresponding to the sine and cosine of the quadratic term of the phase modulated by the sine and cosine of the linear term. The linear term of these interferograms may be recovered with similar procedures of fringe analysis from open fringe interferograms. Once the linear term is retrieved the quadratic phase of the interferogram being analysed can also be calculated. The present approach is also being investigated for interferograms with nearly circularly symmetry given that the phase contains some tilt. The described procedure of Fourier analysis from quadratic phase interferograms of nearly symmetric interferograms could be used instead of complex and time consuming algorithms for phase recovery from fringe patterns with closed fringes. Finally, the method is tested in simulated and real data.

  19. Robust linear quadratic designs with respect to parameter uncertainty

    NASA Technical Reports Server (NTRS)

    Douglas, Joel; Athans, Michael

    1992-01-01

    The authors derive a linear quadratic regulator (LQR) which is robust to parametric uncertainty by using the overbounding method of I. R. Petersen and C. V. Hollot (1986). The resulting controller is determined from the solution of a single modified Riccati equation. It is shown that, when applied to a structural system, the controller gains add robustness by minimizing the potential energy of uncertain stiffness elements, and minimizing the rate of dissipation of energy through uncertain damping elements. A worst-case disturbance in the direction of the uncertainty is also considered. It is proved that performance robustness has been increased with the robust LQR when compared to a mismatched LQR design where the controller is designed on the nominal system, but applied to the actual uncertain system.

  20. The Space-Time Conservative Schemes for Large-Scale, Time-Accurate Flow Simulations with Tetrahedral Meshes

    NASA Technical Reports Server (NTRS)

    Venkatachari, Balaji Shankar; Streett, Craig L.; Chang, Chau-Lyan; Friedlander, David J.; Wang, Xiao-Yen; Chang, Sin-Chung

    2016-01-01

    Despite decades of development of unstructured mesh methods, high-fidelity time-accurate simulations are still predominantly carried out on structured, or unstructured hexahedral meshes by using high-order finite-difference, weighted essentially non-oscillatory (WENO), or hybrid schemes formed by their combinations. In this work, the space-time conservation element solution element (CESE) method is used to simulate several flow problems including supersonic jet/shock interaction and its impact on launch vehicle acoustics, and direct numerical simulations of turbulent flows using tetrahedral meshes. This paper provides a status report for the continuing development of the space-time conservation element solution element (CESE) numerical and software framework under the Revolutionary Computational Aerosciences (RCA) project. Solution accuracy and large-scale parallel performance of the numerical framework is assessed with the goal of providing a viable paradigm for future high-fidelity flow physics simulations.

  1. Factorization using the quadratic sieve algorithm

    SciTech Connect

    Davis, J.A.; Holdridge, D.B.

    1983-12-01

    Since the cryptosecurity of the RSA two key cryptoalgorithm is no greater than the difficulty of factoring the modulus (product of two secret primes), a code that implements the Quadratic Sieve factorization algorithm on the CRAY I computer has been developed at the Sandia National Laboratories to determine as sharply as possible the current state-of-the-art in factoring. Because all viable attacks on RSA thus far proposed are equivalent to factorization of the modulus, sharper bounds on the computational difficulty of factoring permit improved estimates for the size of RSA parameters needed for given levels of cryptosecurity. Analysis of the Quadratic Sieve indicates that it may be faster than any previously published general purpose algorithm for factoring large integers. The high speed of the CRAY I coupled with the capability of the CRAY to pipeline certain vectorized operations make this algorithm (and code) the front runner in current factoring techniques.

  2. Factorization using the quadratic sieve algorithm

    SciTech Connect

    Davis, J.A.; Holdridge, D.B.

    1983-01-01

    Since the cryptosecurity of the RSA two key cryptoalgorithm is no greater than the difficulty of factoring the modulus (product of two secret primes), a code that implements the Quadratic Sieve factorization algorithm on the CRAY I computer has been developed at the Sandia National Laboratories to determine as sharply as possible the current state-of-the-art in factoring. Because all viable attacks on RSA thus far proposed are equivalent to factorization of the modulus, sharper bounds on the computational difficulty of factoring permit improved estimates for the size of RSA parameters needed for given levels of cryptosecurity. Analysis of the Quadratic Sieve indicates that it may be faster than any previously published general purpose algorithm for factoring large integers. The high speed of the CRAY I coupled with the capability of the CRAY to pipeline certain vectorized operations make this algorithm (and code) the front runner in current factoring techniques.

  3. Characterization of a Quadratic Function in Rn

    ERIC Educational Resources Information Center

    Xu, Conway

    2010-01-01

    It is proved that a scalar-valued function "f"(x) defined in "n"-dimensional space must be quadratic, if the intersection of tangent planes at x[subscript 1] and x[subscript 2] always contains the midpoint of the line joining x[subscript 1] and x[subscript 2]. This is the converse of a result of Stenlund proved in this JOURNAL in 2001.

  4. Communications circuit including a linear quadratic estimator

    DOEpatents

    Ferguson, Dennis D.

    2015-07-07

    A circuit includes a linear quadratic estimator (LQE) configured to receive a plurality of measurements a signal. The LQE is configured to weight the measurements based on their respective uncertainties to produce weighted averages. The circuit further includes a controller coupled to the LQE and configured to selectively adjust at least one data link parameter associated with a communication channel in response to receiving the weighted averages.

  5. An alternative method on quadratic programming problems

    NASA Astrophysics Data System (ADS)

    Dasril, Y.; Mohd, I. B.; Mustaffa, I.; Aminuddin, MMM.

    2015-05-01

    In this paper we proposed an alternative approach to find the optimum solution of quadratic programming problems (QPP) in its original form without additional information such as slack variable, surplus variable or artificial variable as done in other favourite methods. This approached is based on the violated constraints by the unconstrained optimum. The optimal solution of QPP obtained by searching from initial point to another point alongside of feasible region.

  6. Extended Decentralized Linear-Quadratic-Gaussian Control

    NASA Technical Reports Server (NTRS)

    Carpenter, J. Russell

    2000-01-01

    A straightforward extension of a solution to the decentralized linear-Quadratic-Gaussian problem is proposed that allows its use for commonly encountered classes of problems that are currently solved with the extended Kalman filter. This extension allows the system to be partitioned in such a way as to exclude the nonlinearities from the essential algebraic relationships that allow the estimation and control to be optimally decentralized.

  7. Quadratic Programming for Allocating Control Effort

    NASA Technical Reports Server (NTRS)

    Singh, Gurkirpal

    2005-01-01

    A computer program calculates an optimal allocation of control effort in a system that includes redundant control actuators. The program implements an iterative (but otherwise single-stage) algorithm of the quadratic-programming type. In general, in the quadratic-programming problem, one seeks the values of a set of variables that minimize a quadratic cost function, subject to a set of linear equality and inequality constraints. In this program, the cost function combines control effort (typically quantified in terms of energy or fuel consumed) and control residuals (differences between commanded and sensed values of variables to be controlled). In comparison with prior control-allocation software, this program offers approximately equal accuracy but much greater computational efficiency. In addition, this program offers flexibility, robustness to actuation failures, and a capability for selective enforcement of control requirements. The computational efficiency of this program makes it suitable for such complex, real-time applications as controlling redundant aircraft actuators or redundant spacecraft thrusters. The program is written in the C language for execution in a UNIX operating system.

  8. Optimal Approximation of Quadratic Interval Functions

    NASA Technical Reports Server (NTRS)

    Koshelev, Misha; Taillibert, Patrick

    1997-01-01

    Measurements are never absolutely accurate, as a result, after each measurement, we do not get the exact value of the measured quantity; at best, we get an interval of its possible values, For dynamically changing quantities x, the additional problem is that we cannot measure them continuously; we can only measure them at certain discrete moments of time t(sub 1), t(sub 2), ... If we know that the value x(t(sub j)) at a moment t(sub j) of the last measurement was in the interval [x-(t(sub j)), x + (t(sub j))], and if we know the upper bound D on the rate with which x changes, then, for any given moment of time t, we can conclude that x(t) belongs to the interval [x-(t(sub j)) - D (t - t(sub j)), x + (t(sub j)) + D (t - t(sub j))]. This interval changes linearly with time, an is, therefore, called a linear interval function. When we process these intervals, we get an expression that is quadratic and higher order w.r.t. time t, Such "quadratic" intervals are difficult to process and therefore, it is necessary to approximate them by linear ones. In this paper, we describe an algorithm that gives the optimal approximation of quadratic interval functions by linear ones.

  9. THE EFFECTIVENESS OF QUADRATS FOR MEASURING VASCULAR PLANT DIVERSITY

    EPA Science Inventory

    Quadrats are widely used for measuring characteristics of vascular plant communities. It is well recognized that quadrat size affects measurements of frequency and cover. The ability of quadrats of varying sizes to adequately measure diversity has not been established. An exha...

  10. Graphical Solution of the Monic Quadratic Equation with Complex Coefficients

    ERIC Educational Resources Information Center

    Laine, A. D.

    2015-01-01

    There are many geometrical approaches to the solution of the quadratic equation with real coefficients. In this article it is shown that the monic quadratic equation with complex coefficients can also be solved graphically, by the intersection of two hyperbolas; one hyperbola being derived from the real part of the quadratic equation and one from…

  11. Photoconductive Detection of Tetrahedrally Coordinated Hydrogen in ZnO

    NASA Astrophysics Data System (ADS)

    Koch, S. G.; Lavrov, E. V.; Weber, J.

    2012-04-01

    In this Letter we apply an innovative experimental approach, which allows us to improve the sensitivity of detecting local vibrational modes (LVMs) even in highly absorbing spectral regions. This photoconductive technique allowed us to confirm a recent suggestion of a new multicenter bond for hydrogen in ZnO [A. Janotti and C. G. Van de Walle, Nature Mater.NMAACR1476-1122 6, 44 (2007)10.1038/nmat1795]. The two LVMs of the hydrogen substituting oxygen in ZnO are identified at 742 and 792cm-1. The modes belong to a nondegenerated A1 and a twofold degenerated E representations of the C3v point group. The tetrahedral coordination of the hydrogen atom is the result of a newly detected multicenter bond for defects in solids.

  12. Search for Tetrahedral Symmetry in Nuclei:. a Short Overview

    NASA Astrophysics Data System (ADS)

    Curien, D.; Dudek, J.; Molique, H.; Sengele, L.; Góźdź, A.; Mazurek, K.

    Following a series of experiments launched by the TetraNuc collaboration to possibly demonstrate the existence of high-rank symmetries in subatomic physics, the first experimental results on the Rare Earth region appear in publications. Meanwhile an important progress has been made on the theory side strongly suggesting that the original criterion of the static tetrahedral symmetry in the form of vanishing quadrupole moments may need to be revised to include explicitly the vibrational motion. The Actinide region seems of particular interest because of the extra stability provided by the octahedral symmetry. In this article a summary of the current experimental efforts on both the Rare-Earth and Actinide regions is given. Finally the ELMA project addressing the experimental search for the symmetries in the Actinides is briefly discussed.

  13. Optimization of Time-Dependent Particle Tracing Using Tetrahedral Decomposition

    NASA Technical Reports Server (NTRS)

    Kenwright, David; Lane, David

    1995-01-01

    An efficient algorithm is presented for computing particle paths, streak lines and time lines in time-dependent flows with moving curvilinear grids. The integration, velocity interpolation and step-size control are all performed in physical space which avoids the need to transform the velocity field into computational space. This leads to higher accuracy because there are no Jacobian matrix approximations or expensive matrix inversions. Integration accuracy is maintained using an adaptive step-size control scheme which is regulated by the path line curvature. The problem of cell-searching, point location and interpolation in physical space is simplified by decomposing hexahedral cells into tetrahedral cells. This enables the point location to be done analytically and substantially faster than with a Newton-Raphson iterative method. Results presented show this algorithm is up to six times faster than particle tracers which operate on hexahedral cells yet produces almost identical particle trajectories.

  14. Why Is the Tetrahedral Bond Angle 109 Degrees? The Tetrahedron-in-a-Cube

    ERIC Educational Resources Information Center

    Lim, Kieran F.

    2012-01-01

    The common question of why the tetrahedral angle is 109.471 degrees can be answered using a tetrahedron-in-a-cube, along with some Year 10 level mathematics. The tetrahedron-in-a-cube can also be used to demonstrate the non-polarity of tetrahedral molecules, the relationship between different types of lattice structures, and to demonstrate that…

  15. High order Godunov mixed methods on tetrahedral meshes for density driven flow simulations in porous media

    NASA Astrophysics Data System (ADS)

    Mazzia, Annamaria; Putti, Mario

    2005-09-01

    Two-dimensional Godunov mixed methods have been shown to be effective for the numerical solution of density-dependent flow and transport problems in groundwater even when concentration gradients are high and the process is dominated by density effects. This class of discretization approaches solves the flow equation by means of the mixed finite element method, thus guaranteeing mass conserving velocity fields, and discretizes the transport equation by mixed finite element and finite volumes techniques combined together via appropriate time splitting. In this paper, we extend this approach to three dimensions employing tetrahedral meshes and introduce a spatially variable time stepping procedure that improves computational efficiency while preserving accuracy by adapting the time step size according to the local Courant-Friedrichs-Lewy (CFL) constraint. Careful attention is devoted to the choice of a truly three-dimensional limiter for the advection equation in the time-splitting technique, so that to preserve second order accuracy in space (in the sense that linear functions are exactly interpolated). The three-dimensional Elder problem and the saltpool problem, recently introduced as a new benchmark for testing three-dimensional density models, provide assessments with respect to accuracy and reliability of this numerical approach.

  16. Holographic entropy increases in quadratic curvature gravity

    NASA Astrophysics Data System (ADS)

    Bhattacharjee, Srijit; Sarkar, Sudipta; Wall, Aron C.

    2015-09-01

    Standard methods for calculating the black hole entropy beyond general relativity are ambiguous when the horizon is nonstationary. We fix these ambiguities in all quadratic curvature gravity theories, by demanding that the entropy be increasing at every time, for linear perturbations to a stationary black hole. Our result matches with the entropy formula found previously in holographic entanglement entropy calculations. We explicitly calculate the entropy increase for Vaidya-like solutions in Ricci-tensor gravity to show that (unlike the Wald entropy) the holographic entropy obeys a second law.

  17. Theoretical Studies of Routes to Synthesis of Tetrahedral N4

    NASA Technical Reports Server (NTRS)

    Lee, Timothy J.; Dateo, Christopher E.

    2007-01-01

    A paper [Chem. Phys. Lett. 345, 295 (2001)] describes theoretical studies of excited electronic states of nitrogen molecules, with a view toward utilizing those states in synthesizing tetrahedral N4, or Td N4 a metastable substance under consideration as a high-energy-density rocket fuel. Several ab initio theoretical approaches were followed in these studies, including complete active space self-consistent field (CASSCF), state-averaged CASSCF (SA-CASSCF), singles configuration interaction (CIS), CIS with second-order and third-order correlation corrections [CIS(D) and CIS(3)], and linear response singles and doubles coupled-cluster (LRCCSD). Standard double zeta polarized and triple zeta double polarized one-particle basis sets were used. The CASSCF calculations overestimated the excitation energies, while SACASSCF calculations partly corrected these overestimates. The accuracy of the CIS calculations varied, depending on the particular state, while the CIS(D), CIS(3), and LRCCSD results were in generally good agreement. The energies of the lowest six excited singlet states of Td N4 as calculated by the LRCCSD were compared with the energies of possible excited states of N2 + N2 fragments, leading to the conclusion that the most likely route for synthesis of Td N4 would involve a combination of two bound quintet states of N2.

  18. Slow dynamics in a primitive tetrahedral network model.

    PubMed

    De Michele, Cristiano; Tartaglia, Piero; Sciortino, Francesco

    2006-11-28

    We report extensive Monte Carlo and event-driven molecular dynamics simulations of the fluid and liquid phase of a primitive model for silica recently introduced by Ford et al. [J. Chem. Phys. 121, 8415 (2004)]. We evaluate the isodiffusivity lines in the temperature-density plane to provide an indication of the shape of the glass transition line. Except for large densities, arrest is driven by the onset of the tetrahedral bonding pattern and the resulting dynamics is strong in Angell's classification scheme [J. Non-Cryst. Solids 131-133, 13 (1991)]. We compare structural and dynamic properties with corresponding results of two recently studied primitive models of network forming liquids-a primitive model for water and an angular-constraint-free model of four-coordinated particles-to pin down the role of the geometric constraints associated with bonding. Eventually we discuss the similarities between "glass" formation in network forming liquids and "gel" formation in colloidal dispersions of patchy particles. PMID:17144726

  19. Polyamorphism in tetrahedral substances: Similarities between silicon and ice.

    PubMed

    Garcez, K M S; Antonelli, A

    2015-07-21

    Tetrahedral substances, such as silicon, water, germanium, and silica, share various unusual phase behaviors. Among them, the so-called polyamorphism, i.e., the existence of more than one amorphous form, has been intensively investigated in the last three decades. In this work, we study the metastable relations between amorphous states of silicon in a wide range of pressures, using Monte Carlo simulations. Our results indicate that the two amorphous forms of silicon at high pressures, the high density amorphous (HDA) and the very high density amorphous (VHDA), can be decompressed from high pressure (∼20 GPa) down to the tensile regime, where both convert into the same low density amorphous. Such behavior is also observed in ice. While at high pressure (∼20 GPa), HDA is less stable than VHDA, at the pressure of 10 GPa both forms exhibit similar stability. On the other hand, at much lower pressure (∼5 GPa), HDA and VHDA are no longer the most stable forms, and, upon isobaric annealing, an even less dense form of amorphous silicon emerges, the expanded high density amorphous, again in close similarity to what occurs in ice. PMID:26203030

  20. Natural frequency of uniform and optimized tetrahedral truss platforms

    NASA Technical Reports Server (NTRS)

    Wu, K. Chauncey; Lake, Mark S.

    1994-01-01

    Qualitative and quantitative estimates for the fundamental frequency of uniform and optimized tetrahedral truss platforms are determined. A semiempirical equation is developed for the frequency of free-free uniform trusses as a function of member material properties, truss dimensions, and parasitic (nonstructural) mass fraction Mp/Mt. Optimized trusses with frequencies approximately two times those of uniform trusses are determined by varying the cross-sectional areas of member groups. Trusses with 3 to 8 rings, no parasitic mass, and member areas up to 25 times the minimum area are optimized. Frequencies computed for ranges of both Mp/Mt and the ratio of maximum area to minimum area are normalized to the frequency of a uniform truss with no parasitic mass. The normalized frequency increases with the number of rings, and both frequency and the ratio of maximum area to minimum area decrease with increasing Mp/Mt. Frequency improvements that are achievable with a limited number of member areas are estimated for a 3-ring truss by using Taguchi methods. Joint stiffness knockdown effects are also considered. Comparison of optimized and baseline uniform truss frequencies indicates that tailoring can significantly increase structural frequency; maximum gains occur for trusses with low values of Mp/Mt. This study examines frequency trends for ranges of structural parameters and may be used as a preliminary design guide.

  1. Glutamate detection by amino functionalized tetrahedral amorphous carbon surfaces.

    PubMed

    Kaivosoja, Emilia; Tujunen, Noora; Jokinen, Ville; Protopopova, Vera; Heinilehto, Santtu; Koskinen, Jari; Laurila, Tomi

    2015-08-15

    In this paper, a novel amperometric glutamate biosensor with glutamate oxidase (GlOx) immobilized directly on NH2 functionalized, platinum doped tetrahedral amorphous carbon (ta-C) film, has been successfully developed. First, we demonstrate that direct GlOx immobilization is more effective on amino-groups than on carboxyl- or hydroxyl-groups. Second, we show that anodizing and plasma treatments increase the amount of nitrogen and the proportion of protonated amino groups relative to amino groups on the aminosilane coating, which subsequently results in an increased amount of active GlOx on the surface. This effect, however, is found to be unstable due to unstable electrostatic interactions between GlOx and NH3(+). We demonstrate the detection of glutamate in the concentration range of 10µM-1mM using the NH2 functionalized Pt doped ta-C surface. The biosensor showed high sensitivity (2.9nA μM(-1)cm(-2)), low detection limit (10μM) and good storage stability. The electrode response to glutamate was linear in the concentrations ranging from 10µM to 500µM. In conclusion, the study shows that GlOx immobilization is most effective on aminosilane treated ta-C surface without any pre-treatments and the fabricated sensor structure is able to detect glutamate in the micromolar range. PMID:25966399

  2. Multi-Criterion Preliminary Design of a Tetrahedral Truss Platform

    NASA Technical Reports Server (NTRS)

    Wu, K. Chauncey

    1995-01-01

    An efficient method is presented for multi-criterion preliminary design and demonstrated for a tetrahedral truss platform. The present method requires minimal analysis effort and permits rapid estimation of optimized truss behavior for preliminary design. A 14-m-diameter, 3-ring truss platform represents a candidate reflector support structure for space-based science spacecraft. The truss members are divided into 9 groups by truss ring and position. Design variables are the cross-sectional area of all members in a group, and are either 1, 3 or 5 times the minimum member area. Non-structural mass represents the node and joint hardware used to assemble the truss structure. Taguchi methods are used to efficiently identify key points in the set of Pareto-optimal truss designs. Key points identified using Taguchi methods are the maximum frequency, minimum mass, and maximum frequency-to-mass ratio truss designs. Low-order polynomial curve fits through these points are used to approximate the behavior of the full set of Pareto-optimal designs. The resulting Pareto-optimal design curve is used to predict frequency and mass for optimized trusses. Performance improvements are plotted in frequency-mass (criterion) space and compared to results for uniform trusses. Application of constraints to frequency and mass and sensitivity to constraint variation are demonstrated.

  3. Galactic chemical evolution and nucleocosmochronology - Analytic quadratic models

    NASA Technical Reports Server (NTRS)

    Clayton, D. D.

    1985-01-01

    Quadratic models of the chemical evolution of the Galaxy for a star formation rate proportional to the square of the gas mass are studied. The search for analytic solutions to the gas mass and star mass for time-dependent rates of gaseous infall onto the disk is examined. The quadratic models are compared to models having linear star formation rates. The mass, metallicity, number of stars, and U-235/U-238 isotopic ratio for the models which are subjected to the same infall rate, the same initial disk mass, and the same final gas fraction are compared. The results of the comparison indicate that: (1) the average dwarf age is greater in the quadratic model, (2) the metallicity grows initially faster in the quadratic model, (3) the quadratic model has a smaller percentage of low-Z dwarfs, and (4) the U-235/U-238 isotopic ratio indicates a younger quadratic model.

  4. Numerical solution of quadratic matrix equations for free vibration analysis of structures

    NASA Technical Reports Server (NTRS)

    Gupta, K. K.

    1975-01-01

    This paper is concerned with the efficient and accurate solution of the eigenvalue problem represented by quadratic matrix equations. Such matrix forms are obtained in connection with the free vibration analysis of structures, discretized by finite 'dynamic' elements, resulting in frequency-dependent stiffness and inertia matrices. The paper presents a new numerical solution procedure of the quadratic matrix equations, based on a combined Sturm sequence and inverse iteration technique enabling economical and accurate determination of a few required eigenvalues and associated vectors. An alternative procedure based on a simultaneous iteration procedure is also described when only the first few modes are the usual requirement. The employment of finite dynamic elements in conjunction with the presently developed eigenvalue routines results in a most significant economy in the dynamic analysis of structures.

  5. Quadratic relations in continuous and discrete Painlevé equations

    NASA Astrophysics Data System (ADS)

    Ramani, A.; Grammaticos, B.; Tamizhmani, T.

    2000-04-01

    The quadratic relations between the solutions of a Painlevé equation and that of a different one, or the same one with a different set of parameters, are investigated in the continuous and discrete cases. We show that the quadratic relations existing for the continuous PII , PIII , PV and PVI have analogues as well as consequences in the discrete case. Moreover, the discrete Painlevé equations have quadratic relations of their own without any reference to the continuous case.

  6. A point-centered arbitrary Lagrangian Eulerian hydrodynamic approach for tetrahedral meshes

    SciTech Connect

    Morgan, Nathaniel R.; Waltz, Jacob I.; Burton, Donald E.; Charest, Marc R.; Canfield, Thomas R.; Wohlbier, John G.

    2015-02-24

    We present a three dimensional (3D) arbitrary Lagrangian Eulerian (ALE) hydrodynamic scheme suitable for modeling complex compressible flows on tetrahedral meshes. The new approach stores the conserved variables (mass, momentum, and total energy) at the nodes of the mesh and solves the conservation equations on a control volume surrounding the point. This type of an approach is termed a point-centered hydrodynamic (PCH) method. The conservation equations are discretized using an edge-based finite element (FE) approach with linear basis functions. All fluxes in the new approach are calculated at the center of each tetrahedron. A multidirectional Riemann-like problem is solved at the center of the tetrahedron. The advective fluxes are calculated by solving a 1D Riemann problem on each face of the nodal control volume. A 2-stage Runge–Kutta method is used to evolve the solution forward in time, where the advective fluxes are part of the temporal integration. The mesh velocity is smoothed by solving a Laplacian equation. The details of the new ALE hydrodynamic scheme are discussed. Results from a range of numerical test problems are presented.

  7. A point-centered arbitrary Lagrangian Eulerian hydrodynamic approach for tetrahedral meshes

    DOE PAGESBeta

    Morgan, Nathaniel R.; Waltz, Jacob I.; Burton, Donald E.; Charest, Marc R.; Canfield, Thomas R.; Wohlbier, John G.

    2015-02-24

    We present a three dimensional (3D) arbitrary Lagrangian Eulerian (ALE) hydrodynamic scheme suitable for modeling complex compressible flows on tetrahedral meshes. The new approach stores the conserved variables (mass, momentum, and total energy) at the nodes of the mesh and solves the conservation equations on a control volume surrounding the point. This type of an approach is termed a point-centered hydrodynamic (PCH) method. The conservation equations are discretized using an edge-based finite element (FE) approach with linear basis functions. All fluxes in the new approach are calculated at the center of each tetrahedron. A multidirectional Riemann-like problem is solved atmore » the center of the tetrahedron. The advective fluxes are calculated by solving a 1D Riemann problem on each face of the nodal control volume. A 2-stage Runge–Kutta method is used to evolve the solution forward in time, where the advective fluxes are part of the temporal integration. The mesh velocity is smoothed by solving a Laplacian equation. The details of the new ALE hydrodynamic scheme are discussed. Results from a range of numerical test problems are presented.« less

  8. Compact stars with quadratic equation of state

    NASA Astrophysics Data System (ADS)

    Ngubelanga, Sifiso A.; Maharaj, Sunil D.; Ray, Subharthi

    2015-05-01

    We provide new exact solutions to the Einstein-Maxwell system of equations for matter configurations with anisotropy and charge. The spacetime is static and spherically symmetric. A quadratic equation of state is utilised for the matter distribution. By specifying a particular form for one of the gravitational potentials and the electric field intensity we obtain new exact solutions in isotropic coordinates. In our general class of models, an earlier model with a linear equation of state is regained. For particular choices of parameters we regain the masses of the stars PSR J1614-2230, 4U 1608-52, PSR J1903+0327, EXO 1745-248 and SAX J1808.4-3658. A comprehensive physical analysis for the star PSR J1903+0327 reveals that our model is reasonable.

  9. Quadratic dynamical decoupling with nonuniform error suppression

    SciTech Connect

    Quiroz, Gregory; Lidar, Daniel A.

    2011-10-15

    We analyze numerically the performance of the near-optimal quadratic dynamical decoupling (QDD) single-qubit decoherence errors suppression method [J. West et al., Phys. Rev. Lett. 104, 130501 (2010)]. The QDD sequence is formed by nesting two optimal Uhrig dynamical decoupling sequences for two orthogonal axes, comprising N{sub 1} and N{sub 2} pulses, respectively. Varying these numbers, we study the decoherence suppression properties of QDD directly by isolating the errors associated with each system basis operator present in the system-bath interaction Hamiltonian. Each individual error scales with the lowest order of the Dyson series, therefore immediately yielding the order of decoherence suppression. We show that the error suppression properties of QDD are dependent upon the parities of N{sub 1} and N{sub 2}, and near-optimal performance is achieved for general single-qubit interactions when N{sub 1}=N{sub 2}.

  10. Large-scale sequential quadratic programming algorithms

    SciTech Connect

    Eldersveld, S.K.

    1992-09-01

    The problem addressed is the general nonlinear programming problem: finding a local minimizer for a nonlinear function subject to a mixture of nonlinear equality and inequality constraints. The methods studied are in the class of sequential quadratic programming (SQP) algorithms, which have previously proved successful for problems of moderate size. Our goal is to devise an SQP algorithm that is applicable to large-scale optimization problems, using sparse data structures and storing less curvature information but maintaining the property of superlinear convergence. The main features are: 1. The use of a quasi-Newton approximation to the reduced Hessian of the Lagrangian function. Only an estimate of the reduced Hessian matrix is required by our algorithm. The impact of not having available the full Hessian approximation is studied and alternative estimates are constructed. 2. The use of a transformation matrix Q. This allows the QP gradient to be computed easily when only the reduced Hessian approximation is maintained. 3. The use of a reduced-gradient form of the basis for the null space of the working set. This choice of basis is more practical than an orthogonal null-space basis for large-scale problems. The continuity condition for this choice is proven. 4. The use of incomplete solutions of quadratic programming subproblems. Certain iterates generated by an active-set method for the QP subproblem are used in place of the QP minimizer to define the search direction for the nonlinear problem. An implementation of the new algorithm has been obtained by modifying the code MINOS. Results and comparisons with MINOS and NPSOL are given for the new algorithm on a set of 92 test problems.

  11. Interpenetrating metal-organic frameworks formed by self-assembly of tetrahedral and octahedral building blocks

    SciTech Connect

    Lu Yongming; Lan Yaqian; Xu Yanhong; Su Zhongmin; Li Shunli; Zang Hongying; Xu Guangjuan

    2009-11-15

    To investigate the relationship between topological types and molecular building blocks (MBBs), we have designed and synthesized a series of three-dimensional (3D) interpenetrating metal-organic frameworks based on different polygons or polyhedra under hydrothermal conditions, namely [Cd(bpib){sub 0.5}(L{sup 1})] (1), [Cd(bpib){sub 0.5}(L{sup 2})].H{sub 2}O (2), [Cd(bpib){sub 0.5}(L{sup 3})] (3) and [Cd(bib){sub 0.5}(L{sup 1})] (4), where bpib=1,4-bis(2-(pyridin-2-yl)-1H-imidazol-1-yl)butane, bib=1,4-bis(1H-imidazol-1-yl)butane, H{sub 2}L{sup 1}=4-(4-carboxybenzyloxy)benzoic acid, H{sub 2}L{sup 2}=4,4'-(ethane-1,2-diylbis(oxy))dibenzoic acid and H{sub 2}L{sup 3}=4,4'-(1,4-phenylenebis(methylene))bis(oxy)dibenzoic acid, respectively. Their structures have been determined by single crystal X-ray diffraction analyses and further characterized by elemental analyses, IR spectra, and thermogravimetric (TG) analyses. Compounds 1-3 display alpha-Po topological nets with different degrees of interpenetration based on the similar octahedral [Cd{sub 2}(-COO){sub 4}] building blocks. Compound 4 is a six-fold interpenetrating diamondoid net based on tetrahedral MBBs. By careful inspection of these structures, we find that various carboxylic ligands and N-donor ligands with different coordination modes and conformations, and metal centers with different geometries are important for the formation of the different MBBs. It is believed that different topological types lie on different MBBs with various polygons or polyhedra. Such as four- and six-connected topologies are formed by tetrahedral and octahedral building blocks. In addition, with the increase of carboxylic ligands' length, the degrees of interpenetration have been changed in the alpha-Po topological nets. And the luminescent properties of these compounds have been investigated in detail. - Graphical abstract: A series of three-dimensional interpenetrating metal-organic frameworks based on different polygons or polyhedra

  12. Linear state feedback, quadratic weights, and closed loop eigenstructures. M.S. Thesis. Final Report

    NASA Technical Reports Server (NTRS)

    Thompson, P. M.

    1980-01-01

    Equations are derived for the angles of general multivariable root loci and linear quadratic optimal root loci, including angles of departure and approach. The generalized eigenvalue problem is used to compute angles of approach. Equations are also derived to find the sensitivity of closed loop eigenvalue and the directional derivatives of closed loop eigenvectors. An equivalence class of quadratic weights that produce the same asymptotic eigenstructure is defined, a canonical element is defined, and an algorithm to find it is given. The behavior of the optimal root locus in the nonasymptotic region is shown to be different for quadratic weights with the same asymptotic properties. An algorithm is presented that can be used to select a feedback gain matrix for the linear state feedback problem which produces a specified asymptotic eigenstructure. Another algorithm is given to compute the asymptotic eigenstructure properties inherent in a given set of quadratic weights. Finally, it is shown that optimal root loci for nongeneric problems can be approximated by generic ones in the nonasymptotic region.

  13. Tetrahedrality and structural order for hydrophobic interactions in a coarse-grained water model

    NASA Astrophysics Data System (ADS)

    Chaimovich, Aviel; Shell, M. Scott

    2014-02-01

    The hydrophobic interaction manifests two separate regimes in terms of size: Small nonpolar bodies exhibit a weak oscillatory force (versus distance) while large nonpolar surfaces exhibit a strong monotonic one. This crossover in hydrophobic behavior is typically explained in terms of water's tetrahedral structure: Its tetrahedrality is enhanced near small solutes and diminished near large planar ones. Here, we demonstrate that water's tetrahedral correlations signal this switch even in a highly simplified, isotropic, "core-softened" water model. For this task, we introduce measures of tetrahedrality based on the angular distribution of water's nearest neighbors. On a quantitative basis, the coarse-grained model of course is only approximate: (1) While greater than simple Lennard-Jones liquids, its bulk tetrahedrality remains lower than that of fully atomic models; and (2) the decay length of the large-scale hydrophobic interaction is less than has been found in experiments. Even so, the qualitative behavior of the model is surprisingly rich and exhibits numerous waterlike hydrophobic behaviors, despite its simplicity. We offer several arguments for the manner in which it should be able to (at least partially) reproduce tetrahedral correlations underlying these effects.

  14. Analysis of Students' Error in Learning of Quadratic Equations

    ERIC Educational Resources Information Center

    Zakaria, Effandi; Ibrahim; Maat, Siti Mistima

    2010-01-01

    The purpose of the study was to determine the students' error in learning quadratic equation. The samples were 30 form three students from a secondary school in Jambi, Indonesia. Diagnostic test was used as the instrument of this study that included three components: factorization, completing the square and quadratic formula. Diagnostic interview…

  15. Tangent Lines without Derivatives for Quadratic and Cubic Equations

    ERIC Educational Resources Information Center

    Carroll, William J.

    2009-01-01

    In the quadratic equation, y = ax[superscript 2] + bx + c, the equation y = bx + c is identified as the equation of the line tangent to the parabola at its y-intercept. This is extended to give a convenient method of graphing tangent lines at any point on the graph of a quadratic or a cubic equation. (Contains 5 figures.)

  16. Effects of Classroom Instruction on Students' Understanding of Quadratic Equations

    ERIC Educational Resources Information Center

    Vaiyavutjamai, Pongchawee; Clements, M. A.

    2006-01-01

    Two hundred and thirty-one students in six Grade 9 classes in two government secondary schools located near Chiang Mai, Thailand, attempted to solve the same 18 quadratic equations before and after participating in 11 lessons on quadratic equations. Data from the students' written responses to the equations, together with data in the form of…

  17. Visualising the Roots of Quadratic Equations with Complex Coefficients

    ERIC Educational Resources Information Center

    Bardell, Nicholas S.

    2014-01-01

    This paper is a natural extension of the root visualisation techniques first presented by Bardell (2012) for quadratic equations with real coefficients. Consideration is now given to the familiar quadratic equation "y = ax[superscript 2] + bx + c" in which the coefficients "a," "b," "c" are generally…

  18. Geometric quadratic stochastic operator on countable infinite set

    SciTech Connect

    Ganikhodjaev, Nasir; Hamzah, Nur Zatul Akmar

    2015-02-03

    In this paper we construct the family of Geometric quadratic stochastic operators defined on the countable sample space of nonnegative integers and investigate their trajectory behavior. Such operators can be reinterpreted in terms of of evolutionary operator of free population. We show that Geometric quadratic stochastic operators are regular transformations.

  19. Quadratic divergences and quantum gravitational contributions to gauge coupling constants

    NASA Astrophysics Data System (ADS)

    Toms, David J.

    2011-10-01

    The calculation of quadratic divergences in Einstein-Maxwell theory with a possible cosmological constant is considered. We describe a method of calculation, using the background-field method, that is sensitive to quadratic divergences, is respectful of gauge invariance, and is independent of gauge conditions. A standard renormalization group analysis is applied to the result where it is shown that the quadratic divergences do lead to asymptotic freedom as found in the original paper of Robinson and Wilczek. The role and nature of these quadratic divergences is critically evaluated in light of recent criticism. Within the context of the background-field method, it is shown that it is possible to define the charge in a physically motivated way in which the quadratic divergences do not play a role. This latter view is studied in more depth in a toy model described in an appendix.

  20. Quadratic algebras for three-dimensional superintegrable systems

    SciTech Connect

    Daskaloyannis, C. Tanoudis, Y.

    2010-02-15

    The three-dimensional superintegrable systems with quadratic integrals of motion have five functionally independent integrals, one among them is the Hamiltonian. Kalnins, Kress, and Miller have proved that in the case of nondegenerate potentials with quadratic integrals of motion there is a sixth quadratic integral, which is linearly independent of the other integrals. The existence of this sixth integral implies that the integrals of motion form a ternary parafermionic-like quadratic Poisson algebra with five generators. In this contribution we investigate the structure of this algebra. We show that in all the nondegenerate cases there is at least one subalgebra of three integrals having a Poisson quadratic algebra structure, which is similar to the two-dimensional case.

  1. Polychromatic solitons in a quadratic medium.

    PubMed

    Towers, I N; Malomed, B A

    2002-10-01

    We introduce the simplest model to describe parametric interactions in a quadratically nonlinear optical medium with the fundamental harmonic containing two components with (slightly) different carrier frequencies [which is a direct analog of wavelength-division multiplexed models, well known in media with cubic nonlinearity]. The model takes a closed form with three different second-harmonic components, and it is formulated in the spatial domain. We demonstrate that the model supports both polychromatic solitons (PCSs), with all the components present in them, and two types of mutually orthogonal simple solitons, both types being stable in a broad parametric region. An essential peculiarity of PCS is that its power is much smaller than that of a simple (usual) soliton (taken at the same values of control parameters), which may be an advantage for experimental generation of PCSs. Collisions between the orthogonal simple solitons are simulated in detail, leading to the conclusion that the collisions are strongly inelastic, converting the simple solitons into polychromatic ones, and generating one or two additional PCSs. A collision velocity at which the inelastic effects are strongest is identified, and it is demonstrated that the collision may be used as a basis to design a simple all-optical XOR logic gate. PMID:12443362

  2. A Quadratic Closure for Compressible Turbulence

    SciTech Connect

    Futterman, J A

    2008-09-16

    We have investigated a one-point closure model for compressible turbulence based on third- and higher order cumulant discard for systems undergoing rapid deformation, such as might occur downstream of a shock or other discontinuity. In so doing, we find the lowest order contributions of turbulence to the mean flow, which lead to criteria for Adaptive Mesh Refinement. Rapid distortion theory (RDT) as originally applied by Herring closes the turbulence hierarchy of moment equations by discarding third order and higher cumulants. This is similar to the fourth-order cumulant discard hypothesis of Millionshchikov, except that the Millionshchikov hypothesis was taken to apply to incompressible homogeneous isotropic turbulence generally, whereas RDT is applied only to fluids undergoing a distortion that is 'rapid' in the sense that the interaction of the mean flow with the turbulence overwhelms the interaction of the turbulence with itself. It is also similar to Gaussian closure, in which both second and fourth-order cumulants are retained. Motivated by RDT, we develop a quadratic one-point closure for rapidly distorting compressible turbulence, without regard to homogeneity or isotropy, and make contact with two equation turbulence models, especially the K-{var_epsilon} and K-L models, and with linear instability growth. In the end, we arrive at criteria for Adaptive Mesh Refinement in Finite Volume simulations.

  3. Elastic Model Transitions Using Quadratic Inequality Constrained Least Squares

    NASA Technical Reports Server (NTRS)

    Orr, Jeb S.

    2012-01-01

    A technique is presented for initializing multiple discrete finite element model (FEM) mode sets for certain types of flight dynamics formulations that rely on superposition of orthogonal modes for modeling the elastic response. Such approaches are commonly used for modeling launch vehicle dynamics, and challenges arise due to the rapidly time-varying nature of the rigid-body and elastic characteristics. By way of an energy argument, a quadratic inequality constrained least squares (LSQI) algorithm is employed to e ect a smooth transition from one set of FEM eigenvectors to another with no requirement that the models be of similar dimension or that the eigenvectors be correlated in any particular way. The physically unrealistic and controversial method of eigenvector interpolation is completely avoided, and the discrete solution approximates that of the continuously varying system. The real-time computational burden is shown to be negligible due to convenient features of the solution method. Simulation results are presented, and applications to staging and other discontinuous mass changes are discussed

  4. Equation for disentangling time-ordered exponentials with arbitrary quadratic generators

    SciTech Connect

    Budanov, V.G.

    1987-12-01

    In many quantum-mechanical constructions, it is necessary to disentangle an operator-valued time-ordered exponential with time-dependent generators quadratic in the creation and annihilation operators. By disentangling, one understands the finding of the matrix elements of the time-ordered exponential or, in a more general formulation. The solution of the problem can also be reduced to calculation of a matrix time-ordered exponential that solves the corresponding classical problem. However, in either case the evolution equations in their usual form do not enable one to take into account explicitly the symmetry of the system. In this paper the methods of Weyl analysis are used to find an ordinary differential equation on a matrix Lie algebra that is invariant with respect to the adjoint action of the dynamical symmetry group of a quadratic Hamiltonian and replaces the operator evolution equation for the Green's function.

  5. RELATIONSHIP BETWEEN RIGIDITY OF EXTERNAL FIXATOR AND NUMBER OF PINS: COMPUTER ANALYSIS USING FINITE ELEMENTS

    PubMed Central

    Sternick, Marcelo Back; Dallacosta, Darlan; Bento, Daniela Águida; do Reis, Marcelo Lemos

    2015-01-01

    Objective: To analyze the rigidity of a platform-type external fixator assembly, according to different numbers of pins on each clamp. Methods: Computer simulation on a large-sized Cromus dynamic external fixator (Baumer SA) was performed using a finite element method, in accordance with the standard ASTM F1541. The models were generated with approximately 450,000 quadratic tetrahedral elements. Assemblies with two, three and four Schanz pins of 5.5 mm in diameter in each clamp were compared. Every model was subjected to a maximum force of 200 N, divided into 10 sub-steps. For the components, the behavior of the material was assumed to be linear, elastic, isotropic and homogeneous. For each model, the rigidity of the assembly and the Von Mises stress distribution were evaluated. Results: The rigidity of the system was 307.6 N/mm for two pins, 369.0 N/mm for three and 437.9 N/mm for four. Conclusion: The results showed that four Schanz pins in each clamp promoted rigidity that was 19% greater than in the configuration with three pins and 42% greater than with two pins. Higher tension occurred in configurations with fewer pins. In the models analyzed, the maximum tension occurred on the surface of the pin, close to the fixation area. PMID:27047879

  6. A decentralized linear quadratic control design method for flexible structures

    NASA Technical Reports Server (NTRS)

    Su, Tzu-Jeng; Craig, Roy R., Jr.

    1990-01-01

    A decentralized suboptimal linear quadratic control design procedure which combines substructural synthesis, model reduction, decentralized control design, subcontroller synthesis, and controller reduction is proposed for the design of reduced-order controllers for flexible structures. The procedure starts with a definition of the continuum structure to be controlled. An evaluation model of finite dimension is obtained by the finite element method. Then, the finite element model is decomposed into several substructures by using a natural decomposition called substructuring decomposition. Each substructure, at this point, still has too large a dimension and must be reduced to a size that is Riccati-solvable. Model reduction of each substructure can be performed by using any existing model reduction method, e.g., modal truncation, balanced reduction, Krylov model reduction, or mixed-mode method. Then, based on the reduced substructure model, a subcontroller is designed by an LQ optimal control method for each substructure independently. After all subcontrollers are designed, a controller synthesis method called substructural controller synthesis is employed to synthesize all subcontrollers into a global controller. The assembling scheme used is the same as that employed for the structure matrices. Finally, a controller reduction scheme, called the equivalent impulse response energy controller (EIREC) reduction algorithm, is used to reduce the global controller to a reasonable size for implementation. The EIREC reduced controller preserves the impulse response energy of the full-order controller and has the property of matching low-frequency moments and low-frequency power moments. An advantage of the substructural controller synthesis method is that it relieves the computational burden associated with dimensionality. Besides that, the SCS design scheme is also a highly adaptable controller synthesis method for structures with varying configuration, or varying mass

  7. Phase recovery based on quadratic programming

    NASA Astrophysics Data System (ADS)

    Zhang, Quan Bing; Ge, Xiao Juan; Cheng, Ya Dong; Ni, Na

    2014-11-01

    Most of the information of optical wavefront is encoded in the phase which includes more details of the object. Conventional optical measuring apparatus is relatively easy to record the intensity of light, but can not measure the phase of light directly. Thus it is important to recovery the phase from the intensity measurements of the object. In recent years, the methods based on quadratic programming such as PhaseLift and PhaseCut can recover the phase of general signal exactly for overdetermined system. To retrieve the phase of sparse signal, the Compressive Phase Retrieval (CPR) algorithm combines the l1-minimization in Compressive Sensing (CS) with low-rank matrix completion problem in PhaseLift, but the result is unsatisfied. This paper focus on the recovery of the phase of sparse signal and propose a new method called the Compressive Phase Cut Retrieval (CPCR) by combining the CPR algorithm with the PhaseCut algorithm. To ensure the sparsity of the recovered signal, we use CPR method to solve a semi-definite programming problem firstly. Then apply linear transformation to the recovered signal, and set the phase of the result as the initial value of the PhaseCut problem. We use TFOCS (a library of Matlab-files) to implement the proposed CPCR algorithm in order to improve the recovered results of the CPR algorithm. Experimental results show that the proposed method can improve the accuracy of the CPR algorithm, and overcome the shortcoming of the PhaseCut method that it can not recover the sparse signal effectively.

  8. Degenerate nonlinear programming with a quadratic growth condition.

    SciTech Connect

    Anitescu, M.; Mathematics and Computer Science

    2000-01-01

    We show that the quadratic growth condition and the Mangasarian-Fromovitz constraint qualification (MFCQ) imply that local minima of nonlinear programs are isolated stationary points. As a result, when started sufficiently close to such points, an L1 exact penalty sequential quadratic programming algorithm will induce at least R-linear convergence of the iterates to such a local minimum. We construct an example of a degenerate nonlinear program with a unique local minimum satisfying the quadratic growth and the MFCQ but for which no positive semidefinite augmented Lagrangian exists. We present numerical results obtained using several nonlinear programming packages on this example and discuss its implications for some algorithms.

  9. Serine Protease Catalysis: A Computational Study of Tetrahedral Intermediates and Inhibitory Adducts.

    PubMed

    Ngo, Phong D; Mansoorabadi, Steven O; Frey, Perry A

    2016-08-01

    Peptide boronic acids and peptidyl trifluoromethyl ketones (TFKs) inhibit serine proteases by forming monoanionic, tetrahedral adducts to serine in the active sites. Investigators regard these adducts as analogs of monoanionic, tetrahedral intermediates. Density functional theory (DFT) calculations and fractional charge analysis show that tetrahedral adducts of model peptidyl TFKs are structurally and electrostatically very similar to corresponding tetrahedral intermediates. In contrast, the DFT calculations show the structures and electrostatic properties of analogous peptide boronate adducts to be significantly different. The peptide boronates display highly electrostatically positive boron, with correspondingly negative ligands in the tetrahedra. In addition, the computed boron-oxygen and boron-carbon bond lengths in peptide boronates (which are identical or very similar to the corresponding bonds in a peptide boronate adduct of α-lytic protease determined by X-ray crystallography at subangstrom resolution) are significantly longer than the corresponding bond lengths in model tetrahedral intermediates. Since protease-peptidyl TFKs incorporate low-barrier hydrogen bonds (LBHBs) between an active site histidine and aspartate, while the protease-peptide boronates do not, these data complement the spectroscopic and chemical evidence for the participation of LBHBs in catalysis by serine proteases. Moreover, while the potency of these classes of inhibitors can be correlated to the structures of the peptide moieties, the present results indicate that the strength of their bonds to serine contribute significantly to their inhibitory properties. PMID:27387593

  10. On a 'Mysterious' Case of a Quadratic Hamiltonian

    NASA Astrophysics Data System (ADS)

    Sakovich, Sergei

    2006-07-01

    We show that one of the five cases of a quadratic Hamiltonian, which were recently selected by Sokolov and Wolf who used the Kovalevskaya-Lyapunov test, fails to pass the Painlevé test for integrability.

  11. Direct Orthogonal Distance to Quadratic Surfaces in 3D.

    PubMed

    Lott, Gus K

    2014-09-01

    Discovering the orthogonal distance to a quadratic surface is a classic geometric task in vision, modeling, and robotics. I describe a simple, efficient, and stable direct solution for the orthogonal distance (foot-point) to an arbitrary quadratic surface from a general finite 3D point. The problem is expressed as the intersection of three quadratic surfaces, two of which are derived from the requirement of orthogonality of two non-coincident planes with the tangent plane to the quadric. A sixth order single-variable polynomial is directly generated in one coordinate of the surface point. The method detects intersection points at infinity and operates smoothly across all real quadratic surface classes. The method also geometrically detects continuums of orthogonal points (i.e., from the exact center of a sphere). I discuss algorithm performance, compare it to a state-of-the-art estimator, demonstrate the algorithm on synthetic data, and describe extension to arbitrary dimension. PMID:26352239

  12. Observation of all-in type tetrahedral displacements in nonmagnetic pyrochlore niobates

    NASA Astrophysics Data System (ADS)

    Torigoe, S.; Ishimoto, Y.; Aoishi, Y.; Murakawa, H.; Matsumura, D.; Yoshii, K.; Yoneda, Y.; Nishihata, Y.; Kodama, K.; Tomiyasu, K.; Ikeda, K.; Nakao, H.; Nogami, Y.; Ikeda, N.; Otomo, T.; Hanasaki, N.

    2016-02-01

    We observed all-in type Nb tetrahedral displacement in nonmagnetic pyrochlore niobates A2Nb2O7 (A =Nd0.5Ca0.5 and Y0.5Ca0.5 ) through the analysis of the neutron pair distribution function and the extended x-ray absorption function spectroscopy. The all-in type Nb tetrahedral displacement, which has the character of a charge singlet state, is driven by the formation of the bonding orbital. The diffuse scattering in the x-ray diffraction, which has the resonant component in the Nb L3 edge, indicates that the all-in type Nb tetrahedral displacement has the periodicity with its short-range correlation.

  13. LETTER TO THE EDITOR: A new approach to modelling tetrahedral amorphous carbon

    NASA Astrophysics Data System (ADS)

    Walters, J. K.; Gilkes, K. W. R.; Wicks, J. D.; Newport, R. J.

    1997-08-01

    We have generated a new model for the structure of tetrahedral amorphous carbon using a modified reverse Monte Carlo modelling method. The novel feature of this approach is the definition of three different types of carbon atom, corresponding to tetrahedral, planar and linear bonding conformations. The particular strengths of the method are the large model size (3000 atoms), that all the possible arrangements of 0953-8984/9/34/001/img7 and 0953-8984/9/34/001/img8 bonds are allowed, and that no interatomic potential is required. For the first time we have determined the distribution of 0953-8984/9/34/001/img8 bonded sites within the predominantly disordered tetrahedral structure, and we find that they form polymer-like chains and small clusters which connect the 0953-8984/9/34/001/img7 bonded regions.

  14. Haldane-Hubbard Mott Insulator: From Tetrahedral Spin Crystal to Chiral Spin Liquid.

    PubMed

    Hickey, Ciarán; Cincio, Lukasz; Papić, Zlatko; Paramekanti, Arun

    2016-04-01

    Motivated by cold atom experiments on Chern insulators, we study the honeycomb lattice Haldane-Hubbard Mott insulator of spin-1/2 fermions using exact diagonalization and density matrix renormalization group methods. We show that this model exhibits various chiral magnetic orders including a wide regime of triple-Q tetrahedral order. Incorporating third-neighbor hopping frustrates and ultimately melts this tetrahedral spin crystal. From analyzing the low energy spectrum, many-body Chern numbers, entanglement spectra, and modular matrices, we identify the molten state as a chiral spin liquid (CSL) with gapped semion excitations. We formulate and study the Chern-Simons-Higgs field theory of the exotic CSL-to-tetrahedral spin crystallization transition. PMID:27082001

  15. Haldane-Hubbard Mott Insulator: From Tetrahedral Spin Crystal to Chiral Spin Liquid

    NASA Astrophysics Data System (ADS)

    Hickey, Ciarán; Cincio, Lukasz; Papić, Zlatko; Paramekanti, Arun

    2016-04-01

    Motivated by cold atom experiments on Chern insulators, we study the honeycomb lattice Haldane-Hubbard Mott insulator of spin-1 /2 fermions using exact diagonalization and density matrix renormalization group methods. We show that this model exhibits various chiral magnetic orders including a wide regime of triple-Q tetrahedral order. Incorporating third-neighbor hopping frustrates and ultimately melts this tetrahedral spin crystal. From analyzing the low energy spectrum, many-body Chern numbers, entanglement spectra, and modular matrices, we identify the molten state as a chiral spin liquid (CSL) with gapped semion excitations. We formulate and study the Chern-Simons-Higgs field theory of the exotic CSL-to-tetrahedral spin crystallization transition.

  16. Self-equilibrium and stability of regular truncated tetrahedral tensegrity structures

    NASA Astrophysics Data System (ADS)

    Zhang, J. Y.; Ohsaki, M.

    2012-10-01

    This paper presents analytical conditions of self-equilibrium and super-stability for the regular truncated tetrahedral tensegrity structures, nodes of which have one-to-one correspondence to the tetrahedral group. These conditions are presented in terms of force densities, by investigating the block-diagonalized force density matrix. The block-diagonalized force density matrix, with independent sub-matrices lying on its leading diagonal, is derived by making use of the tetrahedral symmetry via group representation theory. The condition for self-equilibrium is found by enforcing the force density matrix to have the necessary number of nullities, which is four for three-dimensional structures. The condition for super-stability is further presented by guaranteeing positive semi-definiteness of the force density matrix.

  17. Tetrahedral-mesh-based computational human phantom for fast Monte Carlo dose calculations

    NASA Astrophysics Data System (ADS)

    Yeom, Yeon Soo; Jeong, Jong Hwi; Han, Min Cheol; Kim, Chan Hyeong

    2014-06-01

    Although polygonal-surface computational human phantoms can address several critical limitations of conventional voxel phantoms, their Monte Carlo simulation speeds are much slower than those of voxel phantoms. In this study, we sought to overcome this problem by developing a new type of computational human phantom, a tetrahedral mesh phantom, by converting a polygonal surface phantom to a tetrahedral mesh geometry. The constructed phantom was implemented in the Geant4 Monte Carlo code to calculate organ doses as well as to measure computation speed, the values were then compared with those for the original polygonal surface phantom. It was found that using the tetrahedral mesh phantom significantly improved the computation speed by factors of between 150 and 832 considering all of the particles and simulated energies other than the low-energy neutrons (0.01 and 1 MeV), for which the improvement was less significant (17.2 and 8.8 times, respectively).

  18. AdS waves as exact solutions to quadratic gravity

    SciTech Connect

    Guellue, Ibrahim; Sisman, Tahsin Cagri; Tekin, Bayram; Guerses, Metin

    2011-04-15

    We give an exact solution of the quadratic gravity in D dimensions. The solution is a plane-fronted wave metric with a cosmological constant. This metric solves not only the full quadratic gravity field equations but also the linearized ones which include the linearized equations of the recently found critical gravity. A subset of the solutions change the asymptotic structure of the anti-de Sitter space due to their logarithmic behavior.

  19. Quadratic function approaching method for magnetotelluric soundingdata inversion

    SciTech Connect

    Liangjun, Yan; Wenbao, Hu; Zhang, Keni

    2004-04-05

    The quadratic function approaching method (QFAM) is introduced for magnetotelluric sounding (MT) data inversion. The method takes the advantage of that quadratic function has single extreme value, which avoids leading to an inversion solution for local minimum and ensures the solution for global minimization of an objective function. The method does not need calculation of sensitivity matrix and not require a strict initial earth model. Examples for synthetic data and field measurement data indicate that the proposed inversion method is effective.

  20. Dynamic mesh adaption for triangular and tetrahedral grids

    NASA Technical Reports Server (NTRS)

    Biswas, Rupak; Strawn, Roger

    1993-01-01

    The following topics are discussed: requirements for dynamic mesh adaption; linked-list data structure; edge-based data structure; adaptive-grid data structure; three types of element subdivision; mesh refinement; mesh coarsening; additional constraints for coarsening; anisotropic error indicator for edges; unstructured-grid Euler solver; inviscid 3-D wing; and mesh quality for solution-adaptive grids. The discussion is presented in viewgraph form.

  1. Dynamic electron scattering distinguishes dodecahedral and tetrahedral crystallographic sites in garnet structure

    NASA Astrophysics Data System (ADS)

    Lábár, János L.

    1999-07-01

    It is shown in this letter that, in contrast to the accepted belief in the literature, it is possible to determine if a minority component is located on the dodecahedral, octahedral, or tetrahedral sites in a garnet single crystal. This prediction of ours is based on dynamical Bloch-wave calculations and proved experimentally with x-ray measurements in a transmission electron microscope. The previous literature seemed to agree on the assumption that the dodecahedral and tetrahedral sites are indistinguishable from each other.

  2. Haldane-Hubbard Mott Insulator: From Tetrahedral Spin Crystal to Chiral Spin Liquid

    NASA Astrophysics Data System (ADS)

    Hickey, Ciaran; Cincio, Lukasz; Papic, Zlatko; Paramekanti, Arun

    Motivated by recent experimental realizations of artificial gauge fields in ultracold atoms, we study the honeycomb lattice Haldane-Hubbard Mott insulator of spin- 1 / 2 fermions using exact diagonalization and density matrix renormalization group methods. We show that this model exhibits various chiral magnetic orders including a wide regime of triple-Q tetrahedral order. Incorporating third-neighbor hopping frustrates and ultimately melts this tetrahedral spin crystal. From analyzing low energy spectra, many-body Chern numbers, entanglement spectra, and modular matrices, we identify the molten state as a chiral spin liquid with gapped semion excitations.

  3. Transformational part-count in layered octahedral-tetrahedral truss configurations

    NASA Technical Reports Server (NTRS)

    Lalvani, Haresh

    1990-01-01

    The number of component part (nodes, struts and panels) termed part count, is an important factor in the design, manufacture, and assembly of modular space structures. Part count expressions are presented for a variety of profiles derived from the layered octahedral-tetrahedral truss configuration. Referred to as the tetrahedral truss in the NASA projects, this specific geometry has been used in several missions. The general expressions presented here transforms to others as one profile changes to another. Such transformational part count relations provide a measure of flexibility and generality, and may be useful when dealing with a wider range of geometric configurations.

  4. Refining quadrilateral and brick element meshes

    SciTech Connect

    Schneiders, R.; Debye, J.

    1995-12-31

    We consider the problem of refining unstructured quadrilateral and brick element meshes. We present an algorithm which is a generalization of an algorithm developed by Cheng et. al. for structured quadrilateral element meshes. The problem is solved for the two-dimensional case. Concerning three dimensions we present a solution for some special cases and a general solution that introduces tetrahedral and pyramidal transition elements.

  5. Use of non-quadratic yield surfaces in design of optimal deep-draw blank geometry

    SciTech Connect

    Logan, R.W.

    1995-12-01

    Planar anisotropy in the deep-drawing of sheet can lead to the formation of ears in cylindrical cups and to undesirable metal flow in the blankholder in the general case. For design analysis purposes in non-linear finite-element codes, this anisotropy is characterized by the use of an appropriate yield surface which is then implemented into codes such as DYNA3D . The quadratic Hill yield surface offers a relatively straightforward implementation and can be formulated to be invariant to the coordinate system. Non-quadratic yield surfaces can provide more realistic strength or strain increment ratios, but they may not provide invariance and thus demand certain approximations. Forms due to Hosford and Badat et al. have been shown to more accurately address the earning phenomenon. in this work, use is made of these non-quadratic yield surfaces in order to determine the optimal blank shape for cups and other shapes using ferrous and other metal blank materials with planar anisotropy. The analyses are compared to previous experimental studies on non-uniform blank motion due to anisotropy and asymmetric geometry.

  6. Lattice thermal expansion for normal tetrahedral compound semiconductors

    SciTech Connect

    Omar, M.S. . E-mail: dr_m_s_omar@yahoo.com

    2007-02-15

    The cubic root of the deviation of the lattice thermal expansion from that of the expected value of diamond for group IV semiconductors, binary compounds of III-V and II-VI, as well as several ternary compounds from groups I-III-VI{sub 2}, II-IV-V{sub 2} and I-IV{sub 2}V{sub 3} semiconductors versus their bonding length are given straight lines. Their slopes were found to be 0.0256, 0.0210, 0.0170, 0.0259, 0.0196, and 0.02840 for the groups above, respectively. Depending on the valence electrons of the elements forming these groups, a formula was found to correlate all the values of the slopes mentioned above to that of group IV. This new formula which depends on the melting point and the bonding length as well as the number of valence electrons for the elements forming the compounds, will gives best calculated values for lattice thermal expansion for all compounds forming the groups mentioned above. An empirical relation is also found between the mean ionicity of the compounds forming the groups and their slopes mentioned above and that gave the mean ionicity for the compound CuGe{sub 2}P{sub 3} in the range of 0.442.

  7. A transient, quadratic nodal method for triangular-Z geometry

    SciTech Connect

    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.

  8. Numerical approximation for the infinite-dimensional discrete-time optimal linear-quadratic regulator problem

    NASA Technical Reports Server (NTRS)

    Gibson, J. S.; Rosen, I. G.

    1986-01-01

    An abstract approximation framework is developed for the finite and infinite time horizon discrete-time linear-quadratic regulator problem for systems whose state dynamics are described by a linear semigroup of operators on an infinite dimensional Hilbert space. The schemes included the framework yield finite dimensional approximations to the linear state feedback gains which determine the optimal control law. Convergence arguments are given. Examples involving hereditary and parabolic systems and the vibration of a flexible beam are considered. Spline-based finite element schemes for these classes of problems, together with numerical results, are presented and discussed.

  9. Automated Tetrahedral Mesh Generation for CFD Analysis of Aircraft in Conceptual Design

    NASA Technical Reports Server (NTRS)

    Ordaz, Irian; Li, Wu; Campbell, Richard L.

    2014-01-01

    The paper introduces an automation process of generating a tetrahedral mesh for computational fluid dynamics (CFD) analysis of aircraft configurations in early conceptual design. The method was developed for CFD-based sonic boom analysis of supersonic configurations, but can be applied to aerodynamic analysis of aircraft configurations in any flight regime.

  10. Regioselective placement of alkanethiolate domains on tetrahedral and octahedral gold nanocrystals.

    PubMed

    Wang, Yifeng; Zeiri, Offer; Meshi, Louisa; Stellacci, Francesco; Weinstock, Ira A

    2012-10-01

    Electrostatically stabilized monolayer shells of metal-oxide cluster anions (polyoxometalates, or POMs) on the surfaces of ca. 8 nm tetrahedral and octahedral gold nanocrystals regioselectively direct water-soluble alkanethiolate ligands to the corners and edges of the gold polyhedra. PMID:22918232

  11. Multi-Dimensional, Inviscid Flux Reconstruction for Simulation of Hypersonic Heating on Tetrahedral Grids

    NASA Technical Reports Server (NTRS)

    Gnoffo, Peter A.

    2009-01-01

    The quality of simulated hypersonic stagnation region heating on tetrahedral meshes is investigated by using a three-dimensional, upwind reconstruction algorithm for the inviscid flux vector. Two test problems are investigated: hypersonic flow over a three-dimensional cylinder with special attention to the uniformity of the solution in the spanwise direction and hypersonic flow over a three-dimensional sphere. The tetrahedral cells used in the simulation are derived from a structured grid where cell faces are bisected across the diagonal resulting in a consistent pattern of diagonals running in a biased direction across the otherwise symmetric domain. This grid is known to accentuate problems in both shock capturing and stagnation region heating encountered with conventional, quasi-one-dimensional inviscid flux reconstruction algorithms. Therefore the test problem provides a sensitive test for algorithmic effects on heating. This investigation is believed to be unique in its focus on three-dimensional, rotated upwind schemes for the simulation of hypersonic heating on tetrahedral grids. This study attempts to fill the void left by the inability of conventional (quasi-one-dimensional) approaches to accurately simulate heating in a tetrahedral grid system. Results show significant improvement in spanwise uniformity of heating with some penalty of ringing at the captured shock. Issues with accuracy near the peak shear location are identified and require further study.

  12. Extended orientational correlation study for molecular liquids containing distorted tetrahedral molecules: application to methylene halides.

    PubMed

    Pothoczki, Szilvia; Temleitner, László; Pusztai, László

    2010-04-28

    The method of Rey [Rey, J. Chem. Phys. 126, 164506 (2007)] for describing how molecules orient toward each other in systems with perfect tetrahedral molecules is extended to the case of distorted tetrahedral molecules of c(2v) symmetry by means of introducing 28 subgroups. Additionally, the original analysis developed for perfect tetrahedral molecules, based on six groups, is adapted for molecules with imperfect tetrahedral shape. Deriving orientational correlation functions have been complemented with detailed analyses of dipole-dipole correlations. This way, (up to now) the most complete structure determination can be carried out for such molecular systems. In the present work, these calculations have been applied for particle configurations resulting from reverse Monte Carlo computer modeling. These particle arrangements are fully consistent with structure factors from neutron and x-ray diffraction measurements. Here we present a complex structural study for methylene halide (chloride, bromide, and iodide) molecular liquids, as possibly the best representative examples. It has been found that the most frequent orientations of molecules are of the 2:2 type over the entire distance range in these liquids. Focusing on the short range orientation, neighboring molecules turn toward each other with there "H,Y"-"H,Y" (Y: Cl, Br, I) edges, apart from CH(2)Cl(2) where the H,H-H,Cl arrangement is the most frequent. In general, the structure of methylene chloride appears to be different from the structure of the other two liquids. PMID:20441292

  13. Fostering Teacher Development to a Tetrahedral Orientation in the Teaching of Chemistry

    ERIC Educational Resources Information Center

    Lewthwaite, Brian; Wiebe, Rick

    2011-01-01

    This paper reports on the initial outcomes from the end of the fourth year of a 5 year research and professional development project to improve chemistry teaching among three cohorts of chemistry teachers in Manitoba, Canada. The project responds to a new curriculum introduction advocating a tetrahedral orientation (Mahaffy, "Journal of Chemical…

  14. An efficient method for energy levels calculation using full symmetry and exact kinetic energy operator: Tetrahedral molecules

    SciTech Connect

    Nikitin, A. V.; Rey, M.; Tyuterev, Vl. G.

    2015-03-07

    A simultaneous use of the full molecular symmetry and of an exact kinetic energy operator (KEO) is of key importance for accurate predictions of vibrational levels at a high energy range from a potential energy surface (PES). An efficient method that permits a fast convergence of variational calculations would allow iterative optimization of the PES parameters using experimental data. In this work, we propose such a method applied to tetrahedral AB{sub 4} molecules for which a use of high symmetry is crucial for vibrational calculations. A symmetry-adapted contracted angular basis set for six redundant angles is introduced. Simple formulas using this basis set for explicit calculation of the angular matrix elements of KEO and PES are reported. The symmetric form (six redundant angles) of vibrational KEO without the sin(q){sup −2} type singularity is derived. The efficient recursive algorithm based on the tensorial formalism is used for the calculation of vibrational matrix elements. A good basis set convergence for the calculations of vibrational levels of the CH{sub 4} molecule is demonstrated.

  15. An efficient method for energy levels calculation using full symmetry and exact kinetic energy operator: Tetrahedral molecules

    NASA Astrophysics Data System (ADS)

    Nikitin, A. V.; Rey, M.; Tyuterev, Vl. G.

    2015-03-01

    A simultaneous use of the full molecular symmetry and of an exact kinetic energy operator (KEO) is of key importance for accurate predictions of vibrational levels at a high energy range from a potential energy surface (PES). An efficient method that permits a fast convergence of variational calculations would allow iterative optimization of the PES parameters using experimental data. In this work, we propose such a method applied to tetrahedral AB4 molecules for which a use of high symmetry is crucial for vibrational calculations. A symmetry-adapted contracted angular basis set for six redundant angles is introduced. Simple formulas using this basis set for explicit calculation of the angular matrix elements of KEO and PES are reported. The symmetric form (six redundant angles) of vibrational KEO without the sin(q)-2 type singularity is derived. The efficient recursive algorithm based on the tensorial formalism is used for the calculation of vibrational matrix elements. A good basis set convergence for the calculations of vibrational levels of the CH4 molecule is demonstrated.

  16. Parallel Higher-order Finite Element Method for Accurate Field Computations in Wakefield and PIC Simulations

    SciTech Connect

    Candel, A.; Kabel, A.; Lee, L.; Li, Z.; Limborg, C.; Ng, C.; Prudencio, E.; Schussman, G.; Uplenchwar, R.; Ko, K.; /SLAC

    2009-06-19

    Over the past years, SLAC's Advanced Computations Department (ACD), under SciDAC sponsorship, has developed a suite of 3D (2D) parallel higher-order finite element (FE) codes, T3P (T2P) and Pic3P (Pic2P), aimed at accurate, large-scale simulation of wakefields and particle-field interactions in radio-frequency (RF) cavities of complex shape. The codes are built on the FE infrastructure that supports SLAC's frequency domain codes, Omega3P and S3P, to utilize conformal tetrahedral (triangular)meshes, higher-order basis functions and quadratic geometry approximation. For time integration, they adopt an unconditionally stable implicit scheme. Pic3P (Pic2P) extends T3P (T2P) to treat charged-particle dynamics self-consistently using the PIC (particle-in-cell) approach, the first such implementation on a conformal, unstructured grid using Whitney basis functions. Examples from applications to the International Linear Collider (ILC), Positron Electron Project-II (PEP-II), Linac Coherent Light Source (LCLS) and other accelerators will be presented to compare the accuracy and computational efficiency of these codes versus their counterparts using structured grids.

  17. On Volterra quadratic stochastic operators with continual state space

    SciTech Connect

    Ganikhodjaev, Nasir; Hamzah, Nur Zatul Akmar

    2015-05-15

    Let (X,F) be a measurable space, and S(X,F) be the set of all probability measures on (X,F) where X is a state space and F is σ - algebraon X. We consider a nonlinear transformation (quadratic stochastic operator) defined by (Vλ)(A) = ∫{sub X}∫{sub X}P(x,y,A)dλ(x)dλ(y), where P(x, y, A) is regarded as a function of two variables x and y with fixed A ∈ F . A quadratic stochastic operator V is called a regular, if for any initial measure the strong limit lim{sub n→∞} V{sup n }(λ) is exists. In this paper, we construct a family of quadratic stochastic operators defined on the segment X = [0,1] with Borel σ - algebra F on X , prove their regularity and show that the limit measure is a Dirac measure.

  18. On Volterra quadratic stochastic operators with continual state space

    NASA Astrophysics Data System (ADS)

    Ganikhodjaev, Nasir; Hamzah, Nur Zatul Akmar

    2015-05-01

    Let (X,F) be a measurable space, and S(X,F) be the set of all probability measures on (X,F) where X is a state space and F is σ - algebraon X. We consider a nonlinear transformation (quadratic stochastic operator) defined by (V λ )(A ) = ∫X ∫X P (x ,y ,A )d λ (x )d λ (y ), where P(x, y, A) is regarded as a function of two variables x and y with fixed A ∈ F . A quadratic stochastic operator V is called a regular, if for any initial measure the strong limit lim n →∞ Vn(λ ) is exists. In this paper, we construct a family of quadratic stochastic operators defined on the segment X = [0,1] with Borel σ - algebra F on X , prove their regularity and show that the limit measure is a Dirac measure.

  19. Symmetric quadratic Hamiltonians with pseudo-Hermitian matrix representation

    NASA Astrophysics Data System (ADS)

    Fernández, Francisco M.

    2016-06-01

    We prove that any symmetric Hamiltonian that is a quadratic function of the coordinates and momenta has a pseudo-Hermitian adjoint or regular matrix representation. The eigenvalues of the latter matrix are the natural frequencies of the Hamiltonian operator. When all the eigenvalues of the matrix are real, then the spectrum of the symmetric Hamiltonian is real and the operator is Hermitian. As illustrative examples we choose the quadratic Hamiltonians that model a pair of coupled resonators with balanced gain and loss, the electromagnetic self-force on an oscillating charged particle and an active LRC circuit.

  20. Analysis of integral controls in linear quadratic regulator design

    NASA Technical Reports Server (NTRS)

    Slater, G. L.

    1979-01-01

    The application of linear optimal control to the design of systems with integral control action on specified outputs is considered. Using integral terms in a quadratic performance index, an asymptotic analysis is used to determine the effect of variable quadratic weights on the eigenvalues and eigenvectors of the closed loop system. It is shown that for small integral terms the placement of integrator poles and gain calculation can be effectively decoupled from placement of the primary system eigenvalues. This technique is applied to the design of integral controls for a STOL aircraft outer loop guidance system.

  1. A Note on the Linearly and Quadratically Weighted Kappa Coefficients.

    PubMed

    Li, Pingke

    2016-09-01

    The linearly and quadratically weighted kappa coefficients are popular statistics in measuring inter-rater agreement on an ordinal scale. It has been recently demonstrated that the linearly weighted kappa is a weighted average of the kappa coefficients of the embedded 2 by 2 agreement matrices, while the quadratically weighted kappa is insensitive to the agreement matrices that are row or column reflection symmetric. A rank-one matrix decomposition approach to the weighting schemes is presented in this note such that these phenomena can be demonstrated in a concise manner. PMID:27246436

  2. Cr(1/3)Zr₂P₃O₁₂ with unusual tetrahedral coordination of Cr(III): peculiarities of the formation, thermal stability and application as a pigment.

    PubMed

    Gorodylova, Nataliia; Kosinová, Veronika; Šulcová, Petra; Bělina, Petr; Vlček, Milan

    2014-11-01

    All the known chromium(III) NASICON-related phosphates are considered to be solid solutions. In these compounds chromium atoms share their position in the basic framework of the crystal lattice with other structure forming elements such as zirconium. In our study, we have hypothesised a completely new way of structural organisation of the chromium(III) zirconium(IV) NASICON framework, consisting in the distribution of chromium over the charge-compensating atom sites with tetrahedral oxygen coordination. The possibility of formation of the corresponding phosphate, Cr(1/3)Zr2P3O12, was studied using a classical ceramic route and a sol-gel method. Structural affiliation of the obtained pure phase product was studied using XRD analysis. The results confirmed that the Cr(1/3)Zr2P3O12 phosphate belongs to monoclinic SW-subtype of the NASICON family. In this structure, chromium atoms occupy charge-compensating sites with a strongly distorted tetrahedral oxygen environment. To the best of our knowledge, it is the first example of tetrahedral coordination of chromium(III) in phosphates. Along with the unusual crystallographic characteristics of chromium, special attention in this paper is devoted to the thermal stability of this phosphate and to its performance as an inorganic pigment. The sample was characterised by heating microscopy and DTA study, particle size distribution analysis, and IR- and VIS-spectroscopy. The stability of the obtained powder in a glaze environment, its colouring performance and lightfastness are discussed as well. PMID:25189199

  3. Formation of pyramid elements for hexahedra to tetrahedra transitions

    SciTech Connect

    OWEN,STEVEN J.; SAIGAL,SUNIL

    2000-02-24

    New algorithms are proposed for the modification of a mixed hexahedra-tetrahedra element mesh to maintain compatibility by the insertion of pyramid elements. Several methods for generation of the pyramids are presented involving local tetrahedral transformations and/or node insertion near the hex/tet interface. Local smoothing and topological operations improve the quality of the transition region. Results show superior performance of the resulting elements in a commercial finite element code over non-conforming interface conditions.

  4. A quadratic weight selection algorithm. [for optimal flight control

    NASA Technical Reports Server (NTRS)

    Broussard, J. R.

    1981-01-01

    A new numerical algorithm is presented which determines a positive semi-definite state weighting matrix in the linear-quadratic optimal control design problem. The algorithm chooses the weighting matrix by placing closed-loop eigenvalues and eigenvectors near desired locations using optimal feedback gains. A simplified flight control design example is used to illustrate the algorithms capabilities.

  5. Tuning a fuzzy controller using quadratic response surfaces

    NASA Technical Reports Server (NTRS)

    Schott, Brian; Whalen, Thomas

    1992-01-01

    Response surface methodology, an alternative method to traditional tuning of a fuzzy controller, is described. An example based on a simulated inverted pendulum 'plant' shows that with (only) 15 trial runs, the controller can be calibrated using a quadratic form to approximate the response surface.

  6. A Version of Quadratic Regression with Interpretable Parameters.

    ERIC Educational Resources Information Center

    Cudeck, Robert; du Toit, Stephen H. C.

    2002-01-01

    Suggests an alternative form of the quadratic model that has the same expectation function of the original model but has the useful feature that its parameters are interpretable. Provides examples of a simple regression problem and a nonlinear mixed-effects model. (SLD)

  7. Solving the Quadratic Capacitated Facilities Location Problem by Computer.

    ERIC Educational Resources Information Center

    Cote, Leon C.; Smith, Wayland P.

    Several computer programs were developed to solve various versions of the quadratic capacitated facilities location problem. Matrices, which represent various business costs, are defined for the factors of sites, facilities, customers, commodities, and production units. The objective of the program is to find an optimization matrix for the lowest…

  8. Confidence set interference with a prior quadratic bound. [in geophysics

    NASA Technical Reports Server (NTRS)

    Backus, George E.

    1989-01-01

    Neyman's (1937) theory of confidence sets is developed as a replacement for Bayesian interference (BI) and stochastic inversion (SI) when the prior information is a hard quadratic bound. It is recommended that BI and SI be replaced by confidence set interference (CSI) only in certain circumstances. The geomagnetic problem is used to illustrate the general theory of CSI.

  9. Finding Optimal Gains In Linear-Quadratic Control Problems

    NASA Technical Reports Server (NTRS)

    Milman, Mark H.; Scheid, Robert E., Jr.

    1990-01-01

    Analytical method based on Volterra factorization leads to new approximations for optimal control gains in finite-time linear-quadratic control problem of system having infinite number of dimensions. Circumvents need to analyze and solve Riccati equations and provides more transparent connection between dynamics of system and optimal gain.

  10. Radar Rainfall Estimation using a Quadratic Z-R equation

    NASA Astrophysics Data System (ADS)

    Hall, Will; Rico-Ramirez, Miguel Angel; Kramer, Stefan

    2016-04-01

    The aim of this work is to test a method that enables the input of event based drop size distributions to alter a quadratic reflectivity (Z) to rainfall (R) equation that is limited by fixed upper and lower points. Results will be compared to the Marshall-Palmer Z-R relation outputs and validated by a network of gauges and a single polarisation weather radar located close to Essen, Germany. The time window over which the drop size distribution measurements will be collected is varied to note any effect on the generated quadratic Z-R relation. The new quadratic algorithm shows some distinct improvement over the Marshall-Palmer relationship through multiple events. The inclusion of a minimum number of Z-R points helped to decrease the associated error by defaulting back to the Marshall-Palmer equation if the limit was not reached. More research will be done to discover why the quadratic performs poorly in some events as there appears to be little correlation between number of drops or mean rainfall amount and the associated error. In some cases it seems the spatial distribution of the disdrometers has a significant effect as a large percentage of the rain bands pass to the north of two of the three disdrometers, frequently in a slightly north-easterly direction. However during widespread precipitation events the new algorithm works very well with reductions compared to the Marshall-Palmer relation.

  11. Analysis of Quadratic Diophantine Equations with Fibonacci Number Solutions

    ERIC Educational Resources Information Center

    Leyendekkers, J. V.; Shannon, A. G.

    2004-01-01

    An analysis is made of the role of Fibonacci numbers in some quadratic Diophantine equations. A general solution is obtained for finding factors in sums of Fibonacci numbers. Interpretation of the results is facilitated by the use of a modular ring which also permits extension of the analysis.

  12. Quadratic invariants for discrete clusters of weakly interacting waves

    NASA Astrophysics Data System (ADS)

    Harper, Katie L.; Bustamante, Miguel D.; Nazarenko, Sergey V.

    2013-06-01

    We consider discrete clusters of quasi-resonant triads arising from a Hamiltonian three-wave equation. A cluster consists of N modes forming a total of M connected triads. We investigate the problem of constructing a functionally independent set of quadratic constants of motion. We show that this problem is equivalent to an underlying basic linear problem, consisting of finding the null space of a rectangular M × N matrix {A} with entries 1, -1 and 0. In particular, we prove that the number of independent quadratic invariants is equal to J ≡ N - M* ⩾ N - M, where M* is the number of linearly independent rows in {A}. Thus, the problem of finding all independent quadratic invariants is reduced to a linear algebra problem in the Hamiltonian case. We establish that the properties of the quadratic invariants (e.g., locality) are related to the topological properties of the clusters (e.g., types of linkage). To do so, we formulate an algorithm for decomposing large clusters into smaller ones and show how various invariants are related to certain parts of a cluster, including the basic structures leading to M* < M. We illustrate our findings by presenting examples from the Charney-Hasegawa-Mima wave model, and by showing a classification of small (up to three-triad) clusters.

  13. Selfsimilarity in a Class of Quadratic-Quasiperiodic Chains

    NASA Astrophysics Data System (ADS)

    Kaneko, Masanobu; Odagaki, Takashi

    1993-04-01

    We prove that quasiperiodic chains associated with a class of quadratic irrational numbers have an inflation symmetry and can be generated from a regular chain by a hyperinflation. We devise the explicit method to find the hyperinflation symmetry and discuss the properties of such a class of quasiperiodic sequences.

  14. Unravelling Student Challenges with Quadratics: A Cognitive Approach

    ERIC Educational Resources Information Center

    Kotsopoulos, Donna

    2007-01-01

    The author's secondary school mathematics students have often reported to her that quadratic relations are one of the most conceptually challenging aspects of the high school curriculum. From her own classroom experiences there seemed to be several aspects to the students' challenges. Many students, even in their early secondary education, have…

  15. Quadratic Expressions by Means of "Summing All the Matchsticks"

    ERIC Educational Resources Information Center

    Gierdien, M. Faaiz

    2012-01-01

    This note presents demonstrations of quadratic expressions that come about when particular problems are posed with respect to matchsticks that form regular triangles, squares, pentagons and so on. Usually when such "matchstick" problems are used as ways to foster algebraic thinking, the expressions for the number of matchstick quantities are…

  16. Visualising the Complex Roots of Quadratic Equations with Real Coefficients

    ERIC Educational Resources Information Center

    Bardell, Nicholas S.

    2012-01-01

    The roots of the general quadratic equation y = ax[superscript 2] + bx + c (real a, b, c) are known to occur in the following sets: (i) real and distinct; (ii) real and coincident; and (iii) a complex conjugate pair. Case (iii), which provides the focus for this investigation, can only occur when the values of the real coefficients a, b, and c are…

  17. A Unified Approach to Teaching Quadratic and Cubic Equations.

    ERIC Educational Resources Information Center

    Ward, A. J. B.

    2003-01-01

    Presents a simple method for teaching the algebraic solution of cubic equations via completion of the cube. Shows that this method is readily accepted by students already familiar with completion of the square as a method for quadratic equations. (Author/KHR)

  18. Exploration of tetrahedral structures in silicate cathodes using a motif-network scheme

    SciTech Connect

    Zhao, Xin; Wu, Shunqing; Lv, Xiaobao; Nguyen, Manh Cuong; Wang, Cai -Zhuang; Lin, Zijing; Zhu, Zi -Zhong; Ho, Kai -Ming

    2015-10-26

    Using a motif-network search scheme, we studied the tetrahedral structures of the dilithium/disodium transition metal orthosilicates A2MSiO4 with A = Li or Na and M = Mn, Fe or Co. In addition to finding all previously reported structures, we discovered many other different tetrahedral-network-based crystal structures which are highly degenerate in energy. In addition, these structures can be classified into structures with 1D, 2D and 3D M-Si-O frameworks. A clear trend of the structural preference in different systems was revealed and possible indicators that affect the structure stabilities were introduced. For the case of Na systems which have been much less investigated in the literature relative to the Li systems, we predicted their ground state structures and found evidence for the existence of new structural motifs.

  19. Flexibility of C3h -Symmetrical Linkers in Tris-oligonucleotide-Based Tetrahedral Scaffolds.

    PubMed

    Panagiotidis, Christos; Kath-Schorr, Stephanie; von Kiedrowski, Günter

    2016-02-01

    Flexibility of tris-oligonucleotides is determined by the length of their connecting hydrocarbon chains. Tris-oligonucleotides are branched DNA building blocks with three oligonucleotide arms attached to a C3h -symmetrical linker core at these chains. Four tris-oligonucleotides hybridise into a tetrahedral nanocage by sequence-determined self-assembly. The influence of methylene, ethylene and propylene chains was studied by synthesising sets of tris-oligonucleotides and analysing the relative stability of the hybridisation products against digestion by mung bean nuclease by using gel electrophoresis. Linkers with ethylene chains showed sufficient flexibility, whereas methylene-chain linkers were too rigid. Tris-oligonucleotides based on the latter still formed tetrahedral scaffolds in intermixing experiments with linkers of higher flexibility. Thus, a new generation of versatile isocyanurate-based linkers was established. PMID:26593127

  20. Exploration of tetrahedral structures in silicate cathodes using a motif-network scheme

    DOE PAGESBeta

    Zhao, Xin; Wu, Shunqing; Lv, Xiaobao; Nguyen, Manh Cuong; Wang, Cai -Zhuang; Lin, Zijing; Zhu, Zi -Zhong; Ho, Kai -Ming

    2015-10-26

    Using a motif-network search scheme, we studied the tetrahedral structures of the dilithium/disodium transition metal orthosilicates A2MSiO4 with A = Li or Na and M = Mn, Fe or Co. In addition to finding all previously reported structures, we discovered many other different tetrahedral-network-based crystal structures which are highly degenerate in energy. In addition, these structures can be classified into structures with 1D, 2D and 3D M-Si-O frameworks. A clear trend of the structural preference in different systems was revealed and possible indicators that affect the structure stabilities were introduced. For the case of Na systems which have been muchmore » less investigated in the literature relative to the Li systems, we predicted their ground state structures and found evidence for the existence of new structural motifs.« less

  1. Atomic Scale Picture of the Ion Conduction Mechanism in Tetrahedral Network of Lanthanum Barium Gallate

    SciTech Connect

    Jalarvo, Niina H; Gourdon, Olivier; Bi, Zhonghe; Gout, Delphine J; Ohl, Michael E; Paranthaman, Mariappan Parans

    2013-01-01

    Combined experimental study of impedance spectroscopy, neutron powder diffraction and quasielastic neutron scattering was performed to shed light into the atomic scale ion migration processes in proton and oxide ion conductor; La0.8Ba1.2GaO3.9 . This material consist of tetrahedral GaO4 units, which are rather flexible and rocking motion of these units promotes the ionic migration process. The oxide ion (vacancy) conduction takes place on channels along c axis, involving a single elementary step, which occurs between adjacent tetrahedron (inter-tetrahedron jump). The proton conduction mechanism consists of intra-tetrahedron and inter-tetrahedron elementary processes. The intra-tetrahedron proton transport is the rate-limiting process, with activation energy of 0.44 eV. The rocking motion of the GaO4 tetrahedron aids the inter-tetrahedral proton transport, which has the activation energy of 0.068 eV.

  2. Near-field testing of the 5-meter model of the tetrahedral truss antenna

    NASA Technical Reports Server (NTRS)

    Kefauver, Neill; Cencich, Tom; Osborn, Jim; Osmanski, J. T.

    1986-01-01

    This report documents the technical results from near-field testing of the General Dynamics 5-meter model of the tetrahedral truss antenna at the Martin Marietta Denver Aerospace facility. A 5-meter square side of the tetrahedral served as the perimeter of the antenna, and a mesh surface and extensive surface contouring cord network was used to create a parabolic aperture shape to within an rms accuracy of 30 mils or better. Pattern measurements were made with offset feed systems radiating at frequencies of 7.73, 11.60, 2.27, and 4.26 (all in GHz). This report discusses the method of collecting the data, system measurement accuracy, the test data compiled, and diagostics and isolation of causes of pattern results. The technique of using near-field phase for measuring surface mechanical tolerances is included. Detailed far field antenna patterns and their implications are provided for all tests conducted.

  3. Exploration of tetrahedral structures in silicate cathodes using a motif-network scheme

    PubMed Central

    Zhao, Xin; Wu, Shunqing; Lv, Xiaobao; Nguyen, Manh Cuong; Wang, Cai-Zhuang; Lin, Zijing; Zhu, Zi-Zhong; Ho, Kai-Ming

    2015-01-01

    Using a motif-network search scheme, we studied the tetrahedral structures of the dilithium/disodium transition metal orthosilicates A2MSiO4 with A = Li or Na and M = Mn, Fe or Co. In addition to finding all previously reported structures, we discovered many other different tetrahedral-network-based crystal structures which are highly degenerate in energy. These structures can be classified into structures with 1D, 2D and 3D M-Si-O frameworks. A clear trend of the structural preference in different systems was revealed and possible indicators that affect the structure stabilities were introduced. For the case of Na systems which have been much less investigated in the literature relative to the Li systems, we predicted their ground state structures and found evidence for the existence of new structural motifs. PMID:26497381

  4. A first collision source method for ATTILA, an unstructured tetrahedral mesh discrete ordinates code

    SciTech Connect

    Wareing, T.A.; Morel, J.E.; Parsons, D.K.

    1998-12-01

    A semi-analytic first collision source method is developed for the transport code, ATTILA, a three-dimensional, unstructured tetrahedral mesh, discrete-ordinates code. This first collision source method is intended to mitigate ray effects due to point sources. The method is third-order accurate, which is the same order of accuracy as the linear-discontinuous spatial differencing scheme used in ATTILA. Numerical results are provided to demonstrate the accuracy and efficiency of the first collision source method.

  5. A Reactor Pressure Vessel Dosimetry Calculation Using ATTILA, An Unstructured Tetrahedral Mesh Discrete-Ordinates Code

    SciTech Connect

    Wareing, T.A.; Parsons, D.K.; Pautz, S.

    1997-12-31

    Recently, a new state-of-the-art discrete-ordinates code, ATTILA, was developed. ATTILA provides the capabilities to solve geometrically complex 3-D transport problems by using an unstructured tetrahedral mesh. In this paper we describe the application of ATTILA to a 3-D reactor pressure vessel dosimetry problem. We provide numerical results from ATTILA and the Monte Carlo code, MCNP. The results demonstrate the effectiveness and efficiency of ATTILA for such calculations.

  6. Theory of Andreev reflection spectroscopy for tetrahedral and non-unitary superconductors.

    PubMed

    Bohloul, S; Curnoe, S H

    2016-02-01

    A general formula for point contact conductance from a normal metal tip into a superconductor is derived using the Blonder-Tinkham-Klapwijk theory of Andreev reflection, with special emphasis on non-unitary superconductors. The results of a comprehensive set of conductance spectrum calculations are presented: all symmetry-allowed gap functions for superconductors with tetrahedral symmetry, such as PrOs4Sb12, are considered, including several non-unitary cases. PMID:26750247

  7. Theory of Andreev reflection spectroscopy for tetrahedral and non-unitary superconductors

    NASA Astrophysics Data System (ADS)

    Bohloul, S.; Curnoe, S. H.

    2016-02-01

    A general formula for point contact conductance from a normal metal tip into a superconductor is derived using the Blonder-Tinkham-Klapwijk theory of Andreev reflection, with special emphasis on non-unitary superconductors. The results of a comprehensive set of conductance spectrum calculations are presented: all symmetry-allowed gap functions for superconductors with tetrahedral symmetry, such as PrOs4Sb12, are considered, including several non-unitary cases.

  8. A FAST ITERATIVE METHOD FOR SOLVING THE EIKONAL EQUATION ON TETRAHEDRAL DOMAINS

    PubMed Central

    Fu, Zhisong; Kirby, Robert M.; Whitaker, Ross T.

    2014-01-01

    Generating numerical solutions to the eikonal equation and its many variations has a broad range of applications in both the natural and computational sciences. Efficient solvers on cutting-edge, parallel architectures require new algorithms that may not be theoretically optimal, but that are designed to allow asynchronous solution updates and have limited memory access patterns. This paper presents a parallel algorithm for solving the eikonal equation on fully unstructured tetrahedral meshes. The method is appropriate for the type of fine-grained parallelism found on modern massively-SIMD architectures such as graphics processors and takes into account the particular constraints and capabilities of these computing platforms. This work builds on previous work for solving these equations on triangle meshes; in this paper we adapt and extend previous two-dimensional strategies to accommodate three-dimensional, unstructured, tetrahedralized domains. These new developments include a local update strategy with data compaction for tetrahedral meshes that provides solutions on both serial and parallel architectures, with a generalization to inhomogeneous, anisotropic speed functions. We also propose two new update schemes, specialized to mitigate the natural data increase observed when moving to three dimensions, and the data structures necessary for efficiently mapping data to parallel SIMD processors in a way that maintains computational density. Finally, we present descriptions of the implementations for a single CPU, as well as multicore CPUs with shared memory and SIMD architectures, with comparative results against state-of-the-art eikonal solvers. PMID:25221418

  9. Tetrahedral Finite-Volume Solutions to the Navier-Stokes Equations on Complex Configurations

    NASA Technical Reports Server (NTRS)

    Frink, Neal T.; Pirzadeh, Shahyar Z.

    1998-01-01

    A review of the algorithmic features and capabilities of the unstructured-grid flow solver USM3Dns is presented. This code, along with the tetrahedral grid generator, VGRIDns, is being extensively used throughout the U.S. for solving the Euler and Navier-Stokes equations on complex aerodynamic problems. Spatial discretization is accomplished by a tetrahedral cell-centered finite-volume formulation using Roe's upwind flux difference splitting. The fluxes are limited by either a Superbee or MinMod limiter. Solution reconstruction within the tetrahedral cells is accomplished with a simple, but novel, multidimensional analytical formula. Time is advanced by an implicit backward-Euler time-stepping scheme. Flow turbulence effects are modeled by the Spalart-Allmaras one-equation model, which is coupled with a wall function to reduce the number of cells in the near-wall region of the boundary layer. The issues of accuracy and robustness of USM3Dns Navier-Stokes capabilities are addressed for a flat-plate boundary layer, and a full F-16 aircraft with external stores at transonic speed.

  10. Simulation of Stagnation Region Heating in Hypersonic Flow on Tetrahedral Grids

    NASA Technical Reports Server (NTRS)

    Gnoffo, Peter A.

    2007-01-01

    Hypersonic flow simulations using the node based, unstructured grid code FUN3D are presented. Applications include simple (cylinder) and complex (towed ballute) configurations. Emphasis throughout is on computation of stagnation region heating in hypersonic flow on tetrahedral grids. Hypersonic flow over a cylinder provides a simple test problem for exposing any flaws in a simulation algorithm with regard to its ability to compute accurate heating on such grids. Such flaws predominantly derive from the quality of the captured shock. The importance of pure tetrahedral formulations are discussed. Algorithm adjustments for the baseline Roe / Symmetric, Total-Variation-Diminishing (STVD) formulation to deal with simulation accuracy are presented. Formulations of surface normal gradients to compute heating and diffusion to the surface as needed for a radiative equilibrium wall boundary condition and finite catalytic wall boundary in the node-based unstructured environment are developed. A satisfactory resolution of the heating problem on tetrahedral grids is not realized here; however, a definition of a test problem, and discussion of observed algorithm behaviors to date are presented in order to promote further research on this important problem.