Approximation theory for LQG (Linear-Quadratic-Gaussian) optimal control of flexible structures

NASA Technical Reports Server (NTRS)

Gibson, J. S.; Adamian, A.

1988-01-01

An approximation theory is presented for the LQG (Linear-Quadratic-Gaussian) optimal control problem for flexible structures whose distributed models have bounded input and output operators. The main purpose of the theory is to guide the design of finite dimensional compensators that approximate closely the optimal compensator. The optimal LQG problem separates into an optimal linear-quadratic regulator problem and an optimal state estimation problem. The solution of the former problem lies in the solution to an infinite dimensional Riccati operator equation. The approximation scheme approximates the infinite dimensional LQG problem with a sequence of finite dimensional LQG problems defined for a sequence of finite dimensional, usually finite element or modal, approximations of the distributed model of the structure. Two Riccati matrix equations determine the solution to each approximating problem. The finite dimensional equations for numerical approximation are developed, including formulas for converting matrix control and estimator gains to their functional representation to allow comparison of gains based on different orders of approximation. Convergence of the approximating control and estimator gains and of the corresponding finite dimensional compensators is studied. Also, convergence and stability of the closed-loop systems produced with the finite dimensional compensators are discussed. The convergence theory is based on the convergence of the solutions of the finite dimensional Riccati equations to the solutions of the infinite dimensional Riccati equations. A numerical example with a flexible beam, a rotating rigid body, and a lumped mass is given.

A Biogeography-Based Optimization Algorithm Hybridized with Tabu Search for the Quadratic Assignment Problem

PubMed Central

Lim, Wee Loon; Wibowo, Antoni; Desa, Mohammad Ishak; Haron, Habibollah

2016-01-01

The quadratic assignment problem (QAP) is an NP-hard combinatorial optimization problem with a wide variety of applications. Biogeography-based optimization (BBO), a relatively new optimization technique based on the biogeography concept, uses the idea of migration strategy of species to derive algorithm for solving optimization problems. It has been shown that BBO provides performance on a par with other optimization methods. A classical BBO algorithm employs the mutation operator as its diversification strategy. However, this process will often ruin the quality of solutions in QAP. In this paper, we propose a hybrid technique to overcome the weakness of classical BBO algorithm to solve QAP, by replacing the mutation operator with a tabu search procedure. Our experiments using the benchmark instances from QAPLIB show that the proposed hybrid method is able to find good solutions for them within reasonable computational times. Out of 61 benchmark instances tested, the proposed method is able to obtain the best known solutions for 57 of them. PMID:26819585

A Biogeography-Based Optimization Algorithm Hybridized with Tabu Search for the Quadratic Assignment Problem.

PubMed

Lim, Wee Loon; Wibowo, Antoni; Desa, Mohammad Ishak; Haron, Habibollah

2016-01-01

The quadratic assignment problem (QAP) is an NP-hard combinatorial optimization problem with a wide variety of applications. Biogeography-based optimization (BBO), a relatively new optimization technique based on the biogeography concept, uses the idea of migration strategy of species to derive algorithm for solving optimization problems. It has been shown that BBO provides performance on a par with other optimization methods. A classical BBO algorithm employs the mutation operator as its diversification strategy. However, this process will often ruin the quality of solutions in QAP. In this paper, we propose a hybrid technique to overcome the weakness of classical BBO algorithm to solve QAP, by replacing the mutation operator with a tabu search procedure. Our experiments using the benchmark instances from QAPLIB show that the proposed hybrid method is able to find good solutions for them within reasonable computational times. Out of 61 benchmark instances tested, the proposed method is able to obtain the best known solutions for 57 of them.

Sustainable Odds

NASA Astrophysics Data System (ADS)

Smith, L. A.

2016-12-01

While probability forecasting has many philosophical and mathematical attractions, it is something of a dishonest nonsense if acting on such forecasts is expected to lead to rapid ruin. Model-based probabilities, when interpreted as actionable, are shown to lead to the rapid ruin of a cooperative entity offering odds interpreting the probability forecasts at face value. Arguably, these odds would not be considered "fair", but inasmuch as some definitions of "fair odds" include this case, this presentation will focus on "sustainable odds": Odds which are not expected to lead to the rapid ruin of the cooperative under the assumption that those placing bets have no information beyond that available to the forecast system. It is argued that sustainable odds will not correspond to probabilities outside the Perfect Model Scenario, that the "implied probabilities" determined from sustainable odds will always sum to more than one, and that the excess of this sum over one reflects the skill of the forecast system, being a quantitative measure of structural model error.

User's manual for two dimensional FDTD version TEA and TMA codes for scattering from frequency-independent dielectic materials

NASA Technical Reports Server (NTRS)

Beggs, John H.; Luebbers, Raymond J.; Kunz, Karl S.

1991-01-01

The Penn State Finite Difference Time Domain Electromagnetic Scattering Code Versions TEA and TMA are two dimensional numerical electromagnetic scattering codes based upon the Finite Difference Time Domain Technique (FDTD) first proposed by Yee in 1966. The supplied version of the codes are two versions of our current two dimensional FDTD code set. This manual provides a description of the codes and corresponding results for the default scattering problem. The manual is organized into eleven sections: introduction, Version TEA and TMA code capabilities, a brief description of the default scattering geometry, a brief description of each subroutine, a description of the include files (TEACOM.FOR TMACOM.FOR), a section briefly discussing scattering width computations, a section discussing the scattering results, a sample problem set section, a new problem checklist, references and figure titles.

Filtering techniques for efficient inversion of two-dimensional Nuclear Magnetic Resonance data

NASA Astrophysics Data System (ADS)

Bortolotti, V.; Brizi, L.; Fantazzini, P.; Landi, G.; Zama, F.

2017-10-01

The inversion of two-dimensional Nuclear Magnetic Resonance (NMR) data requires the solution of a first kind Fredholm integral equation with a two-dimensional tensor product kernel and lower bound constraints. For the solution of this ill-posed inverse problem, the recently presented 2DUPEN algorithm [V. Bortolotti et al., Inverse Problems, 33(1), 2016] uses multiparameter Tikhonov regularization with automatic choice of the regularization parameters. In this work, I2DUPEN, an improved version of 2DUPEN that implements Mean Windowing and Singular Value Decomposition filters, is deeply tested. The reconstruction problem with filtered data is formulated as a compressed weighted least squares problem with multi-parameter Tikhonov regularization. Results on synthetic and real 2D NMR data are presented with the main purpose to deeper analyze the separate and combined effects of these filtering techniques on the reconstructed 2D distribution.

Analysis of the Hessian for Aerodynamic Optimization: Inviscid Flow

NASA Technical Reports Server (NTRS)

Arian, Eyal; Ta'asan, Shlomo

1996-01-01

In this paper we analyze inviscid aerodynamic shape optimization problems governed by the full potential and the Euler equations in two and three dimensions. The analysis indicates that minimization of pressure dependent cost functions results in Hessians whose eigenvalue distributions are identical for the full potential and the Euler equations. However the optimization problems in two and three dimensions are inherently different. While the two dimensional optimization problems are well-posed the three dimensional ones are ill-posed. Oscillations in the shape up to the smallest scale allowed by the design space can develop in the direction perpendicular to the flow, implying that a regularization is required. A natural choice of such a regularization is derived. The analysis also gives an estimate of the Hessian's condition number which implies that the problems at hand are ill-conditioned. Infinite dimensional approximations for the Hessians are constructed and preconditioners for gradient based methods are derived from these approximate Hessians.

Confined One Dimensional Harmonic Oscillator as a Two-Mode System

DOE Office of Scientific and Technical Information (OSTI.GOV)

Gueorguiev, V G; Rau, A P; Draayer, J P

2005-07-11

The one-dimensional harmonic oscillator in a box problem is possibly the simplest example of a two-mode system. This system has two exactly solvable limits, the harmonic oscillator and a particle in a (one-dimensional) box. Each of the two limits has a characteristic spectral structure describing the two different excitation modes of the system. Near each of these limits, one can use perturbation theory to achieve an accurate description of the eigenstates. Away from the exact limits, however, one has to carry out a matrix diagonalization because the basis-state mixing that occurs is typically too large to be reproduced in anymore » other way. An alternative to casting the problem in terms of one or the other basis set consists of using an ''oblique'' basis that uses both sets. Through a study of this alternative in this one-dimensional problem, we are able to illustrate practical solutions and infer the applicability of the concept for more complex systems, such as in the study of complex nuclei where oblique-basis calculations have been successful.« less

Isolation and characterization of melanopsin and pinopsin expression within photoreceptive sites of reptiles

NASA Astrophysics Data System (ADS)

Frigato, Elena; Vallone, Daniela; Bertolucci, Cristiano; Foulkes, Nicholas S.

2006-08-01

Non-mammalian vertebrates have multiple extraocular photoreceptors, mainly localised in the pineal complex and the brain, to mediate irradiance detection. In this study, we report the full-length cDNA cloning of ruin lizard melanopsin and pinopsin. The high level of identity with opsins in both the transmembrane regions, where the chromophore binding site is located, and the intracellular loops, where the G-proteins interact, suggests that both melanopsin and pinopsin should be able to generate a stable photopigment, capable of triggering a transduction cascade mediated by G-proteins. Phylogenetic analysis showed that both opsins are located on the expected branches of the corresponding sequences of ortholog proteins. Subsequently, using RT-PCR and RPA analysis, we verified the expression of ruin lizard melanopsin and pinopsin in directly photosensitive organs, such as the lateral eye, brain, pineal gland and parietal eye. Melanopsin expression was detected in the lateral eye and all major regions of the brain. However, different from the situation in Xenopus and chicken, melanopsin is not expressed in the ruin lizard pineal. Pinopsin mRNA expression was only detected in the pineal complex. As a result of their phylogenetic position and ecology, reptiles provide the circadian field with some of the most interesting models for understanding the evolution of the vertebrate circadian timing system and its response to light. This characterization of melanopsin and pinopsin expression in the ruin lizard will be important for future studies aimed at understanding the molecular basis of circadian light detection in reptiles.

New Methodologies for the Documentation of Fortified Architecture in the State of Ruins

NASA Astrophysics Data System (ADS)

Fallavollita, F.; Ugolini, A.

2017-05-01

Fortresses and castles are important symbols of social and cultural identity providing tangible evidence of cultural unity in Europe. They are items for which it is always difficult to outline a credible prospect of reuse, their old raison d'être- namely the military, political and economic purposes for which they were built- having been lost. In recent years a Research Unit of the University of Bologna composed of architects from different disciplines has conducted a series of studies on fortified heritage in the Emilia Romagna region (and not only) often characterized by buildings in ruins. The purpose of this study is mainly to document a legacy, which has already been studied in depth by historians, and previously lacked reliable architectural surveys for the definition of a credible as well as sustainable conservation project. Our contribution will focus on different techniques and methods used for the survey of these architectures, the characteristics of which- in the past- have made an effective survey of these buildings difficult, if not impossible. The survey of a ruin requires, much more than the evaluation of an intact building, reading skills and an interpretation of architectural spaces to better manage the stages of documentation and data processing. Through a series of case studies of fortified buildings in ruins, we intend to describe the reasons that guided the choice of the methods and tools used and to highlight the potentials and the limits of these choices in financial terms.

Uniform high order spectral methods for one and two dimensional Euler equations

NASA Technical Reports Server (NTRS)

Cai, Wei; Shu, Chi-Wang

1991-01-01

Uniform high order spectral methods to solve multi-dimensional Euler equations for gas dynamics are discussed. Uniform high order spectral approximations with spectral accuracy in smooth regions of solutions are constructed by introducing the idea of the Essentially Non-Oscillatory (ENO) polynomial interpolations into the spectral methods. The authors present numerical results for the inviscid Burgers' equation, and for the one dimensional Euler equations including the interactions between a shock wave and density disturbance, Sod's and Lax's shock tube problems, and the blast wave problem. The interaction between a Mach 3 two dimensional shock wave and a rotating vortex is simulated.

Current status of one- and two-dimensional numerical models: Successes and limitations

NASA Technical Reports Server (NTRS)

Schwartz, R. J.; Gray, J. L.; Lundstrom, M. S.

1985-01-01

The capabilities of one and two-dimensional numerical solar cell modeling programs (SCAP1D and SCAP2D) are described. The occasions when a two-dimensional model is required are discussed. The application of the models to design, analysis, and prediction are presented along with a discussion of problem areas for solar cell modeling.

Pressurized-Flat-Interface Heat Exchanger

NASA Technical Reports Server (NTRS)

Voss, F. E.; Howell, H. R.; Winkler, R. V.

1990-01-01

High thermal conductance obtained without leakage between loops. Heat-exchanger interface enables efficient transfer of heat between two working fluids without allowing fluids to intermingle. Interface thin, flat, and easy to integrate into thermal system. Possible application in chemical or pharmaceutical manufacturing when even trace contamination of process stream with water or other coolant ruins product. Reduces costs when highly corrosive fluids must be cooled or heated.

Robust L1-norm two-dimensional linear discriminant analysis.

PubMed

Li, Chun-Na; Shao, Yuan-Hai; Deng, Nai-Yang

2015-05-01

In this paper, we propose an L1-norm two-dimensional linear discriminant analysis (L1-2DLDA) with robust performance. Different from the conventional two-dimensional linear discriminant analysis with L2-norm (L2-2DLDA), where the optimization problem is transferred to a generalized eigenvalue problem, the optimization problem in our L1-2DLDA is solved by a simple justifiable iterative technique, and its convergence is guaranteed. Compared with L2-2DLDA, our L1-2DLDA is more robust to outliers and noises since the L1-norm is used. This is supported by our preliminary experiments on toy example and face datasets, which show the improvement of our L1-2DLDA over L2-2DLDA. Copyright © 2015 Elsevier Ltd. All rights reserved.

Implementation of pattern generation algorithm in forming Gilmore and Gomory model for two dimensional cutting stock problem

NASA Astrophysics Data System (ADS)

Octarina, Sisca; Radiana, Mutia; Bangun, Putra B. J.

2018-01-01

Two dimensional cutting stock problem (CSP) is a problem in determining the cutting pattern from a set of stock with standard length and width to fulfill the demand of items. Cutting patterns were determined in order to minimize the usage of stock. This research implemented pattern generation algorithm to formulate Gilmore and Gomory model of two dimensional CSP. The constraints of Gilmore and Gomory model was performed to assure the strips which cut in the first stage will be used in the second stage. Branch and Cut method was used to obtain the optimal solution. Based on the results, it found many patterns combination, if the optimal cutting patterns which correspond to the first stage were combined with the second stage.

Applications of FEM and BEM in two-dimensional fracture mechanics problems

NASA Technical Reports Server (NTRS)

Min, J. B.; Steeve, B. E.; Swanson, G. R.

1992-01-01

A comparison of the finite element method (FEM) and boundary element method (BEM) for the solution of two-dimensional plane strain problems in fracture mechanics is presented in this paper. Stress intensity factors (SIF's) were calculated using both methods for elastic plates with either a single-edge crack or an inclined-edge crack. In particular, two currently available programs, ANSYS for finite element analysis and BEASY for boundary element analysis, were used.

Convergence acceleration of the Proteus computer code with multigrid methods

NASA Technical Reports Server (NTRS)

Demuren, A. O.; Ibraheem, S. O.

1995-01-01

This report presents the results of a study to implement convergence acceleration techniques based on the multigrid concept in the two-dimensional and three-dimensional versions of the Proteus computer code. The first section presents a review of the relevant literature on the implementation of the multigrid methods in computer codes for compressible flow analysis. The next two sections present detailed stability analysis of numerical schemes for solving the Euler and Navier-Stokes equations, based on conventional von Neumann analysis and the bi-grid analysis, respectively. The next section presents details of the computational method used in the Proteus computer code. Finally, the multigrid implementation and applications to several two-dimensional and three-dimensional test problems are presented. The results of the present study show that the multigrid method always leads to a reduction in the number of iterations (or time steps) required for convergence. However, there is an overhead associated with the use of multigrid acceleration. The overhead is higher in 2-D problems than in 3-D problems, thus overall multigrid savings in CPU time are in general better in the latter. Savings of about 40-50 percent are typical in 3-D problems, but they are about 20-30 percent in large 2-D problems. The present multigrid method is applicable to steady-state problems and is therefore ineffective in problems with inherently unstable solutions.

Scalable posterior approximations for large-scale Bayesian inverse problems via likelihood-informed parameter and state reduction

NASA Astrophysics Data System (ADS)

Cui, Tiangang; Marzouk, Youssef; Willcox, Karen

2016-06-01

Two major bottlenecks to the solution of large-scale Bayesian inverse problems are the scaling of posterior sampling algorithms to high-dimensional parameter spaces and the computational cost of forward model evaluations. Yet incomplete or noisy data, the state variation and parameter dependence of the forward model, and correlations in the prior collectively provide useful structure that can be exploited for dimension reduction in this setting-both in the parameter space of the inverse problem and in the state space of the forward model. To this end, we show how to jointly construct low-dimensional subspaces of the parameter space and the state space in order to accelerate the Bayesian solution of the inverse problem. As a byproduct of state dimension reduction, we also show how to identify low-dimensional subspaces of the data in problems with high-dimensional observations. These subspaces enable approximation of the posterior as a product of two factors: (i) a projection of the posterior onto a low-dimensional parameter subspace, wherein the original likelihood is replaced by an approximation involving a reduced model; and (ii) the marginal prior distribution on the high-dimensional complement of the parameter subspace. We present and compare several strategies for constructing these subspaces using only a limited number of forward and adjoint model simulations. The resulting posterior approximations can rapidly be characterized using standard sampling techniques, e.g., Markov chain Monte Carlo. Two numerical examples demonstrate the accuracy and efficiency of our approach: inversion of an integral equation in atmospheric remote sensing, where the data dimension is very high; and the inference of a heterogeneous transmissivity field in a groundwater system, which involves a partial differential equation forward model with high dimensional state and parameters.

Squeezing the Efimov effect

NASA Astrophysics Data System (ADS)

Sandoval, J. H.; Bellotti, F. F.; Yamashita, M. T.; Frederico, T.; Fedorov, D. V.; Jensen, A. S.; Zinner, N. T.

2018-03-01

The quantum mechanical three-body problem is a source of continuing interest due to its complexity and not least due to the presence of fascinating solvable cases. The prime example is the Efimov effect where infinitely many bound states of identical bosons can arise at the threshold where the two-body problem has zero binding energy. An important aspect of the Efimov effect is the effect of spatial dimensionality; it has been observed in three dimensional systems, yet it is believed to be impossible in two dimensions. Using modern experimental techniques, it is possible to engineer trap geometry and thus address the intricate nature of quantum few-body physics as function of dimensionality. Here we present a framework for studying the three-body problem as one (continuously) changes the dimensionality of the system all the way from three, through two, and down to a single dimension. This is done by considering the Efimov favorable case of a mass-imbalanced system and with an external confinement provided by a typical experimental case with a (deformed) harmonic trap.

Assessment of numerical techniques for unsteady flow calculations

NASA Technical Reports Server (NTRS)

Hsieh, Kwang-Chung

1989-01-01

The characteristics of unsteady flow motions have long been a serious concern in the study of various fluid dynamic and combustion problems. With the advancement of computer resources, numerical approaches to these problems appear to be feasible. The objective of this paper is to assess the accuracy of several numerical schemes for unsteady flow calculations. In the present study, Fourier error analysis is performed for various numerical schemes based on a two-dimensional wave equation. Four methods sieved from the error analysis are then adopted for further assessment. Model problems include unsteady quasi-one-dimensional inviscid flows, two-dimensional wave propagations, and unsteady two-dimensional inviscid flows. According to the comparison between numerical and exact solutions, although second-order upwind scheme captures the unsteady flow and wave motions quite well, it is relatively more dissipative than sixth-order central difference scheme. Among various numerical approaches tested in this paper, the best performed one is Runge-Kutta method for time integration and six-order central difference for spatial discretization.

A Two-Dimensional Linear Bicharacteristic FDTD Method

NASA Technical Reports Server (NTRS)

Beggs, John H.

2002-01-01

The linear bicharacteristic scheme (LBS) was originally developed to improve unsteady solutions in computational acoustics and aeroacoustics. The LBS has previously been extended to treat lossy materials for one-dimensional problems. It is a classical leapfrog algorithm, but is combined with upwind bias in the spatial derivatives. This approach preserves the time-reversibility of the leapfrog algorithm, which results in no dissipation, and it permits more flexibility by the ability to adopt a characteristic based method. The use of characteristic variables allows the LBS to include the Perfectly Matched Layer boundary condition with no added storage or complexity. The LBS offers a central storage approach with lower dispersion than the Yee algorithm, plus it generalizes much easier to nonuniform grids. It has previously been applied to two and three-dimensional free-space electromagnetic propagation and scattering problems. This paper extends the LBS to the two-dimensional case. Results are presented for point source radiation problems, and the FDTD algorithm is chosen as a convenient reference for comparison.

Plane Poiseuille flow of a rarefied gas in the presence of strong gravitation.

PubMed

Doi, Toshiyuki

2011-02-01

Plane Poiseuille flow of a rarefied gas, which flows horizontally in the presence of strong gravitation, is studied based on the Boltzmann equation. Applying the asymptotic analysis for a small variation in the flow direction [Y. Sone, Molecular Gas Dynamics (Birkhäuser, 2007)], the two-dimensional problem is reduced to a one-dimensional problem, as in the case of a Poiseuille flow in the absence of gravitation, and the solution is obtained in a semianalytical form. The reduced one-dimensional problem is solved numerically for a hard sphere molecular gas over a wide range of the gas-rarefaction degree and the gravitational strength. The presence of gravitation reduces the mass flow rate, and the effect of gravitation is significant for large Knudsen numbers. To verify the validity of the asymptotic solution, a two-dimensional problem of a flow through a long channel is directly solved numerically, and the validity of the asymptotic solution is confirmed. ©2011 American Physical Society

Geomatic Archaeological Reconstruction and a Hybrid Viewer for the Archaelogical Site of CÁPARRA (spain)

NASA Astrophysics Data System (ADS)

Tejeda-Sánchez, C.; Muñoz-Nieto, A.; Rodríguez-Gonzálvez, P.

2018-05-01

Visualization and analysis use to be the final steps in Geomatics. This paper shows the workflow followed to set up a hybrid 3D archaeological viewer. Data acquisition of the site survey was done by means of low-cost close-range photogrammetric methods. With the aim not only to satisfy the general public but also the technicians, a large group of Geomatic products has been obtained (2d plans, 3d models, orthophotos, CAD models coming from vectorization, virtual anastylosis, and cross sections). Finally, all these products have been integrated into a three-dimensional archaeological information system. The hybrid archaeological viewer designed allows a metric and quality approach to the scientific analysis of the ruins, improving, thanks to the implementation of a database, and its potential for queries, the benefits of an ordinary topographic survey.

Numerical two-dimensional calculations of the formation of the solar nebula

NASA Technical Reports Server (NTRS)

Bodenheimer, Peter H.

1991-01-01

Numerical two dimensional calculations of the formation of the solar nebula are presented. The following subject areas are covered: (1) observational constraints of the properties of the initial solar nebula; (2) the physical problem; (3) review if two dimensional calculations of the formation phase; (4) recent models with hydrodynamics and radiative transport; and (5) further evolution of the system.

The Goertler vortex instability mechanism in three-dimensional boundary layers

NASA Technical Reports Server (NTRS)

Hall, P.

1984-01-01

The two dimensional boundary layer on a concave wall is centrifugally unstable with respect to vortices aligned with the basic flow for sufficiently high values of the Goertler number. However, in most situations of practical interest the basic flow is three dimensional and previous theoretical investigations do not apply. The linear stability of the flow over an infinitely long swept wall of variable curvature is considered. If there is no pressure gradient in the boundary layer the instability problem can always be related to an equivalent two dimensional calculation. However, in general, this is not the case and even for small values of the crossflow velocity field dramatic differences between the two and three dimensional problems emerge. When the size of the crossflow is further increased, the vortices in the neutral location have their axes locally perpendicular to the vortex lines of the basic flow.

An equivalent domain integral method in the two-dimensional analysis of mixed mode crack problems

NASA Technical Reports Server (NTRS)

Raju, I. S.; Shivakumar, K. N.

1990-01-01

An equivalent domain integral (EDI) method for calculating J-integrals for two-dimensional cracked elastic bodies is presented. The details of the method and its implementation are presented for isoparametric elements. The EDI method gave accurate values of the J-integrals for two mode I and two mixed mode problems. Numerical studies showed that domains consisting of one layer of elements are sufficient to obtain accurate J-integral values. Two procedures for separating the individual modes from the domain integrals are presented.

Bird-nest puzzle: can the study of unisexual flowers such as cucumber solve the problem of plant sex determination?

PubMed

Bai, Shu-Nong; Xu, Zhi-Hong

2012-06-01

Unisexual flower development has long been used as a model system to understand the mechanism of plant sex determination. However, based on our investigation of the mechanisms regulating the development of unisexual cucumber flowers, we have realized that understanding how organ development is inhibited may not necessarily reveal how an organ is formed. We refer to this problem as a "bird-nest puzzle," meaning one cannot understand how a bird lays and hatches its eggs by understanding how its nest is ruined. To understand the biological significance of unisexual flowers, we reexamine the original meaning of sex and its application in plants. Additionally, we propose that the fundamental biological advantage for the selection and maintenance of unisexual flowers during evolution is to promote cross pollination.

Optimality of the barrier strategy in de Finetti's dividend problem for spectrally negative Lévy processes: An alternative approach

NASA Astrophysics Data System (ADS)

Yin, Chuancun; Wang, Chunwei

2009-11-01

The optimal dividend problem proposed in de Finetti [1] is to find the dividend-payment strategy that maximizes the expected discounted value of dividends which are paid to the shareholders until the company is ruined. Avram et al. [9] studied the case when the risk process is modelled by a general spectrally negative Lévy process and Loeffen [10] gave sufficient conditions under which the optimal strategy is of the barrier type. Recently Kyprianou et al. [11] strengthened the result of Loeffen [10] which established a larger class of Lévy processes for which the barrier strategy is optimal among all admissible ones. In this paper we use an analytical argument to re-investigate the optimality of barrier dividend strategies considered in the three recent papers.

Computer simulation of plasma and N-body problems

NASA Technical Reports Server (NTRS)

Harries, W. L.; Miller, J. B.

1975-01-01

The following FORTRAN language computer codes are presented: (1) efficient two- and three-dimensional central force potential solvers; (2) a three-dimensional simulator of an isolated galaxy which incorporates the potential solver; (3) a two-dimensional particle-in-cell simulator of the Jeans instability in an infinite self-gravitating compressible gas; and (4) a two-dimensional particle-in-cell simulator of a rotating self-gravitating compressible gaseous system of which rectangular coordinate and superior polar coordinate versions were written.

Free boundary problems in shock reflection/diffraction and related transonic flow problems

PubMed Central

Chen, Gui-Qiang; Feldman, Mikhail

2015-01-01

Shock waves are steep wavefronts that are fundamental in nature, especially in high-speed fluid flows. When a shock hits an obstacle, or a flying body meets a shock, shock reflection/diffraction phenomena occur. In this paper, we show how several long-standing shock reflection/diffraction problems can be formulated as free boundary problems, discuss some recent progress in developing mathematical ideas, approaches and techniques for solving these problems, and present some further open problems in this direction. In particular, these shock problems include von Neumann's problem for shock reflection–diffraction by two-dimensional wedges with concave corner, Lighthill's problem for shock diffraction by two-dimensional wedges with convex corner, and Prandtl-Meyer's problem for supersonic flow impinging onto solid wedges, which are also fundamental in the mathematical theory of multidimensional conservation laws. PMID:26261363

Thomas-Fermi model for a bulk self-gravitating stellar object in two dimensions

NASA Astrophysics Data System (ADS)

De, Sanchari; Chakrabarty, Somenath

2015-09-01

In this article we have solved a hypothetical problem related to the stability and gross properties of two-dimensional self-gravitating stellar objects using the Thomas-Fermi model. The formalism presented here is an extension of the standard three-dimensional problem discussed in the book on statistical physics, Part I by Landau and Lifshitz. Further, the formalism presented in this article may be considered a class problem for post-graduate-level students of physics or may be assigned as a part of their dissertation project.

Classification Objects, Ideal Observers & Generative Models

ERIC Educational Resources Information Center

Olman, Cheryl; Kersten, Daniel

2004-01-01

A successful vision system must solve the problem of deriving geometrical information about three-dimensional objects from two-dimensional photometric input. The human visual system solves this problem with remarkable efficiency, and one challenge in vision research is to understand how neural representations of objects are formed and what visual…

A discontinuous Galerkin method for two-dimensional PDE models of Asian options

NASA Astrophysics Data System (ADS)

Hozman, J.; Tichý, T.; Cvejnová, D.

2016-06-01

In our previous research we have focused on the problem of plain vanilla option valuation using discontinuous Galerkin method for numerical PDE solution. Here we extend a simple one-dimensional problem into two-dimensional one and design a scheme for valuation of Asian options, i.e. options with payoff depending on the average of prices collected over prespecified horizon. The algorithm is based on the approach combining the advantages of the finite element methods together with the piecewise polynomial generally discontinuous approximations. Finally, an illustrative example using DAX option market data is provided.

Solving time-dependent two-dimensional eddy current problems

NASA Technical Reports Server (NTRS)

Lee, Min Eig; Hariharan, S. I.; Ida, Nathan

1988-01-01

Results of transient eddy current calculations are reported. For simplicity, a two-dimensional transverse magnetic field which is incident on an infinitely long conductor is considered. The conductor is assumed to be a good but not perfect conductor. The resulting problem is an interface initial boundary value problem with the boundary of the conductor being the interface. A finite difference method is used to march the solution explicitly in time. The method is shown. Treatment of appropriate radiation conditions is given special consideration. Results are validated with approximate analytic solutions. Two stringent test cases of high and low frequency incident waves are considered to validate the results.

Intertwined Hamiltonians in two-dimensional curved spaces

NASA Astrophysics Data System (ADS)

Aghababaei Samani, Keivan; Zarei, Mina

2005-04-01

The problem of intertwined Hamiltonians in two-dimensional curved spaces is investigated. Explicit results are obtained for Euclidean plane, Minkowski plane, Poincaré half plane (AdS2), de Sitter plane (dS2), sphere, and torus. It is shown that the intertwining operator is related to the Killing vector fields and the isometry group of corresponding space. It is shown that the intertwined potentials are closely connected to the integral curves of the Killing vector fields. Two problems are considered as applications of the formalism presented in the paper. The first one is the problem of Hamiltonians with equispaced energy levels and the second one is the problem of Hamiltonians whose spectrum is like the spectrum of a free particle.

Aerodynamics of an airfoil with a jet issuing from its surface

NASA Technical Reports Server (NTRS)

Tavella, D. A.; Karamcheti, K.

1982-01-01

A simple, two dimensional, incompressible and inviscid model for the problem posed by a two dimensional wing with a jet issuing from its lower surface is considered and a parametric analysis is carried out to observe how the aerodynamic characteristics depend on the different parameters. The mathematical problem constitutes a boundary value problem where the position of part of the boundary is not known a priori. A nonlinear optimization approach was used to solve the problem, and the analysis reveals interesting characteristics that may help to better understand the physics involved in more complex situations in connection with high lift systems.

Transonic small disturbances equation applied to the solution of two-dimensional nonsteady flows

NASA Technical Reports Server (NTRS)

Couston, M.; Angelini, J. J.; Mulak, P.

1980-01-01

Transonic nonsteady flows are of large practical interest. Aeroelastic instability prediction, control figured vehicle techniques or rotary wings in forward flight are some examples justifying the effort undertaken to improve knowledge of these problems is described. The numerical solution of these problems under the potential flow hypothesis is described. The use of an alternating direction implicit scheme allows the efficient resolution of the two dimensional transonic small perturbations equation.

Aerodynamics of Engine-Airframe Interaction

NASA Technical Reports Server (NTRS)

Caughey, D. A.

1986-01-01

The report describes progress in research directed towards the efficient solution of the inviscid Euler and Reynolds-averaged Navier-Stokes equations for transonic flows through engine inlets, and past complete aircraft configurations, with emphasis on the flowfields in the vicinity of engine inlets. The research focusses upon the development of solution-adaptive grid procedures for these problems, and the development of multi-grid algorithms in conjunction with both, implicit and explicit time-stepping schemes for the solution of three-dimensional problems. The work includes further development of mesh systems suitable for inlet and wing-fuselage-inlet geometries using a variational approach. Work during this reporting period concentrated upon two-dimensional problems, and has been in two general areas: (1) the development of solution-adaptive procedures to cluster the grid cells in regions of high (truncation) error;and (2) the development of a multigrid scheme for solution of the two-dimensional Euler equations using a diagonalized alternating direction implicit (ADI) smoothing algorithm.

Identification of the heat transfer coefficient in the two-dimensional model of binary alloy solidification

NASA Astrophysics Data System (ADS)

Hetmaniok, Edyta; Hristov, Jordan; Słota, Damian; Zielonka, Adam

2017-05-01

The paper presents the procedure for solving the inverse problem for the binary alloy solidification in a two-dimensional space. This is a continuation of some previous works of the authors investigating a similar problem but in the one-dimensional domain. Goal of the problem consists in identification of the heat transfer coefficient on boundary of the region and in reconstruction of the temperature distribution inside the considered region in case when the temperature measurements in selected points of the alloy are known. Mathematical model of the problem is based on the heat conduction equation with the substitute thermal capacity and with the liquidus and solidus temperatures varying in dependance on the concentration of the alloy component. For describing this concentration the Scheil model is used. Investigated procedure involves also the parallelized Ant Colony Optimization algorithm applied for minimizing a functional expressing the error of approximate solution.

A critical assessment of flux and source term closures in shallow water models with porosity for urban flood simulations

NASA Astrophysics Data System (ADS)

Guinot, Vincent

2017-11-01

The validity of flux and source term formulae used in shallow water models with porosity for urban flood simulations is assessed by solving the two-dimensional shallow water equations over computational domains representing periodic building layouts. The models under assessment are the Single Porosity (SP), the Integral Porosity (IP) and the Dual Integral Porosity (DIP) models. 9 different geometries are considered. 18 two-dimensional initial value problems and 6 two-dimensional boundary value problems are defined. This results in a set of 96 fine grid simulations. Analysing the simulation results leads to the following conclusions: (i) the DIP flux and source term models outperform those of the SP and IP models when the Riemann problem is aligned with the main street directions, (ii) all models give erroneous flux closures when is the Riemann problem is not aligned with one of the main street directions or when the main street directions are not orthogonal, (iii) the solution of the Riemann problem is self-similar in space-time when the street directions are orthogonal and the Riemann problem is aligned with one of them, (iv) a momentum balance confirms the existence of the transient momentum dissipation model presented in the DIP model, (v) none of the source term models presented so far in the literature allows all flow configurations to be accounted for(vi) future laboratory experiments aiming at the validation of flux and source term closures should focus on the high-resolution, two-dimensional monitoring of both water depth and flow velocity fields.

Complete Sets of Radiating and Nonradiating Parts of a Source and Their Fields with Applications in Inverse Scattering Limited-Angle Problems

PubMed Central

Louis, A. K.

2006-01-01

Many algorithms applied in inverse scattering problems use source-field systems instead of the direct computation of the unknown scatterer. It is well known that the resulting source problem does not have a unique solution, since certain parts of the source totally vanish outside of the reconstruction area. This paper provides for the two-dimensional case special sets of functions, which include all radiating and all nonradiating parts of the source. These sets are used to solve an acoustic inverse problem in two steps. The problem under discussion consists of determining an inhomogeneous obstacle supported in a part of a disc, from data, known for a subset of a two-dimensional circle. In a first step, the radiating parts are computed by solving a linear problem. The second step is nonlinear and consists of determining the nonradiating parts. PMID:23165060

In Carthage ruins: the illness of Sir Winston Churchill at Carthage, December 1943.

PubMed

Vale, J A; Scadding, J W

2017-09-01

This paper reviews Churchill's illness in Carthage in December 1943. It was characterised by fever that lasted 6 days, left lower lobe pneumonia and two episodes of atrial fibrillation. He was managed in a private villa by Lord Moran, his personal physician, with the assistance of two nurses and the expert advice of colleagues. Sulphadiazine and digitalis leaf were prescribed and Churchill recovered. It is remarkable that, despite the severity of his illness, he continued to direct the affairs of State from his bed.

Multi-Dimensional, Non-Pyrolyzing Ablation Test Problems

NASA Technical Reports Server (NTRS)

Risch, Tim; Kostyk, Chris

2016-01-01

Non-pyrolyzingcarbonaceous materials represent a class of candidate material for hypersonic vehicle components providing both structural and thermal protection system capabilities. Two problems relevant to this technology are presented. The first considers the one-dimensional ablation of a carbon material subject to convective heating. The second considers two-dimensional conduction in a rectangular block subject to radiative heating. Surface thermochemistry for both problems includes finite-rate surface kinetics at low temperatures, diffusion limited ablation at intermediate temperatures, and vaporization at high temperatures. The first problem requires the solution of both the steady-state thermal profile with respect to the ablating surface and the transient thermal history for a one-dimensional ablating planar slab with temperature-dependent material properties. The slab front face is convectively heated and also reradiates to a room temperature environment. The back face is adiabatic. The steady-state temperature profile and steady-state mass loss rate should be predicted. Time-dependent front and back face temperature, surface recession and recession rate along with the final temperature profile should be predicted for the time-dependent solution. The second problem requires the solution for the transient temperature history for an ablating, two-dimensional rectangular solid with anisotropic, temperature-dependent thermal properties. The front face is radiatively heated, convectively cooled, and also reradiates to a room temperature environment. The back face and sidewalls are adiabatic. The solution should include the following 9 items: final surface recession profile, time-dependent temperature history of both the front face and back face at both the centerline and sidewall, as well as the time-dependent surface recession and recession rate on the front face at both the centerline and sidewall. The results of the problems from all submitters will be collected, summarized, and presented at a later conference.

Gun Violence, Mental Illness, And Laws That Prohibit Gun Possession: Evidence From Two Florida Counties

PubMed Central

Swanson, Jeffrey W.; Easter, Michele M.; Robertson, Allison G.; Swartz, Marvin S.; Alanis-Hirsch, Kelly; Moseley, Daniel; Dion, Charles; Petrila, John

2016-01-01

Gun violence kills about ninety people every day in the United States, a toll measured in wasted and ruined lives and with an annual economic price tag exceeding $200 billion. Some policy makers suggest that reforming mental health care systems and improving point-of-purchase background checks to keep guns from mentally disturbed people will address the problem. Epidemiological research shows that serious mental illness contributes little to the risk of interpersonal violence but is a strong factor in suicide, which accounts for most firearm fatalities. Meanwhile, the effectiveness of gun restrictions focused on mental illness remains poorly understood. This article examines gun-related suicide and violent crime in people with serious mental illnesses, and whether legal restrictions on firearm sales to people with a history of mental health adjudication are effective in preventing gun violence. Among the study population in two large Florida counties, we found that 62 percent of violent gun crime arrests and 28 percent of gun suicides involved individuals not legally permitted to have a gun at the time. Suggested policy reforms include enacting risk-based gun removal laws and prohibiting guns from people involuntarily detained in short-term psychiatric hospitalizations. PMID:27269024

The Particle inside a Ring: A Two-Dimensional Quantum Problem Visualized by Scanning Tunneling Microscopy

ERIC Educational Resources Information Center

Ellison, Mark D.

2008-01-01

The one-dimensional particle-in-a-box model used to introduce quantum mechanics to students suffers from a tenuous connection to a real physical system. This article presents a two-dimensional model, the particle confined within a ring, that directly corresponds to observations of surface electrons in a metal trapped inside a circular barrier.…

Approximation and Numerical Analysis of Nonlinear Equations of Evolution.

DTIC Science & Technology

1980-01-31

dominant convective terms, or Stefan type problems such as the flow of fluids through porous media or the melting and freezing of ice. Such problems...means of formulating time-dependent Stefan problems was initiated. Classes of problems considered here include the one-phase and two-phase Stefan ...some new numerical methods were 2 developed for two dimensional, two-phase Stefan problems with time dependent boundary conditions. A variety of example

A solution for two-dimensional Fredholm integral equations of the second kind with periodic, semiperiodic, or nonperiodic kernels. [integral representation of the stationary Navier-Stokes problem

NASA Technical Reports Server (NTRS)

Gabrielsen, R. E.; Uenal, A.

1981-01-01

A numerical scheme for solving two dimensional Fredholm integral equations of the second kind is developed. The proof of the convergence of the numerical scheme is shown for three cases: the case of periodic kernels, the case of semiperiodic kernels, and the case of nonperiodic kernels. Applications to the incompressible, stationary Navier-Stokes problem are of primary interest.

Developing cross entropy genetic algorithm for solving Two-Dimensional Loading Heterogeneous Fleet Vehicle Routing Problem (2L-HFVRP)

NASA Astrophysics Data System (ADS)

Paramestha, D. L.; Santosa, B.

2018-04-01

Two-dimensional Loading Heterogeneous Fleet Vehicle Routing Problem (2L-HFVRP) is a combination of Heterogeneous Fleet VRP and a packing problem well-known as Two-Dimensional Bin Packing Problem (BPP). 2L-HFVRP is a Heterogeneous Fleet VRP in which these costumer demands are formed by a set of two-dimensional rectangular weighted item. These demands must be served by a heterogeneous fleet of vehicles with a fix and variable cost from the depot. The objective function 2L-HFVRP is to minimize the total transportation cost. All formed routes must be consistent with the capacity and loading process of the vehicle. Sequential and unrestricted scenarios are considered in this paper. We propose a metaheuristic which is a combination of the Genetic Algorithm (GA) and the Cross Entropy (CE) named Cross Entropy Genetic Algorithm (CEGA) to solve the 2L-HFVRP. The mutation concept on GA is used to speed up the algorithm CE to find the optimal solution. The mutation mechanism was based on local improvement (2-opt, 1-1 Exchange, and 1-0 Exchange). The probability transition matrix mechanism on CE is used to avoid getting stuck in the local optimum. The effectiveness of CEGA was tested on benchmark instance based 2L-HFVRP. The result of experiments shows a competitive result compared with the other algorithm.

Learning control system design based on 2-D theory - An application to parallel link manipulator

NASA Technical Reports Server (NTRS)

Geng, Z.; Carroll, R. L.; Lee, J. D.; Haynes, L. H.

1990-01-01

An approach to iterative learning control system design based on two-dimensional system theory is presented. A two-dimensional model for the iterative learning control system which reveals the connections between learning control systems and two-dimensional system theory is established. A learning control algorithm is proposed, and the convergence of learning using this algorithm is guaranteed by two-dimensional stability. The learning algorithm is applied successfully to the trajectory tracking control problem for a parallel link robot manipulator. The excellent performance of this learning algorithm is demonstrated by the computer simulation results.

Dynamics of the cosmological relaxation after reheating

NASA Astrophysics Data System (ADS)

Choi, Kiwoon; Kim, Hyungjin; Sekiguchi, Toyokazu

2017-04-01

We examine if the cosmological relaxation mechanism, which was proposed recently as a new solution to the hierarchy problem, can be compatible with high reheating temperature well above the weak scale. As the barrier potential disappears at high temperature, the relaxion rolls down further after the reheating, which may ruin the successful implementation of the relaxation mechanism. It is noted that if the relaxion is coupled to a dark gauge boson, the new frictional force arising from dark gauge boson production can efficiently slow down the relaxion motion, which allows the relaxion to be stabilized after the electroweak phase transition for a wide range of model parameters, while satisfying the known observational constraints.

An equivalent domain integral for analysis of two-dimensional mixed mode problems

NASA Technical Reports Server (NTRS)

Raju, I. S.; Shivakumar, K. N.

1989-01-01

An equivalent domain integral (EDI) method for calculating J-integrals for two-dimensional cracked elastic bodies subjected to mixed mode loading is presented. The total and product integrals consist of the sum of an area or domain integral and line integrals on the crack faces. The EDI method gave accurate values of the J-integrals for two mode I and two mixed mode problems. Numerical studies showed that domains consisting of one layer of elements are sufficient to obtain accurate J-integral values. Two procedures for separating the individual modes from the domain integrals are presented. The procedure that uses the symmetric and antisymmetric components of the stress and displacement fields to calculate the individual modes gave accurate values of the integrals for all the problems analyzed.

Design of a rotational three-dimensional nonimaging device by a compensated two-dimensional design process.

PubMed

Yang, Yi; Qian, Ke-Yuan; Luo, Yi

2006-07-20

A compensation process has been developed to design rotational three-dimensional (3D) nonimaging devices. By compensating the desired light distribution during a two-dimensional (2D) design process for an extended Lambertian source using a compensation coefficient, the meridian plane of a 3D device with good performance can be obtained. This method is suitable in many cases with fast calculation speed. Solutions to two kinds of optical design problems have been proposed, and the limitation of this compensated 2D design method is discussed.

Two solvable problems of planar geometrical optics.

PubMed

Borghero, Francesco; Bozis, George

2006-12-01

In the framework of geometrical optics we consider a two-dimensional transparent inhomogeneous isotropic medium (dispersive or not). We show that (i) for any family belonging to a certain class of planar monoparametric families of monochromatic light rays given in the form f(x,y)=c of any definite color and satisfying a differential condition, all the refractive index profiles n=n(x,y) allowing for the creation of the given family can be found analytically (inverse problem) and that (ii) for any member of a class of two-dimensional refractive index profiles n=n(x,y) satisfying a differential condition, all the compatible families of light rays can be found analytically (direct problem). We present appropriate examples.

Analysis of a Two-Dimensional Thermal Cloaking Problem on the Basis of Optimization

NASA Astrophysics Data System (ADS)

Alekseev, G. V.

2018-04-01

For a two-dimensional model of thermal scattering, inverse problems arising in the development of tools for cloaking material bodies on the basis of a mixed thermal cloaking strategy are considered. By applying the optimization approach, these problems are reduced to optimization ones in which the role of controls is played by variable parameters of the medium occupying the cloaking shell and by the heat flux through a boundary segment of the basic domain. The solvability of the direct and optimization problems is proved, and an optimality system is derived. Based on its analysis, sufficient conditions on the input data are established that ensure the uniqueness and stability of optimal solutions.

High-frequency modes in a two-dimensional rectangular room with windows

NASA Astrophysics Data System (ADS)

Shabalina, E. D.; Shirgina, N. V.; Shanin, A. V.

2010-07-01

We examine a two-dimensional model problem of architectural acoustics on sound propagation in a rectangular room with windows. It is supposed that the walls are ideally flat and hard; the windows absorb all energy that falls upon them. We search for the modes of such a room having minimal attenuation indices, which have the expressed structure of billiard trajectories. The main attenuation mechanism for such modes is diffraction at the edges of the windows. We construct estimates for the attenuation indices of the given modes based on the solution to the Weinstein problem. We formulate diffraction problems similar to the statement of the Weinstein problem that describe the attenuation of billiard modes in complex situations.

On Born's Conjecture about Optimal Distribution of Charges for an Infinite Ionic Crystal

NASA Astrophysics Data System (ADS)

Bétermin, Laurent; Knüpfer, Hans

2018-04-01

We study the problem for the optimal charge distribution on the sites of a fixed Bravais lattice. In particular, we prove Born's conjecture about the optimality of the rock salt alternate distribution of charges on a cubic lattice (and more generally on a d-dimensional orthorhombic lattice). Furthermore, we study this problem on the two-dimensional triangular lattice and we prove the optimality of a two-component honeycomb distribution of charges. The results hold for a class of completely monotone interaction potentials which includes Coulomb-type interactions for d≥3 . In a more general setting, we derive a connection between the optimal charge problem and a minimization problem for the translated lattice theta function.

Robust Multigrid Smoothers for Three Dimensional Elliptic Equations with Strong Anisotropies

NASA Technical Reports Server (NTRS)

Llorente, Ignacio M.; Melson, N. Duane

1998-01-01

We discuss the behavior of several plane relaxation methods as multigrid smoothers for the solution of a discrete anisotropic elliptic model problem on cell-centered grids. The methods compared are plane Jacobi with damping, plane Jacobi with partial damping, plane Gauss-Seidel, plane zebra Gauss-Seidel, and line Gauss-Seidel. Based on numerical experiments and local mode analysis, we compare the smoothing factor of the different methods in the presence of strong anisotropies. A four-color Gauss-Seidel method is found to have the best numerical and architectural properties of the methods considered in the present work. Although alternating direction plane relaxation schemes are simpler and more robust than other approaches, they are not currently used in industrial and production codes because they require the solution of a two-dimensional problem for each plane in each direction. We verify the theoretical predictions of Thole and Trottenberg that an exact solution of each plane is not necessary and that a single two-dimensional multigrid cycle gives the same result as an exact solution, in much less execution time. Parallelization of the two-dimensional multigrid cycles, the kernel of the three-dimensional implicit solver, is also discussed. Alternating-plane smoothers are found to be highly efficient multigrid smoothers for anisotropic elliptic problems.

A Two-Dimensional Linear Bicharacteristic Scheme for Electromagnetics

NASA Technical Reports Server (NTRS)

Beggs, John H.

2002-01-01

The upwind leapfrog or Linear Bicharacteristic Scheme (LBS) has previously been implemented and demonstrated on one-dimensional electromagnetic wave propagation problems. This memorandum extends the Linear Bicharacteristic Scheme for computational electromagnetics to model lossy dielectric and magnetic materials and perfect electrical conductors in two dimensions. This is accomplished by proper implementation of the LBS for homogeneous lossy dielectric and magnetic media and for perfect electrical conductors. Both the Transverse Electric and Transverse Magnetic polarizations are considered. Computational requirements and a Fourier analysis are also discussed. Heterogeneous media are modeled through implementation of surface boundary conditions and no special extrapolations or interpolations at dielectric material boundaries are required. Results are presented for two-dimensional model problems on uniform grids, and the Finite Difference Time Domain (FDTD) algorithm is chosen as a convenient reference algorithm for comparison. The results demonstrate that the two-dimensional explicit LBS is a dissipation-free, second-order accurate algorithm which uses a smaller stencil than the FDTD algorithm, yet it has less phase velocity error.

Phase-space finite elements in a least-squares solution of the transport equation

DOE Office of Scientific and Technical Information (OSTI.GOV)

Drumm, C.; Fan, W.; Pautz, S.

2013-07-01

The linear Boltzmann transport equation is solved using a least-squares finite element approximation in the space, angular and energy phase-space variables. The method is applied to both neutral particle transport and also to charged particle transport in the presence of an electric field, where the angular and energy derivative terms are handled with the energy/angular finite elements approximation, in a manner analogous to the way the spatial streaming term is handled. For multi-dimensional problems, a novel approach is used for the angular finite elements: mapping the surface of a unit sphere to a two-dimensional planar region and using a meshingmore » tool to generate a mesh. In this manner, much of the spatial finite-elements machinery can be easily adapted to handle the angular variable. The energy variable and the angular variable for one-dimensional problems make use of edge/beam elements, also building upon the spatial finite elements capabilities. The methods described here can make use of either continuous or discontinuous finite elements in space, angle and/or energy, with the use of continuous finite elements resulting in a smaller problem size and the use of discontinuous finite elements resulting in more accurate solutions for certain types of problems. The work described in this paper makes use of continuous finite elements, so that the resulting linear system is symmetric positive definite and can be solved with a highly efficient parallel preconditioned conjugate gradients algorithm. The phase-space finite elements capability has been built into the Sceptre code and applied to several test problems, including a simple one-dimensional problem with an analytic solution available, a two-dimensional problem with an isolated source term, showing how the method essentially eliminates ray effects encountered with discrete ordinates, and a simple one-dimensional charged-particle transport problem in the presence of an electric field. (authors)« less

A variational principle for compressible fluid mechanics: Discussion of the multi-dimensional theory

NASA Technical Reports Server (NTRS)

Prozan, R. J.

1982-01-01

The variational principle for compressible fluid mechanics previously introduced is extended to two dimensional flow. The analysis is stable, exactly conservative, adaptable to coarse or fine grids, and very fast. Solutions for two dimensional problems are included. The excellent behavior and results lend further credence to the variational concept and its applicability to the numerical analysis of complex flow fields.

Building on ruins: Copernicus' defense of ancient astronomers against modern critics.

PubMed

Galle, Karl

2008-09-01

Nicholas Copernicus'De revolutionibus (1543) is a frequent starting point for histories of the Scientific Revolution on account of his dramatic reversal of the cosmic order handed down from antiquity. Nevertheless, Copernicus also mounted a surprisingly sharp attack on one of his contemporaries who tried to correct ancient astronomers. This uneasy balance between respecting and criticizing ancient works was part of broader contemporaneous attempts to grapple with the fragmentary legacies of past generations. Debates over stone ruins as well as manuscript texts shaped the evolution of early printed books, and artists, printers and instrument makers joined natural philosophers in pursuing novelties even while emulating tradition.

A solution procedure for behavior of thick plates on a nonlinear foundation and postbuckling behavior of long plates

NASA Technical Reports Server (NTRS)

Stein, M.; Stein, P. A.

1978-01-01

Approximate solutions for three nonlinear orthotropic plate problems are presented: (1) a thick plate attached to a pad having nonlinear material properties which, in turn, is attached to a substructure which is then deformed; (2) a long plate loaded in inplane longitudinal compression beyond its buckling load; and (3) a long plate loaded in inplane shear beyond its buckling load. For all three problems, the two dimensional plate equations are reduced to one dimensional equations in the y-direction by using a one dimensional trigonometric approximation in the x-direction. Each problem uses different trigonometric terms. Solutions are obtained using an existing algorithm for simultaneous, first order, nonlinear, ordinary differential equations subject to two point boundary conditions. Ordinary differential equations are derived to determine the variable coefficients of the trigonometric terms.

Inference of Vohradský's Models of Genetic Networks by Solving Two-Dimensional Function Optimization Problems

PubMed Central

Kimura, Shuhei; Sato, Masanao; Okada-Hatakeyama, Mariko

2013-01-01

The inference of a genetic network is a problem in which mutual interactions among genes are inferred from time-series of gene expression levels. While a number of models have been proposed to describe genetic networks, this study focuses on a mathematical model proposed by Vohradský. Because of its advantageous features, several researchers have proposed the inference methods based on Vohradský's model. When trying to analyze large-scale networks consisting of dozens of genes, however, these methods must solve high-dimensional non-linear function optimization problems. In order to resolve the difficulty of estimating the parameters of the Vohradský's model, this study proposes a new method that defines the problem as several two-dimensional function optimization problems. Through numerical experiments on artificial genetic network inference problems, we showed that, although the computation time of the proposed method is not the shortest, the method has the ability to estimate parameters of Vohradský's models more effectively with sufficiently short computation times. This study then applied the proposed method to an actual inference problem of the bacterial SOS DNA repair system, and succeeded in finding several reasonable regulations. PMID:24386175

Two-dimensional problem of two Coulomb centers at small intercenter distances

NASA Astrophysics Data System (ADS)

Bondar, D. I.; Hnatich, M.; Lazur, V. Yu.

2006-08-01

We use analytic methods to analyze the discrete spectrum for the problem (Z1eZ2)2 in the united-atom limit ( R ≪ 1) and obtain asymptotic expansions for the quantum defect and energy terms of the system (Z1eZ2)2 at small intercenter distances R up to terms of the order O(R6). We investigate the effect of the dimensionality factor on the energy spectrum of the hydrogen molecular ion H{2/+}.

On Riemann boundary value problems for null solutions of the two dimensional Helmholtz equation

NASA Astrophysics Data System (ADS)

Bory Reyes, Juan; Abreu Blaya, Ricardo; Rodríguez Dagnino, Ramón Martin; Kats, Boris Aleksandrovich

2018-01-01

The Riemann boundary value problem (RBVP to shorten notation) in the complex plane, for different classes of functions and curves, is still widely used in mathematical physics and engineering. For instance, in elasticity theory, hydro and aerodynamics, shell theory, quantum mechanics, theory of orthogonal polynomials, and so on. In this paper, we present an appropriate hyperholomorphic approach to the RBVP associated to the two dimensional Helmholtz equation in R^2 . Our analysis is based on a suitable operator calculus.

Extrapolation techniques applied to matrix methods in neutron diffusion problems

NASA Technical Reports Server (NTRS)

Mccready, Robert R

1956-01-01

A general matrix method is developed for the solution of characteristic-value problems of the type arising in many physical applications. The scheme employed is essentially that of Gauss and Seidel with appropriate modifications needed to make it applicable to characteristic-value problems. An iterative procedure produces a sequence of estimates to the answer; and extrapolation techniques, based upon previous behavior of iterants, are utilized in speeding convergence. Theoretically sound limits are placed on the magnitude of the extrapolation that may be tolerated. This matrix method is applied to the problem of finding criticality and neutron fluxes in a nuclear reactor with control rods. The two-dimensional finite-difference approximation to the two-group neutron fluxes in a nuclear reactor with control rods. The two-dimensional finite-difference approximation to the two-group neutron-diffusion equations is treated. Results for this example are indicated.

User's guide for NASCRIN: A vectorized code for calculating two-dimensional supersonic internal flow fields

NASA Technical Reports Server (NTRS)

Kumar, A.

1984-01-01

A computer program NASCRIN has been developed for analyzing two-dimensional flow fields in high-speed inlets. It solves the two-dimensional Euler or Navier-Stokes equations in conservation form by an explicit, two-step finite-difference method. An explicit-implicit method can also be used at the user's discretion for viscous flow calculations. For turbulent flow, an algebraic, two-layer eddy-viscosity model is used. The code is operational on the CDC CYBER 203 computer system and is highly vectorized to take full advantage of the vector-processing capability of the system. It is highly user oriented and is structured in such a way that for most supersonic flow problems, the user has to make only a few changes. Although the code is primarily written for supersonic internal flow, it can be used with suitable changes in the boundary conditions for a variety of other problems.

Comparison of an algebraic multigrid algorithm to two iterative solvers used for modeling ground water flow and transport

USGS Publications Warehouse

Detwiler, R.L.; Mehl, S.; Rajaram, H.; Cheung, W.W.

2002-01-01

Numerical solution of large-scale ground water flow and transport problems is often constrained by the convergence behavior of the iterative solvers used to solve the resulting systems of equations. We demonstrate the ability of an algebraic multigrid algorithm (AMG) to efficiently solve the large, sparse systems of equations that result from computational models of ground water flow and transport in large and complex domains. Unlike geometric multigrid methods, this algorithm is applicable to problems in complex flow geometries, such as those encountered in pore-scale modeling of two-phase flow and transport. We integrated AMG into MODFLOW 2000 to compare two- and three-dimensional flow simulations using AMG to simulations using PCG2, a preconditioned conjugate gradient solver that uses the modified incomplete Cholesky preconditioner and is included with MODFLOW 2000. CPU times required for convergence with AMG were up to 140 times faster than those for PCG2. The cost of this increased speed was up to a nine-fold increase in required random access memory (RAM) for the three-dimensional problems and up to a four-fold increase in required RAM for the two-dimensional problems. We also compared two-dimensional numerical simulations of steady-state transport using AMG and the generalized minimum residual method with an incomplete LU-decomposition preconditioner. For these transport simulations, AMG yielded increased speeds of up to 17 times with only a 20% increase in required RAM. The ability of AMG to solve flow and transport problems in large, complex flow systems and its ready availability make it an ideal solver for use in both field-scale and pore-scale modeling.

Knowledge and Valorization of Historical Sites Through 3d Documentation and Modeling

NASA Astrophysics Data System (ADS)

Farella, E.; Menna, F.; Nocerino, E.; Morabito, D.; Remondino, F.; Campi, M.

2016-06-01

The paper presents the first results of an interdisciplinary project related to the 3D documentation, dissemination, valorization and digital access of archeological sites. Beside the mere 3D documentation aim, the project has two goals: (i) to easily explore and share via web references and results of the interdisciplinary work, including the interpretative process and the final reconstruction of the remains; (ii) to promote and valorize archaeological areas using reality-based 3D data and Virtual Reality devices. This method has been verified on the ruins of the archeological site of Pausilypon, a maritime villa of Roman period (Naples, Italy). Using Unity3D, the virtual tour of the heritage site was integrated and enriched with the surveyed 3D data, text documents, CAAD reconstruction hypotheses, drawings, photos, etc. In this way, starting from the actual appearance of the ruins (panoramic images), passing through the 3D digital surveying models and several other historical information, the user is able to access virtual contents and reconstructed scenarios, all in a single virtual, interactive and immersive environment. These contents and scenarios allow to derive documentation and geometrical information, understand the site, perform analyses, see interpretative processes, communicate historical information and valorize the heritage location.

Ground Penetrating Radar Survey at Yoros Fortesss,Istanbul

NASA Astrophysics Data System (ADS)

Kucukdemirci, M.; Yalçın, A. B.

2016-12-01

Geophysical methods are effective tool to detect the archaeological remains and materials, which were hidden under the ground. One of the most frequently used methods for archaeological prospection is Ground Penetrating Radar (GPR). This paper illustrates the small scale GPR survey to determine the buried archaeological features around the Yoros Fortress, located on shores of the Bosporus strait in Istanbul, during the archaeological excavations. The survey was carried out with a GSSI SIR 3000 system, using 400 Mhz center frequency bistatic antenna with the configuration of 16 bits dynamic range and 512 samples per scan. The data were collected along parallel profiles with an interval of 0.50 meters with zigzag profile configuration on the survey grids. The GPR data were processed by GPR-Slice V.7 (Ground Penetrating Radar Imaging Software). As a result, in the first shallow depths, some scattered anomalies were detected. These can be related to a small portion of archaeological ruins close to the surface. In the deeper levels, the geometry of the anomalies related to the possible archaeological ruins, looks clearer. Two horizontal and parallel anomalies were detected, with the direction NS in the depth of 1.45 meters, possibly related to the ancient channels.

Three-dimensional finite element analysis for high velocity impact. [of projectiles from space debris

NASA Technical Reports Server (NTRS)

Chan, S. T. K.; Lee, C. H.; Brashears, M. R.

1975-01-01

A finite element algorithm for solving unsteady, three-dimensional high velocity impact problems is presented. A computer program was developed based on the Eulerian hydroelasto-viscoplastic formulation and the utilization of the theorem of weak solutions. The equations solved consist of conservation of mass, momentum, and energy, equation of state, and appropriate constitutive equations. The solution technique is a time-dependent finite element analysis utilizing three-dimensional isoparametric elements, in conjunction with a generalized two-step time integration scheme. The developed code was demonstrated by solving one-dimensional as well as three-dimensional impact problems for both the inviscid hydrodynamic model and the hydroelasto-viscoplastic model.

Arbitrarily high-order time-stepping schemes based on the operator spectrum theory for high-dimensional nonlinear Klein-Gordon equations

NASA Astrophysics Data System (ADS)

Liu, Changying; Wu, Xinyuan

2017-07-01

In this paper we explore arbitrarily high-order Lagrange collocation-type time-stepping schemes for effectively solving high-dimensional nonlinear Klein-Gordon equations with different boundary conditions. We begin with one-dimensional periodic boundary problems and first formulate an abstract ordinary differential equation (ODE) on a suitable infinity-dimensional function space based on the operator spectrum theory. We then introduce an operator-variation-of-constants formula which is essential for the derivation of our arbitrarily high-order Lagrange collocation-type time-stepping schemes for the nonlinear abstract ODE. The nonlinear stability and convergence are rigorously analysed once the spatial differential operator is approximated by an appropriate positive semi-definite matrix under some suitable smoothness assumptions. With regard to the two dimensional Dirichlet or Neumann boundary problems, our new time-stepping schemes coupled with discrete Fast Sine / Cosine Transformation can be applied to simulate the two-dimensional nonlinear Klein-Gordon equations effectively. All essential features of the methodology are present in one-dimensional and two-dimensional cases, although the schemes to be analysed lend themselves with equal to higher-dimensional case. The numerical simulation is implemented and the numerical results clearly demonstrate the advantage and effectiveness of our new schemes in comparison with the existing numerical methods for solving nonlinear Klein-Gordon equations in the literature.

Two-dimensional radiative transfer. I - Planar geometry. [in stellar atmospheres

NASA Technical Reports Server (NTRS)

Mihalas, D.; Auer, L. H.; Mihalas, B. R.

1978-01-01

Differential-equation methods for solving the transfer equation in two-dimensional planar geometries are developed. One method, which uses a Hermitian integration formula on ray segments through grid points, proves to be extremely well suited to velocity-dependent problems. An efficient elimination scheme is developed for which the computing time scales linearly with the number of angles and frequencies; problems with large velocity amplitudes can thus be treated accurately. A very accurate and efficient method for performing a formal solution is also presented. A discussion is given of several examples of periodic media and free-standing slabs, both in static cases and with velocity fields. For the free-standing slabs, two-dimensional transport effects are significant near boundaries, but no important effects were found in any of the periodic cases studied.

Two-dimensional integrating matrices on rectangular grids. [solving differential equations associated with rotating structures

NASA Technical Reports Server (NTRS)

Lakin, W. D.

1981-01-01

The use of integrating matrices in solving differential equations associated with rotating beam configurations is examined. In vibration problems, by expressing the equations of motion of the beam in matrix notation, utilizing the integrating matrix as an operator, and applying the boundary conditions, the spatial dependence is removed from the governing partial differential equations and the resulting ordinary differential equations can be cast into standard eigenvalue form. Integrating matrices are derived based on two dimensional rectangular grids with arbitrary grid spacings allowed in one direction. The derivation of higher dimensional integrating matrices is the initial step in the generalization of the integrating matrix methodology to vibration and stability problems involving plates and shells.

A Multifunctional Interface Method for Coupling Finite Element and Finite Difference Methods: Two-Dimensional Scalar-Field Problems

NASA Technical Reports Server (NTRS)

Ransom, Jonathan B.

2002-01-01

A multifunctional interface method with capabilities for variable-fidelity modeling and multiple method analysis is presented. The methodology provides an effective capability by which domains with diverse idealizations can be modeled independently to exploit the advantages of one approach over another. The multifunctional method is used to couple independently discretized subdomains, and it is used to couple the finite element and the finite difference methods. The method is based on a weighted residual variational method and is presented for two-dimensional scalar-field problems. A verification test problem and a benchmark application are presented, and the computational implications are discussed.

electromagnetics, eddy current, computer codes

DOE Office of Scientific and Technical Information (OSTI.GOV)

Gartling, David

TORO Version 4 is designed for finite element analysis of steady, transient and time-harmonic, multi-dimensional, quasi-static problems in electromagnetics. The code allows simulation of electrostatic fields, steady current flows, magnetostatics and eddy current problems in plane or axisymmetric, two-dimensional geometries. TORO is easily coupled to heat conduction and solid mechanics codes to allow multi-physics simulations to be performed.

2D and 3D Traveling Salesman Problem

ERIC Educational Resources Information Center

Haxhimusa, Yll; Carpenter, Edward; Catrambone, Joseph; Foldes, David; Stefanov, Emil; Arns, Laura; Pizlo, Zygmunt

2011-01-01

When a two-dimensional (2D) traveling salesman problem (TSP) is presented on a computer screen, human subjects can produce near-optimal tours in linear time. In this study we tested human performance on a real and virtual floor, as well as in a three-dimensional (3D) virtual space. Human performance on the real floor is as good as that on a…

COMOC 2: Two-dimensional aerodynamics sequence, computer program user's guide

NASA Technical Reports Server (NTRS)

Manhardt, P. D.; Orzechowski, J. A.; Baker, A. J.

1977-01-01

The COMOC finite element fluid mechanics computer program system is applicable to diverse problem classes. The two dimensional aerodynamics sequence was established for solution of the potential and/or viscous and turbulent flowfields associated with subsonic flight of elementary two dimensional isolated airfoils. The sequence is constituted of three specific flowfield options in COMOC for two dimensional flows. These include the potential flow option, the boundary layer option, and the parabolic Navier-Stokes option. By sequencing through these options, it is possible to computationally construct a weak-interaction model of the aerodynamic flowfield. This report is the user's guide to operation of COMOC for the aerodynamics sequence.

Three-Dimensional Inverse Transport Solver Based on Compressive Sensing Technique

NASA Astrophysics Data System (ADS)

Cheng, Yuxiong; Wu, Hongchun; Cao, Liangzhi; Zheng, Youqi

2013-09-01

According to the direct exposure measurements from flash radiographic image, a compressive sensing-based method for three-dimensional inverse transport problem is presented. The linear absorption coefficients and interface locations of objects are reconstructed directly at the same time. It is always very expensive to obtain enough measurements. With limited measurements, compressive sensing sparse reconstruction technique orthogonal matching pursuit is applied to obtain the sparse coefficients by solving an optimization problem. A three-dimensional inverse transport solver is developed based on a compressive sensing-based technique. There are three features in this solver: (1) AutoCAD is employed as a geometry preprocessor due to its powerful capacity in graphic. (2) The forward projection matrix rather than Gauss matrix is constructed by the visualization tool generator. (3) Fourier transform and Daubechies wavelet transform are adopted to convert an underdetermined system to a well-posed system in the algorithm. Simulations are performed and numerical results in pseudo-sine absorption problem, two-cube problem and two-cylinder problem when using compressive sensing-based solver agree well with the reference value.

Constrained optimization by radial basis function interpolation for high-dimensional expensive black-box problems with infeasible initial points

NASA Astrophysics Data System (ADS)

Regis, Rommel G.

2014-02-01

This article develops two new algorithms for constrained expensive black-box optimization that use radial basis function surrogates for the objective and constraint functions. These algorithms are called COBRA and Extended ConstrLMSRBF and, unlike previous surrogate-based approaches, they can be used for high-dimensional problems where all initial points are infeasible. They both follow a two-phase approach where the first phase finds a feasible point while the second phase improves this feasible point. COBRA and Extended ConstrLMSRBF are compared with alternative methods on 20 test problems and on the MOPTA08 benchmark automotive problem (D.R. Jones, Presented at MOPTA 2008), which has 124 decision variables and 68 black-box inequality constraints. The alternatives include a sequential penalty derivative-free algorithm, a direct search method with kriging surrogates, and two multistart methods. Numerical results show that COBRA algorithms are competitive with Extended ConstrLMSRBF and they generally outperform the alternatives on the MOPTA08 problem and most of the test problems.

Phase unwrapping in three dimensions with application to InSAR time series.

PubMed

Hooper, Andrew; Zebker, Howard A

2007-09-01

The problem of phase unwrapping in two dimensions has been studied extensively in the past two decades, but the three-dimensional (3D) problem has so far received relatively little attention. We develop here a theoretical framework for 3D phase unwrapping and also describe two algorithms for implementation, both of which can be applied to synthetic aperture radar interferometry (InSAR) time series. We test the algorithms on simulated data and find both give more accurate results than a two-dimensional algorithm. When applied to actual InSAR time series, we find good agreement both between the algorithms and with ground truth.

Second order symmetry-preserving conservative Lagrangian scheme for compressible Euler equations in two-dimensional cylindrical coordinates

DOE Office of Scientific and Technical Information (OSTI.GOV)

Cheng, Juan, E-mail: cheng_juan@iapcm.ac.cn; Shu, Chi-Wang, E-mail: shu@dam.brown.edu

In applications such as astrophysics and inertial confinement fusion, there are many three-dimensional cylindrical-symmetric multi-material problems which are usually simulated by Lagrangian schemes in the two-dimensional cylindrical coordinates. For this type of simulation, a critical issue for the schemes is to keep spherical symmetry in the cylindrical coordinate system if the original physical problem has this symmetry. In the past decades, several Lagrangian schemes with such symmetry property have been developed, but all of them are only first order accurate. In this paper, we develop a second order cell-centered Lagrangian scheme for solving compressible Euler equations in cylindrical coordinates, basedmore » on the control volume discretizations, which is designed to have uniformly second order accuracy and capability to preserve one-dimensional spherical symmetry in a two-dimensional cylindrical geometry when computed on an equal-angle-zoned initial grid. The scheme maintains several good properties such as conservation for mass, momentum and total energy, and the geometric conservation law. Several two-dimensional numerical examples in cylindrical coordinates are presented to demonstrate the good performance of the scheme in terms of accuracy, symmetry, non-oscillation and robustness. The advantage of higher order accuracy is demonstrated in these examples.« less

Two-dimensional unsteady lift problems in supersonic flight

NASA Technical Reports Server (NTRS)

Heaslet, Max A; Lomax, Harvard

1949-01-01

The variation of pressure distribution is calculated for a two-dimensional supersonic airfoil either experiencing a sudden angle-of-attack change or entering a sharp-edge gust. From these pressure distributions the indicial lift functions applicable to unsteady lift problems are determined for two cases. Results are presented which permit the determination of maximum increment in lift coefficient attained by an unrestrained airfoil during its flight through a gust. As an application of these results, the minimum altitude for safe flight through a specific gust is calculated for a particular supersonic wing of given strength and wing loading.

High-resolution two dimensional advective transport

USGS Publications Warehouse

Smith, P.E.; Larock, B.E.

1989-01-01

The paper describes a two-dimensional high-resolution scheme for advective transport that is based on a Eulerian-Lagrangian method with a flux limiter. The scheme is applied to the problem of pure-advection of a rotated Gaussian hill and shown to preserve the monotonicity property of the governing conservation law.

An iterative bidirectional heuristic placement algorithm for solving the two-dimensional knapsack packing problem

NASA Astrophysics Data System (ADS)

Shiangjen, Kanokwatt; Chaijaruwanich, Jeerayut; Srisujjalertwaja, Wijak; Unachak, Prakarn; Somhom, Samerkae

2018-02-01

This article presents an efficient heuristic placement algorithm, namely, a bidirectional heuristic placement, for solving the two-dimensional rectangular knapsack packing problem. The heuristic demonstrates ways to maximize space utilization by fitting the appropriate rectangle from both sides of the wall of the current residual space layer by layer. The iterative local search along with a shift strategy is developed and applied to the heuristic to balance the exploitation and exploration tasks in the solution space without the tuning of any parameters. The experimental results on many scales of packing problems show that this approach can produce high-quality solutions for most of the benchmark datasets, especially for large-scale problems, within a reasonable duration of computational time.

Asymptotic analysis of the narrow escape problem in dendritic spine shaped domain: three dimensions

NASA Astrophysics Data System (ADS)

Li, Xiaofei; Lee, Hyundae; Wang, Yuliang

2017-08-01

This paper deals with the three-dimensional narrow escape problem in a dendritic spine shaped domain, which is composed of a relatively big head and a thin neck. The narrow escape problem is to compute the mean first passage time of Brownian particles traveling from inside the head to the end of the neck. The original model is to solve a mixed Dirichlet-Neumann boundary value problem for the Poisson equation in the composite domain, and is computationally challenging. In this paper we seek to transfer the original problem to a mixed Robin-Neumann boundary value problem by dropping the thin neck part, and rigorously derive the asymptotic expansion of the mean first passage time with high order terms. This study is a nontrivial three-dimensional generalization of the work in Li (2014 J. Phys. A: Math. Theor. 47 505202), where a two-dimensional analogue domain is considered.

A Hidden Surface Algorithm for Computer Generated Halftone Pictures

DTIC Science & Technology

converting data describing three-dimensional objects into data that can be used to generate two-dimensional halftone images. It deals with some problems that arise in black and white, and color shading.

Numerical simulation of three-dimensional transonic turbulent projectile aerodynamics by TVD schemes

NASA Technical Reports Server (NTRS)

Shiau, Nae-Haur; Hsu, Chen-Chi; Chyu, Wei-Jao

1989-01-01

The two-dimensional symmetric TVD scheme proposed by Yee has been extended to and investigated for three-dimensional thin-layer Navier-Stokes simulation of complex aerodynamic problems. An existing three-dimensional Navier-stokes code based on the beam and warming algorithm is modified to provide an option of using the TVD algorithm and the flow problem considered is a transonic turbulent flow past a projectile with sting at ten-degree angle of attack. Numerical experiments conducted for three flow cases, free-stream Mach numbers of 0.91, 0.96 and 1.20 show that the symmetric TVD algorithm can provide surface pressure distribution in excellent agreement with measured data; moreover, the rate of convergence to attain a steady state solution is about two times faster than the original beam and warming algorithm.

High Performance Parallel Analysis of Coupled Problems for Aircraft Propulsion

NASA Technical Reports Server (NTRS)

Felippa, C. A.; Farhat, C.; Lanteri, S.; Maman, N.; Piperno, S.; Gumaste, U.

1994-01-01

In order to predict the dynamic response of a flexible structure in a fluid flow, the equations of motion of the structure and the fluid must be solved simultaneously. In this paper, we present several partitioned procedures for time-integrating this focus coupled problem and discuss their merits in terms of accuracy, stability, heterogeneous computing, I/O transfers, subcycling, and parallel processing. All theoretical results are derived for a one-dimensional piston model problem with a compressible flow, because the complete three-dimensional aeroelastic problem is difficult to analyze mathematically. However, the insight gained from the analysis of the coupled piston problem and the conclusions drawn from its numerical investigation are confirmed with the numerical simulation of the two-dimensional transient aeroelastic response of a flexible panel in a transonic nonlinear Euler flow regime.

Multiexponential models of (1+1)-dimensional dilaton gravity and Toda-Liouville integrable models

NASA Astrophysics Data System (ADS)

de Alfaro, V.; Filippov, A. T.

2010-01-01

We study general properties of a class of two-dimensional dilaton gravity (DG) theories with potentials containing several exponential terms. We isolate and thoroughly study a subclass of such theories in which the equations of motion reduce to Toda and Liouville equations. We show that the equation parameters must satisfy a certain constraint, which we find and solve for the most general multiexponential model. It follows from the constraint that integrable Toda equations in DG theories generally cannot appear without accompanying Liouville equations. The most difficult problem in the two-dimensional Toda-Liouville (TL) DG is to solve the energy and momentum constraints. We discuss this problem using the simplest examples and identify the main obstacles to solving it analytically. We then consider a subclass of integrable two-dimensional theories where scalar matter fields satisfy the Toda equations and the two-dimensional metric is trivial. We consider the simplest case in some detail. In this example, we show how to obtain the general solution. We also show how to simply derive wavelike solutions of general TL systems. In the DG theory, these solutions describe nonlinear waves coupled to gravity and also static states and cosmologies. For static states and cosmologies, we propose and study a more general one-dimensional TL model typically emerging in one-dimensional reductions of higher-dimensional gravity and supergravity theories. We especially attend to making the analytic structure of the solutions of the Toda equations as simple and transparent as possible.

Coupling Finite Element and Meshless Local Petrov-Galerkin Methods for Two-Dimensional Potential Problems

NASA Technical Reports Server (NTRS)

Chen, T.; Raju, I. S.

2002-01-01

A coupled finite element (FE) method and meshless local Petrov-Galerkin (MLPG) method for analyzing two-dimensional potential problems is presented in this paper. The analysis domain is subdivided into two regions, a finite element (FE) region and a meshless (MM) region. A single weighted residual form is written for the entire domain. Independent trial and test functions are assumed in the FE and MM regions. A transition region is created between the two regions. The transition region blends the trial and test functions of the FE and MM regions. The trial function blending is achieved using a technique similar to the 'Coons patch' method that is widely used in computer-aided geometric design. The test function blending is achieved by using either FE or MM test functions on the nodes in the transition element. The technique was evaluated by applying the coupled method to two potential problems governed by the Poisson equation. The coupled method passed all the patch test problems and gave accurate solutions for the problems studied.

Iterative spectral methods and spectral solutions to compressible flows

NASA Technical Reports Server (NTRS)

Hussaini, M. Y.; Zang, T. A.

1982-01-01

A spectral multigrid scheme is described which can solve pseudospectral discretizations of self-adjoint elliptic problems in O(N log N) operations. An iterative technique for efficiently implementing semi-implicit time-stepping for pseudospectral discretizations of Navier-Stokes equations is discussed. This approach can handle variable coefficient terms in an effective manner. Pseudospectral solutions of compressible flow problems are presented. These include one dimensional problems and two dimensional Euler solutions. Results are given both for shock-capturing approaches and for shock-fitting ones.

On the Origins of Vortex Shedding in Two-dimensional Incompressible Flows

PubMed Central

Boghosian, M. E.; Cassel, K. W.

2016-01-01

An exegesis of a novel mechanism leading to vortex splitting and subsequent shedding that is valid for two-dimensional incompressible, inviscid or viscous, and external or internal or wall-bounded flows, is detailed in this research. The mechanism, termed the Vortex-Shedding Mechanism (VSM), is simple and intuitive, requiring only two coincident conditions in the flow: (1) the existence of a location with zero momentum and (2) the presence of a net force having a positive divergence. Numerical solutions of several model problems illustrate causality of the VSM. Moreover, the VSM criteria is proved to be a necessary and sufficient condition for a vortex splitting event in any two-dimensional, incompressible flow. The VSM is shown to exist in several canonical problems including the external flow past a circular cylinder. Suppression of the von Kármán vortex street is demonstrated for Reynolds numbers of 100 and 400 by mitigating the VSM. PMID:27795617

On the Origins of Vortex Shedding in Two-dimensional Incompressible Flows.

PubMed

Boghosian, M E; Cassel, K W

2016-12-01

An exegesis of a novel mechanism leading to vortex splitting and subsequent shedding that is valid for two-dimensional incompressible, inviscid or viscous, and external or internal or wall-bounded flows, is detailed in this research. The mechanism, termed the Vortex-Shedding Mechanism (VSM), is simple and intuitive, requiring only two coincident conditions in the flow: (1) the existence of a location with zero momentum and (2) the presence of a net force having a positive divergence. Numerical solutions of several model problems illustrate causality of the VSM. Moreover, the VSM criteria is proved to be a necessary and sufficient condition for a vortex splitting event in any two-dimensional, incompressible flow. The VSM is shown to exist in several canonical problems including the external flow past a circular cylinder. Suppression of the von Kármán vortex street is demonstrated for Reynolds numbers of 100 and 400 by mitigating the VSM.

Apollo 11 - Prime and Backup Crews - Geology Training - TX

NASA Image and Video Library

1969-03-03

S69-25199 (25 Feb. 1969) --- Two Apollo 11 astronauts study a rock specimen during a geological field trip to the Quitman Mountains area near the Fort Quitman ruins in far west Texas. On the left is James A. Lovell Jr., Apollo 11 backup crew commander; and on the right is Fred W. Haise Jr., backup crew lunar module pilot. Lovell holds a camera which was used in simulating taking pictures of actual lunar samples on the surface of the Moon.

The onset of layer undulations in smectic A liquid crystals due to a strong magnetic field

NASA Astrophysics Data System (ADS)

Contreras, A.; Garcia-Azpeitia, C.; García-Cervera, C. J.; Joo, S.

2016-08-01

We investigate the effect of a strong magnetic field on a three dimensional smectic A liquid crystal. We identify a critical field above which the uniform layered state loses stability; this is associated to the onset of layer undulations. In a previous work García-Cervera and Joo (2012 Arch. Ration. Mech. Anal. 203 1-43), García-Cervera and Joo considered the two dimensional case and analyzed the transition to the undulated state via a simple bifurcation. In dimension n  =  3 the situation is more delicate because the first eigenvalue of the corresponding linearized problem is not simple. We overcome the difficulties inherent to this higher dimensional setting by identifying the irreducible representations for natural actions on the functional that take into account the invariances of the problem thus allowing for reducing the bifurcation analysis to a subspace with symmetries. We are able to describe at least two bifurcation branches, highlighting the richer landscape of energy critical states in the three dimensional setting. Finally, we analyze a reduced two dimensional problem, assuming the magnetic field is very strong, and are able to relate this to a model in micromagnetics studied in Alouges et al (2002 ESAIM Control Optim. Calc. Var. 8 31-68), from where we deduce the periodicity property of minimizers.

Aerodynamic Analyses Requiring Advanced Computers, part 2

NASA Technical Reports Server (NTRS)

1975-01-01

Papers given at the conference present the results of theoretical research on aerodynamic flow problems requiring the use of advanced computers. Topics discussed include two-dimensional configurations, three-dimensional configurations, transonic aircraft, and the space shuttle.

An incompressible two-dimensional multiphase particle-in-cell model for dense particle flows

DOE Office of Scientific and Technical Information (OSTI.GOV)

Snider, D.M.; O`Rourke, P.J.; Andrews, M.J.

1997-06-01

A two-dimensional, incompressible, multiphase particle-in-cell (MP-PIC) method is presented for dense particle flows. The numerical technique solves the governing equations of the fluid phase using a continuum model and those of the particle phase using a Lagrangian model. Difficulties associated with calculating interparticle interactions for dense particle flows with volume fractions above 5% have been eliminated by mapping particle properties to a Eulerian grid and then mapping back computed stress tensors to particle positions. This approach utilizes the best of Eulerian/Eulerian continuum models and Eulerian/Lagrangian discrete models. The solution scheme allows for distributions of types, sizes, and density of particles,more » with no numerical diffusion from the Lagrangian particle calculations. The computational method is implicit with respect to pressure, velocity, and volume fraction in the continuum solution thus avoiding courant limits on computational time advancement. MP-PIC simulations are compared with one-dimensional problems that have analytical solutions and with two-dimensional problems for which there are experimental data.« less

Three Dimensional Modeling via Photographs for Documentation of a Village Bath

NASA Astrophysics Data System (ADS)

Balta, H. B.; Hamamcioglu-Turan, M.; Ocali, O.

2013-07-01

The aim of this study is supporting the conceptual discussions of architectural restoration with three dimensional modeling of monuments based on photogrammetric survey. In this study, a 16th century village bath in Ulamış, Seferihisar, and Izmir is modeled for documentation. Ulamış is one of the historical villages within which Turkish population first settled in the region of Seferihisar - Urla. The methodology was tested on an antique monument; a bath with a cubical form. Within the limits of this study, only the exterior of the bath was modeled. The presentation scale for the bath was determined as 1 / 50, considering the necessities of designing structural interventions and architectural ones within the scope of a restoration project. The three dimensional model produced is a realistic document presenting the present situation of the ruin. Traditional plan, elevation and perspective drawings may be produced from the model, in addition to the realistic textured renderings and wireframe representations. The model developed in this study provides opportunity for presenting photorealistic details of historical morphologies in scale. Compared to conventional drawings, the renders based on the 3d models provide an opportunity for conceiving architectural details such as color, material and texture. From these documents, relatively more detailed restitution hypothesis can be developed and intervention decisions can be taken. Finally, the principles derived from the case study can be used for 3d documentation of historical structures with irregular surfaces.

Coupling continuous damage and debris fragmentation for energy absorption prediction by cfrp structures during crushing

NASA Astrophysics Data System (ADS)

Espinosa, Christine; Lachaud, Frédéric; Limido, Jérome; Lacome, Jean-Luc; Bisson, Antoine; Charlotte, Miguel

2015-05-01

Energy absorption during crushing is evaluated using a thermodynamic based continuum damage model inspired from the Matzenmiller-Lubliner-Taylors model. It was found that for crash-worthiness applications, it is necessary to couple the progressive ruin of the material to a representation of the matter openings and debris generation. Element kill technique (erosion) and/or cohesive elements are efficient but not predictive. A technique switching finite elements into discrete particles at rupture is used to create debris and accumulated mater during the crushing of the structure. Switching criteria are evaluated using the contribution of the different ruin modes in the damage evolution, energy absorption, and reaction force generation.

The program FANS-3D (finite analytic numerical simulation 3-dimensional) and its applications

NASA Technical Reports Server (NTRS)

Bravo, Ramiro H.; Chen, Ching-Jen

1992-01-01

In this study, the program named FANS-3D (Finite Analytic Numerical Simulation-3 Dimensional) is presented. FANS-3D was designed to solve problems of incompressible fluid flow and combined modes of heat transfer. It solves problems with conduction and convection modes of heat transfer in laminar flow, with provisions for radiation and turbulent flows. It can solve singular or conjugate modes of heat transfer. It also solves problems in natural convection, using the Boussinesq approximation. FANS-3D was designed to solve heat transfer problems inside one, two and three dimensional geometries that can be represented by orthogonal planes in a Cartesian coordinate system. It can solve internal and external flows using appropriate boundary conditions such as symmetric, periodic and user specified.

The two-dimensional Stefan problem with slightly varying heat flux

DOE Office of Scientific and Technical Information (OSTI.GOV)

Gammon, J.; Howarth, J.A.

1995-09-01

The authors solve the two-dimensional stefan problem of solidification in a half-space, where the heat flux at the wall is a slightly varying function of positioning along the wall, by means of a large Stefan number approximation (which turns out to be equivalent to a small time solution), and then by means of the Heat Balance Integral Method, which is valid for all time, and which agrees with the large Stefan number solution for small times. A representative solution is given for a particular form of the heat flux perturbation.

Solving time-dependent two-dimensional eddy current problems

NASA Technical Reports Server (NTRS)

Lee, Min Eig; Hariharan, S. I.; Ida, Nathan

1990-01-01

Transient eddy current calculations are presented for an EM wave-scattering and field-penetrating case in which a two-dimensional transverse magnetic field is incident on a good (i.e., not perfect) and infinitely long conductor. The problem thus posed is of initial boundary-value interface type, where the boundary of the conductor constitutes the interface. A potential function is used for time-domain modeling of the situation, and finite difference-time domain techniques are used to march the potential function explicitly in time. Attention is given to the case of LF radiation conditions.

Scientific literature: can we keep it from ruining science

DOE Office of Scientific and Technical Information (OSTI.GOV)

Lentzner, H.L.

1978-08-24

The scientific literature is failing us, because today's scientific writing is pompous, murky, and dull. Since the U.S. will spend 12 billion dollars on written scientific communication in 1978, this failure represents a significant waste of resources. Some problems can be attributed to training; university English courses and laboratory reports do not prepare science students for the type of writing required after graduation. The specialized literature has done its part by enshrining bad writing. What is to be done. Authors must carefully define their audience, thoughtfully organize their material, use concrete words, use the active voice, and take time tomore » revise. Teachers and referees should be less tolerant of students and colleagues who write badly. 2 figures.« less

Maximizing kinetic energy transfer in one-dimensional many-body collisions

NASA Astrophysics Data System (ADS)

Ricardo, Bernard; Lee, Paul

2015-03-01

The main problem discussed in this paper involves a simple one-dimensional two-body collision, in which the problem can be extended into a chain of one-dimensional many-body collisions. The result is quite interesting, as it provides us with a thorough mathematical understanding that will help in designing a chain system for maximum energy transfer for a range of collision types. In this paper, we will show that there is a way to improve the kinetic energy transfer between two masses, and the idea can be applied recursively. However, this method only works for a certain range of collision types, which is indicated by a range of coefficients of restitution. Although the concept of momentum, elastic and inelastic collision, as well as Newton’s laws, are taught in junior college physics, especially in Singapore schools, students in this level are not expected to be able to do this problem quantitatively, as it requires rigorous mathematics, including calculus. Nevertheless, this paper provides nice analytical steps that address some common misconceptions in students’ way of thinking about one-dimensional collisions.

Inverse regression-based uncertainty quantification algorithms for high-dimensional models: Theory and practice

NASA Astrophysics Data System (ADS)

Li, Weixuan; Lin, Guang; Li, Bing

2016-09-01

Many uncertainty quantification (UQ) approaches suffer from the curse of dimensionality, that is, their computational costs become intractable for problems involving a large number of uncertainty parameters. In these situations, the classic Monte Carlo often remains the preferred method of choice because its convergence rate O (n - 1 / 2), where n is the required number of model simulations, does not depend on the dimension of the problem. However, many high-dimensional UQ problems are intrinsically low-dimensional, because the variation of the quantity of interest (QoI) is often caused by only a few latent parameters varying within a low-dimensional subspace, known as the sufficient dimension reduction (SDR) subspace in the statistics literature. Motivated by this observation, we propose two inverse regression-based UQ algorithms (IRUQ) for high-dimensional problems. Both algorithms use inverse regression to convert the original high-dimensional problem to a low-dimensional one, which is then efficiently solved by building a response surface for the reduced model, for example via the polynomial chaos expansion. The first algorithm, which is for the situations where an exact SDR subspace exists, is proved to converge at rate O (n-1), hence much faster than MC. The second algorithm, which doesn't require an exact SDR, employs the reduced model as a control variate to reduce the error of the MC estimate. The accuracy gain could still be significant, depending on how well the reduced model approximates the original high-dimensional one. IRUQ also provides several additional practical advantages: it is non-intrusive; it does not require computing the high-dimensional gradient of the QoI; and it reports an error bar so the user knows how reliable the result is.

Riemann–Hilbert problem approach for two-dimensional flow inverse scattering

DOE Office of Scientific and Technical Information (OSTI.GOV)

Agaltsov, A. D., E-mail: agalets@gmail.com; Novikov, R. G., E-mail: novikov@cmap.polytechnique.fr; IEPT RAS, 117997 Moscow

2014-10-15

We consider inverse scattering for the time-harmonic wave equation with first-order perturbation in two dimensions. This problem arises in particular in the acoustic tomography of moving fluid. We consider linearized and nonlinearized reconstruction algorithms for this problem of inverse scattering. Our nonlinearized reconstruction algorithm is based on the non-local Riemann–Hilbert problem approach. Comparisons with preceding results are given.

Well-posedness of the Cauchy problem for models of large amplitude internal waves

NASA Astrophysics Data System (ADS)

Guyenne, Philippe; Lannes, David; Saut, Jean-Claude

2010-02-01

We consider in this paper the 'shallow-water/shallow-water' asymptotic model obtained in Choi and Camassa (1999 J. Fluid Mech. 396 1-36), Craig et al (2005 Commun. Pure. Appl. Math. 58 1587-641) (one-dimensional interface) and Bona et al (2008 J. Math. Pures Appl. 89 538-66) (two-dimensional interface) from the two-layer system with rigid lid, for the description of large amplitude internal waves at the interface of two layers of immiscible fluids of different densities. For one-dimensional interfaces, this system is of hyperbolic type and its local well-posedness does not raise serious difficulties, although other issues (blow-up, loss of hyperbolicity, etc) turn out to be delicate. For two-dimensional interfaces, the system is nonlocal. Nevertheless, we prove that it conserves some properties of 'hyperbolic type' and show that the associated Cauchy problem is locally well posed in suitable Sobolev classes provided some natural restrictions are imposed on the data. These results are illustrated by numerical simulations with emphasis on the formation of shock waves.

Finite dimensional approximation of a class of constrained nonlinear optimal control problems

NASA Technical Reports Server (NTRS)

Gunzburger, Max D.; Hou, L. S.

1994-01-01

An abstract framework for the analysis and approximation of a class of nonlinear optimal control and optimization problems is constructed. Nonlinearities occur in both the objective functional and in the constraints. The framework includes an abstract nonlinear optimization problem posed on infinite dimensional spaces, and approximate problem posed on finite dimensional spaces, together with a number of hypotheses concerning the two problems. The framework is used to show that optimal solutions exist, to show that Lagrange multipliers may be used to enforce the constraints, to derive an optimality system from which optimal states and controls may be deduced, and to derive existence results and error estimates for solutions of the approximate problem. The abstract framework and the results derived from that framework are then applied to three concrete control or optimization problems and their approximation by finite element methods. The first involves the von Karman plate equations of nonlinear elasticity, the second, the Ginzburg-Landau equations of superconductivity, and the third, the Navier-Stokes equations for incompressible, viscous flows.

The application of Green's theorem to the solution of boundary-value problems in linearized supersonic wing theory

NASA Technical Reports Server (NTRS)

Heaslet, Max A; Lomax, Harvard

1950-01-01

Following the introduction of the linearized partial differential equation for nonsteady three-dimensional compressible flow, general methods of solution are given for the two and three-dimensional steady-state and two-dimensional unsteady-state equations. It is also pointed out that, in the absence of thickness effects, linear theory yields solutions consistent with the assumptions made when applied to lifting-surface problems for swept-back plan forms at sonic speeds. The solutions of the particular equations are determined in all cases by means of Green's theorem, and thus depend on the use of Green's equivalent layer of sources, sinks, and doublets. Improper integrals in the supersonic theory are treated by means of Hadamard's "finite part" technique.

How Students Solve Problems in Spatial Geometry while Using a Software Application for Visualizing 3D Geometric Objects

ERIC Educational Resources Information Center

Widder, Mirela; Gorsky, Paul

2013-01-01

In schools, learning spatial geometry is usually dependent upon a student's ability to visualize three dimensional geometric configurations from two dimensional drawings. Such a process, however, often creates visual obstacles which are unique to spatial geometry. Useful software programs which realistically depict three dimensional geometric…

Simplified computational methods for elastic and elastic-plastic fracture problems

NASA Technical Reports Server (NTRS)

Atluri, Satya N.

1992-01-01

An overview is given of some of the recent (1984-1991) developments in computational/analytical methods in the mechanics of fractures. Topics covered include analytical solutions for elliptical or circular cracks embedded in isotropic or transversely isotropic solids, with crack faces being subjected to arbitrary tractions; finite element or boundary element alternating methods for two or three dimensional crack problems; a 'direct stiffness' method for stiffened panels with flexible fasteners and with multiple cracks; multiple site damage near a row of fastener holes; an analysis of cracks with bonded repair patches; methods for the generation of weight functions for two and three dimensional crack problems; and domain-integral methods for elastic-plastic or inelastic crack mechanics.

Difference equation state approximations for nonlinear hereditary control problems

NASA Technical Reports Server (NTRS)

Rosen, I. G.

1982-01-01

Discrete approximation schemes for the solution of nonlinear hereditary control problems are constructed. The methods involve approximation by a sequence of optimal control problems in which the original infinite dimensional state equation has been approximated by a finite dimensional discrete difference equation. Convergence of the state approximations is argued using linear semigroup theory and is then used to demonstrate that solutions to the approximating optimal control problems in some sense approximate solutions to the original control problem. Two schemes, one based upon piecewise constant approximation, and the other involving spline functions are discussed. Numerical results are presented, analyzed and used to compare the schemes to other available approximation methods for the solution of hereditary control problems.

Analytical solutions of the two-dimensional Dirac equation for a topological channel intersection

NASA Astrophysics Data System (ADS)

Anglin, J. R.; Schulz, A.

2017-01-01

Numerical simulations in a tight-binding model have shown that an intersection of topologically protected one-dimensional chiral channels can function as a beam splitter for noninteracting fermions on a two-dimensional lattice [Qiao, Jung, and MacDonald, Nano Lett. 11, 3453 (2011), 10.1021/nl201941f; Qiao et al., Phys. Rev. Lett. 112, 206601 (2014), 10.1103/PhysRevLett.112.206601]. Here we confirm this result analytically in the corresponding continuum k .p model, by solving the associated two-dimensional Dirac equation, in the presence of a "checkerboard" potential that provides a right-angled intersection between two zero-line modes. The method by which we obtain our analytical solutions is systematic and potentially generalizable to similar problems involving intersections of one-dimensional systems.

A geometrical multi-scale numerical method for coupled hygro-thermo-mechanical problems in photovoltaic laminates.

PubMed

Lenarda, P; Paggi, M

A comprehensive computational framework based on the finite element method for the simulation of coupled hygro-thermo-mechanical problems in photovoltaic laminates is herein proposed. While the thermo-mechanical problem takes place in the three-dimensional space of the laminate, moisture diffusion occurs in a two-dimensional domain represented by the polymeric layers and by the vertical channel cracks in the solar cells. Therefore, a geometrical multi-scale solution strategy is pursued by solving the partial differential equations governing heat transfer and thermo-elasticity in the three-dimensional space, and the partial differential equation for moisture diffusion in the two dimensional domains. By exploiting a staggered scheme, the thermo-mechanical problem is solved first via a fully implicit solution scheme in space and time, with a specific treatment of the polymeric layers as zero-thickness interfaces whose constitutive response is governed by a novel thermo-visco-elastic cohesive zone model based on fractional calculus. Temperature and relative displacements along the domains where moisture diffusion takes place are then projected to the finite element model of diffusion, coupled with the thermo-mechanical problem by the temperature and crack opening dependent diffusion coefficient. The application of the proposed method to photovoltaic modules pinpoints two important physical aspects: (i) moisture diffusion in humidity freeze tests with a temperature dependent diffusivity is a much slower process than in the case of a constant diffusion coefficient; (ii) channel cracks through Silicon solar cells significantly enhance moisture diffusion and electric degradation, as confirmed by experimental tests.

Solitary wave solutions of two-dimensional nonlinear Kadomtsev-Petviashvili dynamic equation in dust-acoustic plasmas

NASA Astrophysics Data System (ADS)

Seadawy, Aly R.

2017-09-01

Nonlinear two-dimensional Kadomtsev-Petviashvili (KP) equation governs the behaviour of nonlinear waves in dusty plasmas with variable dust charge and two temperature ions. By using the reductive perturbation method, the two-dimensional dust-acoustic solitary waves (DASWs) in unmagnetized cold plasma consisting of dust fluid, ions and electrons lead to a KP equation. We derived the solitary travelling wave solutions of the two-dimensional nonlinear KP equation by implementing sech-tanh, sinh-cosh, extended direct algebraic and fraction direct algebraic methods. We found the electrostatic field potential and electric field in the form travelling wave solutions for two-dimensional nonlinear KP equation. The solutions for the KP equation obtained by using these methods can be demonstrated precisely and efficiency. As an illustration, we used the readymade package of Mathematica program 10.1 to solve the original problem. These solutions are in good agreement with the analytical one.

Experiences with explicit finite-difference schemes for complex fluid dynamics problems on STAR-100 and CYBER-203 computers

NASA Technical Reports Server (NTRS)

Kumar, A.; Rudy, D. H.; Drummond, J. P.; Harris, J. E.

1982-01-01

Several two- and three-dimensional external and internal flow problems solved on the STAR-100 and CYBER-203 vector processing computers are described. The flow field was described by the full Navier-Stokes equations which were then solved by explicit finite-difference algorithms. Problem results and computer system requirements are presented. Program organization and data base structure for three-dimensional computer codes which will eliminate or improve on page faulting, are discussed. Storage requirements for three-dimensional codes are reduced by calculating transformation metric data in each step. As a result, in-core grid points were increased in number by 50% to 150,000, with a 10% execution time increase. An assessment of current and future machine requirements shows that even on the CYBER-205 computer only a few problems can be solved realistically. Estimates reveal that the present situation is more storage limited than compute rate limited, but advancements in both storage and speed are essential to realistically calculate three-dimensional flow.

Vascular Plant and Vertebrate Inventory of Casa Grande Ruins National Monument

USGS Publications Warehouse

Powell, Brian F.; Albrecht, Eric W.; Schmidt, Cecilia A.; Halvorson, William L.; Anning, Pamela; Docherty, Kathleen

2006-01-01

Executive Summary This report summarizes results of the first comprehensive biological inventory of Casa Grande Ruins National Monument (NM) in southern Arizona. Surveys at the monument were part of a larger effort to inventory vascular plants and vertebrates in eight National Park Service units in Arizona and New Mexico. In 2001 and 2002 we surveyed for vascular plants and vertebrates (amphibians, reptiles, birds, and mammals) at Casa Grande Ruins NM to document the presence, and in some cases relative abundance, of species. By using repeatable study designs and standardized field techniques, which included quantified survey effort, we produced inventories that can serve as the basis for a biological monitoring program. Of the National Park Service units in the region, no other has experienced as much recent ecological change as Casa Grande Ruins NM. Once situated in a large and biologically diverse mesquite bosque near the perennially flowing Gila River, the monument is now a patch of sparse desert vegetation surrounded by urban and commercial development that is rapidly replacing agriculture as the dominant land use in the area. Roads, highways, and canals surround the monument. Development, and its associated impacts, has important implications for the plants and animals that live in the monument. The plant species list is small and the distribution and number of non-native plants appears to be increasing. Terrestrial vertebrates are also being impacted by the changing landscape, which is increasing the isolation of these populations from nearby natural areas and thereby reducing the number of species at the monument. These observations are alarming and are based on our review of previous studies, our research in the monument, and our knowledge of the biogeography and ecology of the Sonoran Desert. Together, these data suggest that the monument has lost a significant portion of its historic complement of species and these changes will likely intensify as urbanization continues. Despite isolation of the monument from nearby natural areas, we recorded noteworthy species or observations for all taxonomic groups: * Plants: night-blooming cereus * Amphibians: high abundance of Couch's spadefoot toads * Reptiles: high abundance of long-nosed snakes * Birds: 10 species of diurnal raptors including 4 species of falcons * Mammals: American badger This study is a first step in the process of compiling information about the biological resources of Casa Grande Ruins NM and surrounding areas. We recommend additional inventory and research studies, and we identify aspects of our effort that could be improved upon through application of new techniques or by extending the temporal (and possibly spatial) scope of our work.

On some structure-turbulence interaction problems

NASA Technical Reports Server (NTRS)

Maekawa, S.; Lin, Y. K.

1976-01-01

The interactions between a turbulent flow structure; responding to its excitation were studied. The turbulence was typical of those associated with a boundary layer, having a cross-spectral density indicative of convection and statistical decay. A number of structural models were considered. Among the one-dimensional models were an unsupported infinite beam and a periodically supported infinite beam. The fuselage construction of an aircraft was then considered. For the two-dimensional case a simple membrane was used to illustrate the type of formulation applicable to most two-dimensional structures. Both the one-dimensional and two-dimensional structures studied were backed by a cavity filled with an initially quiescent fluid to simulate the acoustic environment when the structure forms one side of a cabin of a sea vessel or aircraft.

The resistance of an n-dimensional tetrahedron

NASA Astrophysics Data System (ADS)

Griffiths, Martin

2013-01-01

We consider here a problem that is suitable for introducing high-school students to the notion of generalizing shapes and solids to n dimensions. In particular, we calculate the effective resistance between any two vertices of an n-dimensional tetrahedron whose edges are each 1-Ω resistors. This leads, in a natural way, to more demanding problems, and indeed ideas for more advanced work in this area are also suggested.

Oscillations and stability of numerical solutions of the heat conduction equation

NASA Technical Reports Server (NTRS)

Kozdoba, L. A.; Levi, E. V.

1976-01-01

The mathematical model and results of numerical solutions are given for the one dimensional problem when the linear equations are written in a rectangular coordinate system. All the computations are easily realizable for two and three dimensional problems when the equations are written in any coordinate system. Explicit and implicit schemes are shown in tabular form for stability and oscillations criteria; the initial temperature distribution is considered uniform.

Gun Violence, Mental Illness, And Laws That Prohibit Gun Possession: Evidence From Two Florida Counties.

PubMed

Swanson, Jeffrey W; Easter, Michele M; Robertson, Allison G; Swartz, Marvin S; Alanis-Hirsch, Kelly; Moseley, Daniel; Dion, Charles; Petrila, John

2016-06-01

Gun violence kills about ninety people every day in the United States, a toll measured in wasted and ruined lives and with an annual economic price tag exceeding $200 billion. Some policy makers suggest that reforming mental health care systems and improving point-of-purchase background checks to keep guns from mentally disturbed people will address the problem. Epidemiological research shows that serious mental illness contributes little to the risk of interpersonal violence but is a strong factor in suicide, which accounts for most firearm fatalities. Meanwhile, the effectiveness of gun restrictions focused on mental illness remains poorly understood. This article examines gun-related suicide and violent crime in people with serious mental illnesses, and whether legal restrictions on firearm sales to people with a history of mental health adjudication are effective in preventing gun violence. Among the study population in two large Florida counties, we found that 62 percent of violent gun crime arrests and 28 percent of gun suicides involved individuals not legally permitted to have a gun at the time. Suggested policy reforms include enacting risk-based gun removal laws and prohibiting guns from people involuntarily detained in short-term psychiatric hospitalizations. Project HOPE—The People-to-People Health Foundation, Inc.

Computer programs for calculating two-dimensional potential flow in and about propulsion system inlets

NASA Technical Reports Server (NTRS)

Hawk, J. D.; Stockman, N. O.; Farrell, C. A., Jr.

1978-01-01

Incompressible potential flow calculations are presented that were corrected for compressibility in two-dimensional inlets at arbitrary operating conditions. Included are a statement of the problem to be solved, a description of each of the computer programs, and sufficient documentation, including a test case, to enable a user to run the program.

Computer programs for calculating two-dimensional potential flow through deflected nozzles

NASA Technical Reports Server (NTRS)

Hawk, J. D.; Stockman, N. O.

1979-01-01

Computer programs to calculate the incompressible potential flow, corrected for compressibility, in two-dimensional nozzles at arbitrary operating conditions are presented. A statement of the problem to be solved, a description of each of the computer programs, and sufficient documentation, including a test case, to enable a user to run the program are included.

Learning Relative Motion Concepts in Immersive and Non-immersive Virtual Environments

NASA Astrophysics Data System (ADS)

Kozhevnikov, Michael; Gurlitt, Johannes; Kozhevnikov, Maria

2013-12-01

The focus of the current study is to understand which unique features of an immersive virtual reality environment have the potential to improve learning relative motion concepts. Thirty-seven undergraduate students learned relative motion concepts using computer simulation either in immersive virtual environment (IVE) or non-immersive desktop virtual environment (DVE) conditions. Our results show that after the simulation activities, both IVE and DVE groups exhibited a significant shift toward a scientific understanding in their conceptual models and epistemological beliefs about the nature of relative motion, and also a significant improvement on relative motion problem-solving tests. In addition, we analyzed students' performance on one-dimensional and two-dimensional questions in the relative motion problem-solving test separately and found that after training in the simulation, the IVE group performed significantly better than the DVE group on solving two-dimensional relative motion problems. We suggest that egocentric encoding of the scene in IVE (where the learner constitutes a part of a scene they are immersed in), as compared to allocentric encoding on a computer screen in DVE (where the learner is looking at the scene from "outside"), is more beneficial than DVE for studying more complex (two-dimensional) relative motion problems. Overall, our findings suggest that such aspects of virtual realities as immersivity, first-hand experience, and the possibility of changing different frames of reference can facilitate understanding abstract scientific phenomena and help in displacing intuitive misconceptions with more accurate mental models.

Stability of a flow down an incline with respect to two-dimensional and three-dimensional disturbances for Newtonian and non-Newtonian fluids.

PubMed

Allouche, M H; Millet, S; Botton, V; Henry, D; Ben Hadid, H; Rousset, F

2015-12-01

Squire's theorem, which states that the two-dimensional instabilities are more dangerous than the three-dimensional instabilities, is revisited here for a flow down an incline, making use of numerical stability analysis and Squire relationships when available. For flows down inclined planes, one of these Squire relationships involves the slopes of the inclines. This means that the Reynolds number associated with a two-dimensional wave can be shown to be smaller than that for an oblique wave, but this oblique wave being obtained for a larger slope. Physically speaking, this prevents the possibility to directly compare the thresholds at a given slope. The goal of the paper is then to reach a conclusion about the predominance or not of two-dimensional instabilities at a given slope, which is of practical interest for industrial or environmental applications. For a Newtonian fluid, it is shown that, for a given slope, oblique wave instabilities are never the dominant instabilities. Both the Squire relationships and the particular variations of the two-dimensional wave critical curve with regard to the inclination angle are involved in the proof of this result. For a generalized Newtonian fluid, a similar result can only be obtained for a reduced stability problem where some term connected to the perturbation of viscosity is neglected. For the general stability problem, however, no Squire relationships can be derived and the numerical stability results show that the thresholds for oblique waves can be smaller than the thresholds for two-dimensional waves at a given slope, particularly for large obliquity angles and strong shear-thinning behaviors. The conclusion is then completely different in that case: the dominant instability for a generalized Newtonian fluid flowing down an inclined plane with a given slope can be three dimensional.

An adaptive moving mesh method for two-dimensional ideal magnetohydrodynamics

NASA Astrophysics Data System (ADS)

Han, Jianqiang; Tang, Huazhong

2007-01-01

This paper presents an adaptive moving mesh algorithm for two-dimensional (2D) ideal magnetohydrodynamics (MHD) that utilizes a staggered constrained transport technique to keep the magnetic field divergence-free. The algorithm consists of two independent parts: MHD evolution and mesh-redistribution. The first part is a high-resolution, divergence-free, shock-capturing scheme on a fixed quadrangular mesh, while the second part is an iterative procedure. In each iteration, mesh points are first redistributed, and then a conservative-interpolation formula is used to calculate the remapped cell-averages of the mass, momentum, and total energy on the resulting new mesh; the magnetic potential is remapped to the new mesh in a non-conservative way and is reconstructed to give a divergence-free magnetic field on the new mesh. Several numerical examples are given to demonstrate that the proposed method can achieve high numerical accuracy, track and resolve strong shock waves in ideal MHD problems, and preserve divergence-free property of the magnetic field. Numerical examples include the smooth Alfvén wave problem, 2D and 2.5D shock tube problems, two rotor problems, the stringent blast problem, and the cloud-shock interaction problem.

Diffusion Characteristics of Upwind Schemes on Unstructured Triangulations

NASA Technical Reports Server (NTRS)

Wood, William A.; Kleb, William L.

1998-01-01

The diffusive characteristics of two upwind schemes, multi-dimensional fluctuation splitting and dimensionally-split finite volume, are compared for scalar advection-diffusion problems. Algorithms for the two schemes are developed for node-based data representation on median-dual meshes associated with unstructured triangulations in two spatial dimensions. Four model equations are considered: linear advection, non-linear advection, diffusion, and advection-diffusion. Modular coding is employed to isolate the effects of the two approaches for upwind flux evaluation, allowing for head-to-head accuracy and efficiency comparisons. Both the stability of compressive limiters and the amount of artificial diffusion generated by the schemes is found to be grid-orientation dependent, with the fluctuation splitting scheme producing less artificial diffusion than the dimensionally-split finite volume scheme. Convergence rates are compared for the combined advection-diffusion problem, with a speedup of 2-3 seen for fluctuation splitting versus finite volume when solved on the same mesh. However, accurate solutions to problems with small diffusion coefficients can be achieved on coarser meshes using fluctuation splitting rather than finite volume, so that when comparing convergence rates to reach a given accuracy, fluctuation splitting shows a 20-25 speedup over finite volume.

Palmyra as it Once Was: 3d Virtual Reconstruction and Visualization of AN Irreplacable Lost Treasure

NASA Astrophysics Data System (ADS)

Denker, A.

2017-05-01

Palmyra was a mosaic which was composed through its colourful past, by Assyrians, Parthians, Greeks and Romans. For centuries, the spectacular ruins and impressive panorama of the antique city used to captivate and inspire the visitors as the witnesses of its illustrious history. As a grim consequence of the horrific conflict that engulfed Syria, since May 2015 they are no more to be seen. Palmyra has been purposely targeted and obliterated, the ruins have been reduced to rubble. The immense beauty and rich heritage of Palmyra which has been lost forever is reconstructed here as it was once was, at the top of its glory, with the hope of preserving its memory.

Hospitality, culture and regeneration: urban decay, entrepreneurship and the "ruin" bars of Budapest.

PubMed

Lugosi, Peter; Bell, David; Lugosi, Krisztina

2010-01-01

This paper considers the relationships between hospitality, culture and urban regeneration through an examination of rom (ruin) venues, which operate in dilapidated buildings in Budapest, Hungary. The paper reviews previous work on culture and urban regeneration in order to locate the role of hospitality within emerging debates. It subsequently interrogates the evolution of the rom phenomenon and demonstrates how, in this context, hospitality thrives because of social and physical decay in urban locations, how operators and entrepreneurs exploit conflicts among various actors involved in regeneration and how hospitality may be mobilised purposefully in the regeneration process. The paper demonstrates how networked entrepreneurship maintains these operations and how various forms of cultural production are entangled and mobilised in the venues' hospitality propositions.

Computer analysis of sound recordings from two Anasazi sites in northwestern New Mexico

NASA Astrophysics Data System (ADS)

Loose, Richard

2002-11-01

Sound recordings were made at a natural outdoor amphitheater in Chaco Canyon and in a reconstructed great kiva at Aztec Ruins. Recordings included computer-generated tones and swept sine waves, classical concert flute, Native American flute, conch shell trumpet, and prerecorded music. Recording equipment included analog tape deck, digital minidisk recorder, and direct digital recording to a laptop computer disk. Microphones and geophones were used as transducers. The natural amphitheater lies between the ruins of Pueblo Bonito and Chetro Ketl. It is a semicircular arc in a sandstone cliff measuring 500 ft. wide and 75 ft. high. The radius of the arc was verified with aerial photography, and an acoustic ray trace was generated using cad software. The arc is in an overhanging cliff face and brings distant sounds to a line focus. Along this line, there are unusual acoustic effects at conjugate foci. Time history analysis of recordings from both sites showed that a 60-dB reverb decay lasted from 1.8 to 2.0 s, nearly ideal for public performances of music. Echoes from the amphitheater were perceived to be upshifted in pitch, but this was not seen in FFT analysis. Geophones placed on the floor of the great kiva showed a resonance at 95 Hz.

Low frequency acoustic and electromagnetic scattering

NASA Technical Reports Server (NTRS)

Hariharan, S. I.; Maccamy, R. C.

1986-01-01

This paper deals with two classes of problems arising from acoustics and electromagnetics scattering in the low frequency stations. The first class of problem is solving Helmholtz equation with Dirichlet boundary conditions on an arbitrary two dimensional body while the second one is an interior-exterior interface problem with Helmholtz equation in the exterior. Low frequency analysis show that there are two intermediate problems which solve the above problems accurate to 0(k/2/ log k) where k is the frequency. These solutions greatly differ from the zero frequency approximations. For the Dirichlet problem numerical examples are shown to verify the theoretical estimates.

Diffraction of a plane wave on two-dimensional conductive structures and a surface wave

NASA Astrophysics Data System (ADS)

Davidovich, Mikhael V.

2018-04-01

We consider the structures type of two-dimensional electron gas in the form of a thin conductive, in particular, graphene films described by tensor conductivity, which are isolated or located on the dielectric layers. The dispersion equation for hybrid modes, as well as scattering parameters. We show that free wave (eigenwaves) problem follow from the problem of diffraction when linking the amplitude of the current of the linear equations are unsolvable, i.e., the determinant of this system is zero. As a particular case the dispersion equation follow from the conditions of matching (with zero reflection coefficient).

A Three-Dimensional Finite-Element Model for Simulating Water Flow in Variably Saturated Porous Media

NASA Astrophysics Data System (ADS)

Huyakorn, Peter S.; Springer, Everett P.; Guvanasen, Varut; Wadsworth, Terry D.

1986-12-01

A three-dimensional finite-element model for simulating water flow in variably saturated porous media is presented. The model formulation is general and capable of accommodating complex boundary conditions associated with seepage faces and infiltration or evaporation on the soil surface. Included in this formulation is an improved Picard algorithm designed to cope with severely nonlinear soil moisture relations. The algorithm is formulated for both rectangular and triangular prism elements. The element matrices are evaluated using an "influence coefficient" technique that avoids costly numerical integration. Spatial discretization of a three-dimensional region is performed using a vertical slicing approach designed to accommodate complex geometry with irregular boundaries, layering, and/or lateral discontinuities. Matrix solution is achieved using a slice successive overrelaxation scheme that permits a fairly large number of nodal unknowns (on the order of several thousand) to be handled efficiently on small minicomputers. Six examples are presented to verify and demonstrate the utility of the proposed finite-element model. The first four examples concern one- and two-dimensional flow problems used as sample problems to benchmark the code. The remaining examples concern three-dimensional problems. These problems are used to illustrate the performance of the proposed algorithm in three-dimensional situations involving seepage faces and anisotropic soil media.

Linearized compressible-flow theory for sonic flight speeds

NASA Technical Reports Server (NTRS)

Heaslet, Max A; Lomax, Harvard; Spreiter, John R

1950-01-01

The partial differential equation for the perturbation velocity potential is examined for free-stream Mach numbers close to and equal to one. It is found that, under the assumptions of linearized theory, solutions can be found consistent with the theory for lifting-surface problems both in stationary three-dimensional flow and in unsteady two-dimensional flow. Several examples are solved including a three dimensional swept-back wing and two dimensional harmonically-oscillating wing, both for a free stream Mach number equal to one. Momentum relations for the evaluation of wave and vortex drag are also discussed. (author)

User's manual for two dimensional FDTD version TEA and TMA codes for scattering from frequency-independent dielectric materials

NASA Technical Reports Server (NTRS)

Beggs, John H.; Luebbers, Raymond J.; Kunz, Karl S.

1991-01-01

The Penn State Finite Difference Time Domain Electromagnetic Scattering Code Versions TEA and TMA are two dimensional electromagnetic scattering codes based on the Finite Difference Time Domain Technique (FDTD) first proposed by Yee in 1966. The supplied version of the codes are two versions of our current FDTD code set. This manual provides a description of the codes and corresponding results for the default scattering problem. The manual is organized into eleven sections: introduction, Version TEA and TMA code capabilities, a brief description of the default scattering geometry, a brief description of each subroutine, a description of the include files (TEACOM.FOR TMACOM.FOR), a section briefly discussing scattering width computations, a section discussing the scattering results, a sample problem setup section, a new problem checklist, references, and figure titles.

Difference equation state approximations for nonlinear hereditary control problems

NASA Technical Reports Server (NTRS)

Rosen, I. G.

1984-01-01

Discrete approximation schemes for the solution of nonlinear hereditary control problems are constructed. The methods involve approximation by a sequence of optimal control problems in which the original infinite dimensional state equation has been approximated by a finite dimensional discrete difference equation. Convergence of the state approximations is argued using linear semigroup theory and is then used to demonstrate that solutions to the approximating optimal control problems in some sense approximate solutions to the original control problem. Two schemes, one based upon piecewise constant approximation, and the other involving spline functions are discussed. Numerical results are presented, analyzed and used to compare the schemes to other available approximation methods for the solution of hereditary control problems. Previously announced in STAR as N83-33589

Creation of problem-dependent Doppler-broadened cross sections in the KENO Monte Carlo code

DOE Office of Scientific and Technical Information (OSTI.GOV)

Hart, Shane W. D.; Celik, Cihangir; Maldonado, G. Ivan

2015-11-06

In this paper, we introduce a quick method for improving the accuracy of Monte Carlo simulations by generating one- and two-dimensional cross sections at a user-defined temperature before performing transport calculations. A finite difference method is used to Doppler-broaden cross sections to the desired temperature, and unit-base interpolation is done to generate the probability distributions for double differential two-dimensional thermal moderator cross sections at any arbitrarily user-defined temperature. The accuracy of these methods is tested using a variety of contrived problems. In addition, various benchmarks at elevated temperatures are modeled, and results are compared with benchmark results. Lastly, the problem-dependentmore » cross sections are observed to produce eigenvalue estimates that are closer to the benchmark results than those without the problem-dependent cross sections.« less

Convergence of an hp-Adaptive Finite Element Strategy in Two and Three Space-Dimensions

NASA Astrophysics Data System (ADS)

Bürg, Markus; Dörfler, Willy

2010-09-01

We show convergence of an automatic hp-adaptive refinement strategy for the finite element method on the elliptic boundary value problem. The strategy is a generalization of a refinement strategy proposed for one-dimensional situations to problems in two and three space-dimensions.

The P1-RKDG method for two-dimensional Euler equations of gas dynamics

NASA Technical Reports Server (NTRS)

Cockburn, Bernardo; Shu, Chi-Wang

1991-01-01

A class of nonlinearly stable Runge-Kutta local projection discontinuous Galerkin (RKDG) finite element methods for conservation laws is investigated. Two dimensional Euler equations for gas dynamics are solved using P1 elements. The generalization of the local projections, which for scalar nonlinear conservation laws was designed to satisfy a local maximum principle, to systems of conservation laws such as the Euler equations of gas dynamics using local characteristic decompositions is discussed. Numerical examples include the standard regular shock reflection problem, the forward facing step problem, and the double Mach reflection problem. These preliminary numerical examples are chosen to show the capacity of the approach to obtain nonlinearly stable results comparable with the modern nonoscillatory finite difference methods.

A collection of edge-based elements

NASA Technical Reports Server (NTRS)

Kempel, Leo C.; Volakis, John L.

1992-01-01

Edge-based elements have proved useful in solving electromagnetic problems since they are nondivergent. Previous authors have presented several two and three dimensional elements. Herein, we present four types of elements which are suitable for modeling several types of three dimensional geometries. Distorted brick and triangular prism elements are given in cartesian coordinates as well as the specialized cylindrical shell and pie-shaped prism elements which are suitable for problems best described in polar cylindrical coordinates.

Travelling-wave solutions of a weakly nonlinear two-dimensional higher-order Kadomtsev-Petviashvili dynamical equation for dispersive shallow-water waves

NASA Astrophysics Data System (ADS)

Seadawy, Aly R.

2017-01-01

The propagation of three-dimensional nonlinear irrotational flow of an inviscid and incompressible fluid of the long waves in dispersive shallow-water approximation is analyzed. The problem formulation of the long waves in dispersive shallow-water approximation lead to fifth-order Kadomtsev-Petviashvili (KP) dynamical equation by applying the reductive perturbation theory. By using an extended auxiliary equation method, the solitary travelling-wave solutions of the two-dimensional nonlinear fifth-order KP dynamical equation are derived. An analytical as well as a numerical solution of the two-dimensional nonlinear KP equation are obtained and analyzed with the effects of external pressure flow.

Optimal Control for Stochastic Delay Evolution Equations

DOE Office of Scientific and Technical Information (OSTI.GOV)

Meng, Qingxin, E-mail: mqx@hutc.zj.cn; Shen, Yang, E-mail: skyshen87@gmail.com

2016-08-15

In this paper, we investigate a class of infinite-dimensional optimal control problems, where the state equation is given by a stochastic delay evolution equation with random coefficients, and the corresponding adjoint equation is given by an anticipated backward stochastic evolution equation. We first prove the continuous dependence theorems for stochastic delay evolution equations and anticipated backward stochastic evolution equations, and show the existence and uniqueness of solutions to anticipated backward stochastic evolution equations. Then we establish necessary and sufficient conditions for optimality of the control problem in the form of Pontryagin’s maximum principles. To illustrate the theoretical results, we applymore » stochastic maximum principles to study two examples, an infinite-dimensional linear-quadratic control problem with delay and an optimal control of a Dirichlet problem for a stochastic partial differential equation with delay. Further applications of the two examples to a Cauchy problem for a controlled linear stochastic partial differential equation and an optimal harvesting problem are also considered.« less

Whitham modulation theory for (2  +  1)-dimensional equations of Kadomtsev–Petviashvili type

NASA Astrophysics Data System (ADS)

Ablowitz, Mark J.; Biondini, Gino; Rumanov, Igor

2018-05-01

Whitham modulation theory for certain two-dimensional evolution equations of Kadomtsev–Petviashvili (KP) type is presented. Three specific examples are considered in detail: the KP equation, the two-dimensional Benjamin–Ono (2DBO) equation and a modified KP (m2KP) equation. A unified derivation is also provided. In the case of the m2KP equation, the corresponding Whitham modulation system exhibits features different from the other two. The approach presented here does not require integrability of the original evolution equation. Indeed, while the KP equation is known to be a completely integrable equation, the 2DBO equation and the m2KP equation are not known to be integrable. In each of the cases considered, the Whitham modulation system obtained consists of five first-order quasilinear partial differential equations. The Riemann problem (i.e. the analogue of the Gurevich–Pitaevskii problem) for the one-dimensional reduction of the m2KP equation is studied. For the m2KP equation, the system of modulation equations is used to analyze the linear stability of traveling wave solutions.

Lax-Wendroff and TVD finite volume methods for unidimensional thermomechanical numerical simulations of impacts on elastic-plastic solids

NASA Astrophysics Data System (ADS)

Heuzé, Thomas

2017-10-01

We present in this work two finite volume methods for the simulation of unidimensional impact problems, both for bars and plane waves, on elastic-plastic solid media within the small strain framework. First, an extension of Lax-Wendroff to elastic-plastic constitutive models with linear and nonlinear hardenings is presented. Second, a high order TVD method based on flux-difference splitting [1] and Superbee flux limiter [2] is coupled with an approximate elastic-plastic Riemann solver for nonlinear hardenings, and follows that of Fogarty [3] for linear ones. Thermomechanical coupling is accounted for through dissipation heating and thermal softening, and adiabatic conditions are assumed. This paper essentially focuses on one-dimensional problems since analytical solutions exist or can easily be developed. Accordingly, these two numerical methods are compared to analytical solutions and to the explicit finite element method on test cases involving discontinuous and continuous solutions. This allows to study in more details their respective performance during the loading, unloading and reloading stages. Particular emphasis is also paid to the accuracy of the computed plastic strains, some differences being found according to the numerical method used. Lax-Wendoff two-dimensional discretization of a one-dimensional problem is also appended at the end to demonstrate the extensibility of such numerical scheme to multidimensional problems.

Problems of Conducting Research in Organizations: The Case of Police Departments.

ERIC Educational Resources Information Center

Lefkowitz, Joel

This paper presents a description of police research problems in such fashion that it could be generalized to other types of organizations. A two-dimensional taxonomy of problems in conducting psychological research in police departments is discussed. The first dimension concerns generality-uniqueness of the problem, relative to formal…

Active Subspaces of Airfoil Shape Parameterizations

NASA Astrophysics Data System (ADS)

Grey, Zachary J.; Constantine, Paul G.

2018-05-01

Design and optimization benefit from understanding the dependence of a quantity of interest (e.g., a design objective or constraint function) on the design variables. A low-dimensional active subspace, when present, identifies important directions in the space of design variables; perturbing a design along the active subspace associated with a particular quantity of interest changes that quantity more, on average, than perturbing the design orthogonally to the active subspace. This low-dimensional structure provides insights that characterize the dependence of quantities of interest on design variables. Airfoil design in a transonic flow field with a parameterized geometry is a popular test problem for design methodologies. We examine two particular airfoil shape parameterizations, PARSEC and CST, and study the active subspaces present in two common design quantities of interest, transonic lift and drag coefficients, under each shape parameterization. We mathematically relate the two parameterizations with a common polynomial series. The active subspaces enable low-dimensional approximations of lift and drag that relate to physical airfoil properties. In particular, we obtain and interpret a two-dimensional approximation of both transonic lift and drag, and we show how these approximation inform a multi-objective design problem.

Investigation of the relative orientation of the system of optical sensors to monitor the technosphere objects

NASA Astrophysics Data System (ADS)

Petrochenko, Andrey; Konyakhin, Igor

2017-06-01

In connection with the development of robotics have become increasingly popular variety of three-dimensional reconstruction of the system mapping and image-set received from the optical sensors. The main objective of technical and robot vision is the detection, tracking and classification of objects of the space in which these systems and robots operate [15,16,18]. Two-dimensional images sometimes don't contain sufficient information to address those or other problems: the construction of the map of the surrounding area for a route; object identification, tracking their relative position and movement; selection of objects and their attributes to complement the knowledge base. Three-dimensional reconstruction of the surrounding space allows you to obtain information on the relative positions of objects, their shape, surface texture. Systems, providing training on the basis of three-dimensional reconstruction of the results of the comparison can produce two-dimensional images of three-dimensional model that allows for the recognition of volume objects on flat images. The problem of the relative orientation of industrial robots with the ability to build threedimensional scenes of controlled surfaces is becoming actual nowadays.

Some problems of the calculation of three-dimensional boundary layer flows on general configurations

NASA Technical Reports Server (NTRS)

Cebeci, T.; Kaups, K.; Mosinskis, G. J.; Rehn, J. A.

1973-01-01

An accurate solution of the three-dimensional boundary layer equations over general configurations such as those encountered in aircraft and space shuttle design requires a very efficient, fast, and accurate numerical method with suitable turbulence models for the Reynolds stresses. The efficiency, speed, and accuracy of a three-dimensional numerical method together with the turbulence models for the Reynolds stresses are examined. The numerical method is the implicit two-point finite difference approach (Box Method) developed by Keller and applied to the boundary layer equations by Keller and Cebeci. In addition, a study of some of the problems that may arise in the solution of these equations for three-dimensional boundary layer flows over general configurations.

A new Lagrangian method for three-dimensional steady supersonic flows

NASA Technical Reports Server (NTRS)

Loh, Ching-Yuen; Liou, Meng-Sing

1993-01-01

In this report, the new Lagrangian method introduced by Loh and Hui is extended for three-dimensional, steady supersonic flow computation. The derivation of the conservation form and the solution of the local Riemann solver using the Godunov and the high-resolution TVD (total variation diminished) scheme is presented. This new approach is accurate and robust, capable of handling complicated geometry and interactions between discontinuous waves. Test problems show that the extended Lagrangian method retains all the advantages of the two-dimensional method (e.g., crisp resolution of a slip-surface (contact discontinuity) and automatic grid generation). In this report, we also suggest a novel three dimensional Riemann problem in which interesting and intricate flow features are present.

The development of laser speckle velocimetry for the study of vortical flows

NASA Technical Reports Server (NTRS)

Krothapalli, A.

1991-01-01

A new experimental technique commonly known as PIDV (particle image displacement velocity) was developed to measure an instantaneous two dimensional velocity fluid in a selected plane of the flow field. This technique was successfully applied to the study of several problems: (1) unsteady flows with large scale vortical structures; (2) the instantaneous two dimensional flow in the transition region of a rectangular air jet; and (3) the instantaneous flow over a circular bump in a transonic flow. In several other experiments PIDV is routinely used as a non-intrusive measurement technique to obtain instantaneous two dimensional velocity fields.

Optimal one-dimensional inversion and bounding of magnetotelluric apparent resistivity and phase measurements

NASA Astrophysics Data System (ADS)

Parker, Robert L.; Booker, John R.

1996-12-01

The properties of the log of the admittance in the complex frequency plane lead to an integral representation for one-dimensional magnetotelluric (MT) apparent resistivity and impedance phase similar to that found previously for complex admittance. The inverse problem of finding a one-dimensional model for MT data can then be solved using the same techniques as for complex admittance, with similar results. For instance, the one-dimensional conductivity model that minimizes the χ2 misfit statistic for noisy apparent resistivity and phase is a series of delta functions. One of the most important applications of the delta function solution to the inverse problem for complex admittance has been answering the question of whether or not a given set of measurements is consistent with the modeling assumption of one-dimensionality. The new solution allows this test to be performed directly on standard MT data. Recently, it has been shown that induction data must pass the same one-dimensional consistency test if they correspond to the polarization in which the electric field is perpendicular to the strike of two-dimensional structure. This greatly magnifies the utility of the consistency test. The new solution also allows one to compute the upper and lower bounds permitted on phase or apparent resistivity at any frequency given a collection of MT data. Applications include testing the mutual consistency of apparent resistivity and phase data and placing bounds on missing phase or resistivity data. Examples presented demonstrate detection and correction of equipment and processing problems and verification of compatibility with two-dimensional B-polarization for MT data after impedance tensor decomposition and for continuous electromagnetic profiling data.

Possible sources of archaeological maize found in Chaco Canyon and Aztec Ruin, New Mexico

USGS Publications Warehouse

Benson, L.V.; Stein, J.R.; Taylor, Howard E.

2009-01-01

Maize played a major role in Chaco's interaction with outlying communities in the southern Colorado Plateau. This paper seeks to determine where archaeological corn cobs brought to Chaco Canyon were grown. Strontium-isotope and trace-metal ratios of 180 soil-water and 18 surface-water sites in the Southern Colorado Plateau have revealed possible source areas for some of 37 archaeological corn cobs from Chaco Canyon and 10 archaeological corn cobs from Aztec Ruin, New Mexico. The most probable source areas for cobs that predate the middle-12th-century drought include several Upper Rio Chaco sites (not including Chaco Canyon). There are many potential source areas for cobs that date to the late A.D. 1100s and early 1200s, all of which lie in the eastern part of the study area. Some Athapascan-age cobs have potential source areas in the Totah, Lobo Mesa, and Dinetah regions. One Gallo Cliff Dwelling cob has a strontium-isotope ratio that exceeds all measured soil-water values. Field sites for this cob may exist in association with Paleozoic and Precambrian rocks found 80-90 km from Chaco Canyon. Potential source areas for most Aztec Ruin cobs (many of which were found in rooms dating to the first half of the 13th-century) appear to be associated with a loess deposit that blankets the Mesa Verde and McElmo Dome regions.

Advanced numerical methods for three dimensional two-phase flow calculations

DOE Office of Scientific and Technical Information (OSTI.GOV)

Toumi, I.; Caruge, D.

1997-07-01

This paper is devoted to new numerical methods developed for both one and three dimensional two-phase flow calculations. These methods are finite volume numerical methods and are based on the use of Approximate Riemann Solvers concepts to define convective fluxes versus mean cell quantities. The first part of the paper presents the numerical method for a one dimensional hyperbolic two-fluid model including differential terms as added mass and interface pressure. This numerical solution scheme makes use of the Riemann problem solution to define backward and forward differencing to approximate spatial derivatives. The construction of this approximate Riemann solver uses anmore » extension of Roe`s method that has been successfully used to solve gas dynamic equations. As far as the two-fluid model is hyperbolic, this numerical method seems very efficient for the numerical solution of two-phase flow problems. The scheme was applied both to shock tube problems and to standard tests for two-fluid computer codes. The second part describes the numerical method in the three dimensional case. The authors discuss also some improvements performed to obtain a fully implicit solution method that provides fast running steady state calculations. Such a scheme is not implemented in a thermal-hydraulic computer code devoted to 3-D steady-state and transient computations. Some results obtained for Pressurised Water Reactors concerning upper plenum calculations and a steady state flow in the core with rod bow effect evaluation are presented. In practice these new numerical methods have proved to be stable on non staggered grids and capable of generating accurate non oscillating solutions for two-phase flow calculations.« less

The simulation of a two-dimensional (2D) transport problem in a rectangular region with Lattice Boltzmann method with two-relaxation-time

NASA Astrophysics Data System (ADS)

Sugiyanto, S.; Hardyanto, W.; Marwoto, P.

2018-03-01

Transport phenomena are found in many problems in many engineering and industrial sectors. We analyzed a Lattice Boltzmann method with Two-Relaxation Time (LTRT) collision operators for simulation of pollutant moving through the medium as a two-dimensional (2D) transport problem in a rectangular region model. This model consists of a 2D rectangular region with 54 length (x), 27 width (y), and it has isotropic homogeneous medium. Initially, the concentration is zero and is distributed evenly throughout the region of interest. A concentration of 1 is maintained at 9 < y < 18, whereas the concentration of zero is maintained at 0 < y < 9 and 18 < y < 27. A specific discharge (Darcy velocity) of 1.006 is assumed. A diffusion coefficient of 0.8333 is distributed uniformly with a uniform porosity of 0.35. A computer program is written in MATLAB to compute the concentration of pollutant at any specified place and time. The program shows that LTRT solution with quadratic equilibrium distribution functions (EDFs) and relaxation time τa=1.0 are in good agreement result with other numerical solutions methods such as 3DLEWASTE (Hybrid Three-dimensional Lagrangian-Eulerian Finite Element Model of Waste Transport Through Saturated-Unsaturated Media) obtained by Yeh and 3DFEMWATER-LHS (Three-dimensional Finite Element Model of Water Flow Through Saturated-Unsaturated Media with Latin Hypercube Sampling) obtained by Hardyanto.

Two-Dimensional Failure Waves and Ignition Fronts in Premixed Combustion

NASA Technical Reports Server (NTRS)

Vedarajan, T. G.; Buckmaster J.; Ronney, P.

1998-01-01

This paper is a continuation of our work on edge-flames in premixed combustion. An edge-flame is a two-dimensional structure constructed from a one-dimensional configuration that has two stable solutions (bistable equilibrium). Edge-flames can display wavelike behavior, advancing as ignition fronts or retreating as failure waves. Here we consider two one-dimensional configurations: twin deflagrations in a straining flow generated by the counterflow of fresh streams of mixture: and a single deflagration subject to radiation losses. The edge-flames constructed from the first configuration have positive or negative speeds, according to the value of the strain rate. But our numerical solutions strongly suggest that only positive speeds (corresponding to ignition fronts) can exist for the second configuration. We show that this phenomenon can also occur in diffusion flames when the Lewis numbers are small. And we discuss the asymptotics of the one-dimensional twin deflagration configuration. an overlooked problem from the 70s.

Solution of axisymmetric and two-dimensional inviscid flow over blunt bodies by the method of lines

NASA Technical Reports Server (NTRS)

Hamilton, H. H., II

1978-01-01

Comparisons with experimental data and the results of other computational methods demonstrated that very accurate solutions can be obtained by using relatively few lines with the method of lines approach. This method is semidiscrete and has relatively low core storage requirements as compared with fully discrete methods since very little data were stored across the shock layer. This feature is very attractive for three dimensional problems because it enables computer storage requirements to be reduced by approximately an order of magnitude. In the present study it was found that nine lines was a practical upper limit for two dimensional and axisymmetric problems. This condition limits application of the method to smooth body geometries where relatively few lines would be adequate to describe changes in the flow variables around the body. Extension of the method to three dimensions was conceptually straightforward; however, three dimensional applications would also be limited to smooth body geometries although not necessarily to total of nine lines.

Multitasking a three-dimensional Navier-Stokes algorithm on the Cray-2

NASA Technical Reports Server (NTRS)

Swisshelm, Julie M.

1989-01-01

A three-dimensional computational aerodynamics algorithm has been multitasked for efficient parallel execution on the Cray-2. It provides a means for examining the multitasking performance of a complete CFD application code. An embedded zonal multigrid scheme is used to solve the Reynolds-averaged Navier-Stokes equations for an internal flow model problem. The explicit nature of each component of the method allows a spatial partitioning of the computational domain to achieve a well-balanced task load for MIMD computers with vector-processing capability. Experiments have been conducted with both two- and three-dimensional multitasked cases. The best speedup attained by an individual task group was 3.54 on four processors of the Cray-2, while the entire solver yielded a speedup of 2.67 on four processors for the three-dimensional case. The multiprocessing efficiency of various types of computational tasks is examined, performance on two Cray-2s with different memory access speeds is compared, and extrapolation to larger problems is discussed.

Multi-GPU hybrid programming accelerated three-dimensional phase-field model in binary alloy

NASA Astrophysics Data System (ADS)

Zhu, Changsheng; Liu, Jieqiong; Zhu, Mingfang; Feng, Li

2018-03-01

In the process of dendritic growth simulation, the computational efficiency and the problem scales have extremely important influence on simulation efficiency of three-dimensional phase-field model. Thus, seeking for high performance calculation method to improve the computational efficiency and to expand the problem scales has a great significance to the research of microstructure of the material. A high performance calculation method based on MPI+CUDA hybrid programming model is introduced. Multi-GPU is used to implement quantitative numerical simulations of three-dimensional phase-field model in binary alloy under the condition of multi-physical processes coupling. The acceleration effect of different GPU nodes on different calculation scales is explored. On the foundation of multi-GPU calculation model that has been introduced, two optimization schemes, Non-blocking communication optimization and overlap of MPI and GPU computing optimization, are proposed. The results of two optimization schemes and basic multi-GPU model are compared. The calculation results show that the use of multi-GPU calculation model can improve the computational efficiency of three-dimensional phase-field obviously, which is 13 times to single GPU, and the problem scales have been expanded to 8193. The feasibility of two optimization schemes is shown, and the overlap of MPI and GPU computing optimization has better performance, which is 1.7 times to basic multi-GPU model, when 21 GPUs are used.

Algebraic multigrid methods applied to problems in computational structural mechanics

NASA Technical Reports Server (NTRS)

Mccormick, Steve; Ruge, John

1989-01-01

The development of algebraic multigrid (AMG) methods and their application to certain problems in structural mechanics are described with emphasis on two- and three-dimensional linear elasticity equations and the 'jacket problems' (three-dimensional beam structures). Various possible extensions of AMG are also described. The basic idea of AMG is to develop the discretization sequence based on the target matrix and not the differential equation. Therefore, the matrix is analyzed for certain dependencies that permit the proper construction of coarser matrices and attendant transfer operators. In this manner, AMG appears to be adaptable to structural analysis applications.

Simulation and Analysis of Converging Shock Wave Test Problems

DOE Office of Scientific and Technical Information (OSTI.GOV)

Ramsey, Scott D.; Shashkov, Mikhail J.

2012-06-21

Results and analysis pertaining to the simulation of the Guderley converging shock wave test problem (and associated code verification hydrodynamics test problems involving converging shock waves) in the LANL ASC radiation-hydrodynamics code xRAGE are presented. One-dimensional (1D) spherical and two-dimensional (2D) axi-symmetric geometric setups are utilized and evaluated in this study, as is an instantiation of the xRAGE adaptive mesh refinement capability. For the 2D simulations, a 'Surrogate Guderley' test problem is developed and used to obviate subtleties inherent to the true Guderley solution's initialization on a square grid, while still maintaining a high degree of fidelity to the originalmore » problem, and minimally straining the general credibility of associated analysis and conclusions.« less

Pattern-set generation algorithm for the one-dimensional multiple stock sizes cutting stock problem

NASA Astrophysics Data System (ADS)

Cui, Yaodong; Cui, Yi-Ping; Zhao, Zhigang

2015-09-01

A pattern-set generation algorithm (PSG) for the one-dimensional multiple stock sizes cutting stock problem (1DMSSCSP) is presented. The solution process contains two stages. In the first stage, the PSG solves the residual problems repeatedly to generate the patterns in the pattern set, where each residual problem is solved by the column-generation approach, and each pattern is generated by solving a single large object placement problem. In the second stage, the integer linear programming model of the 1DMSSCSP is solved using a commercial solver, where only the patterns in the pattern set are considered. The computational results of benchmark instances indicate that the PSG outperforms existing heuristic algorithms and rivals the exact algorithm in solution quality.

Fast Optimization for Aircraft Descent and Approach Trajectory

NASA Technical Reports Server (NTRS)

Luchinsky, Dmitry G.; Schuet, Stefan; Brenton, J.; Timucin, Dogan; Smith, David; Kaneshige, John

2017-01-01

We address problem of on-line scheduling of the aircraft descent and approach trajectory. We formulate a general multiphase optimal control problem for optimization of the descent trajectory and review available methods of its solution. We develop a fast algorithm for solution of this problem using two key components: (i) fast inference of the dynamical and control variables of the descending trajectory from the low dimensional flight profile data and (ii) efficient local search for the resulting reduced dimensionality non-linear optimization problem. We compare the performance of the proposed algorithm with numerical solution obtained using optimal control toolbox General Pseudospectral Optimal Control Software. We present results of the solution of the scheduling problem for aircraft descent using novel fast algorithm and discuss its future applications.

Plane Poiseuille Flow of a Rarefied Gas in the Presence of a Strong Gravitation

NASA Astrophysics Data System (ADS)

Doi, Toshiyuki

2010-11-01

Poiseuille flow of a rarefied gas between two horizontal planes in the presence of a strong gravitation is considered, where the gravity is so strong that the path of a molecule is curved considerably as it ascends or descends the distance of the planes. The gas behavior is studied based on the Boltzmann equation. An asymptotic analysis for a slow variation in the longitudinal direction is carried out and the problem is reduced to a spatially one dimensional problem, as was in the Poiseuille flow problem in the absence of the gravitation. The mass flow rate as well as the macroscopic variables is obtained for a wide range of the mean free path of the gas and the gravity. A numerical analysis of a two dimensional problem is also carried out and the result of the asymptotic analysis is verified.

A Conserving Discretization for the Free Boundary in a Two-Dimensional Stefan Problem

NASA Astrophysics Data System (ADS)

Segal, Guus; Vuik, Kees; Vermolen, Fred

1998-03-01

The dissolution of a disk-likeAl2Cuparticle is considered. A characteristic property is that initially the particle has a nonsmooth boundary. The mathematical model of this dissolution process contains a description of the particle interface, of which the position varies in time. Such a model is called a Stefan problem. It is impossible to obtain an analytical solution for a general two-dimensional Stefan problem, so we use the finite element method to solve this problem numerically. First, we apply a classical moving mesh method. Computations show that after some time steps the predicted particle interface becomes very unrealistic. Therefore, we derive a new method for the displacement of the free boundary based on the balance of atoms. This method leads to good results, also, for nonsmooth boundaries. Some numerical experiments are given for the dissolution of anAl2Cuparticle in anAl-Cualloy.

Spatial Visualization in Physics Problem Solving

ERIC Educational Resources Information Center

Kozhevnikov, Maria; Motes, Michael A.; Hegarty, Mary

2007-01-01

Three studies were conducted to examine the relation of spatial visualization to solving kinematics problems that involved either predicting the two-dimensional motion of an object, translating from one frame of reference to another, or interpreting kinematics graphs. In Study 1, 60 physics-naive students were administered kinematics problems and…

A three-dimensional Dirichlet-to-Neumann operator for water waves over topography

NASA Astrophysics Data System (ADS)

Andrade, D.; Nachbin, A.

2018-06-01

Surface water waves are considered propagating over highly variable non-smooth topographies. For this three dimensional problem a Dirichlet-to-Neumann (DtN) operator is constructed reducing the numerical modeling and evolution to the two dimensional free surface. The corresponding Fourier-type operator is defined through a matrix decomposition. The topographic component of the decomposition requires special care and a Galerkin method is provided accordingly. One dimensional numerical simulations, along the free surface, validate the DtN formulation in the presence of a large amplitude, rapidly varying topography. An alternative, conformal mapping based, method is used for benchmarking. A two dimensional simulation in the presence of a Luneburg lens (a particular submerged mound) illustrates the accurate performance of the three dimensional DtN operator.

A Multi-Resolution Data Structure for Two-Dimensional Morse Functions

DOE Office of Scientific and Technical Information (OSTI.GOV)

Bremer, P-T; Edelsbrunner, H; Hamann, B

2003-07-30

The efficient construction of simplified models is a central problem in the field of visualization. We combine topological and geometric methods to construct a multi-resolution data structure for functions over two-dimensional domains. Starting with the Morse-Smale complex we build a hierarchy by progressively canceling critical points in pairs. The data structure supports mesh traversal operations similar to traditional multi-resolution representations.

Terahertz signal detection in a short gate length field-effect transistor with a two-dimensional electron gas

DOE Office of Scientific and Technical Information (OSTI.GOV)

Vostokov, N. V., E-mail: vostokov@ipm.sci-nnov.ru; Shashkin, V. I.

2015-11-28

We consider the problem of non-resonant detection of terahertz signals in a short gate length field-effect transistor having a two-dimensional electron channel with zero external bias between the source and the drain. The channel resistance, gate-channel capacitance, and quadratic nonlinearity parameter of the transistor during detection as a function of the gate bias voltage are studied. Characteristics of detection of the transistor connected in an antenna with real impedance are analyzed. The consideration is based on both a simple one-dimensional model of the transistor and allowance for the two-dimensional distribution of the electric field in the transistor structure. The resultsmore » given by the different models are discussed.« less

Unsteady transonic flow calculations for two-dimensional canard-wing configurations with aeroelastic applications

NASA Technical Reports Server (NTRS)

Batina, J. T.

1985-01-01

Unsteady transonic flow calculations for aerodynamically interfering airfoil configurations are performed as a first step toward solving the three dimensional canard wing interaction problem. These calculations are performed by extending the XTRAN2L two dimensional unsteady transonic small disturbance code to include an additional airfoil. Unsteady transonic forces due to plunge and pitch motions of a two dimensional canard and wing are presented. Results for a variety of canard wing separation distances reveal the effects of aerodynamic interference on unsteady transonic airloads. Aeroelastic analyses employing these unsteady airloads demonstrate the effects of aerodynamic interference on aeroelastic stability and flutter. For the configurations studied, increases in wing flutter speed result with the inclusion of the aerodynamically interfering canard.

A generalized rotationally symmetric case of the centroaffine Minkowski problem

NASA Astrophysics Data System (ADS)

Lu, Jian

2018-05-01

In this paper the centroaffine Minkowski problem, a critical case of the Lp-Minkowski problem in the n + 1 dimensional Euclidean space, is studied. By its variational structure and the method of blow-up analyses, we obtain two sufficient conditions for the existence of solutions, for a generalized rotationally symmetric case of the problem.

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.

Two-dimensional wavefront reconstruction based on double-shearing and least squares fitting

NASA Astrophysics Data System (ADS)

Liang, Peiying; Ding, Jianping; Zhu, Yangqing; Dong, Qian; Huang, Yuhua; Zhu, Zhen

2017-06-01

The two-dimensional wavefront reconstruction method based on double-shearing and least squares fitting is proposed in this paper. Four one-dimensional phase estimates of the measured wavefront, which correspond to the two shears and the two orthogonal directions, could be calculated from the differential phase, which solves the problem of the missing spectrum, and then by using the least squares method the two-dimensional wavefront reconstruction could be done. The numerical simulations of the proposed algorithm are carried out to verify the feasibility of this method. The influence of noise generated from different shear amount and different intensity on the accuracy of the reconstruction is studied and compared with the results from the algorithm based on single-shearing and least squares fitting. Finally, a two-grating lateral shearing interference experiment is carried out to verify the wavefront reconstruction algorithm based on doubleshearing and least squares fitting.

A numerical study of blood flow using mixture theory

PubMed Central

Wu, Wei-Tao; Aubry, Nadine; Massoudi, Mehrdad; Kim, Jeongho; Antaki, James F.

2014-01-01

In this paper, we consider the two dimensional flow of blood in a rectangular microfluidic channel. We use Mixture Theory to treat this problem as a two-component system: One component is the red blood cells (RBCs) modeled as a generalized Reiner–Rivlin type fluid, which considers the effects of volume fraction (hematocrit) and influence of shear rate upon viscosity. The other component, plasma, is assumed to behave as a linear viscous fluid. A CFD solver based on OpenFOAM® was developed and employed to simulate a specific problem, namely blood flow in a two dimensional micro-channel, is studied. Finally to better understand this two-component flow system and the effects of the different parameters, the equations are made dimensionless and a parametric study is performed. PMID:24791016

A numerical study of blood flow using mixture theory.

PubMed

Wu, Wei-Tao; Aubry, Nadine; Massoudi, Mehrdad; Kim, Jeongho; Antaki, James F

2014-03-01

In this paper, we consider the two dimensional flow of blood in a rectangular microfluidic channel. We use Mixture Theory to treat this problem as a two-component system: One component is the red blood cells (RBCs) modeled as a generalized Reiner-Rivlin type fluid, which considers the effects of volume fraction (hematocrit) and influence of shear rate upon viscosity. The other component, plasma, is assumed to behave as a linear viscous fluid. A CFD solver based on OpenFOAM ® was developed and employed to simulate a specific problem, namely blood flow in a two dimensional micro-channel, is studied. Finally to better understand this two-component flow system and the effects of the different parameters, the equations are made dimensionless and a parametric study is performed.

General design method for 3-dimensional, potential flow fields. Part 2: Computer program DIN3D1 for simple, unbranched ducts

NASA Technical Reports Server (NTRS)

Stanitz, J. D.

1985-01-01

The general design method for three-dimensional, potential, incompressible or subsonic-compressible flow developed in part 1 of this report is applied to the design of simple, unbranched ducts. A computer program, DIN3D1, is developed and five numerical examples are presented: a nozzle, two elbows, an S-duct, and the preliminary design of a side inlet for turbomachines. The two major inputs to the program are the upstream boundary shape and the lateral velocity distribution on the duct wall. As a result of these inputs, boundary conditions are overprescribed and the problem is ill posed. However, it appears that there are degrees of compatibility between these two major inputs and that, for reasonably compatible inputs, satisfactory solutions can be obtained. By not prescribing the shape of the upstream boundary, the problem presumably becomes well posed, but it is not clear how to formulate a practical design method under this circumstance. Nor does it appear desirable, because the designer usually needs to retain control over the upstream (or downstream) boundary shape. The problem is further complicated by the fact that, unlike the two-dimensional case, and irrespective of the upstream boundary shape, some prescribed lateral velocity distributions do not have proper solutions.

A semi-implicit level set method for multiphase flows and fluid-structure interaction problems

NASA Astrophysics Data System (ADS)

Cottet, Georges-Henri; Maitre, Emmanuel

2016-06-01

In this paper we present a novel semi-implicit time-discretization of the level set method introduced in [8] for fluid-structure interaction problems. The idea stems from a linear stability analysis derived on a simplified one-dimensional problem. The semi-implicit scheme relies on a simple filter operating as a pre-processing on the level set function. It applies to multiphase flows driven by surface tension as well as to fluid-structure interaction problems. The semi-implicit scheme avoids the stability constraints that explicit scheme need to satisfy and reduces significantly the computational cost. It is validated through comparisons with the original explicit scheme and refinement studies on two-dimensional benchmarks.

A 3-dimensional mass conserving element for compressible flows

NASA Technical Reports Server (NTRS)

Fix, G.; Suri, M.

1985-01-01

A variety of finite element schemes has been used in the numerical approximation of compressible flows particularly in underwater acoustics. In many instances instabilities have been generated due to the lack of mass conservation. Two- and three-dimensional elements are developed which avoid these problems.

Implementation of Finite Volume based Navier Stokes Algorithm Within General Purpose Flow Network Code

NASA Technical Reports Server (NTRS)

Schallhorn, Paul; Majumdar, Alok

2012-01-01

This paper describes a finite volume based numerical algorithm that allows multi-dimensional computation of fluid flow within a system level network flow analysis. There are several thermo-fluid engineering problems where higher fidelity solutions are needed that are not within the capacity of system level codes. The proposed algorithm will allow NASA's Generalized Fluid System Simulation Program (GFSSP) to perform multi-dimensional flow calculation within the framework of GFSSP s typical system level flow network consisting of fluid nodes and branches. The paper presents several classical two-dimensional fluid dynamics problems that have been solved by GFSSP's multi-dimensional flow solver. The numerical solutions are compared with the analytical and benchmark solution of Poiseulle, Couette and flow in a driven cavity.

Retraining Eastern Europe.

ERIC Educational Resources Information Center

O'Connor, Robert

1992-01-01

If market economies are to emerge from the ruins of the Soviet empire, people need training in almost everything. Efforts are being made in the areas of management, service, accounting, and police work. (Author/JOW)

Metal Oxide Gas Sensor Drift Compensation Using a Two-Dimensional Classifier Ensemble

PubMed Central

Liu, Hang; Chu, Renzhi; Tang, Zhenan

2015-01-01

Sensor drift is the most challenging problem in gas sensing at present. We propose a novel two-dimensional classifier ensemble strategy to solve the gas discrimination problem, regardless of the gas concentration, with high accuracy over extended periods of time. This strategy is appropriate for multi-class classifiers that consist of combinations of pairwise classifiers, such as support vector machines. We compare the performance of the strategy with those of competing methods in an experiment based on a public dataset that was compiled over a period of three years. The experimental results demonstrate that the two-dimensional ensemble outperforms the other methods considered. Furthermore, we propose a pre-aging process inspired by that applied to the sensors to improve the stability of the classifier ensemble. The experimental results demonstrate that the weight of each multi-class classifier model in the ensemble remains fairly static before and after the addition of new classifier models to the ensemble, when a pre-aging procedure is applied. PMID:25942640

Three-Particle Complexes in Two-Dimensional Semiconductors

NASA Astrophysics Data System (ADS)

Ganchev, Bogdan; Drummond, Neil; Aleiner, Igor; Fal'ko, Vladimir

2015-03-01

We evaluate binding energies of trions X±, excitons bound by a donor or acceptor charge XD (A ) , and overcharged acceptors or donors in two-dimensional atomic crystals by mapping the three-body problem in two dimensions onto one particle in a three-dimensional potential treatable by a purposely developed boundary-matching-matrix method. We find that in monolayers of transition metal dichalcogenides the dissociation energy of X± is typically much larger than that of localized exciton complexes, so that trions are more resilient to heating, despite the fact that their recombination line in optics is less redshifted from the exciton line than the line of XD (A ) .

Scaling between Wind Tunnels-Results Accuracy in Two-Dimensional Testing

NASA Astrophysics Data System (ADS)

Rasuo, Bosko

The establishment of exact two-dimensional flow conditions in wind tunnels is a very difficult problem. This has been evident for wind tunnels of all types and scales. In this paper, the principal factors that influence the accuracy of two-dimensional wind tunnel test results are analyzed. The influences of the Reynolds number, Mach number and wall interference with reference to solid and flow blockage (blockage of wake) as well as the influence of side-wall boundary layer control are analyzed. Interesting results are brought to light regarding the Reynolds number effects of the test model versus the Reynolds number effects of the facility in subsonic and transonic flow.

Two-dimensional supersonic nonlinear Schrödinger flow past an extended obstacle

NASA Astrophysics Data System (ADS)

El, G. A.; Kamchatnov, A. M.; Khodorovskii, V. V.; Annibale, E. S.; Gammal, A.

2009-10-01

Supersonic flow of a superfluid past a slender impenetrable macroscopic obstacle is studied in the framework of the two-dimensional (2D) defocusing nonlinear Schrödinger (NLS) equation. This problem is of fundamental importance as a dispersive analog of the corresponding classical gas-dynamics problem. Assuming the oncoming flow speed is sufficiently high, we asymptotically reduce the original boundary-value problem for a steady flow past a slender body to the one-dimensional dispersive piston problem described by the nonstationary NLS equation, in which the role of time is played by the stretched x coordinate and the piston motion curve is defined by the spatial body profile. Two steady oblique spatial dispersive shock waves (DSWs) spreading from the pointed ends of the body are generated in both half planes. These are described analytically by constructing appropriate exact solutions of the Whitham modulation equations for the front DSW and by using a generalized Bohr-Sommerfeld quantization rule for the oblique dark soliton fan in the rear DSW. We propose an extension of the traditional modulation description of DSWs to include the linear “ship-wave” pattern forming outside the nonlinear modulation region of the front DSW. Our analytic results are supported by direct 2D unsteady numerical simulations and are relevant to recent experiments on Bose-Einstein condensates freely expanding past obstacles.

Improved finite element methodology for integrated thermal structural analysis

NASA Technical Reports Server (NTRS)

Dechaumphai, P.; Thornton, E. A.

1982-01-01

An integrated thermal-structural finite element approach for efficient coupling of thermal and structural analysis is presented. New thermal finite elements which yield exact nodal and element temperatures for one dimensional linear steady state heat transfer problems are developed. A nodeless variable formulation is used to establish improved thermal finite elements for one dimensional nonlinear transient and two dimensional linear transient heat transfer problems. The thermal finite elements provide detailed temperature distributions without using additional element nodes and permit a common discretization with lower order congruent structural finite elements. The accuracy of the integrated approach is evaluated by comparisons with analytical solutions and conventional finite element thermal structural analyses for a number of academic and more realistic problems. Results indicate that the approach provides a significant improvement in the accuracy and efficiency of thermal stress analysis for structures with complex temperature distributions.

Using Three-Dimensional Printing to Fabricate a Tubing Connector for Dilation and Evacuation.

PubMed

Stitely, Michael L; Paterson, Helen

2016-02-01

This is a proof-of-concept study to show that simple instrumentation problems encountered in surgery can be solved by fabricating devices using a three-dimensional printer. The device used in the study is a simple tubing connector fashioned to connect two segments of suction tubing used in a surgical procedure where no commercially available product for this use is available through our usual suppliers in New Zealand. A cylindrical tubing connector was designed using three-dimensional printing design software. The tubing connector was fabricated using the Makerbot Replicator 2X three-dimensional printer. The connector was used in 15 second-trimester dilation and evacuation procedures. Data forms were completed by the primary operating surgeon. Descriptive statistics were used with the expectation that the device would function as intended in all cases. The three-dimensional printed tubing connector functioned as intended in all 15 instances. Commercially available three-dimensional printing technology can be used to overcome simple instrumentation problems encountered during gynecologic surgical procedures.

Desertification in Russia: Problems and solutions (an example in the Republic of Kalmykia-Khalmg Tangch).

PubMed

Zonn, I S

1995-01-01

During the second half of the 20th century Kalmykia has undergone severe desertification. Under Soviet rule, rangelands were increasingly devoted to animal production, and pastures were converted to cropland in a campaign to increase crops. Pastures were grazed two to three times their sustainable production, saiga populations and habitat greatly decreased, more than 17 million ha were subjected to wind erosion, 380,000 ha were transformed into moving sands, and 106,000 ha were ruined by secondary salinization and waterlogging. By the 1990s almost 80% of the Republic had undergone desertification, and 13% had been transformed into a true desert. In 1986 the General Scheme of Desertification Control was formulated. The scheme called for rotating pastures, reclaiming blown sand using silviculture, tilling overgrazed pastures and sowing fodder plants, and developing water supplies for pastures. In its early years the scheme has been successful. But the management of restored pastures usually reverts to the same farms responsible for the poor conditions, and there is great apprehension that degradation could reoccur. This case study concludes that the general cattle and agriculture development in Kalmykia is unviable for ecological and economic reasons, that Kalmykia should implement an adaptive policy oriented toward conservation and accommodating the interrelation and variability of land resources, that the desertification problem can be solved only by changing agrarian policy as a whole, and that a desertification control program must become an integral part of economic and social development of the Republic.

HUFF, a One-Dimensional Hydrodynamics Code for Strong Shocks

DTIC Science & Technology

1978-12-01

results for two sample problems. The first problem discussed is a one-kiloton nuclear burst in infinite sea level air. The second problem is the one...of HUFF as an effective first order hydro- dynamic computer code. 1 KT Explosion The one-kiloton nuclear explosion in infinite sea level air was

Three-dimensional Finite Element Formulation and Scalable Domain Decomposition for High Fidelity Rotor Dynamic Analysis

NASA Technical Reports Server (NTRS)

Datta, Anubhav; Johnson, Wayne R.

2009-01-01

This paper has two objectives. The first objective is to formulate a 3-dimensional Finite Element Model for the dynamic analysis of helicopter rotor blades. The second objective is to implement and analyze a dual-primal iterative substructuring based Krylov solver, that is parallel and scalable, for the solution of the 3-D FEM analysis. The numerical and parallel scalability of the solver is studied using two prototype problems - one for ideal hover (symmetric) and one for a transient forward flight (non-symmetric) - both carried out on up to 48 processors. In both hover and forward flight conditions, a perfect linear speed-up is observed, for a given problem size, up to the point of substructure optimality. Substructure optimality and the linear parallel speed-up range are both shown to depend on the problem size as well as on the selection of the coarse problem. With a larger problem size, linear speed-up is restored up to the new substructure optimality. The solver also scales with problem size - even though this conclusion is premature given the small prototype grids considered in this study.

Ermakov's Superintegrable Toy and Nonlocal Symmetries

NASA Astrophysics Data System (ADS)

Leach, P. G. L.; Karasu Kalkanli, A.; Nucci, M. C.; Andriopoulos, K.

2005-11-01

We investigate the symmetry properties of a pair of Ermakov equations. The system is superintegrable and yet possesses only three Lie point symmetries with the algebra sl(2, R). The number of point symmetries is insufficient and the algebra unsuitable for the complete specification of the system. We use the method of reduction of order to reduce the nonlinear fourth-order system to a third-order system comprising a linear second-order equation and a conservation law. We obtain the representation of the complete symmetry group from this system. Four of the required symmetries are nonlocal and the algebra is the direct sum of a one-dimensional Abelian algebra with the semidirect sum of a two-dimensional solvable algebra with a two-dimensional Abelian algebra. The problem illustrates the difficulties which can arise in very elementary systems. Our treatment demonstrates the existence of possible routes to overcome these problems in a systematic fashion.

Computer analysis of multicircuit shells of revolution by the field method

NASA Technical Reports Server (NTRS)

Cohen, G. A.

1975-01-01

The field method, presented previously for the solution of even-order linear boundary value problems defined on one-dimensional open branch domains, is extended to boundary value problems defined on one-dimensional domains containing circuits. This method converts the boundary value problem into two successive numerically stable initial value problems, which may be solved by standard forward integration techniques. In addition, a new method for the treatment of singular boundary conditions is presented. This method, which amounts to a partial interchange of the roles of force and displacement variables, is problem independent with respect to both accuracy and speed of execution. This method was implemented in a computer program to calculate the static response of ring stiffened orthotropic multicircuit shells of revolution to asymmetric loads. Solutions are presented for sample problems which illustrate the accuracy and efficiency of the method.

Implementation of equivalent domain integral method in the two-dimensional analysis of mixed mode problems

NASA Technical Reports Server (NTRS)

Raju, I. S.; Shivakumar, K. N.

1989-01-01

An equivalent domain integral (EDI) method for calculating J-intergrals for two-dimensional cracked elastic bodies is presented. The details of the method and its implementation are presented for isoparametric elements. The total and product integrals consist of the sum of an area of domain integral and line integrals on the crack faces. The line integrals vanish only when the crack faces are traction free and the loading is either pure mode 1 or pure mode 2 or a combination of both with only the square-root singular term in the stress field. The EDI method gave accurate values of the J-integrals for two mode I and two mixed mode problems. Numerical studies showed that domains consisting of one layer of elements are sufficient to obtain accurate J-integral values. Two procedures for separating the individual modes from the domain integrals are presented. The procedure that uses the symmetric and antisymmetric components of the stress and displacement fields to calculate the individual modes gave accurate values of the integrals for all problems analyzed. The EDI method when applied to a problem of an interface crack in two different materials showed that the mode 1 and mode 2 components are domain dependent while the total integral is not. This behavior is caused by the presence of the oscillatory part of the singularity in bimaterial crack problems. The EDI method, thus, shows behavior similar to the virtual crack closure method for bimaterial problems.

Numerical computations on one-dimensional inverse scattering problems

NASA Technical Reports Server (NTRS)

Dunn, M. H.; Hariharan, S. I.

1983-01-01

An approximate method to determine the index of refraction of a dielectric obstacle is presented. For simplicity one dimensional models of electromagnetic scattering are treated. The governing equations yield a second order boundary value problem, in which the index of refraction appears as a functional parameter. The availability of reflection coefficients yield two additional boundary conditions. The index of refraction by a k-th order spline which can be written as a linear combination of B-splines is approximated. For N distinct reflection coefficients, the resulting N boundary value problems yield a system of N nonlinear equations in N unknowns which are the coefficients of the B-splines.

53. Historic American Buildings Survey Photo by Slevin Photo Taken: ...

Library of Congress Historic Buildings Survey, Historic Engineering Record, Historic Landscapes Survey

53. Historic American Buildings Survey Photo by Slevin Photo Taken: April 1908 RUINS OF MISSION QUANDRANGLE - Mission San Carlos Borromeo, Rio Road & Lausen Drive, Carmel-by-the-Sea, Monterey County, CA

2. Historic American Buildings Survey Frederick D. Nichols, Photographer August ...

Library of Congress Historic Buildings Survey, Historic Engineering Record, Historic Landscapes Survey

2. Historic American Buildings Survey Frederick D. Nichols, Photographer August 1937 VIEW OF CHURCH, LOOKING SOUTHEAST - San Cayetano de Calabasas (Mission, Ruins), Santa Cruz River Vicinity, Nogales, Santa Cruz County, AZ

3. Historic American Buildings Survey Frederick D. Nichols, Photographer September ...

Library of Congress Historic Buildings Survey, Historic Engineering Record, Historic Landscapes Survey

3. Historic American Buildings Survey Frederick D. Nichols, Photographer September 1937 VIEW OF CHURCH - LOOKING NORTHWEST - San Cayetano de Calabasas (Mission, Ruins), Santa Cruz River Vicinity, Nogales, Santa Cruz County, AZ

28. Photocopy of October 25, 1926, photograph in San Francisco ...

Library of Congress Historic Buildings Survey, Historic Engineering Record, Historic Landscapes Survey

28. Photocopy of October 25, 1926, photograph in San Francisco Chronicle Library. RUINS OF TOWERS, LOOKING NORTH - Mission Santa Clara de Asis, Franklin & Grant Streets, Santa Clara, Santa Clara County, CA

29. Photocopy of October 25, 1926, photograph in San Francisco ...

Library of Congress Historic Buildings Survey, Historic Engineering Record, Historic Landscapes Survey

29. Photocopy of October 25, 1926, photograph in San Francisco Chronicle Library. VIEW INSIDE RUINED TOWERS - Mission Santa Clara de Asis, Franklin & Grant Streets, Santa Clara, Santa Clara County, CA

Application of TVD schemes for the Euler equations of gas dynamics. [total variation diminishing for nonlinear hyperbolic systems

NASA Technical Reports Server (NTRS)

Yee, H. C.; Warming, R. F.; Harten, A.

1985-01-01

First-order, second-order, and implicit total variation diminishing (TVD) schemes are reviewed using the modified flux approach. Some transient and steady-state calculations are then carried out to illustrate the applicability of these schemes to the Euler equations. It is shown that the second-order explicit TVD schemes generate good shock resolution for both transient and steady-state one-dimensional and two-dimensional problems. Numerical experiments for a quasi-one-dimensional nozzle problem show that the second-order implicit TVD scheme produces a fairly rapid convergence rate and remains stable even when running with a Courant number of 10 to the 6th.

Open Quantum Random Walks on the Half-Line: The Karlin-McGregor Formula, Path Counting and Foster's Theorem

NASA Astrophysics Data System (ADS)

Jacq, Thomas S.; Lardizabal, Carlos F.

2017-11-01

In this work we consider open quantum random walks on the non-negative integers. By considering orthogonal matrix polynomials we are able to describe transition probability expressions for classes of walks via a matrix version of the Karlin-McGregor formula. We focus on absorbing boundary conditions and, for simpler classes of examples, we consider path counting and the corresponding combinatorial tools. A non-commutative version of the gambler's ruin is studied by obtaining the probability of reaching a certain fortune and the mean time to reach a fortune or ruin in terms of generating functions. In the case of the Hadamard coin, a counting technique for boundary restricted paths in a lattice is also presented. We discuss an open quantum version of Foster's Theorem for the expected return time together with applications.

Romantic versus scientific perspective: the ruins of Radlin palace in Wielkopolska region in the light of remote sensing techniques

NASA Astrophysics Data System (ADS)

Wilgocka, Aleksandra; Ruciński, Dominik; RÄ czkowski, Włodzimierz

2015-06-01

Although ruins of palace in Radlin, localized in Wielkopolska Region (Poland), could have been a great inspiration for romantic landscape painters, they were hardly considered as the subject of artistic interest. Nevertheless they stand as a marker in a landscape as a romantic background for the village on one hand and a memento for the neighbouring graveyard on another. Small scale excavations carried out in late 1950s with historical maps and analysis of still standing remains gave a general idea about wings order, localisation of main entrance and communication routs inside courtyard. Those early research thereby were the first step to change the meaning of this place from romantic to more scientific. New remote sensing technology allows move even further into scientific direction. The ruins in Radlin have been included into project ArchEO - archaeological applications of Earth Observation techniques. The main aim of the project in case of Radlin is an attempt to answer the question to what extent very high resolution optical satellite imagery might allow better understanding the spatial structure of the place. The various processing techniques were applied to facilitate the detection of archaeological features' impact on the vegetation condition. It enabled to assess the usefulness of satellite based data in recognizing specific archaeological remains. Thus, potential and limitations of satellite imagery versus other sources of spatial information like historic maps, excavation results, aerial photographs and Lidar will be discussed.

Optimal swimming of a sheet.

PubMed

Montenegro-Johnson, Thomas D; Lauga, Eric

2014-06-01

Propulsion at microscopic scales is often achieved through propagating traveling waves along hairlike organelles called flagella. Taylor's two-dimensional swimming sheet model is frequently used to provide insight into problems of flagellar propulsion. We derive numerically the large-amplitude wave form of the two-dimensional swimming sheet that yields optimum hydrodynamic efficiency: the ratio of the squared swimming speed to the rate-of-working of the sheet against the fluid. Using the boundary element method, we show that the optimal wave form is a front-back symmetric regularized cusp that is 25% more efficient than the optimal sine wave. This optimal two-dimensional shape is smooth, qualitatively different from the kinked form of Lighthill's optimal three-dimensional flagellum, not predicted by small-amplitude theory, and different from the smooth circular-arc-like shape of active elastic filaments.

Existence of solutions of a two-dimensional boundary value problem for a system of nonlinear equations arising in growing cell populations.

PubMed

Jeribi, Aref; Krichen, Bilel; Mefteh, Bilel

2013-01-01

In the paper [A. Ben Amar, A. Jeribi, and B. Krichen, Fixed point theorems for block operator matrix and an application to a structured problem under boundary conditions of Rotenberg's model type, to appear in Math. Slovaca. (2014)], the existence of solutions of the two-dimensional boundary value problem (1) and (2) was discussed in the product Banach space L(p)×L(p) for p∈(1, ∞). Due to the lack of compactness on L1 spaces, the analysis did not cover the case p=1. The purpose of this work is to extend the results of Ben Amar et al. to the case p=1 by establishing new variants of fixed-point theorems for a 2×2 operator matrix, involving weakly compact operators.

Two-Dimensional and Three-Dimensional Ultrasound of Artificial Skin.

PubMed

Wortsman, Ximena; Navarrete, Nelson

2017-01-01

Wound healing may be a difficult problem, and variable types of artificial skin prototypes have been developed for supporting this process. Using ultrasound, we studied 4 cellulose-derived artificial skin prototypes and assessed their two-dimensional and three-dimensional morphology. These prototypes were identified on ultrasound both on in vitro and in vivo studies. They allowed the sonographic observation of deeper layers on different types of surfaces of the body with good definition on the in vivo examinations performed on healthy skin and cutaneous ulcers. The ultrasound detection of these artificial biomaterials may potentially support the noninvasive monitoring of wound healing. © 2016 by the American Institute of Ultrasound in Medicine.

Direct and inverse problems of studying the properties of multilayer nanostructures based on a two-dimensional model of X-ray reflection and scattering

NASA Astrophysics Data System (ADS)

Khachaturov, R. V.

2014-06-01

A mathematical model of X-ray reflection and scattering by multilayered nanostructures in the quasi-optical approximation is proposed. X-ray propagation and the electric field distribution inside the multilayered structure are considered with allowance for refraction, which is taken into account via the second derivative with respect to the depth of the structure. This model is used to demonstrate the possibility of solving inverse problems in order to determine the characteristics of irregularities not only over the depth (as in the one-dimensional problem) but also over the length of the structure. An approximate combinatorial method for system decomposition and composition is proposed for solving the inverse problems.

Progress on a Taylor weak statement finite element algorithm for high-speed aerodynamic flows

NASA Technical Reports Server (NTRS)

Baker, A. J.; Freels, J. D.

1989-01-01

A new finite element numerical Computational Fluid Dynamics (CFD) algorithm has matured to the point of efficiently solving two-dimensional high speed real-gas compressible flow problems in generalized coordinates on modern vector computer systems. The algorithm employs a Taylor Weak Statement classical Galerkin formulation, a variably implicit Newton iteration, and a tensor matrix product factorization of the linear algebra Jacobian under a generalized coordinate transformation. Allowing for a general two-dimensional conservation law system, the algorithm has been exercised on the Euler and laminar forms of the Navier-Stokes equations. Real-gas fluid properties are admitted, and numerical results verify solution accuracy, efficiency, and stability over a range of test problem parameters.

Initial values for the integration scheme to compute the eigenvalues for propagation in ducts

NASA Technical Reports Server (NTRS)

Eversman, W.

1977-01-01

A scheme for the calculation of eigenvalues in the problem of acoustic propagation in a two-dimensional duct is described. The computation method involves changing the coupled transcendental nonlinear algebraic equations into an initial value problem involving a nonlinear ordinary differential equation. The simplest approach is to use as initial values the hardwall eigenvalues and to integrate away from these values as the admittance varies from zero to its actual value with a linear variation. The approach leads to a powerful root finding routine capable of computing the transverse and axial wave numbers for two-dimensional ducts for any frequency, lining, admittance and Mach number without requiring initial guesses or starting points.

A second-order accurate kinetic-theory-based method for inviscid compressible flows

NASA Technical Reports Server (NTRS)

Deshpande, Suresh M.

1986-01-01

An upwind method for the numerical solution of the Euler equations is presented. This method, called the kinetic numerical method (KNM), is based on the fact that the Euler equations are moments of the Boltzmann equation of the kinetic theory of gases when the distribution function is Maxwellian. The KNM consists of two phases, the convection phase and the collision phase. The method is unconditionally stable and explicit. It is highly vectorizable and can be easily made total variation diminishing for the distribution function by a suitable choice of the interpolation strategy. The method is applied to a one-dimensional shock-propagation problem and to a two-dimensional shock-reflection problem.

Pecos National Monument, New Mexico: Its Geologic Setting

USGS Publications Warehouse

Johnson, Ross Byron

1969-01-01

The ruins of the pueblos and missions of Pecos lie on the east bank of Glorieta Creek near its junction with the Pecos River at the south end of the Sangre de Cristo Mountains in north-central New Mexico. Here the Pecos River and Glorieta Creek have formed a broad rolling valley in which the red adobe walls of the mission church stand as a striking monument to a historic past. This is beautiful country; the bright hues of red rocks are complemented by the varied greens of the junipers, pi?ons, and ponderosa pines. Northward the Sangre de Cristo Mountains stretch for miles in a blue mist toward the Truchas Peaks and forests of the Pecos Wilderness. A few miles south of the ruins the steep high escarpment of Glorieta Mesa marks, in a general way, the southern termination of the Rocky Mountain System, which here is represented by the Sangre de Cristos. The escarpment of Glorieta Mesa has been formed largely by the Pecos River and its tributaries eroding the soft sedimentary layers. The Pecos flows southward from the high mountains in the north, parallels the mesa escarpment for 15 miles, and breaches the mesa near San Jose. About 1-1/2 miles southwest of the Pecos ruins at Cerro de Escobas is the highest point on Glorieta Mesa. It is the most conspicuous feature of the local landscape and rises to an elevation of 8,212 feet - 1,270 feet above the ruins. The slope of the escarpment here is very steep, rising 6 feet in every 10 horizontal feet. Along the north side of the Glorieta Mesa escarpment is a 30-mile-long natural pass around the south end of the Sangre de Cristos that extends from Canoncito on the west to Starvation Peak on the east (fig. 1). The elevation of the pass is greater than 6,000 feet at all places, and it reaches its summit of 7,432 feet near the village of Glorieta near the west end of the pass. This pass has been used as a major travel route for more than 800 years by the Indians, Spanish, and Americans. The famous Santa Fe Trail passed through here and was superseded by the railroad, whose main-line tracks closely parallel the traces of the old wagon ruts. The modern four-lane divided highway, Interstate Highway I-25, carries high-speed automotive traffic through Glorieta Pass alongside the Santa Fe Railway. Glorieta Pass has been the locale of many important historical events, including the passage of Coronado's expedition in 1540-41; the construction of the two large mission churches at Pecos Pueblo; the capture and imprisonment of the men of the Texas Expedition in 1841; the passage of the American Army under General Kearny on its way to Santa Fe, Chihuahua, and California in 1846; and the Civil War Battle of Glorieta Pass in 1862.

Numerical applications of the advective-diffusive codes for the inner magnetosphere

NASA Astrophysics Data System (ADS)

Aseev, N. A.; Shprits, Y. Y.; Drozdov, A. Y.; Kellerman, A. C.

2016-11-01

In this study we present analytical solutions for convection and diffusion equations. We gather here the analytical solutions for the one-dimensional convection equation, the two-dimensional convection problem, and the one- and two-dimensional diffusion equations. Using obtained analytical solutions, we test the four-dimensional Versatile Electron Radiation Belt code (the VERB-4D code), which solves the modified Fokker-Planck equation with additional convection terms. The ninth-order upwind numerical scheme for the one-dimensional convection equation shows much more accurate results than the results obtained with the third-order scheme. The universal limiter eliminates unphysical oscillations generated by high-order linear upwind schemes. Decrease in the space step leads to convergence of a numerical solution of the two-dimensional diffusion equation with mixed terms to the analytical solution. We compare the results of the third- and ninth-order schemes applied to magnetospheric convection modeling. The results show significant differences in electron fluxes near geostationary orbit when different numerical schemes are used.

A multidomain spectral collocation method for the Stokes problem

NASA Technical Reports Server (NTRS)

Landriani, G. Sacchi; Vandeven, H.

1989-01-01

A multidomain spectral collocation scheme is proposed for the approximation of the two-dimensional Stokes problem. It is shown that the discrete velocity vector field is exactly divergence-free and we prove error estimates both for the velocity and the pressure.

Highly Parallel Alternating Directions Algorithm for Time Dependent Problems

NASA Astrophysics Data System (ADS)

Ganzha, M.; Georgiev, K.; Lirkov, I.; Margenov, S.; Paprzycki, M.

2011-11-01

In our work, we consider the time dependent Stokes equation on a finite time interval and on a uniform rectangular mesh, written in terms of velocity and pressure. For this problem, a parallel algorithm based on a novel direction splitting approach is developed. Here, the pressure equation is derived from a perturbed form of the continuity equation, in which the incompressibility constraint is penalized in a negative norm induced by the direction splitting. The scheme used in the algorithm is composed of two parts: (i) velocity prediction, and (ii) pressure correction. This is a Crank-Nicolson-type two-stage time integration scheme for two and three dimensional parabolic problems in which the second-order derivative, with respect to each space variable, is treated implicitly while the other variable is made explicit at each time sub-step. In order to achieve a good parallel performance the solution of the Poison problem for the pressure correction is replaced by solving a sequence of one-dimensional second order elliptic boundary value problems in each spatial direction. The parallel code is implemented using the standard MPI functions and tested on two modern parallel computer systems. The performed numerical tests demonstrate good level of parallel efficiency and scalability of the studied direction-splitting-based algorithm.

Bay of Naples, Italy

NASA Image and Video Library

1981-04-14

STS001-13-442 (14 April 1981) --- This photograph showing much of Italy was taken with a handheld 70mm camera from 276 kilometers above Earth as the NASA space shuttle Columbia and its crew were marking their last few hours in space on the historic first space mission utilizing a reusable vehicle. Included in the area of the frame are Golfo de Napoli, Napoli (Naples), Castellammare, Amalfi, Capri, Sorrento, Mt. Vesuvius and the ruins of Pompei. Astronauts John W. Young and Robert L. Crippen exposed eight magazines of color 70mm film during their two and one-third days in Earth orbit. Photo credit: NASA

Stationary Wavelet-based Two-directional Two-dimensional Principal Component Analysis for EMG Signal Classification

NASA Astrophysics Data System (ADS)

Ji, Yi; Sun, Shanlin; Xie, Hong-Bo

2017-06-01

Discrete wavelet transform (WT) followed by principal component analysis (PCA) has been a powerful approach for the analysis of biomedical signals. Wavelet coefficients at various scales and channels were usually transformed into a one-dimensional array, causing issues such as the curse of dimensionality dilemma and small sample size problem. In addition, lack of time-shift invariance of WT coefficients can be modeled as noise and degrades the classifier performance. In this study, we present a stationary wavelet-based two-directional two-dimensional principal component analysis (SW2D2PCA) method for the efficient and effective extraction of essential feature information from signals. Time-invariant multi-scale matrices are constructed in the first step. The two-directional two-dimensional principal component analysis then operates on the multi-scale matrices to reduce the dimension, rather than vectors in conventional PCA. Results are presented from an experiment to classify eight hand motions using 4-channel electromyographic (EMG) signals recorded in healthy subjects and amputees, which illustrates the efficiency and effectiveness of the proposed method for biomedical signal analysis.

Topics in Two-Dimensional Quantum Gravity and Chern-Simons Gauge Theories

NASA Astrophysics Data System (ADS)

Zemba, Guillermo Raul

A series of studies in two and three dimensional theories is presented. The two dimensional problems are considered in the framework of String Theory. The first one determines the region of integration in the space of inequivalent tori of a tadpole diagram in Closed String Field Theory, using the naive Witten three-string vertex. It is shown that every surface is counted an infinite number of times and the source of this behavior is identified. The second study analyzes the behavior of the discrete matrix model of two dimensional gravity without matter using a mathematically well-defined construction, confirming several conjectures and partial results from the literature. The studies in three dimensions are based on Chern Simons pure gauge theory. The first one deals with the projection of the theory onto a two-dimensional surface of constant time, whereas the second analyzes the large N behavior of the SU(N) theory and makes evident a duality symmetry between the only two parameters of the theory. (Copies available exclusively from MIT Libraries, Rm. 14-0551, Cambridge, MA 02139-4307. Ph. 617-253-5668; Fax 617-253 -1690.).

Problems of interaction longitudinal shear waves with V-shape tunnels defect

NASA Astrophysics Data System (ADS)

Popov, V. G.

2018-04-01

The problem of determining the two-dimensional dynamic stress state near a tunnel defect of V-shaped cross-section is solved. The defect is located in an infinite elastic medium, where harmonic longitudinal shear waves are propagating. The initial problem is reduced to a system of two singular integral or integro-differential equations with fixed singularities. A numerical method for solving these systems with regard to the true asymptotics of the unknown functions is developed.

Computer aided photographic engineering

NASA Technical Reports Server (NTRS)

Hixson, Jeffrey A.; Rieckhoff, Tom

1988-01-01

High speed photography is an excellent source of engineering data but only provides a two-dimensional representation of a three-dimensional event. Multiple cameras can be used to provide data for the third dimension but camera locations are not always available. A solution to this problem is to overlay three-dimensional CAD/CAM models of the hardware being tested onto a film or photographic image, allowing the engineer to measure surface distances, relative motions between components, and surface variations.

On the dynamics of the Ising model of cooperative phenomena

PubMed Central

Montroll, Elliott W.

1981-01-01

A two-dimensional (and to some degree three-dimensional) version of Glauber's one-dimensional spin relaxation model is described. The model is constructed to yield the Ising model of cooperative phenomena at equilibrium. A complete hierarchy of differential equations for multispin correlation functions is constructed. Some remarks are made concerning the solution of them for the initial value problem of determining the relaxation of an initial set of spin distributions. PMID:16592955

Boundary condition computational procedures for inviscid, supersonic steady flow field calculations

NASA Technical Reports Server (NTRS)

Abbett, M. J.

1971-01-01

Results are given of a comparative study of numerical procedures for computing solid wall boundary points in supersonic inviscid flow calculatons. Twenty five different calculation procedures were tested on two sample problems: a simple expansion wave and a simple compression (two-dimensional steady flow). A simple calculation procedure was developed. The merits and shortcomings of the various procedures are discussed, along with complications for three-dimensional and time-dependent flows.

Partial spline models for the inclusion of tropopause and frontal boundary information in otherwise smooth two- and three-dimensional objective analysis

NASA Technical Reports Server (NTRS)

Shiau, Jyh-Jen; Wahba, Grace; Johnson, Donald R.

1986-01-01

A new method, based on partial spline models, is developed for including specified discontinuities in otherwise smooth two- and three-dimensional objective analyses. The method is appropriate for including tropopause height information in two- and three-dimensinal temperature analyses, using the O'Sullivan-Wahba physical variational method for analysis of satellite radiance data, and may in principle be used in a combined variational analysis of observed, forecast, and climate information. A numerical method for its implementation is described and a prototype two-dimensional analysis based on simulated radiosonde and tropopause height data is shown. The method may also be appropriate for other geophysical problems, such as modeling the ocean thermocline, fronts, discontinuities, etc.

1. Historic American Buildings Survey Frederick D. Nichols, Photographer August ...

Library of Congress Historic Buildings Survey, Historic Engineering Record, Historic Landscapes Survey

1. Historic American Buildings Survey Frederick D. Nichols, Photographer August 1937 VIEW OF CHURCH - LOOKING NORTH EAST - San Cayetano de Calabasas (Mission, Ruins), Santa Cruz River Vicinity, Nogales, Santa Cruz County, AZ

4. Historic American Buildings Survey Frederick D. Nichols, Photographer August ...

Library of Congress Historic Buildings Survey, Historic Engineering Record, Historic Landscapes Survey

4. Historic American Buildings Survey Frederick D. Nichols, Photographer August 1937 CANAL IN EAST WALL OF CHURCH - San Cayetano de Calabasas (Mission, Ruins), Santa Cruz River Vicinity, Nogales, Santa Cruz County, AZ

3. Historic American Buildings Survey S.F. Chronicle Library Circa 1906 ...

Library of Congress Historic Buildings Survey, Historic Engineering Record, Historic Landscapes Survey

3. Historic American Buildings Survey S.F. Chronicle Library Circa 1906 RUINS, AFTER EARTHQUAKE AND FIRE (Looking North) - St. Francis Church, Historic View, 610 Vallejo Street, San Francisco, San Francisco County, CA

Guide to the Revised Ground-Water Flow and Heat Transport Simulator: HYDROTHERM - Version 3

USGS Publications Warehouse

Kipp, Kenneth L.; Hsieh, Paul A.; Charlton, Scott R.

2008-01-01

The HYDROTHERM computer program simulates multi-phase ground-water flow and associated thermal energy transport in three dimensions. It can handle high fluid pressures, up to 1 ? 109 pascals (104 atmospheres), and high temperatures, up to 1,200 degrees Celsius. This report documents the release of Version 3, which includes various additions, modifications, and corrections that have been made to the original simulator. Primary changes to the simulator include: (1) the ability to simulate unconfined ground-water flow, (2) a precipitation-recharge boundary condition, (3) a seepage-surface boundary condition at the land surface, (4) the removal of the limitation that a specified-pressure boundary also have a specified temperature, (5) a new iterative solver for the linear equations based on a generalized minimum-residual method, (6) the ability to use time- or depth-dependent functions for permeability, (7) the conversion of the program code to Fortran 90 to employ dynamic allocation of arrays, and (8) the incorporation of a graphical user interface (GUI) for input and output. The graphical user interface has been developed for defining a simulation, running the HYDROTHERM simulator interactively, and displaying the results. The combination of the graphical user interface and the HYDROTHERM simulator forms the HYDROTHERM INTERACTIVE (HTI) program. HTI can be used for two-dimensional simulations only. New features in Version 3 of the HYDROTHERM simulator have been verified using four test problems. Three problems come from the published literature and one problem was simulated by another partially saturated flow and thermal transport simulator. The test problems include: transient partially saturated vertical infiltration, transient one-dimensional horizontal infiltration, two-dimensional steady-state drainage with a seepage surface, and two-dimensional drainage with coupled heat transport. An example application to a hypothetical stratovolcano system with unconfined ground-water flow is presented in detail. It illustrates the use of HTI with the combination precipitation-recharge and seepage-surface boundary condition, and functions as a tutorial example problem for the new user.

Ultrahigh-Dimensional Multiclass Linear Discriminant Analysis by Pairwise Sure Independence Screening

PubMed Central

Pan, Rui; Wang, Hansheng; Li, Runze

2016-01-01

This paper is concerned with the problem of feature screening for multi-class linear discriminant analysis under ultrahigh dimensional setting. We allow the number of classes to be relatively large. As a result, the total number of relevant features is larger than usual. This makes the related classification problem much more challenging than the conventional one, where the number of classes is small (very often two). To solve the problem, we propose a novel pairwise sure independence screening method for linear discriminant analysis with an ultrahigh dimensional predictor. The proposed procedure is directly applicable to the situation with many classes. We further prove that the proposed method is screening consistent. Simulation studies are conducted to assess the finite sample performance of the new procedure. We also demonstrate the proposed methodology via an empirical analysis of a real life example on handwritten Chinese character recognition. PMID:28127109

Solution-adaptive finite element method in computational fracture mechanics

NASA Technical Reports Server (NTRS)

Min, J. B.; Bass, J. M.; Spradley, L. W.

1993-01-01

Some recent results obtained using solution-adaptive finite element method in linear elastic two-dimensional fracture mechanics problems are presented. The focus is on the basic issue of adaptive finite element method for validating the applications of new methodology to fracture mechanics problems by computing demonstration problems and comparing the stress intensity factors to analytical results.

Random Walk Method for Potential Problems

NASA Technical Reports Server (NTRS)

Krishnamurthy, T.; Raju, I. S.

2002-01-01

A local Random Walk Method (RWM) for potential problems governed by Lapalace's and Paragon's equations is developed for two- and three-dimensional problems. The RWM is implemented and demonstrated in a multiprocessor parallel environment on a Beowulf cluster of computers. A speed gain of 16 is achieved as the number of processors is increased from 1 to 23.

Parity and cobordism of free knots

DOE Office of Scientific and Technical Information (OSTI.GOV)

Manturov, Vassily O

2012-02-28

A simple invariant is constructed which obstructs a free knot to be truncated. In particular, this invariant provides an obstruction to the truncatedness of curves immersed in two-dimensional surfaces. A curve on an oriented two-dimensional surface S{sub g} is referred to as truncated (null-cobordant) if there exists a three-dimensional manifold M with boundary S{sub g} and a smooth proper map of a two-disc to M such that the image of the boundary of the disc coincides with the curve. The problem of truncatedness for free knots is solved in this paper using the notion of parity recently introduced by themore » author. Bibliography: 12 titles.« less

Modal Ring Method for the Scattering of Electromagnetic Waves

NASA Technical Reports Server (NTRS)

Baumeister, Kenneth J.; Kreider, Kevin L.

1993-01-01

The modal ring method for electromagnetic scattering from perfectly electric conducting (PEC) symmetrical bodies is presented. The scattering body is represented by a line of finite elements (triangular) on its outer surface. The infinite computational region surrounding the body is represented analytically by an eigenfunction expansion. The modal ring method effectively reduces the two dimensional scattering problem to a one-dimensional problem similar to the method of moments. The modal element method is capable of handling very high frequency scattering because it has a highly banded solution matrix.

Bounded solutions in a T-shaped waveguide and the spectral properties of the Dirichlet ladder

NASA Astrophysics Data System (ADS)

Nazarov, S. A.

2014-08-01

The Dirichlet problem is considered on the junction of thin quantum waveguides (of thickness h ≪ 1) in the shape of an infinite two-dimensional ladder. Passage to the limit as h → +0 is discussed. It is shown that the asymptotically correct transmission conditions at nodes of the corresponding one-dimensional quantum graph are Dirichlet conditions rather than the conventional Kirchhoff transmission conditions. The result is obtained by analyzing bounded solutions of a problem in the T-shaped waveguide that the boundary layer phenomenon.

Inverse energy cascades in three-dimensional turbulence

NASA Technical Reports Server (NTRS)

Hossain, Murshed

1991-01-01

Fully three-dimensional magnetohydrodynamic (MHD) turbulence at large kinetic and low magnetic Reynolds numbers is considered in the presence of a strong uniform magnetic field. It is shown by numerical simulation of a model of MHD that the energy inverse cascades to longer length scales when the interaction parameter is large. While the steady-state dynamics of the driven problem is three-dimensional in character, the behavior has resemblance to two-dimensional hydrodynamics. These results have implications in turbulence theory, MHD power generator, planetary dynamos, and fusion reactor blanket design.

Pseudopotential Method for Higher Partial Wave Scattering

DOE Office of Scientific and Technical Information (OSTI.GOV)

Idziaszek, Zbigniew; Centrum Fizyki Teoretycznej, Polska Akademia Nauk, 02-668 Warsaw; Calarco, Tommaso

2006-01-13

We present a zero-range pseudopotential applicable for all partial wave interactions between neutral atoms. For p and d waves, we derive effective pseudopotentials, which are useful for problems involving anisotropic external potentials. Finally, we consider two nontrivial applications of the p-wave pseudopotential: we solve analytically the problem of two interacting spin-polarized fermions confined in a harmonic trap, and we analyze the scattering of p-wave interacting particles in a quasi-two-dimensional system.

Modeling of thin-walled structures interacting with acoustic media as constrained two-dimensional continua

NASA Astrophysics Data System (ADS)

Rabinskiy, L. N.; Zhavoronok, S. I.

2018-04-01

The transient interaction of acoustic media and elastic shells is considered on the basis of the transition function approach. The three-dimensional hyperbolic initial boundary-value problem is reduced to a two-dimensional problem of shell theory with integral operators approximating the acoustic medium effect on the shell dynamics. The kernels of these integral operators are determined by the elementary solution of the problem of acoustic waves diffraction at a rigid obstacle with the same boundary shape as the wetted shell surface. The closed-form elementary solution for arbitrary convex obstacles can be obtained at the initial interaction stages on the background of the so-called “thin layer hypothesis”. Thus, the shell–wave interaction model defined by integro-differential dynamic equations with analytically determined kernels of integral operators becomes hence two-dimensional but nonlocal in time. On the other hand, the initial interaction stage results in localized dynamic loadings and consequently in complex strain and stress states that require higher-order shell theories. Here the modified theory of I.N.Vekua–A.A.Amosov-type is formulated in terms of analytical continuum dynamics. The shell model is constructed on a two-dimensional manifold within a set of field variables, Lagrangian density, and constraint equations following from the boundary conditions “shifted” from the shell faces to its base surface. Such an approach allows one to construct consistent low-order shell models within a unified formal hierarchy. The equations of the N th-order shell theory are singularly perturbed and contain second-order partial derivatives with respect to time and surface coordinates whereas the numerical integration of systems of first-order equations is more efficient. Such systems can be obtained as Hamilton–de Donder–Weyl-type equations for the Lagrangian dynamical system. The Hamiltonian formulation of the elementary N th-order shell theory is here briefly described.

A BDDC Algorithm with Deluxe Scaling for Three-Dimensional H (curl) Problems

DOE PAGES

Dohrmann, Clark R.; Widlund, Olof B.

2015-04-28

In our paper, we present and analyze a BDDC algorithm for a class of elliptic problems in the three-dimensional H(curl) space. Compared with existing results, our condition number estimate requires fewer assumptions and also involves two fewer powers of log(H/h), making it consistent with optimal estimates for other elliptic problems. Here, H/his the maximum of Hi/hi over all subdomains, where Hi and hi are the diameter and the smallest element diameter for the subdomain Ωi. The analysis makes use of two recent developments. The first is our new approach to averaging across the subdomain interfaces, while the second is amore » new technical tool that allows arguments involving trace classes to be avoided. Furthermore, numerical examples are presented to confirm the theory and demonstrate the importance of the new averaging approach in certain cases.« less

Design of supercritical swept wings

NASA Technical Reports Server (NTRS)

Garabedian, P.; Mcfadden, G.

1982-01-01

Computational fluid dynamics are used to discuss problems inherent to transonic three-dimensional flow past supercritical swept wings. The formulation for a boundary value problem for the flow past the wing is provided, including consideration of weak shock waves and the use of parabolic coordinates. A swept wing code is developed which requires a mesh of 152 x 10 x 12 points and 200 time cycles. A formula for wave drag is calculated, based on the idea that the conservation form of the momentum equation becomes an entropy inequality measuring the drag, expressible in terms of a small-disturbance equation for a potential function in two dimensions. The entropy inequality has been incorporated in a two-dimensional code for the analysis of transonic flow over airfoils. A method of artificial viscosity is explored for optimum pressure distributions with design, and involves a free boundary problem considering speed over only a portion of the wing.

Stress Recovery and Error Estimation for Shell Structures

NASA Technical Reports Server (NTRS)

Yazdani, A. A.; Riggs, H. R.; Tessler, A.

2000-01-01

The Penalized Discrete Least-Squares (PDLS) stress recovery (smoothing) technique developed for two dimensional linear elliptic problems is adapted here to three-dimensional shell structures. The surfaces are restricted to those which have a 2-D parametric representation, or which can be built-up of such surfaces. The proposed strategy involves mapping the finite element results to the 2-D parametric space which describes the geometry, and smoothing is carried out in the parametric space using the PDLS-based Smoothing Element Analysis (SEA). Numerical results for two well-known shell problems are presented to illustrate the performance of SEA/PDLS for these problems. The recovered stresses are used in the Zienkiewicz-Zhu a posteriori error estimator. The estimated errors are used to demonstrate the performance of SEA-recovered stresses in automated adaptive mesh refinement of shell structures. The numerical results are encouraging. Further testing involving more complex, practical structures is necessary.

A 3-D turbulent flow analysis using finite elements with k-ɛ model

NASA Astrophysics Data System (ADS)

Okuda, H.; Yagawa, G.; Eguchi, Y.

1989-03-01

This paper describes the finite element turbulent flow analysis, which is suitable for three-dimensional large scale problems. The k-ɛ turbulence model as well as the conservation equations of mass and momentum are discretized in space using rather low order elements. Resulting coefficient matrices are evaluated by one-point quadrature in order to reduce the computational storage and the CPU cost. The time integration scheme based on the velocity correction method is employed to obtain steady state solutions. For the verification of this FEM program, two-dimensional plenum flow is simulated and compared with experiment. As the application to three-dimensional practical problems, the turbulent flows in the upper plenum of the fast breeder reactor are calculated for various boundary conditions.

Relevance feedback-based building recognition

NASA Astrophysics Data System (ADS)

Li, Jing; Allinson, Nigel M.

2010-07-01

Building recognition is a nontrivial task in computer vision research which can be utilized in robot localization, mobile navigation, etc. However, existing building recognition systems usually encounter the following two problems: 1) extracted low level features cannot reveal the true semantic concepts; and 2) they usually involve high dimensional data which require heavy computational costs and memory. Relevance feedback (RF), widely applied in multimedia information retrieval, is able to bridge the gap between the low level visual features and high level concepts; while dimensionality reduction methods can mitigate the high-dimensional problem. In this paper, we propose a building recognition scheme which integrates the RF and subspace learning algorithms. Experimental results undertaken on our own building database show that the newly proposed scheme appreciably enhances the recognition accuracy.

2. Historic American Buildings Survey Branan Sanders, Photographer March 1934 ...

Library of Congress Historic Buildings Survey, Historic Engineering Record, Historic Landscapes Survey

2. Historic American Buildings Survey Branan Sanders, Photographer March 1934 FRONT This picture also reproduce in 'Garden History of Georgia', opp. p. 1. - Santo Domingo Mission (Ruins), Near Altamaha Canal (River), Broadfield, Glynn County, GA

5. Detail view of masonry foundations and flue openings for ...

Library of Congress Historic Buildings Survey, Historic Engineering Record, Historic Landscapes Survey

5. Detail view of masonry foundations and flue openings for Jamaican Train. - Hacienda Azucarera El Coto, Sugar Mill Ruins, .5 Mi. SW of Rt. 347 Bridge Over Guanajibo River, San German, San German Municipio, PR

What is the latent structure of alcohol use disorders? A taxometric analysis of the Personality Assessment Inventory Alcohol Problems Scale in male and female prison inmates.

PubMed

Walters, Glenn D; Diamond, Pamela M; Magaletta, Philip R

2010-03-01

Three indicators derived from the Personality Assessment Inventory (PAI) Alcohol Problems scale (ALC)-tolerance/high consumption, loss of control, and negative social and psychological consequences-were subjected to taxometric analysis-mean above minus below a cut (MAMBAC), maximum covariance (MAXCOV), and latent mode factor analysis (L-Mode)-in 1,374 federal prison inmates (905 males, 469 females). Whereas the total sample yielded ambiguous results, the male subsample produced dimensional results, and the female subsample produced taxonic results. Interpreting these findings in light of previous taxometric research on alcohol abuse and dependence it is speculated that while alcohol use disorders may be taxonic in female offenders, they are probably both taxonic and dimensional in male offenders. Two models of male alcohol use disorder in males are considered, one in which the diagnostic features are categorical and the severity of symptomatology is dimensional, and one in which some diagnostic features (e.g., withdrawal) are taxonic and other features (e.g., social problems) are dimensional.

Conformal mapping and bound states in bent waveguides

NASA Astrophysics Data System (ADS)

Sadurní, E.; Schleich, W. P.

2010-12-01

Is it possible to trap a quantum particle in an open geometry? In this work we deal with the boundary value problem of the stationary Schroedinger (or Helmholtz) equation within a waveguide with straight segments and a rectangular bending. The problem can be reduced to a one-dimensional matrix Schroedinger equation using two descriptions: oblique modes and conformal coordinates. We use a corner-corrected WKB formalism to find the energies of the one-dimensional problem. It is shown that the presence of bound states is an effect due to the boundary alone, with no classical counterpart for this geometry. The conformal description proves to be simpler, as the coupling of transversal modes is not essential in this case.

Three dimensional elements with Lagrange multipliers for the modified couple stress theory

NASA Astrophysics Data System (ADS)

Kwon, Young-Rok; Lee, Byung-Chai

2018-07-01

Three dimensional mixed elements for the modified couple stress theory are proposed. The C1 continuity for the displacement field, which is required because of the curvature term in the variational form of the theory, is satisfied weakly by introducing a supplementary rotation as an independent variable and constraining the relation between the rotation and the displacement with a Lagrange multiplier vector. An additional constraint about the deviatoric curvature is also considered for three dimensional problems. Weak forms with one constraint and two constraints are derived, and four elements satisfying convergence criteria are developed by applying different approximations to each field of independent variables. The elements pass a [InlineEquation not available: see fulltext.] patch test for three dimensional problems. Numerical examples show that the additional constraint could be considered essential for the three dimensional elements, and one of the elements is recommended for practical applications via the comparison of the performances of the elements. In addition, all the proposed elements can represent the size effect well.

Numerical aerodynamic simulation facility. [for flows about three-dimensional configurations

NASA Technical Reports Server (NTRS)

Bailey, F. R.; Hathaway, A. W.

1978-01-01

Critical to the advancement of computational aerodynamics capability is the ability to simulate flows about three-dimensional configurations that contain both compressible and viscous effects, including turbulence and flow separation at high Reynolds numbers. Analyses were conducted of two solution techniques for solving the Reynolds averaged Navier-Stokes equations describing the mean motion of a turbulent flow with certain terms involving the transport of turbulent momentum and energy modeled by auxiliary equations. The first solution technique is an implicit approximate factorization finite-difference scheme applied to three-dimensional flows that avoids the restrictive stability conditions when small grid spacing is used. The approximate factorization reduces the solution process to a sequence of three one-dimensional problems with easily inverted matrices. The second technique is a hybrid explicit/implicit finite-difference scheme which is also factored and applied to three-dimensional flows. Both methods are applicable to problems with highly distorted grids and a variety of boundary conditions and turbulence models.

Monolithic multigrid methods for two-dimensional resistive magnetohydrodynamics

DOE PAGES

Adler, James H.; Benson, Thomas R.; Cyr, Eric C.; ...

2016-01-06

Magnetohydrodynamic (MHD) representations are used to model a wide range of plasma physics applications and are characterized by a nonlinear system of partial differential equations that strongly couples a charged fluid with the evolution of electromagnetic fields. The resulting linear systems that arise from discretization and linearization of the nonlinear problem are generally difficult to solve. In this paper, we investigate multigrid preconditioners for this system. We consider two well-known multigrid relaxation methods for incompressible fluid dynamics: Braess--Sarazin relaxation and Vanka relaxation. We first extend these to the context of steady-state one-fluid viscoresistive MHD. Then we compare the two relaxationmore » procedures within a multigrid-preconditioned GMRES method employed within Newton's method. To isolate the effects of the different relaxation methods, we use structured grids, inf-sup stable finite elements, and geometric interpolation. Furthermore, we present convergence and timing results for a two-dimensional, steady-state test problem.« less

Chaos and Robustness in a Single Family of Genetic Oscillatory Networks

PubMed Central

Fu, Daniel; Tan, Patrick; Kuznetsov, Alexey; Molkov, Yaroslav I.

2014-01-01

Genetic oscillatory networks can be mathematically modeled with delay differential equations (DDEs). Interpreting genetic networks with DDEs gives a more intuitive understanding from a biological standpoint. However, it presents a problem mathematically, for DDEs are by construction infinitely-dimensional and thus cannot be analyzed using methods common for systems of ordinary differential equations (ODEs). In our study, we address this problem by developing a method for reducing infinitely-dimensional DDEs to two- and three-dimensional systems of ODEs. We find that the three-dimensional reductions provide qualitative improvements over the two-dimensional reductions. We find that the reducibility of a DDE corresponds to its robustness. For non-robust DDEs that exhibit high-dimensional dynamics, we calculate analytic dimension lines to predict the dependence of the DDEs’ correlation dimension on parameters. From these lines, we deduce that the correlation dimension of non-robust DDEs grows linearly with the delay. On the other hand, for robust DDEs, we find that the period of oscillation grows linearly with delay. We find that DDEs with exclusively negative feedback are robust, whereas DDEs with feedback that changes its sign are not robust. We find that non-saturable degradation damps oscillations and narrows the range of parameter values for which oscillations exist. Finally, we deduce that natural genetic oscillators with highly-regular periods likely have solely negative feedback. PMID:24667178

Analysis of two dimensional signals via curvelet transform

NASA Astrophysics Data System (ADS)

Lech, W.; Wójcik, W.; Kotyra, A.; Popiel, P.; Duk, M.

2007-04-01

This paper describes an application of curvelet transform analysis problem of interferometric images. Comparing to two-dimensional wavelet transform, curvelet transform has higher time-frequency resolution. This article includes numerical experiments, which were executed on random interferometric image. In the result of nonlinear approximations, curvelet transform obtains matrix with smaller number of coefficients than is guaranteed by wavelet transform. Additionally, denoising simulations show that curvelet could be a very good tool to remove noise from images.

A Method to Formulate the Unit Cell for Density Functional Theory (DFT) Calculations of the Electronic Band Structure of Heterostructures of Two-dimensional Nanosheets

DTIC Science & Technology

2015-04-01

distribution is unlimited. i CONTENTS Page Introduction 1 Two-dimensional Material Geometry and Analogs with Close-packed Systems 1 Matching... System Lattice Vectors: An Optimization Problem 1 Generating the System Unit Cell 3 Transition Metal Dichalcogenides (TMDCS) with Mismatched... system being analyzed. The creation of a unit cell that accurately describes the system remains one of the largest challenges for DFT calculations

Ancient Maya astronomical tables from Xultun, Guatemala.

PubMed

Saturno, William A; Stuart, David; Aveni, Anthony F; Rossi, Franco

2012-05-11

Maya astronomical tables are recognized in bark-paper books from the Late Postclassic period (1300 to 1521 C.E.), but Classic period (200 to 900 C.E.) precursors have not been found. In 2011, a small painted room was excavated at the extensive ancient Maya ruins of Xultun, Guatemala, dating to the early 9th century C.E. The walls and ceiling of the room are painted with several human figures. Two walls also display a large number of delicate black, red, and incised hieroglyphs. Many of these hieroglyphs are calendrical in nature and relate astronomical computations, including at least two tables concerning the movement of the Moon, and perhaps Mars and Venus. These apparently represent early astronomical tables and may shed light on the later books.

Vertical drying of a suspension of sticks: Monte Carlo simulation for continuous two-dimensional problem

NASA Astrophysics Data System (ADS)

Lebovka, Nikolai I.; Tarasevich, Yuri Yu.; Vygornitskii, Nikolai V.

2018-02-01

The vertical drying of a two-dimensional colloidal film containing zero-thickness sticks (lines) was studied by means of kinetic Monte Carlo (MC) simulations. The continuous two-dimensional problem for both the positions and orientations was considered. The initial state before drying was produced using a model of random sequential adsorption with isotropic orientations of the sticks. During the evaporation, an upper interface falls with a linear velocity in the vertical direction, and the sticks undergo translational and rotational Brownian motions. The MC simulations were run at different initial number concentrations (the numbers of sticks per unit area), pi, and solvent evaporation rates, u . For completely dried films, the spatial distributions of the sticks, the order parameters, and the electrical conductivities of the films in both the horizontal, x , and vertical, y , directions were examined. Significant evaporation-driven self-assembly and stratification of the sticks in the vertical direction was observed. The extent of stratification increased with increasing values of u . The anisotropy of the electrical conductivity of the film can be finely regulated by changes in the values of pi and u .

Intersection of three-dimensional geometric surfaces

NASA Technical Reports Server (NTRS)

Crisp, V. K.; Rehder, J. J.; Schwing, J. L.

1985-01-01

Calculating the line of intersection between two three-dimensional objects and using the information to generate a third object is a key element in a geometry development system. Techniques are presented for the generation of three-dimensional objects, the calculation of a line of intersection between two objects, and the construction of a resultant third object. The objects are closed surfaces consisting of adjacent bicubic parametric patches using Bezier basis functions. The intersection determination involves subdividing the patches that make up the objects until they are approximately planar and then calculating the intersection between planes. The resulting straight-line segments are connected to form the curve of intersection. The polygons in the neighborhood of the intersection are reconstructed and put back into the Bezier representation. A third object can be generated using various combinations of the original two. Several examples are presented. Special cases and problems were encountered, and the method for handling them is discussed. The special cases and problems included intersection of patch edges, gaps between adjacent patches because of unequal subdivision, holes, or islands within patches, and computer round-off error.

Application of the coherent anomaly method to percolation

NASA Astrophysics Data System (ADS)

Takayasu, Misako; Takayasu, Hideki

1988-03-01

Applying the coherent anomaly method (CAM) to site percolation problems, we estimate the percolation threshold pc and critical exponents. We obtain pc=0.589, β=0.140, γ=2.426 on the two-dimensional square lattice. These values are in good agreement with the values already known. We also investigate higher-dimensional cases by this method.

Application of the Coherent Anomaly Method to Percolation

NASA Astrophysics Data System (ADS)

Takayasu, Misako; Takayasu, Hideki

Applying the coherent anomaly method (CAM) to site percolation problems, we estimate the percolation threshold ϱc and critical exponents. We obtain pc = 0.589, Β=0.140, Γ = 2.426 on the two-dimensional square lattice. These values are in good agreement with the values already known. We also investigate higher-dimensional cases by this method.

Moving boundary problems for a rarefied gas: Spatially one-dimensional case

NASA Astrophysics Data System (ADS)

Tsuji, Tetsuro; Aoki, Kazuo

2013-10-01

Unsteady flows of a rarefied gas in a full space caused by an oscillation of an infinitely wide plate in its normal direction are investigated numerically on the basis of the Bhatnagar-Gross-Krook (BGK) model of the Boltzmann equation. The paper aims at showing properties and difficulties inherent to moving boundary problems in kinetic theory of gases using a simple one-dimensional setting. More specifically, the following two problems are considered: (Problem I) the plate starts a forced harmonic oscillation (forced motion); (Problem II) the plate, which is subject to an external restoring force obeying Hooke’s law, is displaced from its equilibrium position and released (free motion). The physical interest in Problem I lies in the propagation of nonlinear acoustic waves in a rarefied gas, whereas that in Problem II in the decay rate of the oscillation of the plate. An accurate numerical method, which is capable of describing singularities caused by the oscillating plate, is developed on the basis of the method of characteristics and is applied to the two problems mentioned above. As a result, the unsteady behavior of the solution, such as the propagation of discontinuities and some weaker singularities in the molecular velocity distribution function, are clarified. Some results are also compared with those based on the existing method.

The quantum-field renormalization group in the problem of a growing phase boundary

DOE Office of Scientific and Technical Information (OSTI.GOV)

Antonov, N.V.; Vasil`ev, A.N.

1995-09-01

Within the quantum-field renormalization-group approach we examine the stochastic equation discussed by S.I. Pavlik in describing a randomly growing phase boundary. We show that, in contrast to Pavlik`s assertion, the model is not multiplicatively renormalizable and that its consistent renormalization-group analysis requires introducing an infinite number of counterterms and the respective coupling constants ({open_quotes}charge{close_quotes}). An explicit calculation in the one-loop approximation shows that a two-dimensional surface of renormalization-group points exits in the infinite-dimensional charge space. If the surface contains an infrared stability region, the problem allows for scaling with the nonuniversal critical dimensionalities of the height of the phase boundarymore » and time, {delta}{sub h} and {delta}{sub t}, which satisfy the exact relationship 2 {delta}{sub h}= {delta}{sub t} + d, where d is the dimensionality of the phase boundary. 23 refs., 1 tab.« less

Quantum key distribution session with 16-dimensional photonic states.

PubMed

Etcheverry, S; Cañas, G; Gómez, E S; Nogueira, W A T; Saavedra, C; Xavier, G B; Lima, G

2013-01-01

The secure transfer of information is an important problem in modern telecommunications. Quantum key distribution (QKD) provides a solution to this problem by using individual quantum systems to generate correlated bits between remote parties, that can be used to extract a secret key. QKD with D-dimensional quantum channels provides security advantages that grow with increasing D. However, the vast majority of QKD implementations has been restricted to two dimensions. Here we demonstrate the feasibility of using higher dimensions for real-world quantum cryptography by performing, for the first time, a fully automated QKD session based on the BB84 protocol with 16-dimensional quantum states. Information is encoded in the single-photon transverse momentum and the required states are dynamically generated with programmable spatial light modulators. Our setup paves the way for future developments in the field of experimental high-dimensional QKD.

Quantum key distribution session with 16-dimensional photonic states

NASA Astrophysics Data System (ADS)

Etcheverry, S.; Cañas, G.; Gómez, E. S.; Nogueira, W. A. T.; Saavedra, C.; Xavier, G. B.; Lima, G.

2013-07-01

The secure transfer of information is an important problem in modern telecommunications. Quantum key distribution (QKD) provides a solution to this problem by using individual quantum systems to generate correlated bits between remote parties, that can be used to extract a secret key. QKD with D-dimensional quantum channels provides security advantages that grow with increasing D. However, the vast majority of QKD implementations has been restricted to two dimensions. Here we demonstrate the feasibility of using higher dimensions for real-world quantum cryptography by performing, for the first time, a fully automated QKD session based on the BB84 protocol with 16-dimensional quantum states. Information is encoded in the single-photon transverse momentum and the required states are dynamically generated with programmable spatial light modulators. Our setup paves the way for future developments in the field of experimental high-dimensional QKD.

Numerical method for solution of systems of non-stationary spatially one-dimensional nonlinear differential equations

NASA Technical Reports Server (NTRS)

Morozov, S. K.; Krasitskiy, O. P.

1978-01-01

A computational scheme and a standard program is proposed for solving systems of nonstationary spatially one-dimensional nonlinear differential equations using Newton's method. The proposed scheme is universal in its applicability and its reduces to a minimum the work of programming. The program is written in the FORTRAN language and can be used without change on electronic computers of type YeS and BESM-6. The standard program described permits the identification of nonstationary (or stationary) solutions to systems of spatially one-dimensional nonlinear (or linear) partial differential equations. The proposed method may be used to solve a series of geophysical problems which take chemical reactions, diffusion, and heat conductivity into account, to evaluate nonstationary thermal fields in two-dimensional structures when in one of the geometrical directions it can take a small number of discrete levels, and to solve problems in nonstationary gas dynamics.

Quantum key distribution session with 16-dimensional photonic states

PubMed Central

Etcheverry, S.; Cañas, G.; Gómez, E. S.; Nogueira, W. A. T.; Saavedra, C.; Xavier, G. B.; Lima, G.

2013-01-01

The secure transfer of information is an important problem in modern telecommunications. Quantum key distribution (QKD) provides a solution to this problem by using individual quantum systems to generate correlated bits between remote parties, that can be used to extract a secret key. QKD with D-dimensional quantum channels provides security advantages that grow with increasing D. However, the vast majority of QKD implementations has been restricted to two dimensions. Here we demonstrate the feasibility of using higher dimensions for real-world quantum cryptography by performing, for the first time, a fully automated QKD session based on the BB84 protocol with 16-dimensional quantum states. Information is encoded in the single-photon transverse momentum and the required states are dynamically generated with programmable spatial light modulators. Our setup paves the way for future developments in the field of experimental high-dimensional QKD. PMID:23897033

Ceramic component reliability with the restructured NASA/CARES computer program

NASA Technical Reports Server (NTRS)

Powers, Lynn M.; Starlinger, Alois; Gyekenyesi, John P.

1992-01-01

The Ceramics Analysis and Reliability Evaluation of Structures (CARES) integrated design program on statistical fast fracture reliability and monolithic ceramic components is enhanced to include the use of a neutral data base, two-dimensional modeling, and variable problem size. The data base allows for the efficient transfer of element stresses, temperatures, and volumes/areas from the finite element output to the reliability analysis program. Elements are divided to insure a direct correspondence between the subelements and the Gaussian integration points. Two-dimensional modeling is accomplished by assessing the volume flaw reliability with shell elements. To demonstrate the improvements in the algorithm, example problems are selected from a round-robin conducted by WELFEP (WEakest Link failure probability prediction by Finite Element Postprocessors).

Stiffener-skin interactions in pressure-loaded composite panels

NASA Technical Reports Server (NTRS)

Loup, D. C.; Hyer, M. W.; Starnes, J. H., Jr.

1986-01-01

The effects of flange thickness, web height, and skin stiffness on the strain distributions in the skin-stiffener interface region of pressure-loaded graphite-epoxy panels, stiffened by the type-T stiffener, were examined at pressure levels up to one atmosphere. The results indicate that at these pressures geometric nonlinearities are important, and that the overall stiffener stiffness has a significant effect on panel response, particularly on the out-of-plane deformation or pillowing of the skin. The strain gradients indicated that the interface between the skin and the stiffener experiences two components of shear stress, in addition to a normal (peel) stress. Thus, the skin-stiffener interface problem is a three-dimensional problem rather than a two-dimensional one, as is often assumed.

SEMI-ANALYTIC CALCULATION OF THE TEMPERATURE DISTRIBUTION IN A PERFORATED CIRCLE

DOE Office of Scientific and Technical Information (OSTI.GOV)

Kennedy, J.M.; Fowler, J.K.

The flow of heat in a tube-in-shell fuel element is closely related to the two-dimensional heat flow in a circular region perforated by a number of circular holes. Mathematical expressions for the two-dimensional temperature distribution were obtained in terms of sources and sinks of increasing complexity located within the holes and beyond the outer circle. A computer program, TINS, which solves the temperature problem for an array of one or two rings of holes, with or without a center hole, is also described. (auth)

A Fast Estimation Algorithm for Two-Dimensional Gravity Data (GEOFAST),

DTIC Science & Technology

1979-11-15

to a wide class of problems (Refs. 9 and 17). The major inhibitor to the widespread appli- ( cation of optimal gravity data processing is the severe...extends directly to two dimensions. Define the nln 2xn1 n2 diagonal window matrix W as the Kronecker product of two one-dimensional windows W = W1 0 W2 (B...Inversion of Separable Matrices Consider the linear system y = T x (B.3-1) where T is block Toeplitz of dimension nln 2xnIn 2 . Its fre- quency domain

A Simple Algebraic Grid Adaptation Scheme with Applications to Two- and Three-dimensional Flow Problems

NASA Technical Reports Server (NTRS)

Hsu, Andrew T.; Lytle, John K.

1989-01-01

An algebraic adaptive grid scheme based on the concept of arc equidistribution is presented. The scheme locally adjusts the grid density based on gradients of selected flow variables from either finite difference or finite volume calculations. A user-prescribed grid stretching can be specified such that control of the grid spacing can be maintained in areas of known flowfield behavior. For example, the grid can be clustered near a wall for boundary layer resolution and made coarse near the outer boundary of an external flow. A grid smoothing technique is incorporated into the adaptive grid routine, which is found to be more robust and efficient than the weight function filtering technique employed by other researchers. Since the present algebraic scheme requires no iteration or solution of differential equations, the computer time needed for grid adaptation is trivial, making the scheme useful for three-dimensional flow problems. Applications to two- and three-dimensional flow problems show that a considerable improvement in flowfield resolution can be achieved by using the proposed adaptive grid scheme. Although the scheme was developed with steady flow in mind, it is a good candidate for unsteady flow computations because of its efficiency.

"The Poison That Ruined the Nation": Native American Men-Alcohol, Identity, and Traditional Healing.

PubMed

Matamonasa-Bennett, Arieahn

2017-07-01

Alcoholism and destructive drinking patterns are serious social problems in many Native American reservation and urban communities. This qualitative study of men from a single Great Lakes reservation community examined the social, cultural, and psychological aspects of their alcohol problems through their life stories. The men were in various stages of recovery and sobriety, and data collection consisted of open-ended interviews and analysis utilizing principles and techniques from grounded theory and ethnographic content analysis. Alcoholism and other serious social problems facing Native American communities need to be understood in the sociocultural and historical contexts of colonization and historical grief and trauma. This study suggests that for Native American men, there are culturally specific perspectives on alcohol that have important implications for prevention and treatment of alcohol abuse. The participants' narratives provided insight into the ways reconnecting with traditional cultural values (retraditionalization) helped them achieve sobriety. For these men, alcohol was highly symbolic of colonization as well as a protest to it. Alcohol was a means for affirming "Indian" identity and sobriety a means for reaffirming traditional tribal identity. Their narratives suggested the ways in which elements of traditional cultural values and practices facilitate healing in syncretic models and Nativized treatment. Understanding the ways in which specific Native cultural groups perceive their problems with drinking and sobriety can create more culturally congruent, culturally sensitive, and effective treatment approaches and inform future research.

Hyperspectral Super-Resolution of Locally Low Rank Images From Complementary Multisource Data.

PubMed

Veganzones, Miguel A; Simoes, Miguel; Licciardi, Giorgio; Yokoya, Naoto; Bioucas-Dias, Jose M; Chanussot, Jocelyn

2016-01-01

Remote sensing hyperspectral images (HSIs) are quite often low rank, in the sense that the data belong to a low dimensional subspace/manifold. This has been recently exploited for the fusion of low spatial resolution HSI with high spatial resolution multispectral images in order to obtain super-resolution HSI. Most approaches adopt an unmixing or a matrix factorization perspective. The derived methods have led to state-of-the-art results when the spectral information lies in a low-dimensional subspace/manifold. However, if the subspace/manifold dimensionality spanned by the complete data set is large, i.e., larger than the number of multispectral bands, the performance of these methods mainly decreases because the underlying sparse regression problem is severely ill-posed. In this paper, we propose a local approach to cope with this difficulty. Fundamentally, we exploit the fact that real world HSIs are locally low rank, that is, pixels acquired from a given spatial neighborhood span a very low-dimensional subspace/manifold, i.e., lower or equal than the number of multispectral bands. Thus, we propose to partition the image into patches and solve the data fusion problem independently for each patch. This way, in each patch the subspace/manifold dimensionality is low enough, such that the problem is not ill-posed anymore. We propose two alternative approaches to define the hyperspectral super-resolution through local dictionary learning using endmember induction algorithms. We also explore two alternatives to define the local regions, using sliding windows and binary partition trees. The effectiveness of the proposed approaches is illustrated with synthetic and semi real data.

Classification of symmetry-protected phases for interacting fermions in two dimensions

NASA Astrophysics Data System (ADS)

Cheng, Meng; Bi, Zhen; You, Yi-Zhuang; Gu, Zheng-Cheng

2018-05-01

Recently, it has been established that two-dimensional bosonic symmetry-protected topological (SPT) phases with on-site unitary symmetry G can be completely classified by the group cohomology H3( G ,U (1 ) ) . Later, group supercohomology was proposed as a partial classification for SPT phases of interacting fermions. In this work, we revisit this problem based on the algebraic theory of symmetry and defects in two-dimensional topological phases. We reproduce the partial classifications given by group supercohomology, and we also show that with an additional H1(G ,Z2) structure, a complete classification of SPT phases for two-dimensional interacting fermion systems with a total symmetry group G ×Z2f is obtained. We also discuss the classification of interacting fermionic SPT phases protected by time-reversal symmetry.

Self-Paced Physics, Segments 11-14.

ERIC Educational Resources Information Center

New York Inst. of Tech., Old Westbury.

Four segments of the Self-Paced Physics Course materials are presented in this problems and solutions book for use as the third part of student course work. The subject-matter topics are related to impulses, inelastic and elastic collisions, two-dimensional collision problems, universal constant of gravitation, gravitational acceleration and…

3. Historic American Buildings Survey, Thomas T. Waterman, Photographer March, ...

Library of Congress Historic Buildings Survey, Historic Engineering Record, Historic Landscapes Survey

3. Historic American Buildings Survey, Thomas T. Waterman, Photographer March, 1936 NORTH ELEVATION. 4. Historic American Buildings Survey, Thomas T. Waterman, Photographer March, 1936 EAST ELEVATION. - Paul Hamilton House (Ruins), Russell Creek, Edisto Island, Charleston County, SC

23. Photocopy of photograph (from California State Library, Sacramento, California, ...

Library of Congress Historic Buildings Survey, Historic Engineering Record, Historic Landscapes Survey

23. Photocopy of photograph (from California State Library, Sacramento, California, C. 1909) EXTERIOR, VIEW OF SOUTH FRONT OF CONVENTO IN RUINS, C. 1909 - Mission San Francisco Solano de Sonoma, First & Spain Streets, Sonoma, Sonoma County, CA

8. Detail view of steam dome attached to top of ...

Library of Congress Historic Buildings Survey, Historic Engineering Record, Historic Landscapes Survey

8. Detail view of steam dome attached to top of Lancashire double flue boiler. - Hacienda Azucarera El Coto, Sugar Mill Ruins, .5 Mi. SW of Rt. 347 Bridge Over Guanajibo River, San German, San German Municipio, PR

Variational asymptotic modeling of composite dimensionally reducible structures

NASA Astrophysics Data System (ADS)

Yu, Wenbin

A general framework to construct accurate reduced models for composite dimensionally reducible structures (beams, plates and shells) was formulated based on two theoretical foundations: decomposition of the rotation tensor and the variational asymptotic method. Two engineering software systems, Variational Asymptotic Beam Sectional Analysis (VABS, new version) and Variational Asymptotic Plate and Shell Analysis (VAPAS), were developed. Several restrictions found in previous work on beam modeling were removed in the present effort. A general formulation of Timoshenko-like cross-sectional analysis was developed, through which the shear center coordinates and a consistent Vlasov model can be obtained. Recovery relations are given to recover the asymptotic approximations for the three-dimensional field variables. A new version of VABS has been developed, which is a much improved program in comparison to the old one. Numerous examples are given for validation. A Reissner-like model being as asymptotically correct as possible was obtained for composite plates and shells. After formulating the three-dimensional elasticity problem in intrinsic form, the variational asymptotic method was used to systematically reduce the dimensionality of the problem by taking advantage of the smallness of the thickness. The through-the-thickness analysis is solved by a one-dimensional finite element method to provide the stiffnesses as input for the two-dimensional nonlinear plate or shell analysis as well as recovery relations to approximately express the three-dimensional results. The known fact that there exists more than one theory that is asymptotically correct to a given order is adopted to cast the refined energy into a Reissner-like form. A two-dimensional nonlinear shell theory consistent with the present modeling process was developed. The engineering computer code VAPAS was developed and inserted into DYMORE to provide an efficient and accurate analysis of composite plates and shells. Numerical results are compared with the exact solutions, and the excellent agreement proves that one can use VAPAS to analyze composite plates and shells efficiently and accurately. In conclusion, rigorous modeling approaches were developed for composite beams, plates and shells within a general framework. No such consistent and general treatment is found in the literature. The associated computer programs VABS and VAPAS are envisioned to have many applications in industry.

Reconstructing high-dimensional two-photon entangled states via compressive sensing

PubMed Central

Tonolini, Francesco; Chan, Susan; Agnew, Megan; Lindsay, Alan; Leach, Jonathan

2014-01-01

Accurately establishing the state of large-scale quantum systems is an important tool in quantum information science; however, the large number of unknown parameters hinders the rapid characterisation of such states, and reconstruction procedures can become prohibitively time-consuming. Compressive sensing, a procedure for solving inverse problems by incorporating prior knowledge about the form of the solution, provides an attractive alternative to the problem of high-dimensional quantum state characterisation. Using a modified version of compressive sensing that incorporates the principles of singular value thresholding, we reconstruct the density matrix of a high-dimensional two-photon entangled system. The dimension of each photon is equal to d = 17, corresponding to a system of 83521 unknown real parameters. Accurate reconstruction is achieved with approximately 2500 measurements, only 3% of the total number of unknown parameters in the state. The algorithm we develop is fast, computationally inexpensive, and applicable to a wide range of quantum states, thus demonstrating compressive sensing as an effective technique for measuring the state of large-scale quantum systems. PMID:25306850

Three-dimensional couette flow of dusty fluid with heat transfer in the presence of magnetic field

NASA Astrophysics Data System (ADS)

Gayathri, R.; Govindarajan, A.; Sasikala, R.

2018-04-01

This paper is focused on the mathematical modelling of three-dimensional couette flow and heat transfer of a dusty fluid between two infinite horizontal parallel porous flat plates in the presence of an induced magnetic field. The problem is formulated using a continuum two-phase model and the resulting equations are solved analytically. The lower plate is stationary while the upper plate is undergoing uniform motion in its plane. These plates are, respectively subjected to transverse exponential injection and its corresponding removal by constant suction. Due to this type of injection velocity, the flow becomes three dimensional. The closed-form expressions for velocity and temperature fields of both the fluid and dust phase are obtained by solving the governing partial differentiation equations using the perturbation method. A selective set of graphical results is presented and discussed to show interesting features of the problem. It is found that the velocity profiles of both fluid and dust particles decrease due to the increase of (magnetic parameter) Hartmann number.

Three-dimensional thermocapillary flow regimes with evaporation

NASA Astrophysics Data System (ADS)

Bekezhanova, V. B.; Goncharova, O. N.

2017-10-01

A three-dimensional problem of evaporative convection in a system of the immiscible media with a common thermocapillary interface is studied. New exact solution, which is a generalization of the Ostroumov - Birikh solution of the Navier - Stokes equations in the Oberbeck - Boussinesq approximation, is presented in order to describe the joint flows of the liquid and gas - vapor mixture in an infinite channel with a rectangular cross-section. The motion occurs in the bulk force field under action of a constant longitudinal temperature gradient. The velocity components depend only on the transverse coordinates. The functions of pressure, temperature and concentration of vapor in the gas are characterized by the linear dependence on the longitudinal coordinate. In the framework of the problem statement, which takes into account diffusive mass flux through the interface and zero vapor flux at the upper boundary of the channel, the influence of the gravity and intensity of the thermal action on flow structure is studied. The original three-dimensional problem is reduced to a chain of two-dimensional problems which are solved numerically with help of modification of the method of alternating directions. Arising flows can be characterized as a translational-rotational motion, under that the symmetrical double, quadruple or sextuple vortex structures are formed. Quantity, shape and structure of the vortexes also depend on properties of the working media.

Internet Addiction: A Brief Summary of Research and Practice

PubMed Central

Cash, Hilarie; Rae, Cosette D; Steel, Ann H; Winkler, Alexander

2012-01-01

Problematic computer use is a growing social issue which is being debated worldwide. Internet Addiction Disorder (IAD) ruins lives by causing neurological complications, psychological disturbances, and social problems. Surveys in the United States and Europe have indicated alarming prevalence rates between 1.5 and 8.2% [1]. There are several reviews addressing the definition, classification, assessment, epidemiology, and co-morbidity of IAD [2-5], and some reviews [6-8] addressing the treatment of IAD. The aim of this paper is to give a preferably brief overview of research on IAD and theoretical considerations from a practical perspective based on years of daily work with clients suffering from Internet addiction. Furthermore, with this paper we intend to bring in practical experience in the debate about the eventual inclusion of IAD in the next version of the Diagnostic and Statistical Manual of Mental Disorders (DSM). PMID:23125561

Evidence is good for your health system: policy reform to remedy catastrophic and impoverishing health spending in Mexico.

PubMed

Knaul, Felicia Marie; Arreola-Ornelas, Héctor; Méndez-Carniado, Oscar; Bryson-Cahn, Chloe; Barofsky, Jeremy; Maguire, Rachel; Miranda, Martha; Sesma, Sergio

2006-11-18

Absence of financial protection in health is a recently diagnosed "disease" of health systems. The most obvious symptom is that families face economic ruin and poverty as a consequence of financing their health care. Mexico was one of the first countries to diagnose the problem, attribute it to lack of financial protection, and propose systemic therapy through health reform. In this article we assess how Mexico turned evidence on catastrophic and impoverishing health spending into a catalyst for institutional renovation through the reform that created Seguro Popular (Popular Health Insurance). We present 15-year trends on the evolution of catastrophic and impoverishing health spending, including evidence on how the situation is improving. The results of the Mexican experience suggest an important role for the organisation and financing of the health system in reducing impoverishment and protecting households during periods of individual and collective financial crisis.

[Evidence is good for your health system: policy reform to remedy catastrophic and impoverishing health spending in Mexico].

PubMed

Knaul, Felicia Marie; Arreola-Ornelas, Héctor; Méndez-Carniado, Oscar; Bryson-Cahn, Chloe; Barofsky, Jeremy; Maguire, Rachel; Miranda, Martha; Sesma, Sergio

2007-01-01

Absence of financial protection in health is a recently diagnosed "disease" of health systems. The most obvious symptom is that families face economic ruin and poverty as a consequence of financing their health care. Mexico was one of the first countries to diagnose the problem, attribute it to lack of financial protection, and propose systemic therapy through health reform. In this article we assess how Mexico turned evidence on catastrophic and impoverishing health spending into a catalyst for institutional renovation through the reform that created Seguro Popular de Salud (Popular Health Insurance). We present 15-year trends on the evolution of catastrophic and impoverishing health spending, including evidence on how the situation is improving. The results of the Mexican experience suggest an important role for the organisation and financing of the health system in reducing impoverishment and protecting households during periods of individual and collective financial crisis.

A 2-dimensional optical architecture for solving Hamiltonian path problem based on micro ring resonators

NASA Astrophysics Data System (ADS)

Shakeri, Nadim; Jalili, Saeed; Ahmadi, Vahid; Rasoulzadeh Zali, Aref; Goliaei, Sama

2015-01-01

The problem of finding the Hamiltonian path in a graph, or deciding whether a graph has a Hamiltonian path or not, is an NP-complete problem. No exact solution has been found yet, to solve this problem using polynomial amount of time and space. In this paper, we propose a two dimensional (2-D) optical architecture based on optical electronic devices such as micro ring resonators, optical circulators and MEMS based mirror (MEMS-M) to solve the Hamiltonian Path Problem, for undirected graphs in linear time. It uses a heuristic algorithm and employs n+1 different wavelengths of a light ray, to check whether a Hamiltonian path exists or not on a graph with n vertices. Then if a Hamiltonian path exists, it reports the path. The device complexity of the proposed architecture is O(n2).

Comments on the Diffusive Behavior of Two Upwind Schemes

NASA Technical Reports Server (NTRS)

Wood, William A.; Kleb, William L.

1998-01-01

The diffusive characteristics of two upwind schemes, multi-dimensional fluctuation splitting and locally one-dimensional finite volume, are compared for scalar advection-diffusion problems. Algorithms for the two schemes are developed for node-based data representation on median-dual meshes associated with unstructured triangulations in two spatial dimensions. Four model equations are considered: linear advection, non-linear advection, diffusion, and advection-diffusion. Modular coding is employed to isolate the effects of the two approaches for upwind flux evaluation, allowing for head-to-head accuracy and efficiency comparisons. Both the stability of compressive limiters and the amount of artificial diffusion generated by the schemes is found to be grid-orientation dependent, with the fluctuation splitting scheme producing less artificial diffusion than the finite volume scheme. Convergence rates are compared for the combined advection-diffusion problem, with a speedup of 2.5 seen for fluctuation splitting versus finite volume when solved on the same mesh. However, accurate solutions to problems with small diffusion coefficients can be achieved on coarser meshes using fluctuation splitting rather than finite volume, so that when comparing convergence rates to reach a given accuracy, fluctuation splitting shows a speedup of 29 over finite volume.

On the Transition from Two-Dimensional to Three-Dimensional MHD Turbulence

NASA Technical Reports Server (NTRS)

Thess, A.; Zikanov, Oleg

2004-01-01

We report a theoretical investigation of the robustness of two-dimensional inviscid MHD flows at low magnetic Reynolds numbers with respect to three-dimensional perturbations. We analyze three model problems, namely flow in the interior of a triaxial ellipsoid, an unbounded vortex with elliptical streamlines, and a vortex sheet parallel to the magnetic field. We demonstrate that motion perpendicular to the magnetic field with elliptical streamlines becomes unstable with respect to the elliptical instability once the velocity has reached a critical magnitude whose value tends to zero as the eccentricity of the streamlines becomes large. Furthermore, vortex sheets parallel to the magnetic field, which are unstable for any velocity and any magnetic field, are found to emit eddies with vorticity perpendicular to the magnetic field and with an aspect ratio proportional to N(sup 1/2). The results suggest that purely two-dimensional motion without Joule energy dissipation is a singular type of flow which does not represent the asymptotic behaviour of three-dimensional MHD turbulence in the limit of infinitely strong magnetic fields.

On the three-dimensional instability of strained vortices

NASA Technical Reports Server (NTRS)

Waleffe, Fabian

1990-01-01

The three-dimensional (3-D) instability of a two-dimensional (2-D) flow with elliptical streamlines has been proposed as a generic mechanism for the breakdown of many 2-D flows. A physical interpretation for the mechanism is presented together with an analytical treatment of the problem. It is shown that the stability of an elliptical flow is governed by an Ince equation. An analytical representation for a localized solution is given and establishes a direct link with previous computations and experiments.

The TORSED method for construction of TORT boundary sources from external DORT flux files

DOE Office of Scientific and Technical Information (OSTI.GOV)

Rhoades, W.A.

1993-08-01

The TORSED method provides a means of coupling cylindrical two-dimensional DORT fluxes or fluences to a three-dimensional TORT calculation in Cartesian geometry through construction of external boundary sources for TORT. This can be important for several reasons. The two-dimensional environment may be too large for TORT simulation. The two-dimensional environment may be truly cylindrical in nature, and thus, better treated in that geometry. It may be desired to use a single environment calculation to study numerous local perturbations. In Section I the TORSED code is described in detail and the diverse demonstration problems that accompany the code distribution are discussed.more » In Section II, an updated discussion of the VISA code is given. VISA is required to preprocess the DORT files for use in TORSED. In Section III, the references are listed.« less

A two dimensional finite difference time domain analysis of the quiet zone fields of an anechoic chamber

NASA Technical Reports Server (NTRS)

Ryan, Deirdre A.; Luebbers, Raymond J.; Nguyen, Truong X.; Kunz, Karl S.; Steich, David J.

1992-01-01

Prediction of anechoic chamber performance is a difficult problem. Electromagnetic anechoic chambers exist for a wide range of frequencies but are typically very large when measured in wavelengths. Three dimensional finite difference time domain (FDTD) modeling of anechoic chambers is possible with current computers but at frequencies lower than most chamber design frequencies. However, two dimensional FDTD (2D-FTD) modeling enables much greater detail at higher frequencies and offers significant insight into compact anechoic chamber design and performance. A major subsystem of an anechoic chamber for which computational electromagnetic analyses exist is the reflector. First, an analysis of the quiet zone fields of a low frequency anechoic chamber produced by a uniform source and a reflector in two dimensions using the FDTD method is presented. The 2D-FDTD results are compared with results from a three dimensional corrected physical optics calculation and show good agreement. Next, a directional source is substituted for the uniform radiator. Finally, a two dimensional anechoic chamber geometry, including absorbing materials, is considered, and the 2D-FDTD results for these geometries appear reasonable.

Phase retrieval from local measurements in two dimensions

NASA Astrophysics Data System (ADS)

Iwen, Mark; Preskitt, Brian; Saab, Rayan; Viswanathan, Aditya

2017-08-01

The phase retrieval problem has appeared in a multitude of applications for decades. While ad hoc solutions have existed since the early 1970s, recent developments have provided algorithms that offer promising theoretical guarantees under increasingly realistic assumptions. Motivated by ptychographic imaging, we generalize a recent result on phase retrieval of a one dimensional objective vector x ∈ ℂd to recover a two dimensional sample Q ∈ ℂd x d from phaseless measurements, using a tensor product formulation to extend the previous work.

Observation of two-dimensional Faraday waves in extremely shallow depth.

PubMed

Li, Xiaochen; Yu, Zhengyue; Liao, Shijun

2015-09-01

A family of two-dimensional Faraday waves in extremely shallow depth (1 mm to 2 mm) of absolute ethanol are observed experimentally using a Hele-Shaw cell that vibrates vertically. The same phenomena are not observed by means of water, ethanol solution, and silicone oil. These Faraday waves are quite different from the traditional ones. These phenomena are helpful to deepen and enrich our understandings about Faraday waves, and besides provide a challenging problem for computational fluid dynamics.

Multi-dimensional Fokker-Planck equation analysis using the modified finite element method

NASA Astrophysics Data System (ADS)

Náprstek, J.; Král, R.

2016-09-01

The Fokker-Planck equation (FPE) is a frequently used tool for the solution of cross probability density function (PDF) of a dynamic system response excited by a vector of random processes. FEM represents a very effective solution possibility, particularly when transition processes are investigated or a more detailed solution is needed. Actual papers deal with single degree of freedom (SDOF) systems only. So the respective FPE includes two independent space variables only. Stepping over this limit into MDOF systems a number of specific problems related to a true multi-dimensionality must be overcome. Unlike earlier studies, multi-dimensional simplex elements in any arbitrary dimension should be deployed and rectangular (multi-brick) elements abandoned. Simple closed formulae of integration in multi-dimension domain have been derived. Another specific problem represents the generation of multi-dimensional finite element mesh. Assembling of system global matrices should be subjected to newly composed algorithms due to multi-dimensionality. The system matrices are quite full and no advantages following from their sparse character can be profited from, as is commonly used in conventional FEM applications in 2D/3D problems. After verification of partial algorithms, an illustrative example dealing with a 2DOF non-linear aeroelastic system in combination with random and deterministic excitations is discussed.

SUSY’s Ladder: Reframing sequestering at Large Volume

DOE PAGES

Reece, Matthew; Xue, Wei

2016-04-07

Theories with approximate no-scale structure, such as the Large Volume Scenario, have a distinctive hierarchy of multiple mass scales in between TeV gaugino masses and the Planck scale, which we call SUSY's Ladder. This is a particular realization of Split Supersymmetry in which the same small parameter suppresses gaugino masses relative to scalar soft masses, scalar soft masses relative to the gravitino mass, and the UV cutoff or string scale relative to the Planck scale. This scenario has many phenomenologically interesting properties, and can avoid dangers including the gravitino problem, flavor problems, and the moduli-induced LSP problem that plague othermore » supersymmetric theories. We study SUSY's Ladder using a superspace formalism that makes the mysterious cancelations in previous computations manifest. This opens the possibility of a consistent effective field theory understanding of the phenomenology of these scenarios, based on power-counting in the small ratio of string to Planck scales. We also show that four-dimensional theories with approximate no-scale structure enforced by a single volume modulus arise only from two special higher-dimensional theories: five-dimensional supergravity and ten-dimensional type IIB supergravity. As a result, this gives a phenomenological argument in favor of ten dimensional ultraviolet physics which is different from standard arguments based on the consistency of superstring theory.« less

A Finite Element Projection Method for the Solution of Particle Transport Problems with Anisotropic Scattering.

DTIC Science & Technology

1984-07-01

piecewise constant energy dependence. This is a seven-dimensional problem with time dependence, three spatial and two angular or directional variables and...in extending the computer implementation of the method to time and energy dependent problems, and to solving and validating this technique on a...problems they have severe limitations. The Monte Carlo method, usually requires the use of many hours of expensive computer time , and for deep

Solution methods for one-dimensional viscoelastic problems

NASA Technical Reports Server (NTRS)

Stubstad, John M.; Simitses, George J.

1987-01-01

A recently developed differential methodology for solution of one-dimensional nonlinear viscoelastic problems is presented. Using the example of an eccentrically loaded cantilever beam-column, the results from the differential formulation are compared to results generated using a previously published integral solution technique. It is shown that the results obtained from these distinct methodologies exhibit a surprisingly high degree of correlation with one another. A discussion of the various factors affecting the numerical accuracy and rate of convergence of these two procedures is also included. Finally, the influences of some 'higher order' effects, such as straining along the centroidal axis are discussed.

Description of a highly symmetric polytope observed in Thomson's problem of charges on a hypersphere

NASA Astrophysics Data System (ADS)

Roth, J.

2007-10-01

In a recent paper, Altschuler and Pérez-Garrido [Phys. Rev. E 76, 016705 (2007)] have presented a four-dimensional polytope with 80 vertices. We demonstrate how this polytope can be derived from the regular four-dimensional 600-cell with 120 vertices if two orthogonal positive disclinations are created. Some related polytopes are also described.

Earth observations taken during STS-1 mission

NASA Image and Video Library

2009-06-24

STS001-13-443 (12-14 April 1981) --- This photograph showing much of Italy was taken with a hand-held 70mm camera from 276 kilometers above Earth as the NASA space shuttle Columbia and its crew were marking their last few hours in space on the historic first space mission utilizing a reusable vehicle. Included in the area of the frame are Golfo de Napoli, Napoli (Naples), Castellammare, Amalfi, Capri, Sorrento, Mt. Vesuvius and the ruins of Pompei. Astronauts John W. Young and Robert L. Crippen exposed eight magazines of color 70mm film during their two-and-one-third days in Earth orbit. Photo credit: NASA

Development of a linearized unsteady Euler analysis for turbomachinery blade rows

NASA Technical Reports Server (NTRS)

Verdon, Joseph M.; Montgomery, Matthew D.; Kousen, Kenneth A.

1995-01-01

A linearized unsteady aerodynamic analysis for axial-flow turbomachinery blading is described in this report. The linearization is based on the Euler equations of fluid motion and is motivated by the need for an efficient aerodynamic analysis that can be used in predicting the aeroelastic and aeroacoustic responses of blade rows. The field equations and surface conditions required for inviscid, nonlinear and linearized, unsteady aerodynamic analyses of three-dimensional flow through a single, blade row operating within a cylindrical duct, are derived. An existing numerical algorithm for determining time-accurate solutions of the nonlinear unsteady flow problem is described, and a numerical model, based upon this nonlinear flow solver, is formulated for the first-harmonic linear unsteady problem. The linearized aerodynamic and numerical models have been implemented into a first-harmonic unsteady flow code, called LINFLUX. At present this code applies only to two-dimensional flows, but an extension to three-dimensions is planned as future work. The three-dimensional aerodynamic and numerical formulations are described in this report. Numerical results for two-dimensional unsteady cascade flows, excited by prescribed blade motions and prescribed aerodynamic disturbances at inlet and exit, are also provided to illustrate the present capabilities of the LINFLUX analysis.

New evidence for the catastrophic demise of a prehistoric settlement (the Lajia Ruins) in the Guanting Basin, upper Yellow River, NW China

NASA Astrophysics Data System (ADS)

Zhao, Hui; Huang, Chun Chang; Zheng, Zixing; Hu, Ying; Zhang, Yuzhu; Guo, Yongqiang; Zhou, Qiang

2017-09-01

The Lajia Ruins in the Guanting Basin, NW China, are a product of the prehistoric Qijia Culture. Like Pompeii, they are a rare example of an archaeological site preserved by a natural disaster and are therefore important in archaeology, anthropology and geology. However, the nature of the disaster(s) responsible for the destruction of the site remains controversial. Most studies have focused on an earthquake and a red clay layer directly overlying the site and a detailed stratigraphic study of the mid-Holocene sedimentary strata combined with other intervals of red clay deposition (hence possible disasters) is lacking. We identified a mid-Holocene paleosol sequence (the Shanglajia section) at the site which contains two layers of red clay, dated to 3950 a BP and 3500 a BP, intercalated within the mid-Holocene paleosol (S0). Subsequent multi-proxy analysis indicated that the characteristics of the two red clay layers resemble those of typical Tertiary red clay deposits and the modern gully deposit at the foot of the Great Red Hills, but are distinctly different from those of the slackwater deposits of the Yellow River and the mid-Holocene paleosol. Our results suggest that, at 3950 a BP and 3500 a BP, two large-scale rainstorm-induced mudflow events, originating from the gullies to the north, flooded the Lajia area on the second terrace of the Yellow River, devastating and burying the human settlements. We infer that the intensified erosion and mass wasting were caused by human activity; in addition, natural factors such as rainstorms and earthquakes, may also have played an important role in triggering catastrophic mudflow events in the Tertiary Red Clay deposits. Overall, our results provide further insights into prehistoric man-land relationships in this environmentally sensitive region which may have implications for modern land use in this region of China and elsewhere.

Multiple-Event, Single-Photon Counting Imaging Sensor

NASA Technical Reports Server (NTRS)

Zheng, Xinyu; Cunningham, Thomas J.; Sun, Chao; Wang, Kang L.

2011-01-01

The single-photon counting imaging sensor is typically an array of silicon Geiger-mode avalanche photodiodes that are monolithically integrated with CMOS (complementary metal oxide semiconductor) readout, signal processing, and addressing circuits located in each pixel and the peripheral area of the chip. The major problem is its single-event method for photon count number registration. A single-event single-photon counting imaging array only allows registration of up to one photon count in each of its pixels during a frame time, i.e., the interval between two successive pixel reset operations. Since the frame time can t be too short, this will lead to very low dynamic range and make the sensor merely useful for very low flux environments. The second problem of the prior technique is a limited fill factor resulting from consumption of chip area by the monolithically integrated CMOS readout in pixels. The resulting low photon collection efficiency will substantially ruin any benefit gained from the very sensitive single-photon counting detection. The single-photon counting imaging sensor developed in this work has a novel multiple-event architecture, which allows each of its pixels to register as more than one million (or more) photon-counting events during a frame time. Because of a consequently boosted dynamic range, the imaging array of the invention is capable of performing single-photon counting under ultra-low light through high-flux environments. On the other hand, since the multiple-event architecture is implemented in a hybrid structure, back-illumination and close-to-unity fill factor can be realized, and maximized quantum efficiency can also be achieved in the detector array.

OPTOTRAK: at last a system with resolution of 10 μm (Abstract Only)

NASA Astrophysics Data System (ADS)

Crouch, David G.; Kehl, L.; Krist, J. R.

1990-08-01

Northern Digital's first active marker point measurement system, the WATSMART, was begun in 1983. Development ended in 1985 with the manufacture of a highly accurate system, which achieved .15 to .25 mm accuracies in three dimensions within a .75-meter cube. Further improvements in accuracy were rendered meaningless, and a great obstacle to usability was presented by a surplus light problem somewhat incorrectly known as "the reflection problem". In 1985, development of a new system to overcome "the reflection problem" was begun. The advantages and disadvantages involved in the use of active versus passive markers were considered. The implications of using a CCD device as the imaging element in a precision measurement device were analyzed, as were device characteristics such as dynamic range, peak readout noise and charge transfer efficiency. A new type of lens was also designed The end result, in 1988, was the first OPTOTRAK system. This system produces three-dimensional data in real-time and is not at all affected by reflections. Accuracies of 30 microns have been achieved in a 1-meter volume. Each two-dimensional camera actually has two separate, one-dimensional, CCD elements and two separate anamorphic lenses. It can locate a point from 1-8 meters away with a resolution of 1 part in 64,000 and an accuracy of 1 part in 20,000 over the field of view.

Propagation in and scattering from a matched metamaterial having a zero index of refraction.

PubMed

Ziolkowski, Richard W

2004-10-01

Planar metamaterials that exhibit a zero index of refraction have been realized experimentally by several research groups. Their existence stimulated the present investigation, which details the properties of a passive, dispersive metamaterial that is matched to free space and has an index of refraction equal to zero. Thus, unlike previous zero-index investigations, both the permittivity and permeability are zero here at a specified frequency. One-, two-, and three-dimensional source problems are treated analytically. The one- and two-dimensional source problem results are confirmed numerically with finite difference time domain (FDTD) simulations. The FDTD simulator is also used to treat the corresponding one- and two-dimensional scattering problems. It is shown that in both the source and scattering configurations the electromagnetic fields in a matched zero-index medium take on a static character in space, yet remain dynamic in time, in such a manner that the underlying physics remains associated with propagating fields. Zero phase variation at various points in the zero-index medium is demonstrated once steady-state conditions are obtained. These behaviors are used to illustrate why a zero-index metamaterial, such as a zero-index electromagnetic band-gap structured medium, significantly narrows the far-field pattern associated with an antenna located within it. They are also used to show how a matched zero-index slab could be used to transform curved wave fronts into planar ones.

AQMAN; linear and quadratic programming matrix generator using two-dimensional ground-water flow simulation for aquifer management modeling

USGS Publications Warehouse

Lefkoff, L.J.; Gorelick, S.M.

1987-01-01

A FORTRAN-77 computer program code that helps solve a variety of aquifer management problems involving the control of groundwater hydraulics. It is intended for use with any standard mathematical programming package that uses Mathematical Programming System input format. The computer program creates the input files to be used by the optimization program. These files contain all the hydrologic information and management objectives needed to solve the management problem. Used in conjunction with a mathematical programming code, the computer program identifies the pumping or recharge strategy that achieves a user 's management objective while maintaining groundwater hydraulic conditions within desired limits. The objective may be linear or quadratic, and may involve the minimization of pumping and recharge rates or of variable pumping costs. The problem may contain constraints on groundwater heads, gradients, and velocities for a complex, transient hydrologic system. Linear superposition of solutions to the transient, two-dimensional groundwater flow equation is used by the computer program in conjunction with the response matrix optimization method. A unit stress is applied at each decision well and transient responses at all control locations are computed using a modified version of the U.S. Geological Survey two dimensional aquifer simulation model. The program also computes discounted cost coefficients for the objective function and accounts for transient aquifer conditions. (Author 's abstract)

Unsteady free surface flow in porous media: One-dimensional model equations including vertical effects and seepage face

NASA Astrophysics Data System (ADS)

Di Nucci, Carmine

2018-05-01

This note examines the two-dimensional unsteady isothermal free surface flow of an incompressible fluid in a non-deformable, homogeneous, isotropic, and saturated porous medium (with zero recharge and neglecting capillary effects). Coupling a Boussinesq-type model for nonlinear water waves with Darcy's law, the two-dimensional flow problem is solved using one-dimensional model equations including vertical effects and seepage face. In order to take into account the seepage face development, the system equations (given by the continuity and momentum equations) are completed by an integral relation (deduced from the Cauchy theorem). After testing the model against data sets available in the literature, some numerical simulations, concerning the unsteady flow through a rectangular dam (with an impermeable horizontal bottom), are presented and discussed.

A Navier-Stokes solution of the three-dimensional viscous compressible flow in a centrifugal compressor impeller

NASA Technical Reports Server (NTRS)

Harp, J. L., Jr.

1977-01-01

A two-dimensional time-dependent computer code was utilized to calculate the three-dimensional steady flow within the impeller blading. The numerical method is an explicit time marching scheme in two spatial dimensions. Initially, an inviscid solution is generated on the hub blade-to-blade surface by the method of Katsanis and McNally (1973). Starting with the known inviscid solution, the viscous effects are calculated through iteration. The approach makes it possible to take into account principal impeller fluid-mechanical effects. It is pointed out that the second iterate provides a complete solution to the three-dimensional, compressible, Navier-Stokes equations for flow in a centrifugal impeller. The problems investigated are related to the study of a radial impeller and a backswept impeller.

Multi-robot task allocation based on two dimensional artificial fish swarm algorithm

NASA Astrophysics Data System (ADS)

Zheng, Taixiong; Li, Xueqin; Yang, Liangyi

2007-12-01

The problem of task allocation for multiple robots is to allocate more relative-tasks to less relative-robots so as to minimize the processing time of these tasks. In order to get optimal multi-robot task allocation scheme, a twodimensional artificial swarm algorithm based approach is proposed in this paper. In this approach, the normal artificial fish is extended to be two dimension artificial fish. In the two dimension artificial fish, each vector of primary artificial fish is extended to be an m-dimensional vector. Thus, each vector can express a group of tasks. By redefining the distance between artificial fish and the center of artificial fish, the behavior of two dimension fish is designed and the task allocation algorithm based on two dimension artificial swarm algorithm is put forward. At last, the proposed algorithm is applied to the problem of multi-robot task allocation and comparer with GA and SA based algorithm is done. Simulation and compare result shows the proposed algorithm is effective.

Misuse of Prescription Pain Relievers: The Buzz Takes Your Breath Away. Permanently.

MedlinePlus

... parents, doctors, relatives, teachers, or school guidance counselors. Substance abuse ruins lives. Don't let it happen to ... is hooked on prescription pain relievers, call the substance abuse treatment 24-hour helpline: 1.800.662.HELP ...

Astronaut Neil Armstrong - Rock Sample Study - Geological Field Trip - TX

NASA Image and Video Library

1969-03-03

S69-25198 (25 Feb. 1969) --- Astronaut Neil A. Armstrong, commander of the Apollo 11 prime crew, studies rock sample during a geological field trip to the Quitman Mountains area near the Fort Quitman ruins in far west Texas.

Casa Grande Ruins National Monument acoustical monitoring 2010

DOT National Transportation Integrated Search

2014-11-01

During September 2010, The Volpe Center collected baseline acoustical data at Casa Grande National Monument (CAGR), at one site for 28 days. The baseline data collected during this period will help park managers and planners estimate the effects of f...

3. Historic American Buildings Survey Philip E. Gardner, Photographer April ...

Library of Congress Historic Buildings Survey, Historic Engineering Record, Historic Landscapes Survey

3. Historic American Buildings Survey Philip E. Gardner, Photographer April 1958 FRONT THRESHOLD FROM N. BEFORE OLD BRICK CAPPED BY NEW - Fort Frederica, Captain John Mackay House (Ruins), Lot No. 6, North Ward, Saint Simons Island, Glynn County, GA

Thermally induced rarefied gas flow in a three-dimensional enclosure with square cross-section

NASA Astrophysics Data System (ADS)

Zhu, Lianhua; Yang, Xiaofan; Guo, Zhaoli

2017-12-01

Rarefied gas flow in a three-dimensional enclosure induced by nonuniform temperature distribution is numerically investigated. The enclosure has a square channel-like geometry with alternatively heated closed ends and lateral walls with a linear temperature distribution. A recently proposed implicit discrete velocity method with a memory reduction technique is used to numerically simulate the problem based on the nonlinear Shakhov kinetic equation. The Knudsen number dependencies of the vortices pattern, slip velocity at the planar walls and edges, and heat transfer are investigated. The influences of the temperature ratio imposed at the ends of the enclosure and the geometric aspect ratio are also evaluated. The overall flow pattern shows similarities with those observed in two-dimensional configurations in literature. However, features due to the three-dimensionality are observed with vortices that are not identified in previous studies on similar two-dimensional enclosures at high Knudsen and small aspect ratios.

An adaptive front tracking technique for three-dimensional transient flows

NASA Astrophysics Data System (ADS)

Galaktionov, O. S.; Anderson, P. D.; Peters, G. W. M.; van de Vosse, F. N.

2000-01-01

An adaptive technique, based on both surface stretching and surface curvature analysis for tracking strongly deforming fluid volumes in three-dimensional flows is presented. The efficiency and accuracy of the technique are demonstrated for two- and three-dimensional flow simulations. For the two-dimensional test example, the results are compared with results obtained using a different tracking approach based on the advection of a passive scalar. Although for both techniques roughly the same structures are found, the resolution for the front tracking technique is much higher. In the three-dimensional test example, a spherical blob is tracked in a chaotic mixing flow. For this problem, the accuracy of the adaptive tracking is demonstrated by the volume conservation for the advected blob. Adaptive front tracking is suitable for simulation of the initial stages of fluid mixing, where the interfacial area can grow exponentially with time. The efficiency of the algorithm significantly benefits from parallelization of the code. Copyright

Exact solutions and conservation laws of the system of two-dimensional viscous Burgers equations

NASA Astrophysics Data System (ADS)

Abdulwahhab, Muhammad Alim

2016-10-01

Fluid turbulence is one of the phenomena that has been studied extensively for many decades. Due to its huge practical importance in fluid dynamics, various models have been developed to capture both the indispensable physical quality and the mathematical structure of turbulent fluid flow. Among the prominent equations used for gaining in-depth insight of fluid turbulence is the two-dimensional Burgers equations. Its solutions have been studied by researchers through various methods, most of which are numerical. Being a simplified form of the two-dimensional Navier-Stokes equations and its wide range of applicability in various fields of science and engineering, development of computationally efficient methods for the solution of the two-dimensional Burgers equations is still an active field of research. In this study, Lie symmetry method is used to perform detailed analysis on the system of two-dimensional Burgers equations. Optimal system of one-dimensional subalgebras up to conjugacy is derived and used to obtain distinct exact solutions. These solutions not only help in understanding the physical effects of the model problem but also, can serve as benchmarks for constructing algorithms and validation of numerical solutions of the system of Burgers equations under consideration at finite Reynolds numbers. Independent and nontrivial conserved vectors are also constructed.

Protein structure-structure alignment with discrete Fréchet distance.

PubMed

Jiang, Minghui; Xu, Ying; Zhu, Binhai

2008-02-01

Matching two geometric objects in two-dimensional (2D) and three-dimensional (3D) spaces is a central problem in computer vision, pattern recognition, and protein structure prediction. In particular, the problem of aligning two polygonal chains under translation and rotation to minimize their distance has been studied using various distance measures. It is well known that the Hausdorff distance is useful for matching two point sets, and that the Fréchet distance is a superior measure for matching two polygonal chains. The discrete Fréchet distance closely approximates the (continuous) Fréchet distance, and is a natural measure for the geometric similarity of the folded 3D structures of biomolecules such as proteins. In this paper, we present new algorithms for matching two polygonal chains in two dimensions to minimize their discrete Fréchet distance under translation and rotation, and an effective heuristic for matching two polygonal chains in three dimensions. We also describe our empirical results on the application of the discrete Fréchet distance to protein structure-structure alignment.

Eigenvalue problems for Beltrami fields arising in a three-dimensional toroidal magnetohydrodynamic equilibrium problem

DOE Office of Scientific and Technical Information (OSTI.GOV)

Hudson, S. R.; Hole, M. J.; Dewar, R. L.

2007-05-15

A generalized energy principle for finite-pressure, toroidal magnetohydrodynamic (MHD) equilibria in general three-dimensional configurations is proposed. The full set of ideal-MHD constraints is applied only on a discrete set of toroidal magnetic surfaces (invariant tori), which act as barriers against leakage of magnetic flux, helicity, and pressure through chaotic field-line transport. It is argued that a necessary condition for such invariant tori to exist is that they have fixed, irrational rotational transforms. In the toroidal domains bounded by these surfaces, full Taylor relaxation is assumed, thus leading to Beltrami fields {nabla}xB={lambda}B, where {lambda} is constant within each domain. Two distinctmore » eigenvalue problems for {lambda} arise in this formulation, depending on whether fluxes and helicity are fixed, or boundary rotational transforms. These are studied in cylindrical geometry and in a three-dimensional toroidal region of annular cross section. In the latter case, an application of a residue criterion is used to determine the threshold for connected chaos.« less

Dynamics and Stability of Rolling Viscoelastic Tires

DOE Office of Scientific and Technical Information (OSTI.GOV)

Potter, Trevor

2013-04-30

Current steady state rolling tire calculations often do not include treads because treads destroy the rotational symmetry of the tire. We describe two methodologies to compute time periodic solutions of a two-dimensional viscoelastic tire with treads: solving a minimization problem and solving a system of equations. We also expand on work by Oden and Lin on free spinning rolling elastic tires in which they disovered a hierachy of N-peak steady state standing wave solutions. In addition to discovering a two-dimensional hierarchy of standing wave solutions that includes their N-peak hiearchy, we consider the eects of viscoelasticity on the standing wavemore » solutions. Finally, a commonplace model of viscoelasticity used in our numerical experiments led to non-physical elastic energy growth for large tire speeds. We show that a viscoelastic model of Govindjee and Reese remedies the problem.« less

Two-dimensional Anderson-Hubbard model in the DMFT + {Sigma} approximation

DOE Office of Scientific and Technical Information (OSTI.GOV)

Kuchinskii, E. Z., E-mail: kuchinsk@iep.uran.ru; Kuleeva, N. A.; Nekrasov, I. A.

The density of states, the dynamic (optical) conductivity, and the phase diagram of the paramagnetic two-dimensional Anderson-Hubbard model with strong correlations and disorder are analyzed within the generalized dynamical mean field theory (DMFT + {Sigma} approximation). Strong correlations are accounted by the DMFT, while disorder is taken into account via the appropriate generalization of the self-consistent theory of localization. We consider the two-dimensional system with the rectangular 'bare' density of states (DOS). The DMFT effective single-impurity problem is solved by numerical renormalization group (NRG). The 'correlated metal,' Mott insulator, and correlated Anderson insulator phases are identified from the evolution ofmore » the density of states, optical conductivity, and localization length, demonstrating both Mott-Hubbard and Anderson metal-insulator transitions in two-dimensional systems of finite size, allowing us to construct the complete zero-temperature phase diagram of the paramagnetic Anderson-Hubbard model. The localization length in our approximation is practically independent of the strength of Hubbard correlations. But the divergence of the localization length in a finite-size two-dimensional system at small disorder signifies the existence of an effective Anderson transition.« less

Secondary School Students' Reasoning about Conditional Probability, Samples, and Sampling Procedures

ERIC Educational Resources Information Center

Prodromou, Theodosia

2016-01-01

In the Australian mathematics curriculum, Year 12 students (aged 16-17) are asked to solve conditional probability problems that involve the representation of the problem situation with two-way tables or three-dimensional diagrams and consider sampling procedures that result in different correct answers. In a small exploratory study, we…

Two dimensional finite element heat transfer models for softwood

Treesearch

Hongmei Gu; John F. Hunt

2004-01-01

The anisotropy of wood creates a complex problem for solving heat and mass transfer problems that require analyses be based on fundamental material properties of the wood structure. Most heat transfer models use average thermal properties across either the radial or tangential directions and have not differentiated the effects of cellular alignment, earlywood/latewood...

Two-Dimensional Crystallography Introduced by the Sprinkler Watering Problem

ERIC Educational Resources Information Center

De Toro, Jose A.; Calvo, Gabriel F.; Muniz, Pablo

2012-01-01

The problem of optimizing the number of circular sprinklers watering large fields is used to introduce, from a purely elementary geometrical perspective, some basic concepts in crystallography and comment on a few size effects in condensed matter physics. We examine square and hexagonal lattices to build a function describing the, so-called, dry…

Two fast approximate wavelet algorithms for image processing, classification, and recognition

NASA Astrophysics Data System (ADS)

Wickerhauser, Mladen V.

1994-07-01

We use large libraries of template waveforms with remarkable orthogonality properties to recast the relatively complex principal orthogonal decomposition (POD) into an optimization problem with a fast solution algorithm. Then it becomes practical to use POD to solve two related problems: recognizing or classifying images, and inverting a complicated map from a low-dimensional configuration space to a high-dimensional measurement space. In the case where the number N of pixels or measurements is more than 1000 or so, the classical O(N3) POD algorithms becomes very costly, but it can be replaced with an approximate best-basis method that has complexity O(N2logN). A variation of POD can also be used to compute an approximate Jacobian for the complicated map.

Flow simulations about steady-complex and unsteady moving configurations using structured-overlapped and unstructured grids

NASA Technical Reports Server (NTRS)

Newman, James C., III

1995-01-01

The limiting factor in simulating flows past realistic configurations of interest has been the discretization of the physical domain on which the governing equations of fluid flow may be solved. In an attempt to circumvent this problem, many Computational Fluid Dynamic (CFD) methodologies that are based on different grid generation and domain decomposition techniques have been developed. However, due to the costs involved and expertise required, very few comparative studies between these methods have been performed. In the present work, the two CFD methodologies which show the most promise for treating complex three-dimensional configurations as well as unsteady moving boundary problems are evaluated. These are namely the structured-overlapped and the unstructured grid schemes. Both methods use a cell centered, finite volume, upwind approach. The structured-overlapped algorithm uses an approximately factored, alternating direction implicit scheme to perform the time integration, whereas, the unstructured algorithm uses an explicit Runge-Kutta method. To examine the accuracy, efficiency, and limitations of each scheme, they are applied to the same steady complex multicomponent configurations and unsteady moving boundary problems. The steady complex cases consist of computing the subsonic flow about a two-dimensional high-lift multielement airfoil and the transonic flow about a three-dimensional wing/pylon/finned store assembly. The unsteady moving boundary problems are a forced pitching oscillation of an airfoil in a transonic freestream and a two-dimensional, subsonic airfoil/store separation sequence. Accuracy was accessed through the comparison of computed and experimentally measured pressure coefficient data on several of the wing/pylon/finned store assembly's components and at numerous angles-of-attack for the pitching airfoil. From this study, it was found that both the structured-overlapped and the unstructured grid schemes yielded flow solutions of comparable accuracy for these simulations. This study also indicated that, overall, the structured-overlapped scheme was slightly more CPU efficient than the unstructured approach.

Engineering two-photon high-dimensional states through quantum interference

PubMed Central

Zhang, Yingwen; Roux, Filippus S.; Konrad, Thomas; Agnew, Megan; Leach, Jonathan; Forbes, Andrew

2016-01-01

Many protocols in quantum science, for example, linear optical quantum computing, require access to large-scale entangled quantum states. Such systems can be realized through many-particle qubits, but this approach often suffers from scalability problems. An alternative strategy is to consider a lesser number of particles that exist in high-dimensional states. The spatial modes of light are one such candidate that provides access to high-dimensional quantum states, and thus they increase the storage and processing potential of quantum information systems. We demonstrate the controlled engineering of two-photon high-dimensional states entangled in their orbital angular momentum through Hong-Ou-Mandel interference. We prepare a large range of high-dimensional entangled states and implement precise quantum state filtering. We characterize the full quantum state before and after the filter, and are thus able to determine that only the antisymmetric component of the initial state remains. This work paves the way for high-dimensional processing and communication of multiphoton quantum states, for example, in teleportation beyond qubits. PMID:26933685

Multigrid one shot methods for optimal control problems: Infinite dimensional control

NASA Technical Reports Server (NTRS)

Arian, Eyal; Taasan, Shlomo

1994-01-01

The multigrid one shot method for optimal control problems, governed by elliptic systems, is introduced for the infinite dimensional control space. ln this case, the control variable is a function whose discrete representation involves_an increasing number of variables with grid refinement. The minimization algorithm uses Lagrange multipliers to calculate sensitivity gradients. A preconditioned gradient descent algorithm is accelerated by a set of coarse grids. It optimizes for different scales in the representation of the control variable on different discretization levels. An analysis which reduces the problem to the boundary is introduced. It is used to approximate the two level asymptotic convergence rate, to determine the amplitude of the minimization steps, and the choice of a high pass filter to be used when necessary. The effectiveness of the method is demonstrated on a series of test problems. The new method enables the solutions of optimal control problems at the same cost of solving the corresponding analysis problems just a few times.

High-resolution Self-Organizing Maps for advanced visualization and dimension reduction.

PubMed

Saraswati, Ayu; Nguyen, Van Tuc; Hagenbuchner, Markus; Tsoi, Ah Chung

2018-05-04

Kohonen's Self Organizing feature Map (SOM) provides an effective way to project high dimensional input features onto a low dimensional display space while preserving the topological relationships among the input features. Recent advances in algorithms that take advantages of modern computing hardware introduced the concept of high resolution SOMs (HRSOMs). This paper investigates the capabilities and applicability of the HRSOM as a visualization tool for cluster analysis and its suitabilities to serve as a pre-processor in ensemble learning models. The evaluation is conducted on a number of established benchmarks and real-world learning problems, namely, the policeman benchmark, two web spam detection problems, a network intrusion detection problem, and a malware detection problem. It is found that the visualization resulted from an HRSOM provides new insights concerning these learning problems. It is furthermore shown empirically that broad benefits from the use of HRSOMs in both clustering and classification problems can be expected. Copyright © 2018 Elsevier Ltd. All rights reserved.

Computation of three-dimensional shock wave and boundary-layer interactions

NASA Technical Reports Server (NTRS)

Hung, C. M.

1985-01-01

Computations of the impingement of an oblique shock wave on a cylinder and a supersonic flow past a blunt fin mounted on a plate are used to study three dimensional shock wave and boundary layer interaction. In the impingement case, the problem of imposing a planar impinging shock as an outer boundary condition is discussed and the details of particle traces in windward and leeward symmetry planes and near the body surface are presented. In the blunt fin case, differences between two dimensional and three dimensional separation are discussed, and the existence of an unique high speed, low pressure region under the separated spiral vortex core is demonstrated. The accessibility of three dimensional separation is discussed.

One-dimensional Gromov minimal filling problem

NASA Astrophysics Data System (ADS)

Ivanov, Alexandr O.; Tuzhilin, Alexey A.

2012-05-01

The paper is devoted to a new branch in the theory of one-dimensional variational problems with branching extremals, the investigation of one-dimensional minimal fillings introduced by the authors. On the one hand, this problem is a one-dimensional version of a generalization of Gromov's minimal fillings problem to the case of stratified manifolds. On the other hand, this problem is interesting in itself and also can be considered as a generalization of another classical problem, the Steiner problem on the construction of a shortest network connecting a given set of terminals. Besides the statement of the problem, we discuss several properties of the minimal fillings and state several conjectures. Bibliography: 38 titles.

Triangular node for Transmission-Line Modeling (TLM) applied to bio-heat transfer.

PubMed

Milan, Hugo F M; Gebremedhin, Kifle G

2016-12-01

Transmission-Line Modeling (TLM) is a numerical method used to solve complex and time-domain bio-heat transfer problems. In TLM, rectangles are used to discretize two-dimensional problems. The drawback in using rectangular shapes is that instead of refining only the domain of interest, a large additional domain will also be refined in the x and y axes, which results in increased computational time and memory space. In this paper, we developed a triangular node for TLM applied to bio-heat transfer that does not have the drawback associated with the rectangular nodes. The model includes heat source, blood perfusion (advection), boundary conditions and initial conditions. The boundary conditions could be adiabatic, temperature, heat flux, or convection. A matrix equation for TLM, which simplifies the solution of time-domain problems or solves steady-state problems, was also developed. The predicted results were compared against results obtained from the solution of a simplified two-dimensional problem, and they agreed within 1% for a mesh length of triangular faces of 59µm±9µm (mean±standard deviation) and a time step of 1ms. Copyright © 2016 Elsevier Ltd. All rights reserved.

An investigation on a two-dimensional problem of Mode-I crack in a thermoelastic medium

NASA Astrophysics Data System (ADS)

Kant, Shashi; Gupta, Manushi; Shivay, Om Namha; Mukhopadhyay, Santwana

2018-04-01

In this work, we consider a two-dimensional dynamical problem of an infinite space with finite linear Mode-I crack and employ a recently proposed heat conduction model: an exact heat conduction with a single delay term. The thermoelastic medium is taken to be homogeneous and isotropic. However, the boundary of the crack is subjected to a prescribed temperature and stress distributions. The Fourier and Laplace transform techniques are used to solve the problem. Mathematical modeling of the present problem reduces the solution of the problem into the solution of a system of four dual integral equations. The solution of these equations is equivalent to the solution of the Fredholm's integral equation of the first kind which has been solved by using the regularization method. Inverse Laplace transform is carried out by using the Bellman method, and we obtain the numerical solution for all the physical field variables in the physical domain. Results are shown graphically, and we highlight the effects of the presence of crack in the behavior of thermoelastic interactions inside the medium in the present context, and its results are compared with the results of the thermoelasticity of type-III.

Heat transfer of phase-change materials in two-dimensional cylindrical coordinates

NASA Technical Reports Server (NTRS)

Labdon, M. B.; Guceri, S. I.

1981-01-01

Two-dimensional phase-change problem is numerically solved in cylindrical coordinates (r and z) by utilizing two Taylor series expansions for the temperature distributions in the neighborhood of the interface location. These two expansions form two polynomials in r and z directions. For the regions sufficiently away from the interface the temperature field equations are numerically solved in the usual way and the results are coupled with the polynomials. The main advantages of this efficient approach include ability to accept arbitrarily time dependent boundary conditions of all types and arbitrarily specified initial temperature distributions. A modified approach using a single Taylor series expansion in two variables is also suggested.

A general method for constructing multidimensional molecular potential energy surfaces from {ital ab} {ital initio} calculations

DOE Office of Scientific and Technical Information (OSTI.GOV)

Ho, T.; Rabitz, H.

1996-02-01

A general interpolation method for constructing smooth molecular potential energy surfaces (PES{close_quote}s) from {ital ab} {ital initio} data are proposed within the framework of the reproducing kernel Hilbert space and the inverse problem theory. The general expression for an {ital a} {ital posteriori} error bound of the constructed PES is derived. It is shown that the method yields globally smooth potential energy surfaces that are continuous and possess derivatives up to second order or higher. Moreover, the method is amenable to correct symmetry properties and asymptotic behavior of the molecular system. Finally, the method is generic and can be easilymore » extended from low dimensional problems involving two and three atoms to high dimensional problems involving four or more atoms. Basic properties of the method are illustrated by the construction of a one-dimensional potential energy curve of the He{endash}He van der Waals dimer using the exact quantum Monte Carlo calculations of Anderson {ital et} {ital al}. [J. Chem. Phys. {bold 99}, 345 (1993)], a two-dimensional potential energy surface of the HeCO van der Waals molecule using recent {ital ab} {ital initio} calculations by Tao {ital et} {ital al}. [J. Chem. Phys. {bold 101}, 8680 (1994)], and a three-dimensional potential energy surface of the H{sup +}{sub 3} molecular ion using highly accurate {ital ab} {ital initio} calculations of R{umlt o}hse {ital et} {ital al}. [J. Chem. Phys. {bold 101}, 2231 (1994)]. In the first two cases the constructed potentials clearly exhibit the correct asymptotic forms, while in the last case the constructed potential energy surface is in excellent agreement with that constructed by R{umlt o}hse {ital et} {ital al}. using a low order polynomial fitting procedure. {copyright} {ital 1996 American Institute of Physics.}« less

Economics and Migration: NAFTA’s Impact on Mexico

DTIC Science & Technology

2011-10-28

immigrants to the U.S. 1 No nation was ever ruined by trade, even seemingly the most disadvantageous . -Benjamin...NAFTA forced rural farmers to cultivate marginal land due to lost income caused by falling commodity prices, leading to extensive deforestation in

EXTERIOR VIEW, LOOKING WEST, AT FORMER TCIUS STEEL COMPANY BEE ...

Library of Congress Historic Buildings Survey, Historic Engineering Record, Historic Landscapes Survey

EXTERIOR VIEW, LOOKING WEST, AT FORMER TCI-US STEEL COMPANY BEE HIVE OVEN RUINS WITH FORESTED OVERGROWTH. - Pratt Coal & Coke Company, Pratt Mines, Coke Ovens & Railroad, Bounded by First Street, Avenue G, Third Place, Birmingham Southern Railroad, Birmingham, Jefferson County, AL

How We Snatched Our Schools from the Brink of Financial Ruin.

ERIC Educational Resources Information Center

Jones, Mike M.

1981-01-01

The Waycross (Georgia) school system retired its debt and reduced its tax rates, without reducing educational quality, by budgeting for an operating reserve, carefully reviewing programs, reducing its fiscal responsibilities for special programs, and tightening its budget management techniques. (RW)

Mending Time.

ERIC Educational Resources Information Center

Flanagan, Joseph

1998-01-01

Hispanic and Native-American seasonal Park Service employees in the Southwest are training younger local workers in masonry preservation in a National Park Service effort to establish a ruins preservation program. These new technicians will also learn to document, photograph, and report on preservation activities and follow legal administrative…

MayaQuest: A Student-Directed Expedition.

ERIC Educational Resources Information Center

Hefte, Rachel

1995-01-01

Describes an educational project linking classrooms using telecommunications with a four-person bicycling team exploring Mayan ruins in Central America. Provides a historical overview of the Mayan civilization. Includes suggested activities and provides information on how to obtain lesson plans on the project. (CFR)

Advancing Humanity through Honors

ERIC Educational Resources Information Center

Rao, Michael

2015-01-01

In February of 1967, newly elected Governor Ronald Reagan delivered a speech about California's austere budget in which he chided the "intellectual luxury" of higher education as something that could lead to "economic ruin" if unchecked. "Taxpayers shouldn't be subsidizing intellectual curiosity," Gov. Reagan decreed.…

1. PERSPECTIVE VIEW, FROM THE NORTHEAST ATOP EAST WING (FOREGROUND), ...

Library of Congress Historic Buildings Survey, Historic Engineering Record, Historic Landscapes Survey

1. PERSPECTIVE VIEW, FROM THE NORTHEAST ATOP EAST WING (FOREGROUND), OF GREAT KIVA (RECONSTUCTED IN BACKGROUND) AND A SMALLER KIVA (MIDDLE GROUND) - Aztec Ruins, Great Kiva, New Mexico 44 near junction of U.S. 550, Aztec, San Juan County, NM

7. Historic American Buildings Survey Photographed by Henry F. Withey, ...

Library of Congress Historic Buildings Survey, Historic Engineering Record, Historic Landscapes Survey

7. Historic American Buildings Survey Photographed by Henry F. Withey, MaY 1937 SHOP RUINS (DETAIL OF SOAP OVEN) - Mission San Juan Capistrano, Industrial Shops, Olive Street, between U.S. Highway 101 & Main Street, San Juan Capistrano, Orange County, CA

2. Historic American Buildings Survey Copied by E. W. Russell, ...

Library of Congress Historic Buildings Survey, Historic Engineering Record, Historic Landscapes Survey

2. Historic American Buildings Survey Copied by E. W. Russell, Photographer, August 31, 1936 75TH ANNIVERSARY YEARBOOK (NOT COPYRIGHT) - RUINS OF OLD COLLEGE AFTER FIRE OF 1869 - Spring Hill College, Original Building, Old Shell Road, Spring Hill, Mobile County, AL

1. Dyea Dock looking south. Note faint evenly spaced circular ...

Library of Congress Historic Buildings Survey, Historic Engineering Record, Historic Landscapes Survey

1. Dyea Dock looking south. Note faint evenly spaced circular dark pieces of grass up through the middle of the picture indicating posts making up the pier. Photograph made from park service cherry picker. - Dyea Dock & Association (Ruins), Skagway, Skagway, AK

Multi-dimensional rheology-based two-phase model for sediment transport and applications to sheet flow and pipeline scour

DOE Office of Scientific and Technical Information (OSTI.GOV)

Lee, Cheng-Hsien; Department of Water Resources and Environmental Engineering, Tamkang University, New Taipei City 25137, Taiwan; Low, Ying Min, E-mail: ceelowym@nus.edu.sg

2016-05-15

Sediment transport is fundamentally a two-phase phenomenon involving fluid and sediments; however, many existing numerical models are one-phase approaches, which are unable to capture the complex fluid-particle and inter-particle interactions. In the last decade, two-phase models have gained traction; however, there are still many limitations in these models. For example, several existing two-phase models are confined to one-dimensional problems; in addition, the existing two-dimensional models simulate only the region outside the sand bed. This paper develops a new three-dimensional two-phase model for simulating sediment transport in the sheet flow condition, incorporating recently published rheological characteristics of sediments. The enduring-contact, inertial,more » and fluid viscosity effects are considered in determining sediment pressure and stresses, enabling the model to be applicable to a wide range of particle Reynolds number. A k − ε turbulence model is adopted to compute the Reynolds stresses. In addition, a novel numerical scheme is proposed, thus avoiding numerical instability caused by high sediment concentration and allowing the sediment dynamics to be computed both within and outside the sand bed. The present model is applied to two classical problems, namely, sheet flow and scour under a pipeline with favorable results. For sheet flow, the computed velocity is consistent with measured data reported in the literature. For pipeline scour, the computed scour rate beneath the pipeline agrees with previous experimental observations. However, the present model is unable to capture vortex shedding; consequently, the sediment deposition behind the pipeline is overestimated. Sensitivity analyses reveal that model parameters associated with turbulence have strong influence on the computed results.« less

Two-dimensional character of internal rotation of furfural and other five-member heterocyclic aromatic aldehydes

NASA Astrophysics Data System (ADS)

Bataev, Vadim A.; Pupyshev, Vladimir I.; Godunov, Igor A.

2016-05-01

The features of nuclear motion corresponding to the rotation of the formyl group (CHO) are studied for the molecules of furfural and some other five-member heterocyclic aromatic aldehydes by the use of MP2/6-311G** quantum chemical approximation. It is demonstrated that the traditional one-dimensional models of internal rotation for the molecules studied have only limited applicability. The reason is the strong kinematic interaction of the rotation of the CHO group and out-of-plane CHO deformation that is realized for the molecules under consideration. The computational procedure based on the two-dimensional approximation is considered for low lying vibrational states as more adequate to the problem.

Face recognition based on two-dimensional discriminant sparse preserving projection

NASA Astrophysics Data System (ADS)

Zhang, Dawei; Zhu, Shanan

2018-04-01

In this paper, a supervised dimensionality reduction algorithm named two-dimensional discriminant sparse preserving projection (2DDSPP) is proposed for face recognition. In order to accurately model manifold structure of data, 2DDSPP constructs within-class affinity graph and between-class affinity graph by the constrained least squares (LS) and l1 norm minimization problem, respectively. Based on directly operating on image matrix, 2DDSPP integrates graph embedding (GE) with Fisher criterion. The obtained projection subspace preserves within-class neighborhood geometry structure of samples, while keeping away samples from different classes. The experimental results on the PIE and AR face databases show that 2DDSPP can achieve better recognition performance.

Finite element meshing approached as a global minimization process

DOE Office of Scientific and Technical Information (OSTI.GOV)

WITKOWSKI,WALTER R.; JUNG,JOSEPH; DOHRMANN,CLARK R.

2000-03-01

The ability to generate a suitable finite element mesh in an automatic fashion is becoming the key to being able to automate the entire engineering analysis process. However, placing an all-hexahedron mesh in a general three-dimensional body continues to be an elusive goal. The approach investigated in this research is fundamentally different from any other that is known of by the authors. A physical analogy viewpoint is used to formulate the actual meshing problem which constructs a global mathematical description of the problem. The analogy used was that of minimizing the electrical potential of a system charged particles within amore » charged domain. The particles in the presented analogy represent duals to mesh elements (i.e., quads or hexes). Particle movement is governed by a mathematical functional which accounts for inter-particles repulsive, attractive and alignment forces. This functional is minimized to find the optimal location and orientation of each particle. After the particles are connected a mesh can be easily resolved. The mathematical description for this problem is as easy to formulate in three-dimensions as it is in two- or one-dimensions. The meshing algorithm was developed within CoMeT. It can solve the two-dimensional meshing problem for convex and concave geometries in a purely automated fashion. Investigation of the robustness of the technique has shown a success rate of approximately 99% for the two-dimensional geometries tested. Run times to mesh a 100 element complex geometry were typically in the 10 minute range. Efficiency of the technique is still an issue that needs to be addressed. Performance is an issue that is critical for most engineers generating meshes. It was not for this project. The primary focus of this work was to investigate and evaluate a meshing algorithm/philosophy with efficiency issues being secondary. The algorithm was also extended to mesh three-dimensional geometries. Unfortunately, only simple geometries were tested before this project ended. The primary complexity in the extension was in the connectivity problem formulation. Defining all of the interparticle interactions that occur in three-dimensions and expressing them in mathematical relationships is very difficult.« less

Augmented design and analysis of computer experiments: a novel tolerance embedded global optimization approach applied to SWIR hyperspectral illumination design.

PubMed

Keresztes, Janos C; John Koshel, R; D'huys, Karlien; De Ketelaere, Bart; Audenaert, Jan; Goos, Peter; Saeys, Wouter

2016-12-26

A novel meta-heuristic approach for minimizing nonlinear constrained problems is proposed, which offers tolerance information during the search for the global optimum. The method is based on the concept of design and analysis of computer experiments combined with a novel two phase design augmentation (DACEDA), which models the entire merit space using a Gaussian process, with iteratively increased resolution around the optimum. The algorithm is introduced through a series of cases studies with increasing complexity for optimizing uniformity of a short-wave infrared (SWIR) hyperspectral imaging (HSI) illumination system (IS). The method is first demonstrated for a two-dimensional problem consisting of the positioning of analytical isotropic point sources. The method is further applied to two-dimensional (2D) and five-dimensional (5D) SWIR HSI IS versions using close- and far-field measured source models applied within the non-sequential ray-tracing software FRED, including inherent stochastic noise. The proposed method is compared to other heuristic approaches such as simplex and simulated annealing (SA). It is shown that DACEDA converges towards a minimum with 1 % improvement compared to simplex and SA, and more importantly requiring only half the number of simulations. Finally, a concurrent tolerance analysis is done within DACEDA for to the five-dimensional case such that further simulations are not required.

Hypergraph-based anomaly detection of high-dimensional co-occurrences.

PubMed

Silva, Jorge; Willett, Rebecca

2009-03-01

This paper addresses the problem of detecting anomalous multivariate co-occurrences using a limited number of unlabeled training observations. A novel method based on using a hypergraph representation of the data is proposed to deal with this very high-dimensional problem. Hypergraphs constitute an important extension of graphs which allow edges to connect more than two vertices simultaneously. A variational Expectation-Maximization algorithm for detecting anomalies directly on the hypergraph domain without any feature selection or dimensionality reduction is presented. The resulting estimate can be used to calculate a measure of anomalousness based on the False Discovery Rate. The algorithm has O(np) computational complexity, where n is the number of training observations and p is the number of potential participants in each co-occurrence event. This efficiency makes the method ideally suited for very high-dimensional settings, and requires no tuning, bandwidth or regularization parameters. The proposed approach is validated on both high-dimensional synthetic data and the Enron email database, where p > 75,000, and it is shown that it can outperform other state-of-the-art methods.

What is integrability of discrete variational systems?

PubMed

Boll, Raphael; Petrera, Matteo; Suris, Yuri B

2014-02-08

We propose a notion of a pluri-Lagrangian problem, which should be understood as an analogue of multi-dimensional consistency for variational systems. This is a development along the line of research of discrete integrable Lagrangian systems initiated in 2009 by Lobb and Nijhoff, however, having its more remote roots in the theory of pluriharmonic functions, in the Z -invariant models of statistical mechanics and their quasiclassical limit, as well as in the theory of variational symmetries going back to Noether. A d -dimensional pluri-Lagrangian problem can be described as follows: given a d -form [Formula: see text] on an m -dimensional space (called multi-time, m > d ), whose coefficients depend on a sought-after function x of m independent variables (called field), find those fields x which deliver critical points to the action functionals [Formula: see text] for any d -dimensional manifold Σ in the multi-time. We derive the main building blocks of the multi-time Euler-Lagrange equations for a discrete pluri-Lagrangian problem with d =2, the so-called corner equations, and discuss the notion of consistency of the system of corner equations. We analyse the system of corner equations for a special class of three-point two-forms, corresponding to integrable quad-equations of the ABS list. This allows us to close a conceptual gap of the work by Lobb and Nijhoff by showing that the corresponding two-forms are closed not only on solutions of (non-variational) quad-equations, but also on general solutions of the corresponding corner equations. We also find an example of a pluri-Lagrangian system not coming from a multi-dimensionally consistent system of quad-equations.

What is integrability of discrete variational systems?

PubMed Central

Boll, Raphael; Petrera, Matteo; Suris, Yuri B.

2014-01-01

We propose a notion of a pluri-Lagrangian problem, which should be understood as an analogue of multi-dimensional consistency for variational systems. This is a development along the line of research of discrete integrable Lagrangian systems initiated in 2009 by Lobb and Nijhoff, however, having its more remote roots in the theory of pluriharmonic functions, in the Z-invariant models of statistical mechanics and their quasiclassical limit, as well as in the theory of variational symmetries going back to Noether. A d-dimensional pluri-Lagrangian problem can be described as follows: given a d-form on an m-dimensional space (called multi-time, m>d), whose coefficients depend on a sought-after function x of m independent variables (called field), find those fields x which deliver critical points to the action functionals for any d-dimensional manifold Σ in the multi-time. We derive the main building blocks of the multi-time Euler–Lagrange equations for a discrete pluri-Lagrangian problem with d=2, the so-called corner equations, and discuss the notion of consistency of the system of corner equations. We analyse the system of corner equations for a special class of three-point two-forms, corresponding to integrable quad-equations of the ABS list. This allows us to close a conceptual gap of the work by Lobb and Nijhoff by showing that the corresponding two-forms are closed not only on solutions of (non-variational) quad-equations, but also on general solutions of the corresponding corner equations. We also find an example of a pluri-Lagrangian system not coming from a multi-dimensionally consistent system of quad-equations. PMID:24511254

On l(1): Optimal decentralized performance

NASA Technical Reports Server (NTRS)

Sourlas, Dennis; Manousiouthakis, Vasilios

1993-01-01

In this paper, the Manousiouthakis parametrization of all decentralized stabilizing controllers is employed in mathematically formulating the l(sup 1) optimal decentralized controller synthesis problem. The resulting optimization problem is infinite dimensional and therefore not directly amenable to computations. It is shown that finite dimensional optimization problems that have value arbitrarily close to the infinite dimensional one can be constructed. Based on this result, an algorithm that solves the l(sup 1) decentralized performance problems is presented. A global optimization approach to the solution of the infinite dimensional approximating problems is also discussed.

Effect of virtual memory on efficient solution of two model problems

NASA Technical Reports Server (NTRS)

Lambiotte, J. J., Jr.

1977-01-01

Computers with virtual memory architecture allow programs to be written as if they were small enough to be contained in memory. Two types of problems are investigated to show that this luxury can lead to quite an inefficient performance if the programmer does not interact strongly with the characteristics of the operating system when developing the program. The two problems considered are the simultaneous solutions of a large linear system of equations by Gaussian elimination and a model three-dimensional finite-difference problem. The Control Data STAR-100 computer runs are made to demonstrate the inefficiencies of programming the problems in the manner one would naturally do if the problems were indeed, small enough to be contained in memory. Program redesigns are presented which achieve large improvements in performance through changes in the computational procedure and the data base arrangement.

The lizard celestial compass detects linearly polarized light in the blue.

PubMed

Beltrami, Giulia; Parretta, Antonio; Petrucci, Ferruccio; Buttini, Paola; Bertolucci, Cristiano; Foà, Augusto

2012-09-15

The present study first examined whether ruin lizards, Podarcis sicula, are able to orientate using plane-polarized light produced by an LCD screen. Ruin lizards were trained and tested indoors, inside a hexagonal Morris water maze positioned under an LCD screen producing white polarized light with a single E-vector, which provided an axial cue. White polarized light did not include wavelengths in the UV. Lizards orientated correctly either when tested with E-vector parallel to the training axis or after 90 deg rotation of the E-vector direction, thus validating the apparatus. Further experiments examined whether there is a preferential region of the light spectrum to perceive the E-vector direction of polarized light. For this purpose, lizards reaching learning criteria under white polarized light were subdivided into four experimental groups. Each group was tested for orientation under a different spectrum of plane-polarized light (red, green, cyan and blue) with equalized photon flux density. Lizards tested under blue polarized light orientated correctly, whereas lizards tested under red polarized light were completely disoriented. Green polarized light was barely discernible by lizards, and thus insufficient for a correct functioning of their compass. When exposed to cyan polarized light, lizard orientation performances were optimal, indistinguishable from lizards detecting blue polarized light. Overall, the present results demonstrate that perception of linear polarization in the blue is necessary - and sufficient - for a proper functioning of the sky polarization compass of ruin lizards. This may be adaptively important, as detection of polarized light in the blue improves functioning of the polarization compass under cloudy skies, i.e. when the alternative celestial compass based on detection of the sun disk is rendered useless because the sun is obscured by clouds.

Is the quality of sushi ruined by freezing raw fish and squid? A randomized double-blind trial with sensory evaluation using discrimination testing.

PubMed

Iwata, Kentaro; Fukuchi, Takahiko; Yoshimura, Kenichi

2015-05-01

Sushi is a traditional Japanese cuisine enjoyed worldwide. However, using raw fish to make sushi may pose risk of certain parasitic infections, such as anisakidosis, which is most reported in Japan. This risk of infection can be eliminated by freezing fish; however, Japanese people are hesitant to freeze fish because it is believed that freezing ruins sushi's taste. A randomized double-blind trial with discrimination testing was conducted to examine the ability of Japanese individuals to distinguish between frozen and unfrozen sushi. A pair of mackerel and squid sushi, one once frozen and the other not, was provided to the participants, and they were asked to answer which one tasted better. Among 120 rounds of discrimination testing involving the consumption of 240 pieces of mackerel sushi, unfrozen sushi was believed to taste better in 42.5% (51 dishes) of cases, frozen sushi was thought to taste better in 49.2% (59 dishes), and the participants felt the taste was the same in 8.3% (10 dishes). The odds ratio for selecting unfrozen sushi as "tastes better" over frozen sushi was 0.86 (95% confidence interval [CI], .59-1.26; P = .45). For squid, unfrozen sushi was believed to be superior 48.3% of the time (58 dishes), and frozen sushi, 35.0% of the time (42 dishes). They were felt to be the same in 16.7% (20 dishes) (odds ratio, 1.38; 95% CI, .93-2.05; P = .11). Freezing raw fish did not ruin sushi's taste. These findings may encourage the practice of freezing fish before using it in sushi, helping to decrease the incidence of anisakidosis. © The Author 2015. Published by Oxford University Press on behalf of the Infectious Diseases Society of America.

Dispersed gold nanoparticles potentially ruin gold barley yellow dwarf virus and eliminate virus infectivity hazards

NASA Astrophysics Data System (ADS)

Alkubaisi, Noorah A.; Aref, Nagwa M. A.

2017-02-01

Gold nanoparticles (AuNPs) application melted barley yellow dwarf virus-PAV (BYDV-PAV) spherical nanoparticle capsids. Synergistic therapeutic effects for plant virus resistance were induced by interaction with binding units of prepared AuNPs in a water solution which was characterized and evaluated by zeta sizer, zeta potential and transmission electron microscopy (TEM). The yield of purified nanoparticles of BYDV-PAV was obtained from Hordeum vulgare (Barley) cultivars, local and Giza 121/Justo. It was 0.62 mg/ml from 27.30 g of infected leaves at an A260/A280 ratio. Virus nanoparticle has a spherical shape 30 nm in size by TEM. BYDV-PAV combined with AuNPs to challenge virus function in vivo and in vitro. Dual AuNPs existence in vivo and in vitro affected compacted configuration of viral capsid protein in the interior surface of capsomers, the outer surface, or between the interface of coat protein subunits for 24 and 48 h incubation period in vitro at room temperature. The sizes of AuNPs that had a potentially dramatic deteriorated effect are 3.151 and 31.67 nm with a different intensity of 75.3% for the former and 24.7% for the latter, which enhances optical sensing applications to eliminate virus infectivity. Damages of capsid protein due to AuNPs on the surface of virus subunits caused variable performance in four different types of TEM named puffed, deteriorated and decorated, ruined and vanished. Viral yield showed remarkably high-intensity degree of particle symmetry and uniformity in the local cultivar greater than in Giza 121/Justo cultivar. A high yield of ruined VLPs in the local cultivar than Justo cultivar was noticed. AuNPs indicated complete lysed VLPs and some deteriorated VLPs at 48 h.

NMR Analysis of Unknowns: An Introduction to 2D NMR Spectroscopy

ERIC Educational Resources Information Center

Alonso, David E.; Warren, Steven E.

2005-01-01

A study combined 1D (one-dimensional) and 2D (two-dimensional) NMR spectroscopy to solve structural organic problems of three unknowns, which include 2-, 3-, and 4-heptanone. Results showed [to the first power]H NMR and [to the thirteenth power]C NMR signal assignments for 2- and 3-heptanone were more challenging than for 4-heptanone owing to the…

Application of an inverse method for calculating three-dimensional fault geometries and clip vectors, Nun River Field, Nigeria

DOE Office of Scientific and Technical Information (OSTI.GOV)

Kerr, H.G.; White, N.

A general, automatic method for determining the three-dimensional geometry of a normal fault of any shape and size is applied to a three-dimensional seismic reflection data set from the Nun River field, Nigeria. In addition to calculating fault geometry, the method also automatically retrieves the extension direction without requiring any previous information about either the fault shape or the extension direction. Solutions are found by minimizing the misfit between sets of faults that are calculated from the observed geometries of two or more hanging-wall beds. In the example discussed here, the predicted fault surface is in excellent agreement with themore » shape of the seismically imaged fault. Although the calculated extension direction is oblique to the average strike of the fault, the value of this parameter is not well resolved. Our approach differs markedly from standard section-balancing models in two important ways. First, we do not assume that the extension direction is known, and second, the use of inverse theory ensures that formal confidence bounds can be determined for calculated fault geometries. This ability has important implications for a range of geological problems encountered at both exploration and production scales. In particular, once the three-dimensional displacement field has been constrained, the difficult but important problem of three-dimensional palinspastic restoration of hanging-wall structures becomes tractable.« less

Some issues in the simulation of two-phase flows: The relative velocity

DOE Office of Scientific and Technical Information (OSTI.GOV)

Gräbel, J.; Hensel, S.; Ueberholz, P.

In this paper we compare numerical approximations for solving the Riemann problem for a hyperbolic two-phase flow model in two-dimensional space. The model is based on mixture parameters of state where the relative velocity between the two-phase systems is taken into account. This relative velocity appears as a main discontinuous flow variable through the complete wave structure and cannot be recovered correctly by some numerical techniques when simulating the associated Riemann problem. Simulations are validated by comparing the results of the numerical calculation qualitatively with OpenFOAM software. Simulations also indicate that OpenFOAM is unable to resolve the relative velocity associatedmore » with the Riemann problem.« less

Exact Solution of the Two-Dimensional Problem on an Impact Ideal-Liquid Jet

NASA Astrophysics Data System (ADS)

Belik, V. D.

2018-05-01

The two-dimensional problem on the collision of a potential ideal-liquid jet, outflowing from a reservoir through a nozzle, with an infinite plane obstacle was considered for the case where the distance between the nozzle exit section and the obstacle is finite. An exact solution of this problem has been found using methods of the complex-variable function theory. Simple analytical expressions for the complex velocity of the liquid, its flow rate, and the force of action of the jet on the obstacle have been obtained. The velocity distributions of the liquid at the nozzle exit section, in the region of spreading of the jet, and at the obstacle have been constructed for different distances between the nozzle exit section and the obstacle. Analytical expressions for the thickness of the boundary layer and the Nusselt number at the point of stagnation of the jet have been obtained. A number of distributions of the local friction coefficient and the Nusselt number of the indicated jet are presented.

Analytic Approximations to the Free Boundary and Multi-dimensional Problems in Financial Derivatives Pricing

NASA Astrophysics Data System (ADS)

Lau, Chun Sing

This thesis studies two types of problems in financial derivatives pricing. The first type is the free boundary problem, which can be formulated as a partial differential equation (PDE) subject to a set of free boundary condition. Although the functional form of the free boundary condition is given explicitly, the location of the free boundary is unknown and can only be determined implicitly by imposing continuity conditions on the solution. Two specific problems are studied in details, namely the valuation of fixed-rate mortgages and CEV American options. The second type is the multi-dimensional problem, which involves multiple correlated stochastic variables and their governing PDE. One typical problem we focus on is the valuation of basket-spread options, whose underlying asset prices are driven by correlated geometric Brownian motions (GBMs). Analytic approximate solutions are derived for each of these three problems. For each of the two free boundary problems, we propose a parametric moving boundary to approximate the unknown free boundary, so that the original problem transforms into a moving boundary problem which can be solved analytically. The governing parameter of the moving boundary is determined by imposing the first derivative continuity condition on the solution. The analytic form of the solution allows the price and the hedging parameters to be computed very efficiently. When compared against the benchmark finite-difference method, the computational time is significantly reduced without compromising the accuracy. The multi-stage scheme further allows the approximate results to systematically converge to the benchmark results as one recasts the moving boundary into a piecewise smooth continuous function. For the multi-dimensional problem, we generalize the Kirk (1995) approximate two-asset spread option formula to the case of multi-asset basket-spread option. Since the final formula is in closed form, all the hedging parameters can also be derived in closed form. Numerical examples demonstrate that the pricing and hedging errors are in general less than 1% relative to the benchmark prices obtained by numerical integration or Monte Carlo simulation. By exploiting an explicit relationship between the option price and the underlying probability distribution, we further derive an approximate distribution function for the general basket-spread variable. It can be used to approximate the transition probability distribution of any linear combination of correlated GBMs. Finally, an implicit perturbation is applied to reduce the pricing errors by factors of up to 100. When compared against the existing methods, the basket-spread option formula coupled with the implicit perturbation turns out to be one of the most robust and accurate approximation methods.

Three New (2+1)-dimensional Integrable Systems and Some Related Darboux Transformations

NASA Astrophysics Data System (ADS)

Guo, Xiu-Rong

2016-06-01

We introduce two operator commutators by using different-degree loop algebras of the Lie algebra A1, then under the framework of zero curvature equations we generate two (2+1)-dimensional integrable hierarchies, including the (2+1)-dimensional shallow water wave (SWW) hierarchy and the (2+1)-dimensional Kaup-Newell (KN) hierarchy. Through reduction of the (2+1)-dimensional hierarchies, we get a (2+1)-dimensional SWW equation and a (2+1)-dimensional KN equation. Furthermore, we obtain two Darboux transformations of the (2+1)-dimensional SWW equation. Similarly, the Darboux transformations of the (2+1)-dimensional KN equation could be deduced. Finally, with the help of the spatial spectral matrix of SWW hierarchy, we generate a (2+1) heat equation and a (2+1) nonlinear generalized SWW system containing inverse operators with respect to the variables x and y by using a reduction spectral problem from the self-dual Yang-Mills equations. Supported by the National Natural Science Foundation of China under Grant No. 11371361, the Shandong Provincial Natural Science Foundation of China under Grant Nos. ZR2012AQ011, ZR2013AL016, ZR2015EM042, National Social Science Foundation of China under Grant No. 13BJY026, the Development of Science and Technology Project under Grant No. 2015NS1048 and A Project of Shandong Province Higher Educational Science and Technology Program under Grant No. J14LI58

A matrix dependent/algebraic multigrid approach for extruded meshes with applications to ice sheet modeling

DOE PAGES

Tuminaro, Raymond S.; Perego, Mauro; Tezaur, Irina Kalashnikova; ...

2016-10-06

A multigrid method is proposed that combines ideas from matrix dependent multigrid for structured grids and algebraic multigrid for unstructured grids. It targets problems where a three-dimensional mesh can be viewed as an extrusion of a two-dimensional, unstructured mesh in a third dimension. Our motivation comes from the modeling of thin structures via finite elements and, more specifically, the modeling of ice sheets. Extruded meshes are relatively common for thin structures and often give rise to anisotropic problems when the thin direction mesh spacing is much smaller than the broad direction mesh spacing. Within our approach, the first few multigridmore » hierarchy levels are obtained by applying matrix dependent multigrid to semicoarsen in a structured thin direction fashion. After sufficient structured coarsening, the resulting mesh contains only a single layer corresponding to a two-dimensional, unstructured mesh. Algebraic multigrid can then be employed in a standard manner to create further coarse levels, as the anisotropic phenomena is no longer present in the single layer problem. The overall approach remains fully algebraic, with the minor exception that some additional information is needed to determine the extruded direction. Furthermore, this facilitates integration of the solver with a variety of different extruded mesh applications.« less

Additive schemes for certain operator-differential equations

NASA Astrophysics Data System (ADS)

Vabishchevich, P. N.

2010-12-01

Unconditionally stable finite difference schemes for the time approximation of first-order operator-differential systems with self-adjoint operators are constructed. Such systems arise in many applied problems, for example, in connection with nonstationary problems for the system of Stokes (Navier-Stokes) equations. Stability conditions in the corresponding Hilbert spaces for two-level weighted operator-difference schemes are obtained. Additive (splitting) schemes are proposed that involve the solution of simple problems at each time step. The results are used to construct splitting schemes with respect to spatial variables for nonstationary Navier-Stokes equations for incompressible fluid. The capabilities of additive schemes are illustrated using a two-dimensional model problem as an example.

The Position Control of the Surface Motor with the Poles Distribution of Triangular Lattice

NASA Astrophysics Data System (ADS)

Watada, Masaya; Katsuyama, Norikazu; Ebihara, Daiki

Recently, as for the machine tools or industrial robots, high performance, accuracy, etc. are demanded. Generally, when drive of many degrees of freedom is required in the machine tools or industrial robots, it has realized by using two or more motors. For example, two-dimensional positioning stages such as the X-Y plotter or the X-Y stage are enabling the two-dimensional drive by using each one motor in the direction of x, y. In order to use plural motors, these, however, have problems that equipment becomes large and complicate control system. From such problems, the Surface Motor (SFM) that can drive two directions by only one motor is researched. Authors have proposed SFM that considered wide range movement and the application to a curved surface. In this paper, the characteristics of the micro step drive by the open loop control are showed. Introduction of closed loop control for highly accurate positioning, moreover, is examined. The drive characteristics by each control are compared.

Analysis of students’ spatial thinking in geometry: 3D object into 2D representation

NASA Astrophysics Data System (ADS)

Fiantika, F. R.; Maknun, C. L.; Budayasa, I. K.; Lukito, A.

2018-05-01

The aim of this study is to find out the spatial thinking process of students in transforming 3-dimensional (3D) object to 2-dimensional (2D) representation. Spatial thinking is helpful in using maps, planning routes, designing floor plans, and creating art. The student can engage geometric ideas by using concrete models and drawing. Spatial thinking in this study is identified through geometrical problems of transforming a 3-dimensional object into a 2-dimensional object image. The problem was resolved by the subject and analyzed by reference to predetermined spatial thinking indicators. Two representative subjects of elementary school were chosen based on mathematical ability and visual learning style. Explorative description through qualitative approach was used in this study. The result of this study are: 1) there are different representations of spatial thinking between a boy and a girl object, 2) the subjects has their own way to invent the fastest way to draw cube net.

Thermomagnetic instabilities in a vertical layer of ferrofluid: nonlinear analysis away from a critical point

NASA Astrophysics Data System (ADS)

Dey, Pinkee; Suslov, Sergey A.

2016-12-01

A finite amplitude instability has been analysed to discover the exact mechanism leading to the appearance of stationary magnetoconvection patterns in a vertical layer of a non-conducting ferrofluid heated from the side and placed in an external magnetic field perpendicular to the walls. The physical results have been obtained using a version of a weakly nonlinear analysis that is based on the disturbance amplitude expansion. It enables a low-dimensional reduction of a full nonlinear problem in supercritical regimes away from a bifurcation point. The details of the reduction are given in comparison with traditional small-parameter expansions. It is also demonstrated that Squire’s transformation can be introduced for higher-order nonlinear terms thus reducing the full three-dimensional problem to its equivalent two-dimensional counterpart and enabling significant computational savings. The full three-dimensional instability patterns are subsequently recovered using the inverse transforms The analysed stationary thermomagnetic instability is shown to occur as a result of a supercritical pitchfork bifurcation.

A Comparative Study on Multifactor Dimensionality Reduction Methods for Detecting Gene-Gene Interactions with the Survival Phenotype

PubMed Central

Lee, Seungyeoun; Kim, Yongkang; Kwon, Min-Seok; Park, Taesung

2015-01-01

Genome-wide association studies (GWAS) have extensively analyzed single SNP effects on a wide variety of common and complex diseases and found many genetic variants associated with diseases. However, there is still a large portion of the genetic variants left unexplained. This missing heritability problem might be due to the analytical strategy that limits analyses to only single SNPs. One of possible approaches to the missing heritability problem is to consider identifying multi-SNP effects or gene-gene interactions. The multifactor dimensionality reduction method has been widely used to detect gene-gene interactions based on the constructive induction by classifying high-dimensional genotype combinations into one-dimensional variable with two attributes of high risk and low risk for the case-control study. Many modifications of MDR have been proposed and also extended to the survival phenotype. In this study, we propose several extensions of MDR for the survival phenotype and compare the proposed extensions with earlier MDR through comprehensive simulation studies. PMID:26339630

Using three dimensional silicone ``boots`` to solve complex remedial design problems in curtain walls

DOE Office of Scientific and Technical Information (OSTI.GOV)

Chen, Y.J.

1998-12-31

Stick system curtain wall leak problems are frequently caused by water entry at the splice joints of the curtain wall frame and failure of the internal metal joinery seals. Remedial solutions involving occupied buildings inevitably face the multiple constraints of existing construction and business operations not present during the original curtain wall construction. In most cases, even partial disassembly of the curtain wall for internal seal repairs is not feasible. Remedial solutions which must be executed from the exterior of the curtain wall often involve wet-applied or preformed sealant tape bridge joints. However, some of the more complex joints cannotmore » be repaired effectively or economically with the conventional bridge joint. Fortunately, custom fabricated three-dimensional preformed sealant boots are becoming available to address these situations. This paper discusses the design considerations and the selective use of three-dimensional preformed boots in sealing complex joint geometry that would not be effective with the conventional two-dimensional bridge joint.« less

A coupled sharp-interface immersed boundary-finite-element method for flow-structure interaction with application to human phonation.

PubMed

Zheng, X; Xue, Q; Mittal, R; Beilamowicz, S

2010-11-01

A new flow-structure interaction method is presented, which couples a sharp-interface immersed boundary method flow solver with a finite-element method based solid dynamics solver. The coupled method provides robust and high-fidelity solution for complex flow-structure interaction (FSI) problems such as those involving three-dimensional flow and viscoelastic solids. The FSI solver is used to simulate flow-induced vibrations of the vocal folds during phonation. Both two- and three-dimensional models have been examined and qualitative, as well as quantitative comparisons, have been made with established results in order to validate the solver. The solver is used to study the onset of phonation in a two-dimensional laryngeal model and the dynamics of the glottal jet in a three-dimensional model and results from these studies are also presented.

Artificial viscosity to cure the carbuncle phenomenon: The three-dimensional case

NASA Astrophysics Data System (ADS)

Rodionov, Alexander V.

2018-05-01

The carbuncle phenomenon (also known as the shock instability) has remained a serious computational challenge since it was first noticed and described [1,2]. In [3] the author presented a summary on this subject and proposed a new technique for curing the problem. Its idea is to introduce some dissipation in the form of right-hand sides of the Navier-Stokes equations into the basic method of solving Euler equations; in so doing, the molecular viscosity coefficient is replaced by the artificial viscosity coefficient. The new cure for the carbuncle flaw was tested and tuned for the case of using first-order schemes in two-dimensional simulations. Its efficiency was demonstrated on several well-known test problems. In this paper we extend the technique of [3] to the case of three-dimensional simulations.

Separation of variables in the special diagonal Hamilton-Jacobi equation: Application to the dynamical problem of a particle constrained on a moving surface

NASA Technical Reports Server (NTRS)

Blanchard, D. L.; Chan, F. K.

1973-01-01

For a time-dependent, n-dimensional, special diagonal Hamilton-Jacobi equation a necessary and sufficient condition for the separation of variables to yield a complete integral of the form was established by specifying the admissible forms in terms of arbitrary functions. A complete integral was then expressed in terms of these arbitrary functions and also the n irreducible constants. As an application of the results obtained for the two-dimensional Hamilton-Jacobi equation, analysis was made for a comparatively wide class of dynamical problems involving a particle moving in Euclidean three-dimensional space under the action of external forces but constrained on a moving surface. All the possible cases in which this equation had a complete integral of the form were obtained and these are tubulated for reference.

General View of Kinckley Knitting Mill site, east end of ...

Library of Congress Historic Buildings Survey, Historic Engineering Record, Historic Landscapes Survey

General View of Kinckley Knitting Mill site, east end of site along Wister Street, looking north, showing ruin of stone mill building and smoke stack of the Germantown Dyeworks in background - Hinckley Knitting Mills, 21-35 East Wister Street, Philadelphia, Philadelphia County, PA

The Phony Funding Crisis

ERIC Educational Resources Information Center

Guthrie, James W.; Peng, Arthur

2010-01-01

If one relies on newspaper headlines for education funding information, one might conclude that America's schools suffer from a perpetual fiscal crisis, every year perched precariously on the brink of financial ruin, never knowing whether there will be sufficient funding to continue operating. Budgetary shortfalls, school district bankruptcies,…

Teaching Archaeology. ERIC Digest.

ERIC Educational Resources Information Center

Brown, Gail William

How could handchipped stones, ancient ruins, old broken dishes, and antiquated garbage help students learn about the world and themselves? Within archaeology, these seemingly irrelevant items can enlighten students about the world around them through science, culture, and history. When teaching archaeology in the classroom, educators can lead…

Ages and Stages: Teen

MedlinePlus

... Syndrome Tanning and Tanning Salon Safety Tips for Young People Teenagers and Gangs Teens and Acne Treatment Teens and Sun: Keeping Them Safe Without Ruining Their Fun The Teenage Brain VIDEO Ways To Build Your Teenager’s Self-Esteem What’s Going On in the Teenage Brain? Zits ...

The Emotionally Sensitive Adolescent.

ERIC Educational Resources Information Center

Shaughnessy, Michael F.; Lehtonen, Kimmo

This paper provides a list of signs, symptoms, and indicators of emotionally sensitive adolescents includes clinging behavior, withdrawn behavior, shy/inhibited behavior, represses anger, poor reaction to criticism, makes self-disparaging statements, low self-esteem, "can't forgive self or others," ruined by a small critical comment, exploding…

19. AERIAL VIEW, LOOKING SOUTHWEST, OF SILVER LAKE WEST SHORE ...

Library of Congress Historic Buildings Survey, Historic Engineering Record, Historic Landscapes Survey

19. AERIAL VIEW, LOOKING SOUTHWEST, OF SILVER LAKE WEST SHORE WITH RUINS OF SILVER LAKE MINE AT BOTTOM CENTER. NOTE CONTINUITY OF GEOLOGICAL FORMATIONS WHERE BOTTOM OF LAKE CAN BE SEEN. - Shenandoah-Dives Mill, 135 County Road 2, Silverton, San Juan County, CO

Prescription for a Revolution.

ERIC Educational Resources Information Center

Olson, Lynn

1990-01-01

Discusses the book "Politics, Markets, and America's Schools," by John Chubb and Terry Moe. They claim the school reform movement has not worked, and a massive, deadening bureaucracy is ruining American education. They suggest a new system for effective schooling based on good organization and parental choice. (SM)

Node-Based Learning of Multiple Gaussian Graphical Models

PubMed Central

Mohan, Karthik; London, Palma; Fazel, Maryam; Witten, Daniela; Lee, Su-In

2014-01-01

We consider the problem of estimating high-dimensional Gaussian graphical models corresponding to a single set of variables under several distinct conditions. This problem is motivated by the task of recovering transcriptional regulatory networks on the basis of gene expression data containing heterogeneous samples, such as different disease states, multiple species, or different developmental stages. We assume that most aspects of the conditional dependence networks are shared, but that there are some structured differences between them. Rather than assuming that similarities and differences between networks are driven by individual edges, we take a node-based approach, which in many cases provides a more intuitive interpretation of the network differences. We consider estimation under two distinct assumptions: (1) differences between the K networks are due to individual nodes that are perturbed across conditions, or (2) similarities among the K networks are due to the presence of common hub nodes that are shared across all K networks. Using a row-column overlap norm penalty function, we formulate two convex optimization problems that correspond to these two assumptions. We solve these problems using an alternating direction method of multipliers algorithm, and we derive a set of necessary and sufficient conditions that allows us to decompose the problem into independent subproblems so that our algorithm can be scaled to high-dimensional settings. Our proposal is illustrated on synthetic data, a webpage data set, and a brain cancer gene expression data set. PMID:25309137

Adjoint shape optimization for fluid-structure interaction of ducted flows

NASA Astrophysics Data System (ADS)

Heners, J. P.; Radtke, L.; Hinze, M.; Düster, A.

2018-03-01

Based on the coupled problem of time-dependent fluid-structure interaction, equations for an appropriate adjoint problem are derived by the consequent use of the formal Lagrange calculus. Solutions of both primal and adjoint equations are computed in a partitioned fashion and enable the formulation of a surface sensitivity. This sensitivity is used in the context of a steepest descent algorithm for the computation of the required gradient of an appropriate cost functional. The efficiency of the developed optimization approach is demonstrated by minimization of the pressure drop in a simple two-dimensional channel flow and in a three-dimensional ducted flow surrounded by a thin-walled structure.

Hybrid-dimensional modelling of two-phase flow through fractured porous media with enhanced matrix fracture transmission conditions

NASA Astrophysics Data System (ADS)

Brenner, Konstantin; Hennicker, Julian; Masson, Roland; Samier, Pierre

2018-03-01

In this work, we extend, to two-phase flow, the single-phase Darcy flow model proposed in [26], [12] in which the (d - 1)-dimensional flow in the fractures is coupled with the d-dimensional flow in the matrix. Three types of so called hybrid-dimensional two-phase Darcy flow models are proposed. They all account for fractures acting either as drains or as barriers, since they allow pressure jumps at the matrix-fracture interfaces. The models also permit to treat gravity dominated flow as well as discontinuous capillary pressure at the material interfaces. The three models differ by their transmission conditions at matrix fracture interfaces: while the first model accounts for the nonlinear two-phase Darcy flux conservations, the second and third ones are based on the linear single phase Darcy flux conservations combined with different approximations of the mobilities. We adapt the Vertex Approximate Gradient (VAG) scheme to this problem, in order to account for anisotropy and heterogeneity aspects as well as for applicability on general meshes. Several test cases are presented to compare our hybrid-dimensional models to the generic equi-dimensional model, in which fractures have the same dimension as the matrix, leading to deep insight about the quality of the proposed reduced models.

Vehicle Routing with Three-dimensional Container Loading Constraints—Comparison of Nested and Joint Algorithms

NASA Astrophysics Data System (ADS)

Koloch, Grzegorz; Kaminski, Bogumil

2010-10-01

In the paper we examine a modification of the classical Vehicle Routing Problem (VRP) in which shapes of transported cargo are accounted for. This problem, known as a three-dimensional VRP with loading constraints (3D-VRP), is appropriate when transported commodities are not perfectly divisible, but they have fixed and heterogeneous dimensions. In the paper restrictions on allowable cargo positionings are also considered. These restrictions are derived from business practice and they extended the baseline 3D-VRP formulation as considered by Koloch and Kaminski (2010). In particular, we investigate how additional restrictions influence relative performance of two proposed optimization algorithms: the nested and the joint one. Performance of both methods is compared on artificial problems and on a big-scale real life case study.

Computer modeling of inversion layer MOS solar cells and arrays

NASA Technical Reports Server (NTRS)

Ho, Fat Duen

1991-01-01

A two dimensional numerical model of the inversion layer metal insulator semiconductor (IL/MIS) solar cell is proposed by using the finite element method. The two-dimensional current flow in the device is taken into account in this model. The electrostatic potential distribution, the electron concentration distribution, and the hole concentration distribution for different terminal voltages are simulated. The results of simple calculation are presented. The existing problems for this model are addressed. Future work is proposed. The MIS structures are studied and some of the results are reported.

Integration of the Two-Dimensional Power Spectral Density into Specifications for the X-ray Domain -- Problems and Opportunities

DOE Office of Scientific and Technical Information (OSTI.GOV)

McKinney, Wayne R.; Howells, M. R.; Yashchuk, V. V.

2008-09-30

An implementation of the two-dimensional statistical scattering theory of Church and Takacs for the prediction of scattering from x-ray mirrors is presented with a graphical user interface. The process of this development has clarified several problems which are of significant interest to the synchrotron community. These problems have been addressed to some extent, for example, for large astronomical telescopes, and at the National Ignition Facility for normal incidence optics, but not in the synchrotron community for grazing incidence optics. Since it is based on the Power Spectral Density (PSD) to provide a description of the deviations from ideal shape ofmore » the surface, accurate prediction of the scattering requires an accurate estimation of the PSD. Specifically, the spatial frequency range of measurement must be the correct one for the geometry of use of the optic--including grazing incidence and coherence effects, and the modifications to the PSD of the Optical Transfer Functions (OTF) of the measuring instruments must be removed. A solution for removal of OTF effects has been presented previously, the Binary Pseudo-Random Grating. Typically, the frequency range of a single instrument does not cover the range of interest, requiring the stitching together of PSD estimations. This combination generates its own set of difficulties in two dimensions. Fitting smooth functions to two dimensional PSDs, particularly in the case of spatial non-isotropy of the surface, which is often the case for optics in synchrotron beam lines, can be difficult. The convenient, and physically accurate fractal for one dimension does not readily transfer to two dimensions. Finally, a completely statistical description of scattering must be integrated with a deterministic low spatial frequency component in order to completely model the intensity near the image. An outline for approaching these problems, and our proposed experimental program is given.« less

A Matlab toolkit for three-dimensional electrical impedance tomography: a contribution to the Electrical Impedance and Diffuse Optical Reconstruction Software project

NASA Astrophysics Data System (ADS)

Polydorides, Nick; Lionheart, William R. B.

2002-12-01

The objective of the Electrical Impedance and Diffuse Optical Reconstruction Software project is to develop freely available software that can be used to reconstruct electrical or optical material properties from boundary measurements. Nonlinear and ill posed problems such as electrical impedance and optical tomography are typically approached using a finite element model for the forward calculations and a regularized nonlinear solver for obtaining a unique and stable inverse solution. Most of the commercially available finite element programs are unsuitable for solving these problems because of their conventional inefficient way of calculating the Jacobian, and their lack of accurate electrode modelling. A complete package for the two-dimensional EIT problem was officially released by Vauhkonen et al at the second half of 2000. However most industrial and medical electrical imaging problems are fundamentally three-dimensional. To assist the development we have developed and released a free toolkit of Matlab routines which can be employed to solve the forward and inverse EIT problems in three dimensions based on the complete electrode model along with some basic visualization utilities, in the hope that it will stimulate further development. We also include a derivation of the formula for the Jacobian (or sensitivity) matrix based on the complete electrode model.

Theoretical thermal conductivity equation for uniform density wood cells

Treesearch

John F. Hunt; Hongmei Gu; Patricia Lebow

2008-01-01

The anisotropy of wood creates a complex problem requiring that analyses be based on fundamental material properties and characteristics of the wood structure to solve heat transfer problems. A two-dimensional finite element model that evaluates the effective thermal conductivity of a wood cell over the full range of moisture contents and porosities was previously...