A numerical analysis on forming limits during spiral and concentric single point incremental forming

Science.gov (United States)

Gipiela, M. L.; Amauri, V.; Nikhare, C.; Marcondes, P. V. P.

2017-01-01

Sheet metal forming is one of the major manufacturing industries, which are building numerous parts for aerospace, automotive and medical industry. Due to the high demand in vehicle industry and environmental regulations on less fuel consumption on other hand, researchers are innovating new methods to build these parts with energy efficient sheet metal forming process instead of conventionally used punch and die to form the parts to achieve the lightweight parts. One of the most recognized manufacturing process in this category is Single Point Incremental Forming (SPIF). SPIF is the die-less sheet metal forming process in which the single point tool incrementally forces any single point of sheet metal at any process time to plastic deformation zone. In the present work, finite element method (FEM) is applied to analyze the forming limits of high strength low alloy steel formed by single point incremental forming (SPIF) by spiral and concentric tool path. SPIF numerical simulations were model with 24 and 29 mm cup depth, and the results were compare with Nakajima results obtained by experiments and FEM. It was found that the cup formed with Nakajima tool failed at 24 mm while cups formed by SPIF surpassed the limit for both depths with both profiles. It was also notice that the strain achieved in concentric profile are lower than that in spiral profile.

Numerical instability of time-discretized one-point kinetic equations

International Nuclear Information System (INIS)

Hashimoto, Kengo; Ikeda, Hideaki; Takeda, Toshikazu

2000-01-01

The one-point kinetic equations with numerical errors induced by the explicit, implicit and Crank-Nicolson integration methods are derived. The zero-power transfer functions based on the present equations are demonstrated to investigate the numerical stability of the discretized systems. These demonstrations indicate unconditional stability for the implicit and Crank-Nicolson methods but present the possibility of numerical instability for the explicit method. An upper limit of time mesh spacing for the stability is formulated and several numerical calculations are made to confirm the validity of this formula

Elastic dipoles of point defects from atomistic simulations

Science.gov (United States)

Varvenne, Céline; Clouet, Emmanuel

2017-12-01

The interaction of point defects with an external stress field or with other structural defects is usually well described within continuum elasticity by the elastic dipole approximation. Extraction of the elastic dipoles from atomistic simulations is therefore a fundamental step to connect an atomistic description of the defect with continuum models. This can be done either by a fitting of the point-defect displacement field, by a summation of the Kanzaki forces, or by a linking equation to the residual stress. We perform here a detailed comparison of these different available methods to extract elastic dipoles, and show that they all lead to the same values when the supercell of the atomistic simulations is large enough and when the anharmonic region around the point defect is correctly handled. But, for small simulation cells compatible with ab initio calculations, only the definition through the residual stress appears tractable. The approach is illustrated by considering various point defects (vacancy, self-interstitial, and hydrogen solute atom) in zirconium, using both empirical potentials and ab initio calculations.

Development of modifications to the material point method for the simulation of thin membranes, compressible fluids, and their interactions

Energy Technology Data Exchange (ETDEWEB)

York, A.R. II [Sandia National Labs., Albuquerque, NM (United States). Engineering and Process Dept.

1997-07-01

The material point method (MPM) is an evolution of the particle in cell method where Lagrangian particles or material points are used to discretize the volume of a material. The particles carry properties such as mass, velocity, stress, and strain and move through a Eulerian or spatial mesh. The momentum equation is solved on the Eulerian mesh. Modifications to the material point method are developed that allow the simulation of thin membranes, compressible fluids, and their dynamic interactions. A single layer of material points through the thickness is used to represent a membrane. The constitutive equation for the membrane is applied in the local coordinate system of each material point. Validation problems are presented and numerical convergence is demonstrated. Fluid simulation is achieved by implementing a constitutive equation for a compressible, viscous, Newtonian fluid and by solution of the energy equation. The fluid formulation is validated by simulating a traveling shock wave in a compressible fluid. Interactions of the fluid and membrane are handled naturally with the method. The fluid and membrane communicate through the Eulerian grid on which forces are calculated due to the fluid and membrane stress states. Validation problems include simulating a projectile impacting an inflated airbag. In some impact simulations with the MPM, bodies may tend to stick together when separating. Several algorithms are proposed and tested that allow bodies to separate from each other after impact. In addition, several methods are investigated to determine the local coordinate system of a membrane material point without relying upon connectivity data.

Evaluation of null-point detection methods on simulation data

Science.gov (United States)

Olshevsky, Vyacheslav; Fu, Huishan; Vaivads, Andris; Khotyaintsev, Yuri; Lapenta, Giovanni; Markidis, Stefano

2014-05-01

We model the measurements of artificial spacecraft that resemble the configuration of CLUSTER propagating in the particle-in-cell simulation of turbulent magnetic reconnection. The simulation domain contains multiple isolated X-type null-points, but the majority are O-type null-points. Simulations show that current pinches surrounded by twisted fields, analogous to laboratory pinches, are formed along the sequences of O-type nulls. In the simulation, the magnetic reconnection is mainly driven by the kinking of the pinches, at spatial scales of several ion inertial lentghs. We compute the locations of magnetic null-points and detect their type. When the satellites are separated by the fractions of ion inertial length, as it is for CLUSTER, they are able to locate both the isolated null-points, and the pinches. We apply the method to the real CLUSTER data and speculate how common are pinches in the magnetosphere, and whether they play a dominant role in the dissipation of magnetic energy.

Numerical analysis on pump turbine runaway points

International Nuclear Information System (INIS)

Guo, L; Liu, J T; Wang, L Q; Jiao, L; Li, Z F

2012-01-01

To research the character of pump turbine runaway points with different guide vane opening, a hydraulic model was established based on a pumped storage power station. The RNG k-ε model and SMPLEC algorithms was used to simulate the internal flow fields. The result of the simulation was compared with the test data and good correspondence was got between experimental data and CFD result. Based on this model, internal flow analysis was carried out. The result show that when the pump turbine ran at the runway speed, lots of vortexes appeared in the flow passage of the runner. These vortexes could always be observed even if the guide vane opening changes. That is an important way of energy loss in the runaway condition. Pressure on two sides of the runner blades were almost the same. So the runner power is very low. High speed induced large centrifugal force and the small guide vane opening gave the water velocity a large tangential component, then an obvious water ring could be observed between the runner blades and guide vanes in small guide vane opening condition. That ring disappeared when the opening bigger than 20°. These conclusions can provide a theory basis for the analysis and simulation of the pump turbine runaway points.

Numerical Treatment of Fixed Point Applied to the Nonlinear Fredholm Integral Equation

Directory of Open Access Journals (Sweden)

Berenguer MI

2009-01-01

Full Text Available The authors present a method of numerical approximation of the fixed point of an operator, specifically the integral one associated with a nonlinear Fredholm integral equation, that uses strongly the properties of a classical Schauder basis in the Banach space .

A full three dimensional Navier-Stokes numerical simulation of flow field inside a power plant Kaplan turbine using some model test turbine hill chart points

Science.gov (United States)

Hosseinalipour, S. M.; Raja, A.; Hajikhani, S.

2012-06-01

A full three dimensional Navier - Stokes numerical simulation has been performed for performance analysis of a Kaplan turbine which is installed in one of the Irans south dams. No simplifications have been enforced in the simulation. The numerical results have been evaluated using some integral parameters such as the turbine efficiency via comparing the results with existing experimental data from the prototype Hill chart. In part of this study the numerical simulations were performed in order to calculate the prototype turbine efficiencies in some specific points which comes from the scaling up of the model efficiency that are available in the model experimental Hill chart. The results are very promising which shows the good ability of the numerical techniques for resolving the flow characteristics in these kind of complex geometries. A parametric study regarding the evaluation of turbine performance in three different runner angles of the prototype is also performed and the results are cited in this paper.

Using a Virtual Experiment to Analyze Infiltration Process from Point to Grid-cell Size Scale

Science.gov (United States)

Barrios, M. I.

2013-12-01

The hydrological science requires the emergence of a consistent theoretical corpus driving the relationships between dominant physical processes at different spatial and temporal scales. However, the strong spatial heterogeneities and non-linearities of these processes make difficult the development of multiscale conceptualizations. Therefore, scaling understanding is a key issue to advance this science. This work is focused on the use of virtual experiments to address the scaling of vertical infiltration from a physically based model at point scale to a simplified physically meaningful modeling approach at grid-cell scale. Numerical simulations have the advantage of deal with a wide range of boundary and initial conditions against field experimentation. The aim of the work was to show the utility of numerical simulations to discover relationships between the hydrological parameters at both scales, and to use this synthetic experience as a media to teach the complex nature of this hydrological process. The Green-Ampt model was used to represent vertical infiltration at point scale; and a conceptual storage model was employed to simulate the infiltration process at the grid-cell scale. Lognormal and beta probability distribution functions were assumed to represent the heterogeneity of soil hydraulic parameters at point scale. The linkages between point scale parameters and the grid-cell scale parameters were established by inverse simulations based on the mass balance equation and the averaging of the flow at the point scale. Results have shown numerical stability issues for particular conditions and have revealed the complex nature of the non-linear relationships between models' parameters at both scales and indicate that the parameterization of point scale processes at the coarser scale is governed by the amplification of non-linear effects. The findings of these simulations have been used by the students to identify potential research questions on scale issues

Downscaling remotely sensed imagery using area-to-point cokriging and multiple-point geostatistical simulation

Science.gov (United States)

Tang, Yunwei; Atkinson, Peter M.; Zhang, Jingxiong

2015-03-01

A cross-scale data integration method was developed and tested based on the theory of geostatistics and multiple-point geostatistics (MPG). The goal was to downscale remotely sensed images while retaining spatial structure by integrating images at different spatial resolutions. During the process of downscaling, a rich spatial correlation model in the form of a training image was incorporated to facilitate reproduction of similar local patterns in the simulated images. Area-to-point cokriging (ATPCK) was used as locally varying mean (LVM) (i.e., soft data) to deal with the change of support problem (COSP) for cross-scale integration, which MPG cannot achieve alone. Several pairs of spectral bands of remotely sensed images were tested for integration within different cross-scale case studies. The experiment shows that MPG can restore the spatial structure of the image at a fine spatial resolution given the training image and conditioning data. The super-resolution image can be predicted using the proposed method, which cannot be realised using most data integration methods. The results show that ATPCK-MPG approach can achieve greater accuracy than methods which do not account for the change of support issue.

Critical point analysis of phase envelope diagram

International Nuclear Information System (INIS)

Soetikno, Darmadi; Siagian, Ucok W. R.; Kusdiantara, Rudy; Puspita, Dila; Sidarto, Kuntjoro A.; Soewono, Edy; Gunawan, Agus Y.

2014-01-01

Phase diagram or phase envelope is a relation between temperature and pressure that shows the condition of equilibria between the different phases of chemical compounds, mixture of compounds, and solutions. Phase diagram is an important issue in chemical thermodynamics and hydrocarbon reservoir. It is very useful for process simulation, hydrocarbon reactor design, and petroleum engineering studies. It is constructed from the bubble line, dew line, and critical point. Bubble line and dew line are composed of bubble points and dew points, respectively. Bubble point is the first point at which the gas is formed when a liquid is heated. Meanwhile, dew point is the first point where the liquid is formed when the gas is cooled. Critical point is the point where all of the properties of gases and liquids are equal, such as temperature, pressure, amount of substance, and others. Critical point is very useful in fuel processing and dissolution of certain chemicals. Here in this paper, we will show the critical point analytically. Then, it will be compared with numerical calculations of Peng-Robinson equation by using Newton-Raphson method. As case studies, several hydrocarbon mixtures are simulated using by Matlab

Critical point analysis of phase envelope diagram

Energy Technology Data Exchange (ETDEWEB)

Soetikno, Darmadi; Siagian, Ucok W. R. [Department of Petroleum Engineering, Institut Teknologi Bandung, Jl. Ganesha 10, Bandung 40132 (Indonesia); Kusdiantara, Rudy, E-mail: rkusdiantara@s.itb.ac.id; Puspita, Dila, E-mail: rkusdiantara@s.itb.ac.id; Sidarto, Kuntjoro A., E-mail: rkusdiantara@s.itb.ac.id; Soewono, Edy; Gunawan, Agus Y. [Department of Mathematics, Institut Teknologi Bandung, Jl. Ganesha 10, Bandung 40132 (Indonesia)

2014-03-24

Phase diagram or phase envelope is a relation between temperature and pressure that shows the condition of equilibria between the different phases of chemical compounds, mixture of compounds, and solutions. Phase diagram is an important issue in chemical thermodynamics and hydrocarbon reservoir. It is very useful for process simulation, hydrocarbon reactor design, and petroleum engineering studies. It is constructed from the bubble line, dew line, and critical point. Bubble line and dew line are composed of bubble points and dew points, respectively. Bubble point is the first point at which the gas is formed when a liquid is heated. Meanwhile, dew point is the first point where the liquid is formed when the gas is cooled. Critical point is the point where all of the properties of gases and liquids are equal, such as temperature, pressure, amount of substance, and others. Critical point is very useful in fuel processing and dissolution of certain chemicals. Here in this paper, we will show the critical point analytically. Then, it will be compared with numerical calculations of Peng-Robinson equation by using Newton-Raphson method. As case studies, several hydrocarbon mixtures are simulated using by Matlab.

Numerical analysis of slowest heating or cooling point in a canned food in oil

Energy Technology Data Exchange (ETDEWEB)

Hanzawa, T.; Wang, Q.; Suzuki, M.; Sakai, N. [Tokyo Univ. of Fisheries (Japan)

1998-06-01

In the sterilizing process of canned food in oil for a fish meat such as tunny, the slowest heating or cooling point is very important for the thermal process determination of the can. To obtain the slowest point, the temperature profiles in solid food are estimated by numerical calculation from the fundamental equations at unsteady state in consideration of a free convection in the space occupied by the oil. The positions of the slowest heating or cooling point in the canned food in oil are obtained accurately, and a correlative equation for the position is obtained numerically under various operating conditions. The calculated temperature profiles and the position of both slowest points are in sufficiently good approximation to the experimental ones. 4 refs., 9 figs.

Experimental research and numerical optimisation of multi-point sheet metal forming implementation using a solid elastic cushion system

Science.gov (United States)

Tolipov, A. A.; Elghawail, A.; Shushing, S.; Pham, D.; Essa, K.

2017-09-01

There is a growing demand for flexible manufacturing techniques that meet the rapid changes in customer needs. A finite element analysis numerical optimisation technique was used to optimise the multi-point sheet forming process. Multi-point forming (MPF) is a flexible sheet metal forming technique where the same tool can be readily changed to produce different parts. The process suffers from some geometrical defects such as wrinkling and dimpling, which have been found to be the cause of the major surface quality problems. This study investigated the influence of parameters such as the elastic cushion hardness, blank holder force, coefficient of friction, cushion thickness and radius of curvature, on the quality of parts formed in a flexible multi-point stamping die. For those reasons, in this investigation, a multipoint forming stamping process using a blank holder was carried out in order to study the effects of the wrinkling, dimpling, thickness variation and forming force. The aim was to determine the optimum values of these parameters. Finite element modelling (FEM) was employed to simulate the multi-point forming of hemispherical shapes. Using the response surface method, the effects of process parameters on wrinkling, maximum deviation from the target shape and thickness variation were investigated. The results show that elastic cushion with proper thickness and polyurethane with the hardness of Shore A90. It has also been found that the application of lubrication cans improve the shape accuracy of the formed workpiece. These final results were compared with the numerical simulation results of the multi-point forming for hemispherical shapes using a blank-holder and it was found that using cushion hardness realistic to reduce wrinkling and maximum deviation.

SU-E-T-154: Calculation of Tissue Dose Point Kernels Using GATE Monte Carlo Simulation Toolkit to Compare with Water Dose Point Kernel

Energy Technology Data Exchange (ETDEWEB)

Khazaee, M [shahid beheshti university, Tehran, Tehran (Iran, Islamic Republic of); Asl, A Kamali [Shahid Beheshti University, Tehran, Iran., Tehran, Tehran (Iran, Islamic Republic of); Geramifar, P [Shariati Hospital, Tehran, Iran., Tehran, Tehran (Iran, Islamic Republic of)

2015-06-15

Purpose: the objective of this study was to assess utilizing water dose point kernel (DPK)instead of tissue dose point kernels in convolution algorithms.to the best of our knowledge, in providing 3D distribution of absorbed dose from a 3D distribution of the activity, the human body is considered equivalent to water. as a Result tissue variations are not considered in patient specific dosimetry. Methods: In this study Gate v7.0 was used to calculate tissue dose point kernel. the beta emitter radionuclides which have taken into consideration in this simulation include Y-90, Lu-177 and P-32 which are commonly used in nuclear medicine. the comparison has been performed for dose point kernels of adipose, bone, breast, heart, intestine, kidney, liver, lung and spleen versus water dose point kernel. Results: In order to validate the simulation the Result of 90Y DPK in water were compared with published results of Papadimitroulas et al (Med. Phys., 2012). The results represented that the mean differences between water DPK and other soft tissues DPKs range between 0.6 % and 1.96% for 90Y, except for lung and bone, where the observed discrepancies are 6.3% and 12.19% respectively. The range of DPK difference for 32P is between 1.74% for breast and 18.85% for bone. For 177Lu, the highest difference belongs to bone which is equal to 16.91%. For other soft tissues the least discrepancy is observed in kidney with 1.68%. Conclusion: In all tissues except for lung and bone, the results of GATE for dose point kernel were comparable to water dose point kernel which demonstrates the appropriateness of applying water dose point kernel instead of soft tissues in the field of nuclear medicine.

SU-E-T-154: Calculation of Tissue Dose Point Kernels Using GATE Monte Carlo Simulation Toolkit to Compare with Water Dose Point Kernel

International Nuclear Information System (INIS)

Khazaee, M; Asl, A Kamali; Geramifar, P

2015-01-01

Purpose: the objective of this study was to assess utilizing water dose point kernel (DPK)instead of tissue dose point kernels in convolution algorithms.to the best of our knowledge, in providing 3D distribution of absorbed dose from a 3D distribution of the activity, the human body is considered equivalent to water. as a Result tissue variations are not considered in patient specific dosimetry. Methods: In this study Gate v7.0 was used to calculate tissue dose point kernel. the beta emitter radionuclides which have taken into consideration in this simulation include Y-90, Lu-177 and P-32 which are commonly used in nuclear medicine. the comparison has been performed for dose point kernels of adipose, bone, breast, heart, intestine, kidney, liver, lung and spleen versus water dose point kernel. Results: In order to validate the simulation the Result of 90Y DPK in water were compared with published results of Papadimitroulas et al (Med. Phys., 2012). The results represented that the mean differences between water DPK and other soft tissues DPKs range between 0.6 % and 1.96% for 90Y, except for lung and bone, where the observed discrepancies are 6.3% and 12.19% respectively. The range of DPK difference for 32P is between 1.74% for breast and 18.85% for bone. For 177Lu, the highest difference belongs to bone which is equal to 16.91%. For other soft tissues the least discrepancy is observed in kidney with 1.68%. Conclusion: In all tissues except for lung and bone, the results of GATE for dose point kernel were comparable to water dose point kernel which demonstrates the appropriateness of applying water dose point kernel instead of soft tissues in the field of nuclear medicine

Generalized transmissibilities for corner point rids in reservoir simulation; Transmissibilidades generalizadas em malhas corner point na simulacao de reservatorios

Energy Technology Data Exchange (ETDEWEB)

Tada, Mauricio P.; Silva, Antonio Fabio C. da; Maliska, Clovis R. [Universidade Federal de Santa Catarina (UFSC), Florianopolis, SC (Brazil). Lab. de Simulacao Numerica em Mecanica dos Fluidos e Transferencia de Calor (SINMEC)

2008-07-01

It is common to use five points schemes in reservoir simulation, since it simplifies the computational implementation and takes the linear resultant system simplest to be solved, giving more process velocity and robust to simulator. However, the use of these schemes may introduce significant errors solutions as function of volume's non orthogonality and medium anisotropy. These errors do not disappear with grid refinement since they are not truncate errors, it is a flux calculus approximation in the control volume faces. In order to get a correct solution and with no errors of these kind, it must be used a nine point scheme based on the correct flux calculus. The objective of this work is to present a new methodology to calculate the transmissibility on simulation reservoir that use a five and nine points scheme with corner-points grids. This mode considers full tensor anisotropy and the heterogeneity. The transmissibility presented are derivative of the discrete flux expression through control volume faces, where a generalized curvilinear coordinate system, located inner to the control volume, it is adopted. The transmissibility is then written on vector form and may to be used for any coordinate system. (author)

An Explicit Finite Difference scheme for numerical solution of fractional neutron point kinetic equation

International Nuclear Information System (INIS)

Saha Ray, S.; Patra, A.

2012-01-01

Highlights: ► In this paper fractional neutron point kinetic equation has been analyzed. ► The numerical solution for fractional neutron point kinetic equation is obtained. ► Explicit Finite Difference Method has been applied. ► Supercritical reactivity, critical reactivity and subcritical reactivity analyzed. ► Comparison between fractional and classical neutron density is presented. - Abstract: In the present article, a numerical procedure to efficiently calculate the solution for fractional point kinetics equation in nuclear reactor dynamics is investigated. The Explicit Finite Difference Method is applied to solve the fractional neutron point kinetic equation with the Grunwald–Letnikov (GL) definition (). Fractional Neutron Point Kinetic Model has been analyzed for the dynamic behavior of the neutron motion in which the relaxation time associated with a variation in the neutron flux involves a fractional order acting as exponent of the relaxation time, to obtain the best operation of a nuclear reactor dynamics. Results for neutron dynamic behavior for subcritical reactivity, supercritical reactivity and critical reactivity and also for different values of fractional order have been presented and compared with the classical neutron point kinetic (NPK) equation as well as the results obtained by the learned researchers .

Evaluation of melting point of UO2 by molecular dynamics simulation

International Nuclear Information System (INIS)

Arima, Tatsumi; Idemitsu, Kazuya; Inagaki, Yaohiro; Tsujita, Yuichi; Kinoshita, Motoyasu; Yakub, Eugene

2009-01-01

The melting point of UO 2 has been evaluated by molecular dynamics simulation (MD) in terms of interatomic potential, pressure and Schottky defect concentration. The Born-Mayer-Huggins potentials with or without a Morse potential were explored in the present study. Two-phase simulation whose supercell at the initial state consisted of solid and liquid phases gave the melting point comparable to the experimental data using the potential proposed by Yakub. The heat of fusion was determined by the difference in enthalpy at the melting point. In addition, MD calculations showed that the melting point increased with pressure applied to the system. Thus, the Clausius-Clapeyron equation was verified. Furthermore, MD calculations clarified that an addition of Schottky defects, which generated the local disorder in the UO 2 crystal, lowered the melting point.

The use of ZFP lossy floating point data compression in tornado-resolving thunderstorm simulations

Science.gov (United States)

Orf, L.

2017-12-01

In the field of atmospheric science, numerical models are used to produce forecasts of weather and climate and serve as virtual laboratories for scientists studying atmospheric phenomena. In both operational and research arenas, atmospheric simulations exploiting modern supercomputing hardware can produce a tremendous amount of data. During model execution, the transfer of floating point data from memory to the file system is often a significant bottleneck where I/O can dominate wallclock time. One way to reduce the I/O footprint is to compress the floating point data, which reduces amount of data saved to the file system. In this presentation we introduce LOFS, a file system developed specifically for use in three-dimensional numerical weather models that are run on massively parallel supercomputers. LOFS utilizes the core (in-memory buffered) HDF5 driver and includes compression options including ZFP, a lossy floating point data compression algorithm. ZFP offers several mechanisms for specifying the amount of lossy compression to be applied to floating point data, including the ability to specify the maximum absolute error allowed in each compressed 3D array. We explore different maximum error tolerances in a tornado-resolving supercell thunderstorm simulation for model variables including cloud and precipitation, temperature, wind velocity and vorticity magnitude. We find that average compression ratios exceeding 20:1 in scientifically interesting regions of the simulation domain produce visually identical results to uncompressed data in visualizations and plots. Since LOFS splits the model domain across many files, compression ratios for a given error tolerance can be compared across different locations within the model domain. We find that regions of high spatial variability (which tend to be where scientifically interesting things are occurring) show the lowest compression ratios, whereas regions of the domain with little spatial variability compress

Multi-scale calculation based on dual domain material point method combined with molecular dynamics

Energy Technology Data Exchange (ETDEWEB)

Dhakal, Tilak Raj [Los Alamos National Lab. (LANL), Los Alamos, NM (United States)

2017-02-27

This dissertation combines the dual domain material point method (DDMP) with molecular dynamics (MD) in an attempt to create a multi-scale numerical method to simulate materials undergoing large deformations with high strain rates. In these types of problems, the material is often in a thermodynamically non-equilibrium state, and conventional constitutive relations are often not available. In this method, the closure quantities, such as stress, at each material point are calculated from a MD simulation of a group of atoms surrounding the material point. Rather than restricting the multi-scale simulation in a small spatial region, such as phase interfaces, or crack tips, this multi-scale method can be used to consider non-equilibrium thermodynamic e ects in a macroscopic domain. This method takes advantage that the material points only communicate with mesh nodes, not among themselves; therefore MD simulations for material points can be performed independently in parallel. First, using a one-dimensional shock problem as an example, the numerical properties of the original material point method (MPM), the generalized interpolation material point (GIMP) method, the convected particle domain interpolation (CPDI) method, and the DDMP method are investigated. Among these methods, only the DDMP method converges as the number of particles increases, but the large number of particles needed for convergence makes the method very expensive especially in our multi-scale method where we calculate stress in each material point using MD simulation. To improve DDMP, the sub-point method is introduced in this dissertation, which provides high quality numerical solutions with a very small number of particles. The multi-scale method based on DDMP with sub-points is successfully implemented for a one dimensional problem of shock wave propagation in a cerium crystal. The MD simulation to calculate stress in each material point is performed in GPU using CUDA to accelerate the

Numerical methods for polyline-to-point-cloud registration with applications to patient-specific stent reconstruction.

Science.gov (United States)

Lin, Claire Yilin; Veneziani, Alessandro; Ruthotto, Lars

2018-03-01

We present novel numerical methods for polyline-to-point-cloud registration and their application to patient-specific modeling of deployed coronary artery stents from image data. Patient-specific coronary stent reconstruction is an important challenge in computational hemodynamics and relevant to the design and improvement of the prostheses. It is an invaluable tool in large-scale clinical trials that computationally investigate the effect of new generations of stents on hemodynamics and eventually tissue remodeling. Given a point cloud of strut positions, which can be extracted from images, our stent reconstruction method aims at finding a geometrical transformation that aligns a model of the undeployed stent to the point cloud. Mathematically, we describe the undeployed stent as a polyline, which is a piecewise linear object defined by its vertices and edges. We formulate the nonlinear registration as an optimization problem whose objective function consists of a similarity measure, quantifying the distance between the polyline and the point cloud, and a regularization functional, penalizing undesired transformations. Using projections of points onto the polyline structure, we derive novel distance measures. Our formulation supports most commonly used transformation models including very flexible nonlinear deformations. We also propose 2 regularization approaches ensuring the smoothness of the estimated nonlinear transformation. We demonstrate the potential of our methods using an academic 2D example and a real-life 3D bioabsorbable stent reconstruction problem. Our results show that the registration problem can be solved to sufficient accuracy within seconds using only a few number of Gauss-Newton iterations. Copyright © 2017 John Wiley & Sons, Ltd.

Dynamical simulation of a linear sigma model near the critical point

Energy Technology Data Exchange (ETDEWEB)

Wesp, Christian; Meistrenko, Alex; Greiner, Carsten [Institut fuer Theoretische Physik, Goethe-Universitaet Frankfurt, Max-von-Laue-Strasse 1, D-60438 Frankfurt (Germany); Hees, Hendrik van [Frankfurt Institute for Advanced Studies, Ruth-Moufang-Strasse 1, D-60438 Frankfurt (Germany)

2014-07-01

The intention of this study is the search for signatures of the chiral phase transition. To investigate the impact of fluctuations, e.g. of the baryon number, on the transition or a critical point, the linear sigma model is treated in a dynamical 3+1D numerical simulation. Chiral fields are approximated as classical fields, quarks are described by quasi particles in a Vlasov equation. Additional dynamic is implemented by quark-quark and quark-sigma-field interaction. For a consistent description of field-particle interactions, a new Monte-Carlo-Langevin-like formalism has been developed and is discussed.

Numerical solutions of a three-point boundary value problem with an ...

African Journals Online (AJOL)

Numerical solutions of a three-point boundary value problem with an integral condition for a third-order partial differential equation by using Laplace transform method Solutions numeriques d'un probleme pour une classe d'equations differentielles d'ordr.

Investigation of a novel dew point indirect evaporative air conditioning system for buildings

OpenAIRE

Duan, Zhiyin

2011-01-01

This study aims to improve the performance of existing indirect evaporative coolers. A new dew point indirect evaporative cooler with counter-current heat/mass exchanger was developed in this research by optimal design, material selection, numerical simulation, experimental investigations and economic, environmental, regional acceptance analysis. A new dew point heat/mass exchanger using a counter-current flow pattern was designed by numerical simulation in terms of material, structure, g...

Simulation of the Nonlinear Dose Dependence of Stabilized Point Defects

International Nuclear Information System (INIS)

Chen, R; Pagonis, V; Lawless, J L

2010-01-01

The dose dependence of the concentration of point defects in alkali-halides as well as other crystals, as exhibited by the dependence of the thermoluminescence (TL), optical absorption and ESR on the dose of non-ionizing UV excitation is studied using numerical simulation. The relevant set of coupled rate equations are first written and plausible sets of trapping parameters are chosen. Instead of using simplifying assumptions previously used for reaching conclusions concerning this dose behavior, exact numerical solutions have now been reached. Depending on the parameters chosen, different dose dependencies are seen. In some cases, linear dose dependence is reached in a broad range. Sublinear dose dependence, close to a D 1/2 dependence when D is the dose of excitation can be reached when retrapping is stronger than trapping in other traps stabilizing the defects. When strong competition between stabilizing traps takes place, an initial linear range is observed followed by strong superlinearity and an approach to saturation. All these behaviors have been observed experimentally in TL measurements as well as ESR and optical absorption in different materials. Similarities and dissimilarities to linear and non-linear dose dependencies obtained experimentally and by simulations when ionizing irradiation is used for excitation are discussed.

Simulation of the accumulation kinetics for radiation point defects in a metals with impurity

International Nuclear Information System (INIS)

Iskakov, B.M.; Nurova, A.B.

2001-01-01

In the work a kinetics of vacancies (V) and interstitial atoms (IA) accumulation for cases when the V and IA are recombining with each other, absorbing by drain and capturing by impurity atoms has been simulated. The differential equations system numerical solution was carried out by the Runge-Kutta method. The dynamical equilibrium time achievement for the point radiation defects accumulation process in the metal with impurity is considered

Theoretical, experimental and numerical diagnose of critical power point of thermoelectric generators

DEFF Research Database (Denmark)

Chen, Min; Gao, Xin

2014-01-01

of the critical power point in the series and parallel TEM arrays. Secondly, experiments of a series-parallel hybrid interconnected TEG are presented to clearly quantify the theoretical analyses. Finally, the hierarchical simulation, based on the SPICE (simulation program with integrated circuit emphasis...

Optimal Point-to-Point Trajectory Tracking of Redundant Manipulators using Generalized Pattern Search

Directory of Open Access Journals (Sweden)

Thi Rein Myo

2008-11-01

Full Text Available Optimal point-to-point trajectory planning for planar redundant manipulator is considered in this study. The main objective is to minimize the sum of the position error of the end-effector at each intermediate point along the trajectory so that the end-effector can track the prescribed trajectory accurately. An algorithm combining Genetic Algorithm and Pattern Search as a Generalized Pattern Search GPS is introduced to design the optimal trajectory. To verify the proposed algorithm, simulations for a 3-D-O-F planar manipulator with different end-effector trajectories have been carried out. A comparison between the Genetic Algorithm and the Generalized Pattern Search shows that The GPS gives excellent tracking performance.

Point vortex description of drift wave vortices: Dynamics and transport

International Nuclear Information System (INIS)

Kono, M.; Horton, W.

1991-05-01

Point-vortex description for drift wave vortices is formulated based on the Hasegawa-Mima equation to study elementary processes for the interactions of vortices as well as statistical properties like vortex diffusion. Dynamical properties of drift wave vortices known by numerical experiments are recovered. Furthermore a vortex diffusion model discussed by Horton based on numerical simulations is shown to be analytically obtained. A variety of phenomena arising from the short-range nature of the interaction force of point vortices are suggested. 12 refs., 10 figs

The Multiscale Material Point Method for Simulating Transient Responses

Science.gov (United States)

Chen, Zhen; Su, Yu-Chen; Zhang, Hetao; Jiang, Shan; Sewell, Thomas

2015-06-01

To effectively simulate multiscale transient responses such as impact and penetration without invoking master/slave treatment, the multiscale material point method (Multi-MPM) is being developed in which molecular dynamics at nanoscale and dissipative particle dynamics at mesoscale might be concurrently handled within the framework of the original MPM at microscale (continuum level). The proposed numerical scheme for concurrently linking different scales is described in this paper with simple examples for demonstration. It is shown from the preliminary study that the mapping and re-mapping procedure used in the original MPM could coarse-grain the information at fine scale and that the proposed interfacial scheme could provide a smooth link between different scales. Since the original MPM is an extension from computational fluid dynamics to solid dynamics, the proposed Multi-MPM might also become robust for dealing with multiphase interactions involving failure evolution. This work is supported in part by DTRA and NSFC.

Throughput analysis of point-to-multi-point hybric FSO/RF network

KAUST Repository

Rakia, Tamer

2017-07-31

This paper presents and analyzes a point-to-multi-point (P2MP) network that uses a number of free-space optical (FSO) links for data transmission from the central node to the different remote nodes. A common backup radio-frequency (RF) link is used by the central node for data transmission to any remote node in case of the failure of any one of FSO links. We develop a cross-layer Markov chain model to study the throughput from central node to a tagged remote node. Numerical examples are presented to compare the performance of the proposed P2MP hybrid FSO/RF network with that of a P2MP FSO-only network and show that the P2MP Hybrid FSO/RF network achieves considerable performance improvement over the P2MP FSO-only network.

Computing bubble-points of CO

NARCIS (Netherlands)

Ramdin, M.; Balaji, S.P.; Vicent Luna, J.M.; Torres-Knoop, A; Chen, Q.; Dubbeldam, D.; Calero, S; de Loos, T.W.; Vlugt, T.J.H.

2016-01-01

Computing bubble-points of multicomponent mixtures using Monte Carlo simulations is a non-trivial task. A new method is used to compute gas compositions from a known temperature, bubble-point pressure, and liquid composition. Monte Carlo simulations are used to calculate the bubble-points of

Experimental and numerical investigation of electrohydrodynamic flow in a point-to-ring corona discharge

Science.gov (United States)

Guan, Yifei; Vaddi, Ravi Sankar; Aliseda, Alberto; Novosselov, Igor

2018-04-01

An electrohydrodynamic (EHD) flow in a point-to-ring corona configuration is investigated experimentally and via a multiphysics computational model. The model couples the ion transport equation and the Navier-Stokes equations (NSE) to solve for the spatiotemporal distribution of electric field, flow field, and charge density. The numerical simulation results are validated against experimental measurements of the cathode voltage, ion concentration, and velocity profiles. The maximum flow velocity is at the centerline, and it decays rapidly with radial distance due to the viscous and electric forces acting on the partially ionized gas. To understand this coupling, a nondimensional parameter, X , is formulated as the ratio of the local electric force to the inertial term in the NSE. In the region of X ≥1 , the electric force dominates the flow dynamics, while in the X ≪1 region, the balance of viscous and inertial terms yields traditional pipe flow characteristics. This approach expands on the analytical model of Guan et al. by adding a description of the developing flow region. The approach allows the model to be used for the entire EHD domain, providing insights into the near-field flow in the corona region.

Statistical properties of several models of fractional random point processes

Science.gov (United States)

Bendjaballah, C.

2011-08-01

Statistical properties of several models of fractional random point processes have been analyzed from the counting and time interval statistics points of view. Based on the criterion of the reduced variance, it is seen that such processes exhibit nonclassical properties. The conditions for these processes to be treated as conditional Poisson processes are examined. Numerical simulations illustrate part of the theoretical calculations.

Comments on numerical simulations

International Nuclear Information System (INIS)

Sato, T.

1984-01-01

The author comments on a couple of things about numerical simulation. One is just about the philosophical discussion that is, spontaneous or driven. The other thing is the numerical or technical one. Frankly, the author didn't want to touch on the technical matter because this should be a common sense one for those who are working at numerical simulation. But since many people take numerical simulation results at their face value, he would like to remind you of the reality hidden behind them. First, he would point out that the meaning of ''driven'' in driven reconnection is different from that defined by Schindler or Akasofu. The author's definition is closer to Axford's definition. In the spontaneous case, for some unpredicted reason an excess energy of the system is suddenly released at a certain point. However, one does not answer how such an unstable state far beyond a stable limit is realized in the magnetotail. In the driven case, there is a definite energy buildup phase starting from a stable state; namely, energy in the black box increases from a stable level subject to an external source. When the state has reached a certain position, the energy is released suddenly. The difference between driven and spontaneous is whether the cause (plasma flow) to trigger reconnection is specified or reconnection is triggered unpredictably. Another difference is that in driven reconnection the reconnection rate is dependent on the speed of the external plasma flow, but in spontaneous reconnection the rate is dependent on the internal condition such as the resistivity

Three-Dimensional Adaptive Mesh Refinement Simulations of Point-Symmetric Nebulae

NARCIS (Netherlands)

Rijkhorst, E.-J.; Icke, V.; Mellema, G.; Meixner, M.; Kastner, J.H.; Balick, B.; Soker, N.

2004-01-01

Previous analytical and numerical work shows that the generalized interacting stellar winds model can explain the observed bipolar shapes of planetary nebulae very well. However, many circumstellar nebulae have a multipolar or point-symmetric shape. With two-dimensional calculations, Icke showed

Multiple point statistical simulation using uncertain (soft) conditional data

Science.gov (United States)

Hansen, Thomas Mejer; Vu, Le Thanh; Mosegaard, Klaus; Cordua, Knud Skou

2018-05-01

Geostatistical simulation methods have been used to quantify spatial variability of reservoir models since the 80s. In the last two decades, state of the art simulation methods have changed from being based on covariance-based 2-point statistics to multiple-point statistics (MPS), that allow simulation of more realistic Earth-structures. In addition, increasing amounts of geo-information (geophysical, geological, etc.) from multiple sources are being collected. This pose the problem of integration of these different sources of information, such that decisions related to reservoir models can be taken on an as informed base as possible. In principle, though difficult in practice, this can be achieved using computationally expensive Monte Carlo methods. Here we investigate the use of sequential simulation based MPS simulation methods conditional to uncertain (soft) data, as a computational efficient alternative. First, it is demonstrated that current implementations of sequential simulation based on MPS (e.g. SNESIM, ENESIM and Direct Sampling) do not account properly for uncertain conditional information, due to a combination of using only co-located information, and a random simulation path. Then, we suggest two approaches that better account for the available uncertain information. The first make use of a preferential simulation path, where more informed model parameters are visited preferentially to less informed ones. The second approach involves using non co-located uncertain information. For different types of available data, these approaches are demonstrated to produce simulation results similar to those obtained by the general Monte Carlo based approach. These methods allow MPS simulation to condition properly to uncertain (soft) data, and hence provides a computationally attractive approach for integration of information about a reservoir model.

Definition of distance for nonlinear time series analysis of marked point process data

Energy Technology Data Exchange (ETDEWEB)

Iwayama, Koji, E-mail: koji@sat.t.u-tokyo.ac.jp [Research Institute for Food and Agriculture, Ryukoku Univeristy, 1-5 Yokotani, Seta Oe-cho, Otsu-Shi, Shiga 520-2194 (Japan); Hirata, Yoshito; Aihara, Kazuyuki [Institute of Industrial Science, The University of Tokyo, 4-6-1 Komaba, Meguro-ku, Tokyo 153-8505 (Japan)

2017-01-30

Marked point process data are time series of discrete events accompanied with some values, such as economic trades, earthquakes, and lightnings. A distance for marked point process data allows us to apply nonlinear time series analysis to such data. We propose a distance for marked point process data which can be calculated much faster than the existing distance when the number of marks is small. Furthermore, under some assumptions, the Kullback–Leibler divergences between posterior distributions for neighbors defined by this distance are small. We performed some numerical simulations showing that analysis based on the proposed distance is effective. - Highlights: • A new distance for marked point process data is proposed. • The distance can be computed fast enough for a small number of marks. • The method to optimize parameter values of the distance is also proposed. • Numerical simulations indicate that the analysis based on the distance is effective.

A Numerical Investigation of CFRP-Steel Interfacial Failure with Material Point Method

International Nuclear Information System (INIS)

Shen Luming; Faleh, Haydar; Al-Mahaidi, Riadh

2010-01-01

The success of retrofitting steel structures by using the Carbon Fibre Reinforced Polymers (CFRP) significantly depends on the performance and integrity of CFRP-steel joint and the effectiveness of the adhesive used. Many of the previous numerical studies focused on the design and structural performance of the CFRP-steel system and neglected the mechanical responses of adhesive layer, which results in the lack of understanding in how the adhesive layer between the CFRP and steel performs during the loading and failure stages. Based on the recent observation on the failure of CFRP-steel bond in the double lap shear tests, a numerical approach is proposed in this study to simulate the delamination process of CFRP sheet from steel plate using the Material Point Method (MPM). In the proposed approach, an elastoplasticity model with a linear hardening and softening law is used to model the epoxy layer. The MPM, which does not employ fixed mesh-connectivity, is employed as a robust spatial discretization method to accommodate the multi-scale discontinuities involved in the CFRP-steel bond failure process. To demonstrate the potential of the proposed approach, a parametric study is conducted to investigate the effects of bond length and loading rates on the capacity and failure modes of CFRP-steel system. The evolution of the CFRP-steel bond failure and the distribution of stress and strain along bond length direction will be presented. The simulation results not only well match the available experimental data but also provide a better understanding on the physics behind the CFRP sheet delamination process.

Numerical simulation of water quality in Yangtze Estuary

Directory of Open Access Journals (Sweden)

Xi Li

2009-12-01

Full Text Available In order to monitor water quality in the Yangtze Estuary, water samples were collected and field observation of current and velocity stratification was carried out using a shipboard acoustic Doppler current profiler (ADCP. Results of two representative variables, the temporal and spatial variation of new point source sewage discharge as manifested by chemical oxygen demand (COD and the initial water quality distribution as manifested by dissolved oxygen (DO, were obtained by application of the Environmental Fluid Dynamics Code (EFDC with solutions for hydrodynamics during tides. The numerical results were compared with field data, and the field data provided verification of numerical application: this numerical model is an effective tool for water quality simulation. For point source discharge, COD concentration was simulated with an initial value in the river of zero. The simulated increments and distribution of COD in the water show acceptable agreement with field data. The concentration of DO is much higher in the North Branch than in the South Branch due to consumption of oxygen in the South Branch resulting from discharge of sewage from Shanghai. The DO concentration is greater in the surface layer than in the bottom layer. The DO concentration is low in areas with a depth of less than 20 m, and high in areas between the 20-m and 30-m isobaths. It is concluded that the numerical model is valuable in simulation of water quality in the case of specific point source pollutant discharge. The EFDC model is also of satisfactory accuracy in water quality simulation of the Yangtze Estuary.

Comparison Between Numerical Modeling and Experimental Testing of a Point Absorber WEC Using Linear Power Take-Off System

DEFF Research Database (Denmark)

Zurkinden, Andrew Stephen; Kramer, Morten; Sichani, Mahdi Teimouri

2012-01-01

systems into a central wave to wire model. The power production then depends on the control strategy which is applied to the device. The objective of this paper is to develop numerical methods for motion analysis of marine structures with a special emphasis on wave energy converters. Two different time...... domain models are applied to a point absorber system working in pitch mode only. The device is similar to the well-known Wavestar prototype located in the Danish North Sea. A laboratory model has been set up in order to validate the numerical simulations of the dynamics. Wave Excitation force...... interaction assumption and linearized equation of motion. The region over which the numerical model is valid will be presented in terms of non-dimensional parameters describing the different wave states....

Numerical analysis for multi-group neutron-diffusion equation using Radial Point Interpolation Method (RPIM)

International Nuclear Information System (INIS)

Kim, Kyung-O; Jeong, Hae Sun; Jo, Daeseong

2017-01-01

Highlights: • Employing the Radial Point Interpolation Method (RPIM) in numerical analysis of multi-group neutron-diffusion equation. • Establishing mathematical formation of modified multi-group neutron-diffusion equation by RPIM. • Performing the numerical analysis for 2D critical problem. - Abstract: A mesh-free method is introduced to overcome the drawbacks (e.g., mesh generation and connectivity definition between the meshes) of mesh-based (nodal) methods such as the finite-element method and finite-difference method. In particular, the Point Interpolation Method (PIM) using a radial basis function is employed in the numerical analysis for the multi-group neutron-diffusion equation. The benchmark calculations are performed for the 2D homogeneous and heterogeneous problems, and the Multiquadrics (MQ) and Gaussian (EXP) functions are employed to analyze the effect of the radial basis function on the numerical solution. Additionally, the effect of the dimensionless shape parameter in those functions on the calculation accuracy is evaluated. According to the results, the radial PIM (RPIM) can provide a highly accurate solution for the multiplication eigenvalue and the neutron flux distribution, and the numerical solution with the MQ radial basis function exhibits the stable accuracy with respect to the reference solutions compared with the other solution. The dimensionless shape parameter directly affects the calculation accuracy and computing time. Values between 1.87 and 3.0 for the benchmark problems considered in this study lead to the most accurate solution. The difference between the analytical and numerical results for the neutron flux is significantly increased in the edge of the problem geometry, even though the maximum difference is lower than 4%. This phenomenon seems to arise from the derivative boundary condition at (x,0) and (0,y) positions, and it may be necessary to introduce additional strategy (e.g., the method using fictitious points and

Point-to-point people with purpose—Exploring the possibility of a commercial traveler market for point-to-point suborbital space transportation

Science.gov (United States)

Webber, Derek

2013-12-01

An argument was made at the First Arcachon Conference on Private Human Access to Space in 2008 [1] that some systematic market research should be conducted into potential market segments for point-to-point suborbital space transportation (PtP), in order to understand whether a commercial market exists which might augment possible government use for such a vehicle. The cargo market potential was subsequently addressed via desk research, and the results, which resulted in a pessimistic business case outlook, were presented in [2]. The same desk research approach is now used in this paper to address the potential business and wealthy individual passenger traveler market segment ("point-to-point people with purpose"). The results, with the assumed ticket pricing, are not encouraging.

Turbulence simulations of X point physics on the L-H transitions

International Nuclear Information System (INIS)

Xu, X.Q.; Cohen, R.H.; Nevins, W.M.

2001-01-01

We show that the resistive X-point mode is the dominant mode in boundary plasmas in X-point divertor geometry. The poloidal fluctuation phase velocity from the simulation results of the resistive X-point turbulence shows experimentally measured structure across separatrix in many fusion devices. The fluctuation phase velocity is larger than ExB velocity both in L and H mode phases. We also demonstrate that there is a strong poloidal asymmetry of particle flux in the proximity of the separatrix. Turbulence suppression in the L-H transition results when sources of energy and particles drive sufficient gradients as in the experiments. (author)

Ewald Electrostatics for Mixtures of Point and Continuous Line Charges.

Science.gov (United States)

Antila, Hanne S; Tassel, Paul R Van; Sammalkorpi, Maria

2015-10-15

Many charged macro- or supramolecular systems, such as DNA, are approximately rod-shaped and, to the lowest order, may be treated as continuous line charges. However, the standard method used to calculate electrostatics in molecular simulation, the Ewald summation, is designed to treat systems of point charges. We extend the Ewald concept to a hybrid system containing both point charges and continuous line charges. We find the calculated force between a point charge and (i) a continuous line charge and (ii) a discrete line charge consisting of uniformly spaced point charges to be numerically equivalent when the separation greatly exceeds the discretization length. At shorter separations, discretization induces deviations in the force and energy, and point charge-point charge correlation effects. Because significant computational savings are also possible, the continuous line charge Ewald method presented here offers the possibility of accurate and efficient electrostatic calculations.

APPLICABILITY ANALYSIS OF CLOTH SIMULATION FILTERING ALGORITHM FOR MOBILE LIDAR POINT CLOUD

Directory of Open Access Journals (Sweden)

S. Cai

2018-04-01

Full Text Available Classifying the original point clouds into ground and non-ground points is a key step in LiDAR (light detection and ranging data post-processing. Cloth simulation filtering (CSF algorithm, which based on a physical process, has been validated to be an accurate, automatic and easy-to-use algorithm for airborne LiDAR point cloud. As a new technique of three-dimensional data collection, the mobile laser scanning (MLS has been gradually applied in various fields, such as reconstruction of digital terrain models (DTM, 3D building modeling and forest inventory and management. Compared with airborne LiDAR point cloud, there are some different features (such as point density feature, distribution feature and complexity feature for mobile LiDAR point cloud. Some filtering algorithms for airborne LiDAR data were directly used in mobile LiDAR point cloud, but it did not give satisfactory results. In this paper, we explore the ability of the CSF algorithm for mobile LiDAR point cloud. Three samples with different shape of the terrain are selected to test the performance of this algorithm, which respectively yields total errors of 0.44 %, 0.77 % and1.20 %. Additionally, large area dataset is also tested to further validate the effectiveness of this algorithm, and results show that it can quickly and accurately separate point clouds into ground and non-ground points. In summary, this algorithm is efficient and reliable for mobile LiDAR point cloud.

Numerical methods in simulation of resistance welding

DEFF Research Database (Denmark)

Nielsen, Chris Valentin; Martins, Paulo A.F.; Zhang, Wenqi

2015-01-01

Finite element simulation of resistance welding requires coupling betweenmechanical, thermal and electrical models. This paper presents the numerical models and theircouplings that are utilized in the computer program SORPAS. A mechanical model based onthe irreducible flow formulation is utilized...... a resistance welding point of view, the most essential coupling between the above mentioned models is the heat generation by electrical current due to Joule heating. The interaction between multiple objects is anothercritical feature of the numerical simulation of resistance welding because it influences...... thecontact area and the distribution of contact pressure. The numerical simulation of resistancewelding is illustrated by a spot welding example that includes subsequent tensile shear testing...

May 2002 Lidar Point Data of Southern California Coastline: Dana Point to Point La Jolla

Data.gov (United States)

National Oceanic and Atmospheric Administration, Department of Commerce — This data set contains lidar point data from a strip of Southern California coastline (including water, beach, cliffs, and top of cliffs) from Dana Point to Point La...

September 2002 Lidar Point Data of Southern California Coastline: Dana Point to Point La Jolla

Data.gov (United States)

National Oceanic and Atmospheric Administration, Department of Commerce — This data set contains lidar point data from a strip of Southern California coastline (including water, beach, cliffs, and top of cliffs) from Dana Point to Point La...

Vortex pinning by point defect in superconductors

International Nuclear Information System (INIS)

Liao Hongyin; Zhou Shiping; Du Haochen

2003-01-01

We apply the periodic time-dependent Ginzburg-Landau model to study vortex distribution in type-II superconductors with a point-like defect and square pinning array. A defect site will pin vortices, and a periodic pinning array with right geometric parameters, which can be any form designed in advance, shapes the vortex pattern as external magnetic field varies. The maximum length over which an attractive interaction between a pinning centre and a vortex extends is estimated to be about 6.0ξ. We also derive spatial distribution expressions for the order parameter, vector potential, magnetic field and supercurrent induced by a point defect. Theoretical results and numerical simulations are compared with each other and they are consistent

Simulation of growing grains under orientation relation - dependent quadruple point dragging

International Nuclear Information System (INIS)

Ito, K

2015-01-01

The growth behaviour of a specified grain embedded in matrix grains, for which the migration mobility of the quadruple points depended on the relation between the orientations of the growing and shrinking grains, was studied using a modified Potts MC-type threedimensional simulation. Large embedded grains continued to grow without being overcome by coarsening matrix grains, whereas small embedded grains disappeared, under the influence of the relative mobilities of the quadruple points, the composition of the matrix grain texture and the width of the grain size distribution of the matrix grains. These results indicate that orientation relation-dependent quadruple point dragging can affect the recrystallization texture during the grain coarsening stage. (paper)

Free Energy, Enthalpy and Entropy from Implicit Solvent End-Point Simulations.

Science.gov (United States)

Fogolari, Federico; Corazza, Alessandra; Esposito, Gennaro

2018-01-01

Free energy is the key quantity to describe the thermodynamics of biological systems. In this perspective we consider the calculation of free energy, enthalpy and entropy from end-point molecular dynamics simulations. Since the enthalpy may be calculated as the ensemble average over equilibrated simulation snapshots the difficulties related to free energy calculation are ultimately related to the calculation of the entropy of the system and in particular of the solvent entropy. In the last two decades implicit solvent models have been used to circumvent the problem and to take into account solvent entropy implicitly in the solvation terms. More recently outstanding advancement in both implicit solvent models and in entropy calculations are making the goal of free energy estimation from end-point simulations more feasible than ever before. We review briefly the basic theory and discuss the advancements in light of practical applications.

Free Energy, Enthalpy and Entropy from Implicit Solvent End-Point Simulations

Directory of Open Access Journals (Sweden)

Federico Fogolari

2018-02-01

Full Text Available Free energy is the key quantity to describe the thermodynamics of biological systems. In this perspective we consider the calculation of free energy, enthalpy and entropy from end-point molecular dynamics simulations. Since the enthalpy may be calculated as the ensemble average over equilibrated simulation snapshots the difficulties related to free energy calculation are ultimately related to the calculation of the entropy of the system and in particular of the solvent entropy. In the last two decades implicit solvent models have been used to circumvent the problem and to take into account solvent entropy implicitly in the solvation terms. More recently outstanding advancement in both implicit solvent models and in entropy calculations are making the goal of free energy estimation from end-point simulations more feasible than ever before. We review briefly the basic theory and discuss the advancements in light of practical applications.

Numerical simulation support to the ESA/THOR mission

Science.gov (United States)

Valentini, F.; Servidio, S.; Perri, S.; Perrone, D.; De Marco, R.; Marcucci, M. F.; Daniele, B.; Bruno, R.; Camporeale, E.

2016-12-01

THOR is a spacecraft concept currently undergoing study phase as acandidate for the next ESA medium size mission M4. THOR has been designedto solve the longstanding physical problems of particle heating andenergization in turbulent plasmas. It will provide high resolutionmeasurements of electromagnetic fields and particle distribution functionswith unprecedented resolution, with the aim of exploring the so-calledkinetic scales. We present the numerical simulation framework which is supporting the THOR mission during the study phase. The THOR teamincludes many scientists developing and running different simulation codes(Eulerian-Vlasov, Particle-In-Cell, Gyrokinetics, Two-fluid, MHD, etc.),addressing the physics of plasma turbulence, shocks, magnetic reconnectionand so on.These numerical codes are being used during the study phase, mainly withthe aim of addressing the following points:(i) to simulate the response of real particle instruments on board THOR, byemploying an electrostatic analyser simulator which mimics the response ofthe CSW, IMS and TEA instruments to the particle velocity distributions ofprotons, alpha particle and electrons, as obtained from kinetic numericalsimulations of plasma turbulence.(ii) to compare multi-spacecraft with single-spacecraft configurations inmeasuring current density, by making use of both numerical models ofsynthetic turbulence and real data from MMS spacecraft.(iii) to investigate the validity of the Taylor hypothesis indifferent configurations of plasma turbulence

Mathematical simulation of point defect interaction with grain boundaries

International Nuclear Information System (INIS)

Bojko, V.S.

1987-01-01

Published works, where the interaction of point defects and grain boundaries was studied by mathematical simulation methods, have been analysed. Energetics of the vacancy formation both in nuclei of large-angle special grain boundaries and in lattice regions adjoining them has been considered. The data obtained permit to explain specific features of grain-boundary diffusion processes. Results of mathematical simulation of the interaction of impurity atoms and boundaries have been considered. Specific features of the helium atom interaction with large-angle grain boundaries are analysed as well

Fixed-Rate Compressed Floating-Point Arrays.

Science.gov (United States)

Lindstrom, Peter

2014-12-01

Current compression schemes for floating-point data commonly take fixed-precision values and compress them to a variable-length bit stream, complicating memory management and random access. We present a fixed-rate, near-lossless compression scheme that maps small blocks of 4(d) values in d dimensions to a fixed, user-specified number of bits per block, thereby allowing read and write random access to compressed floating-point data at block granularity. Our approach is inspired by fixed-rate texture compression methods widely adopted in graphics hardware, but has been tailored to the high dynamic range and precision demands of scientific applications. Our compressor is based on a new, lifted, orthogonal block transform and embedded coding, allowing each per-block bit stream to be truncated at any point if desired, thus facilitating bit rate selection using a single compression scheme. To avoid compression or decompression upon every data access, we employ a software write-back cache of uncompressed blocks. Our compressor has been designed with computational simplicity and speed in mind to allow for the possibility of a hardware implementation, and uses only a small number of fixed-point arithmetic operations per compressed value. We demonstrate the viability and benefits of lossy compression in several applications, including visualization, quantitative data analysis, and numerical simulation.

Point-of-care echocardiography in simulation-based education and assessment

Directory of Open Access Journals (Sweden)

Amini R

2016-05-01

Full Text Available Richard Amini, Lori A Stolz, Parisa P Javedani, Kevin Gaskin, Nicola Baker, Vivienne Ng, Srikar Adhikari Department of Emergency Medicine, University of Arizona Medical Center, Tucson, AZ, USA Background: Emergency medicine milestones released by the Accreditation Council for Graduate Medical Education require residents to demonstrate competency in bedside ultrasound (US. The acquisition of these skills necessitates a combination of exposure to clinical pathology, hands-on US training, and feedback. Objectives: We describe a novel simulation-based educational and assessment tool designed to evaluate emergency medicine residents’ competency in point-of-care echocardiography for evaluation of a hypotensive patient with chest pain using bedside US. Methods: This was a cross-sectional study conducted at an academic medical center. A simulation-based module was developed to teach and assess the use of point-of-care echocardiography in the evaluation of the hypotensive patient. The focus of this module was sonographic imaging of cardiac pathology, and this focus was incorporated in all components of the session: asynchronous learning, didactic lecture, case-based learning, and hands-on stations. Results: A total of 52 residents with varying US experience participated in this study. Questions focused on knowledge assessment demonstrated improvement across the postgraduate year (PGY of training. Objective standardized clinical examination evaluation demonstrated improvement between PGY I and PGY III; however, it was noted that there was a small dip in hands-on scanning skills during the PGY II. Clinical diagnosis and management skills also demonstrated incremental improvement across the PGY of training. Conclusion: The 1-day, simulation-based US workshop was an effective educational and assessment tool at our institution. Keywords: point-of care ultrasound, simulation education

An experimental and numerical study of diffusion flames in cross-flow and quiescent environment at smoke point condition

Science.gov (United States)

Goh, Sien Fong

An experimental and numerical study of a turbulent smoke point diffusion flame in a quiescent and cross-flow condition was performed. The fuel mass flow rate of a turbulent smoke point flame was determined at a quiescent condition and in cross-flow with velocity ranging from 2 to 4 m/s. This fuel mass flow rate is defined as the Critical Fuel Mass Flow Rate (CFMFR). At a fuel mass flow rate below the CFMFR the flame produces smoke. In the dilution study, an amount of inert gas (nitrogen) was added to the fuel stream to achieve the smoke point condition for ten different fractions of CFMFR. From this dilution study, three regions were defined, the chemically-dominated region, transition region, and momentum-dominated region. The first objective of this study was to determine the factors behind the distinction of these three regions. The second objective was to understand the effect of cross-flow velocity on the smoke point flame structure. The flame temperature, radiation, geometrical dimension of flame, velocity, and global emissions and in-flame species concentration were measured. The third objective was to study a numerical model that can simulate the turbulent smoke point flame structure. The dilution study showed that the flames in quiescent condition and in the 3.5 and 4 m/s cross-flow condition had the chemically-dominated region at 5% to 20% CFMFR, the transition region at 20% to 40% CFMFR, and the momentum-dominated region at 40% to 100% CFMFR. On the other hand, the flame in cross-flow of 2 to 3 m/s showed the chemically-dominated region at 5% to 10% CFMFR, the transition region at 10% to 30% CFMFR, and the momentum-dominated region at 30% to 100% CFMFR. The chemically-dominated flame had a sharp dual-peak structure for the flame temperature, CO2 and NO concentration profiles at 25% and 50% flame length. However, the momentum-dominated region flame exhibited a dual peak structure only at 25% flame length. The decrease of flow rate from 30% to 10% CFMFR

Numerical solution of the point reactor kinetics equations with fuel burn-up and temperature feedback

International Nuclear Information System (INIS)

Tashakor, S.; Jahanfarnia, G.; Hashemi-Tilehnoee, M.

2010-01-01

Point reactor kinetics equations are solved numerically using one group of delayed neutrons and with fuel burn-up and temperature feedback included. To calculate the fraction of one-group delayed neutrons, a group of differential equations are solved by an implicit time method. Using point reactor kinetics equations, changes in mean neutrons density, temperature, and reactivity are calculated in different times during the reactor operation. The variation of reactivity, temperature, and maximum power with time are compared with the predictions by other methods.

MAGNETIC NULL POINTS IN KINETIC SIMULATIONS OF SPACE PLASMAS

International Nuclear Information System (INIS)

Olshevsky, Vyacheslav; Innocenti, Maria Elena; Cazzola, Emanuele; Lapenta, Giovanni; Deca, Jan; Divin, Andrey; Peng, Ivy Bo; Markidis, Stefano

2016-01-01

We present a systematic attempt to study magnetic null points and the associated magnetic energy conversion in kinetic particle-in-cell simulations of various plasma configurations. We address three-dimensional simulations performed with the semi-implicit kinetic electromagnetic code iPic3D in different setups: variations of a Harris current sheet, dipolar and quadrupolar magnetospheres interacting with the solar wind, and a relaxing turbulent configuration with multiple null points. Spiral nulls are more likely created in space plasmas: in all our simulations except lunar magnetic anomaly (LMA) and quadrupolar mini-magnetosphere the number of spiral nulls prevails over the number of radial nulls by a factor of 3–9. We show that often magnetic nulls do not indicate the regions of intensive energy dissipation. Energy dissipation events caused by topological bifurcations at radial nulls are rather rare and short-lived. The so-called X-lines formed by the radial nulls in the Harris current sheet and LMA simulations are rather stable and do not exhibit any energy dissipation. Energy dissipation is more powerful in the vicinity of spiral nulls enclosed by magnetic flux ropes with strong currents at their axes (their cross sections resemble 2D magnetic islands). These null lines reminiscent of Z-pinches efficiently dissipate magnetic energy due to secondary instabilities such as the two-stream or kinking instability, accompanied by changes in magnetic topology. Current enhancements accompanied by spiral nulls may signal magnetic energy conversion sites in the observational data

Modelling point patterns with linear structures

DEFF Research Database (Denmark)

Møller, Jesper; Rasmussen, Jakob Gulddahl

2009-01-01

processes whose realizations contain such linear structures. Such a point process is constructed sequentially by placing one point at a time. The points are placed in such a way that new points are often placed close to previously placed points, and the points form roughly line shaped structures. We...... consider simulations of this model and compare with real data....

Modelling point patterns with linear structures

DEFF Research Database (Denmark)

Møller, Jesper; Rasmussen, Jakob Gulddahl

processes whose realizations contain such linear structures. Such a point process is constructed sequentially by placing one point at a time. The points are placed in such a way that new points are often placed close to previously placed points, and the points form roughly line shaped structures. We...... consider simulations of this model and compare with real data....

Assessing accuracy of point fire intervals across landscapes with simulation modelling

Science.gov (United States)

Russell A. Parsons; Emily K. Heyerdahl; Robert E. Keane; Brigitte Dorner; Joseph Fall

2007-01-01

We assessed accuracy in point fire intervals using a simulation model that sampled four spatially explicit simulated fire histories. These histories varied in fire frequency and size and were simulated on a flat landscape with two forest types (dry versus mesic). We used three sampling designs (random, systematic grids, and stratified). We assessed the sensitivity of...

Simulation of the Point Lepreau core-follow history with SORO

Energy Technology Data Exchange (ETDEWEB)

Shanes, F.C.; Olive, C.G.; Cheng, I. [Nuclear Safety Solutions Limited, Toronto, Ontario (Canada); Banica, C. [Ontario Power Generation, Ajax, Ontario (Canada); Newman, C. [NB Power Nuclear, Point Lepreau Generating Station, Lepreau, New Brunswick (Canada); Nainer, O. [Bruce Power, Toronto, Ontario (Canada)

2006-07-01

The core tracking SORO program calculates the flux, burn-ups, bundle and channel powers, and is used for Ontario Power Generation and Bruce Power reactors. This paper presents the results of a comparison of the SORO cell fluxes against the Point Lepreau CANDU-6 measured flux data. A SORO model was created for the Point Lepreau reactor, and simulations were carried out to compare the SORO fluxes to the Travelling Flux Detector (TFD) scans and one year of operating history. Considering the good agreement between measured and computed fluxes, the results provide confidence that SORO accurately calculates the cell fluxes and bundle powers. (author)

Simulation of the Point Lepreau core-follow history with SORO

International Nuclear Information System (INIS)

Shanes, F.C.; Olive, C.G.; Cheng, I.; Banica, C.; Newman, C.; Nainer, O.

2006-01-01

The core tracking SORO program calculates the flux, burn-ups, bundle and channel powers, and is used for Ontario Power Generation and Bruce Power reactors. This paper presents the results of a comparison of the SORO cell fluxes against the Point Lepreau CANDU-6 measured flux data. A SORO model was created for the Point Lepreau reactor, and simulations were carried out to compare the SORO fluxes to the Travelling Flux Detector (TFD) scans and one year of operating history. Considering the good agreement between measured and computed fluxes, the results provide confidence that SORO accurately calculates the cell fluxes and bundle powers. (author)

Research on point source simulating the γ-ray detection efficiencies of stander source

International Nuclear Information System (INIS)

Tian Zining; Jia Mingyan; Shen Maoquan; Yang Xiaoyan; Cheng Zhiwei

2010-01-01

For φ 75 mm x 25 mm sample, the full energy peak efficiencies on different heights of sample radius were obtained using the point sources, and the function parameters about the full energy peak efficiencies of point sources based on radius was fixed. The 59.54 keV γ-ray, 661.66 keV γ-ray, 1173.2 keV γ-ray, 1332.5 keV γ-ray detection efficiencies on different height of samples were obtained, based on the full energy peak efficiencies of point sources and its height, and the function parameters about the full energy peak efficiencies of surface sources based on sample height was fixed. The detection efficiency of (75 mm x 25 mm calibration source can be obtained by integrality, the detection efficiencies simulated by point sources are consistent with the results of stander source in 10%. Therefore, the calibration method of stander source can be substituted by the point source simulation method, and it tis feasible when there is no stander source.) (authors)

Strategies and limits in multi-stage single-point incremental forming

DEFF Research Database (Denmark)

Skjødt, Martin; Silva, M.B.; Martins, P. A. F.

2010-01-01

paths. The results also reveal that the sequence of multi-stage forming has a large effect on the location of strain points in the principal strain space. Strain paths are linear in the first stage and highly non-linear in the subsequent forming stages. The overall results show that the experimentally......Multi-stage single-point incremental forming (SPIF) is a state-of-the-art manufacturing process that allows small-quantity production of complex sheet metal parts with vertical walls. This paper is focused on the application of multi-stage SPIF with the objective of producing cylindrical cups......-limit curves and fracture forming-limit curves (FFLCs), numerical simulation, and experimentation, namely the evaluation of strain paths and fracture strains in actual multi-stage parts. Assessment of numerical simulation with experimentation shows good agreement between computed and measured strain and strain...

Zero-point energy conservation in classical trajectory simulations: Application to H2CO

Science.gov (United States)

Lee, Kin Long Kelvin; Quinn, Mitchell S.; Kolmann, Stephen J.; Kable, Scott H.; Jordan, Meredith J. T.

2018-05-01

A new approach for preventing zero-point energy (ZPE) violation in quasi-classical trajectory (QCT) simulations is presented and applied to H2CO "roaming" reactions. Zero-point energy may be problematic in roaming reactions because they occur at or near bond dissociation thresholds and these channels may be incorrectly open or closed depending on if, or how, ZPE has been treated. Here we run QCT simulations on a "ZPE-corrected" potential energy surface defined as the sum of the molecular potential energy surface (PES) and the global harmonic ZPE surface. Five different harmonic ZPE estimates are examined with four, on average, giving values within 4 kJ/mol—chemical accuracy—for H2CO. The local harmonic ZPE, at arbitrary molecular configurations, is subsequently defined in terms of "projected" Cartesian coordinates and a global ZPE "surface" is constructed using Shepard interpolation. This, combined with a second-order modified Shepard interpolated PES, V, allows us to construct a proof-of-concept ZPE-corrected PES for H2CO, Veff, at no additional computational cost to the PES itself. Both V and Veff are used to model product state distributions from the H + HCO → H2 + CO abstraction reaction, which are shown to reproduce the literature roaming product state distributions. Our ZPE-corrected PES allows all trajectories to be analysed, whereas, in previous simulations, a significant proportion was discarded because of ZPE violation. We find ZPE has little effect on product rotational distributions, validating previous QCT simulations. Running trajectories on V, however, shifts the product kinetic energy release to higher energy than on Veff and classical simulations of kinetic energy release should therefore be viewed with caution.

Numerical Simulation of Missouri River Bed Evolution Downstream of Gavins Point Dam

Science.gov (United States)

Sulaiman, Z. A.; Blum, M. D.; Lephart, G.; Viparelli, E.

2016-12-01

The Missouri River originates in the Rocky Mountains in western Montana and joins the Mississippi River near Saint Louis, Missouri. In the 1900s dam construction and river engineering works, such as river alignment, narrowing and bank protections were performed in the Missouri River basin to control the flood flows, ensure navigation and use the water for agricultural, industrial and municipal needs, for the production of hydroelectric power generation and for recreation. These projects altered the flow and the sediment transport regimes in the river and the exchange of sediment between the river and the adjoining floodplain. Here we focus on the long term effect of dam construction and channel narrowing on the 1200 km long reach of the Missouri River between Gavins Point Dam, Nebraska and South Dakota, and the confluence with the Mississippi River. Field observations show that two downstream migrating waves of channel bed degradation formed in this reach in response to the changes in flow regime, sediment load and channel geometry. We implemented a one dimensional morphodynamic model for large, low slope sand bed rivers, we validated the model at field scale by comparing the numerical results with the available field data and we use the model to 1) predict the magnitude and the migration rate of the waves of degradation at engineering time scales ( 150 years into the future), 2) quantify the changes in the sand load delivered to the Mississippi River, where field observations at Thebes, i.e. downstream of Saint Louis, suggest a decline in the mean annual sand load in the past 50 years, and 3) identify the role of the main tributaries - Little Sioux River, Platte River and Kansas River - on the wave migration speed and the annual sand load in the Missouri River main channel.

Teaching Experience: How to Make and Use PowerPoint-Based Interactive Simulations for Undergraduate IR Teaching

Science.gov (United States)

Meibauer, Gustav; Aagaard Nøhr, Andreas

2018-01-01

This article is about designing and implementing PowerPoint-based interactive simulations for use in International Relations (IR) introductory undergraduate classes based on core pedagogical literature, models of human skill acquisition, and previous research on simulations in IR teaching. We argue that simulations can be usefully employed at the…

Molecular dynamics simulation of a binary mixture near the lower critical point

Energy Technology Data Exchange (ETDEWEB)

Pousaneh, Faezeh; Edholm, Olle, E-mail: oed@kth.se [Theoretical Biological Physics, Department of Theoretical Physics, Royal Institute of Technology (KTH), AlbaNova University Center, SE-106 91 Stockholm (Sweden); Maciołek, Anna [Institute of Physical Chemistry, Polish Academy of Sciences, Kasprzaka 44/52, 01-224 Warsaw (Poland); Max-Planck-Institut für Intelligente Systeme, Heisenbergstrasse 3, D-70569 Stuttgart (Germany)

2016-07-07

2,6-lutidine molecules mix with water at high and low temperatures but in a wide intermediate temperature range a 2,6-lutidine/water mixture exhibits a miscibility gap. We constructed and validated an atomistic model for 2,6-lutidine and performed molecular dynamics simulations of 2,6-lutidine/water mixture at different temperatures. We determined the part of demixing curve with the lower critical point. The lower critical point extracted from our data is located close to the experimental one. The estimates for critical exponents obtained from our simulations are in a good agreement with the values corresponding to the 3D Ising universality class.

A point particle model of lightly bound skyrmions

Directory of Open Access Journals (Sweden)

Mike Gillard

2017-04-01

Full Text Available A simple model of the dynamics of lightly bound skyrmions is developed in which skyrmions are replaced by point particles, each carrying an internal orientation. The model accounts well for the static energy minimizers of baryon number 1≤B≤8 obtained by numerical simulation of the full field theory. For 9≤B≤23, a large number of static solutions of the point particle model are found, all closely resembling size B subsets of a face centred cubic lattice, with the particle orientations dictated by a simple colouring rule. Rigid body quantization of these solutions is performed, and the spin and isospin of the corresponding ground states extracted. As part of the quantization scheme, an algorithm to compute the symmetry group of an oriented point cloud, and to determine its corresponding Finkelstein–Rubinstein constraints, is devised.

Modeling the contribution of point sources and non-point sources to Thachin River water pollution.

Science.gov (United States)

Schaffner, Monika; Bader, Hans-Peter; Scheidegger, Ruth

2009-08-15

Major rivers in developing and emerging countries suffer increasingly of severe degradation of water quality. The current study uses a mathematical Material Flow Analysis (MMFA) as a complementary approach to address the degradation of river water quality due to nutrient pollution in the Thachin River Basin in Central Thailand. This paper gives an overview of the origins and flow paths of the various point- and non-point pollution sources in the Thachin River Basin (in terms of nitrogen and phosphorus) and quantifies their relative importance within the system. The key parameters influencing the main nutrient flows are determined and possible mitigation measures discussed. The results show that aquaculture (as a point source) and rice farming (as a non-point source) are the key nutrient sources in the Thachin River Basin. Other point sources such as pig farms, households and industries, which were previously cited as the most relevant pollution sources in terms of organic pollution, play less significant roles in comparison. This order of importance shifts when considering the model results for the provincial level. Crosschecks with secondary data and field studies confirm the plausibility of our simulations. Specific nutrient loads for the pollution sources are derived; these can be used for a first broad quantification of nutrient pollution in comparable river basins. Based on an identification of the sensitive model parameters, possible mitigation scenarios are determined and their potential to reduce the nutrient load evaluated. A comparison of simulated nutrient loads with measured nutrient concentrations shows that nutrient retention in the river system may be significant. Sedimentation in the slow flowing surface water network as well as nitrogen emission to the air from the warm oxygen deficient waters are certainly partly responsible, but also wetlands along the river banks could play an important role as nutrient sinks.

The effect of drive frequency and set point amplitude on tapping forces in atomic force microscopy: simulation and experiment

International Nuclear Information System (INIS)

Legleiter, Justin

2009-01-01

In tapping mode atomic force microscopy (AFM), a sharp probe tip attached to an oscillating cantilever is allowed to intermittently strike a surface. By raster scanning the probe while monitoring the oscillation amplitude of the cantilever via a feedback loop, topographical maps of surfaces with nanoscale resolution can be acquired. While numerous studies have employed numerical simulations to elucidate the time-resolved tapping force between the probe tip and surface, until recent technique developments, specific read-outs from such models could not be experimentally verified. In this study, we explore, via numerical simulation, the impact of imaging parameters, i.e. set point ratio and drive frequency as a function of resonance, on time-varying tip-sample force interactions, which are directly compared to reconstructed tapping forces from real AFM experiments. As the AFM model contains a feedback loop allowing for the simulation of the entire scanning process, we further explore the impact that various tip-sample force have on the entire imaging process.

Post-Processing in the Material-Point Method

DEFF Research Database (Denmark)

Andersen, Søren; Andersen, Lars Vabbersgaard

The material-point method (MPM) is a numerical method for dynamic or static analysis of solids using a discretization in time and space. The method has shown to be successful in modelling physical problems involving large deformations, which are difficult to model with traditional numerical tools...... such as the finite element method. In the material-point method, a set of material points is utilized to track the problem in time and space, while a computational background grid is utilized to obtain spatial derivatives relevant to the physical problem. Currently, the research within the material-point method......-point method. The first idea involves associating a volume with each material point and displaying the deformation of this volume. In the discretization process, the physical domain is divided into a number of smaller volumes each represented by a simple shape; here quadrilaterals are chosen for the presented...

Multiscale Modeling using Molecular Dynamics and Dual Domain Material Point Method

Energy Technology Data Exchange (ETDEWEB)

Dhakal, Tilak Raj [Los Alamos National Lab. (LANL), Los Alamos, NM (United States). Theoretical Division. Fluid Dynamics and Solid Mechanics Group, T-3; Rice Univ., Houston, TX (United States)

2016-07-07

For problems involving large material deformation rate, the material deformation time scale can be shorter than the material takes to reach a thermodynamical equilibrium. For such problems, it is difficult to obtain a constitutive relation. History dependency become important because of thermodynamic non-equilibrium. Our goal is to build a multi-scale numerical method which can bypass the need for a constitutive relation. In conclusion, multi-scale simulation method is developed based on the dual domain material point (DDMP). Molecular dynamics (MD) simulation is performed to calculate stress. Since the communication among material points is not necessary, the computation can be done embarrassingly parallel in CPU-GPU platform.

ENERGY DISSIPATION IN MAGNETIC NULL POINTS AT KINETIC SCALES

International Nuclear Information System (INIS)

Olshevsky, Vyacheslav; Lapenta, Giovanni; Divin, Andrey; Eriksson, Elin; Markidis, Stefano

2015-01-01

We use kinetic particle-in-cell and MHD simulations supported by an observational data set to investigate magnetic reconnection in clusters of null points in space plasma. The magnetic configuration under investigation is driven by fast adiabatic flux rope compression that dissipates almost half of the initial magnetic field energy. In this phase powerful currents are excited producing secondary instabilities, and the system is brought into a state of “intermittent turbulence” within a few ion gyro-periods. Reconnection events are distributed all over the simulation domain and energy dissipation is rather volume-filling. Numerous spiral null points interconnected via their spines form null lines embedded into magnetic flux ropes; null point pairs demonstrate the signatures of torsional spine reconnection. However, energy dissipation mainly happens in the shear layers formed by adjacent flux ropes with oppositely directed currents. In these regions radial null pairs are spontaneously emerging and vanishing, associated with electron streams and small-scale current sheets. The number of spiral nulls in the simulation outweighs the number of radial nulls by a factor of 5–10, in accordance with Cluster observations in the Earth's magnetosheath. Twisted magnetic fields with embedded spiral null points might indicate the regions of major energy dissipation for future space missions such as the Magnetospheric Multiscale Mission

Numerical simulation of Higgs models

International Nuclear Information System (INIS)

Jaster, A.

1995-10-01

The SU(2) Higgs and the Schwinger model on the lattice were analysed. Numerical simulations of the SU(2) Higgs model were performed to study the finite temperature electroweak phase transition. With the help of the multicanonical method the distribution of an order parameter at the phase transition point was measured. This was used to obtain the order of the phase transition and the value of the interface tension with the histogram method. Numerical simulations were also performed at zero temperature to perform renormalization. The measured values for the Wilson loops were used to determine the static potential and from this the renormalized gauge coupling. The Schwinger model was simulated at different gauge couplings to analyse the properties of the Kaplan-Shamir fermions. The prediction that the mass parameter gets only multiplicative renormalization was tested and verified. (orig.)

Critical points for finite Fibonacci chains of point delta-interactions and orthogonal polynomials

International Nuclear Information System (INIS)

De Prunele, E

2011-01-01

For a one-dimensional Schroedinger operator with a finite number n of point delta-interactions with a common intensity, the parameters are the intensity, the n - 1 intercenter distances and the mass. Critical points are points in the parameters space of the Hamiltonian where one bound state appears or disappears. The study of critical points for Hamiltonians with point delta-interactions arranged along a Fibonacci chain is shown to be closely related to the study of the so-called Fibonacci operator, a discrete one-dimensional Schroedinger-type operator, which occurs in the context of tight binding Hamiltonians. These critical points are the zeros of orthogonal polynomials previously studied in the context of special diatomic linear chains with elastic nearest-neighbor interaction. Properties of the zeros (location, asymptotic behavior, gaps, ...) are investigated. The perturbation series from the solvable periodic case is determined. The measure which yields orthogonality is investigated numerically from the zeros. It is shown that the transmission coefficient at zero energy can be expressed in terms of the orthogonal polynomials and their associated polynomials. In particular, it is shown that when the number of point delta-interactions is equal to a Fibonacci number minus 1, i.e. when the intervals between point delta-interactions form a palindrome, all the Fibonacci chains at critical points are completely transparent at zero energy. (paper)

Numerical study and ex vivo assessment of HIFU treatment time reduction through optimization of focal point trajectory

Science.gov (United States)

Grisey, A.; Yon, S.; Pechoux, T.; Letort, V.; Lafitte, P.

2017-03-01

Treatment time reduction is a key issue to expand the use of high intensity focused ultrasound (HIFU) surgery, especially for benign pathologies. This study aims at quantitatively assessing the potential reduction of the treatment time arising from moving the focal point during long pulses. In this context, the optimization of the focal point trajectory is crucial to achieve a uniform thermal dose repartition and avoid boiling. At first, a numerical optimization algorithm was used to generate efficient trajectories. Thermal conduction was simulated in 3D with a finite difference code and damages to the tissue were modeled using the thermal dose formula. Given an initial trajectory, the thermal dose field was first computed, then, making use of Pontryagin's maximum principle, the trajectory was iteratively refined. Several initial trajectories were tested. Then, an ex vivo study was conducted in order to validate the efficicency of the resulting optimized strategies. Single pulses were performed at 3MHz on fresh veal liver samples with an Echopulse and the size of each unitary lesion was assessed by cutting each sample along three orthogonal planes and measuring the dimension of the whitened area based on photographs. We propose a promising approach to significantly shorten HIFU treatment time: the numerical optimization algorithm was shown to provide a reliable insight on trajectories that can improve treatment strategies. The model must now be improved in order to take in vivo conditions into account and extensively validated.

Flow area optimization in point to area or area to point flows

International Nuclear Information System (INIS)

Ghodoossi, Lotfollah; Egrican, Niluefer

2003-01-01

This paper deals with the constructal theory of generation of shape and structure in flow systems connecting one point to a finite size area. The flow direction may be either from the point to the area or the area to the point. The formulation of the problem remains the same if the flow direction is reversed. Two models are used in optimization of the point to area or area to point flow problem: cost minimization and revenue maximization. The cost minimization model enables one to predict the shape of the optimized flow areas, but the geometric sizes of the flow areas are not predictable. That is, as an example, if the area of flow is a rectangle with a fixed area size, optimization of the point to area or area to point flow problem by using the cost minimization model will only predict the height/length ratio of the rectangle not the height and length itself. By using the revenue maximization model in optimization of the flow problems, all optimized geometric aspects of the interested flow areas will be derived as well. The aim of this paper is to optimize the point to area or area to point flow problems in various elemental flow area shapes and various structures of the flow system (various combinations of elemental flow areas) by using the revenue maximization model. The elemental flow area shapes used in this paper are either rectangular or triangular. The forms of the flow area structure, made up of an assembly of optimized elemental flow areas to obtain bigger flow areas, are rectangle-in-rectangle, rectangle-in-triangle, triangle-in-triangle and triangle-in-rectangle. The global maximum revenue, revenue collected per unit flow area and the shape and sizes of each flow area structure have been derived in optimized conditions. The results for each flow area structure have been compared with the results of the other structures to determine the structure that provides better performance. The conclusion is that the rectangle-in-triangle flow area structure

Logistic simulation of the Cryomagnet Activities at Point 18

CERN Document Server

Lundgren, M

2002-01-01

The purpose of this work is to build a model in a simulation software (AutoMod) that can be used to improve the logistic planning and scheduling of the assembly and installation procedures concerning the Cryomagnets for the LHC project at point 18. The model should provide necessary information that simplifies the scheduling of activities based on time and resource availability. This report contains the results gathered from these simulations. This work been performed in the framework of a Master thesis. It is based on data available in November 2001. It has not been updated following latter decisions but the methodology and main results are still valid.

Bayesian analysis of Markov point processes

DEFF Research Database (Denmark)

Berthelsen, Kasper Klitgaard; Møller, Jesper

2006-01-01

Recently Møller, Pettitt, Berthelsen and Reeves introduced a new MCMC methodology for drawing samples from a posterior distribution when the likelihood function is only specified up to a normalising constant. We illustrate the method in the setting of Bayesian inference for Markov point processes...... a partially ordered Markov point process as the auxiliary variable. As the method requires simulation from the "unknown" likelihood, perfect simulation algorithms for spatial point processes become useful....

Influence of point defects' concentration on the ZnO matrix – a ...

African Journals Online (AJOL)

Numerical simulations of the effect of the presence of point defects in the matrix of ZnO were carried out using MATLAB computation procedure. The variation of impurity atom concentration with the energy gap was simulated. The plot obtained exhibited an exponential increase in energy gap with respect to the concentration ...

ERROR DISTRIBUTION EVALUATION OF THE THIRD VANISHING POINT BASED ON RANDOM STATISTICAL SIMULATION

Directory of Open Access Journals (Sweden)

C. Li

2012-07-01

Full Text Available POS, integrated by GPS / INS (Inertial Navigation Systems, has allowed rapid and accurate determination of position and attitude of remote sensing equipment for MMS (Mobile Mapping Systems. However, not only does INS have system error, but also it is very expensive. Therefore, in this paper error distributions of vanishing points are studied and tested in order to substitute INS for MMS in some special land-based scene, such as ground façade where usually only two vanishing points can be detected. Thus, the traditional calibration approach based on three orthogonal vanishing points is being challenged. In this article, firstly, the line clusters, which parallel to each others in object space and correspond to the vanishing points, are detected based on RANSAC (Random Sample Consensus and parallelism geometric constraint. Secondly, condition adjustment with parameters is utilized to estimate nonlinear error equations of two vanishing points (VX, VY. How to set initial weights for the adjustment solution of single image vanishing points is presented. Solving vanishing points and estimating their error distributions base on iteration method with variable weights, co-factor matrix and error ellipse theory. Thirdly, under the condition of known error ellipses of two vanishing points (VX, VY and on the basis of the triangle geometric relationship of three vanishing points, the error distribution of the third vanishing point (VZ is calculated and evaluated by random statistical simulation with ignoring camera distortion. Moreover, Monte Carlo methods utilized for random statistical estimation are presented. Finally, experimental results of vanishing points coordinate and their error distributions are shown and analyzed.

Error Distribution Evaluation of the Third Vanishing Point Based on Random Statistical Simulation

Science.gov (United States)

Li, C.

2012-07-01

POS, integrated by GPS / INS (Inertial Navigation Systems), has allowed rapid and accurate determination of position and attitude of remote sensing equipment for MMS (Mobile Mapping Systems). However, not only does INS have system error, but also it is very expensive. Therefore, in this paper error distributions of vanishing points are studied and tested in order to substitute INS for MMS in some special land-based scene, such as ground façade where usually only two vanishing points can be detected. Thus, the traditional calibration approach based on three orthogonal vanishing points is being challenged. In this article, firstly, the line clusters, which parallel to each others in object space and correspond to the vanishing points, are detected based on RANSAC (Random Sample Consensus) and parallelism geometric constraint. Secondly, condition adjustment with parameters is utilized to estimate nonlinear error equations of two vanishing points (VX, VY). How to set initial weights for the adjustment solution of single image vanishing points is presented. Solving vanishing points and estimating their error distributions base on iteration method with variable weights, co-factor matrix and error ellipse theory. Thirdly, under the condition of known error ellipses of two vanishing points (VX, VY) and on the basis of the triangle geometric relationship of three vanishing points, the error distribution of the third vanishing point (VZ) is calculated and evaluated by random statistical simulation with ignoring camera distortion. Moreover, Monte Carlo methods utilized for random statistical estimation are presented. Finally, experimental results of vanishing points coordinate and their error distributions are shown and analyzed.

Floating-to-Fixed-Point Conversion for Digital Signal Processors

Directory of Open Access Journals (Sweden)

Menard Daniel

2006-01-01

Full Text Available Digital signal processing applications are specified with floating-point data types but they are usually implemented in embedded systems with fixed-point arithmetic to minimise cost and power consumption. Thus, methodologies which establish automatically the fixed-point specification are required to reduce the application time-to-market. In this paper, a new methodology for the floating-to-fixed point conversion is proposed for software implementations. The aim of our approach is to determine the fixed-point specification which minimises the code execution time for a given accuracy constraint. Compared to previous methodologies, our approach takes into account the DSP architecture to optimise the fixed-point formats and the floating-to-fixed-point conversion process is coupled with the code generation process. The fixed-point data types and the position of the scaling operations are optimised to reduce the code execution time. To evaluate the fixed-point computation accuracy, an analytical approach is used to reduce the optimisation time compared to the existing methods based on simulation. The methodology stages are described and several experiment results are presented to underline the efficiency of this approach.

Floating-to-Fixed-Point Conversion for Digital Signal Processors

Science.gov (United States)

Menard, Daniel; Chillet, Daniel; Sentieys, Olivier

2006-12-01

Digital signal processing applications are specified with floating-point data types but they are usually implemented in embedded systems with fixed-point arithmetic to minimise cost and power consumption. Thus, methodologies which establish automatically the fixed-point specification are required to reduce the application time-to-market. In this paper, a new methodology for the floating-to-fixed point conversion is proposed for software implementations. The aim of our approach is to determine the fixed-point specification which minimises the code execution time for a given accuracy constraint. Compared to previous methodologies, our approach takes into account the DSP architecture to optimise the fixed-point formats and the floating-to-fixed-point conversion process is coupled with the code generation process. The fixed-point data types and the position of the scaling operations are optimised to reduce the code execution time. To evaluate the fixed-point computation accuracy, an analytical approach is used to reduce the optimisation time compared to the existing methods based on simulation. The methodology stages are described and several experiment results are presented to underline the efficiency of this approach.

Clinical value of homodynamic numerical simulation applied in the treatment of cerebral aneurysm.

Science.gov (United States)

Zhang, Hailin; Li, Li; Cheng, Chongjie; Sun, Xiaochuan

2017-12-01

Our objective was to evaluate the clinical value of numerical simulation in diagnosing cerebral aneurysm based on the analysis of numerical simulation of hemodynamic model. The experimental method used was the numerical model of cerebral aneurysm hemodynamic, and the numerical value of blood flow at each point was analyzed. The results showed that, the wall shear stress (WSS) value on the top of CA1 was significantly lower than that of the top (Pvalue of each point on the CA2 tumor was significantly lower than that of tumor neck (Pvalue on the tumor top and tumor neck between CA1 and CA2 had no significant difference (P>0.05); the unsteady index of shear (UIS) value at the points of 20 had distinctly changed, the wave range was 0.6-1.5; the unsteady index of pressure value of every point was significantly lower than UIS value, the wave range was 0.25-0.40. In conclusion, the application of cerebral aneurysm hemodynamic research can help doctors to diagnose cerebral aneurysm more precisely and to grasp the opportunity of treatment during the formulating of the treatment strategies.

An equilibrium point stabilization strategy for the Chen system

International Nuclear Information System (INIS)

Alvarez-Ramirez, Jose; Cevantes, Ilse; Femat, Ricardo

2004-01-01

The aim of this Letter is to address the equilibrium point stabilization problem of the Chen system by employing a simple linear feedback controller derived from time-scaling the dynamics of a single variable. The controller has the advantage of being easy to implement and a rigorous stability proof is provided based on singular perturbation arguments. Results are illustrated via numerical simulations

NUMERICAL ANALYSIS OF MATHEMATICAL MODELS OF THE FACTUAL CONTRIBUTION DISTRIBUTION IN ASYMMETRY AND DEVIATION OF VOLTAGE AT THE COMMON COUPLING POINTS OF ENERGY SUPPLY SYSTEMS

Directory of Open Access Journals (Sweden)

Yu.L. Sayenko

2016-05-01

Full Text Available Purpose. Perform numerical analysis of the distribution of the factual contributions of line sources of distortion in the voltage distortion at the point of common coupling, based on the principles of superposition and exclusions. Methodology. Numerical analysis was performed on the results of the simulation steady state operation of power supply system of seven electricity consumers. Results. Mathematical model for determining the factual contribution of line sources of distortion in the voltage distortion at the point of common coupling, based on the principles of superposition and exclusions, are equivalent. To assess the degree of participation of each source of distortion in the voltage distortion at the point of common coupling and distribution of financial compensation to the injured party by all sources of distortion developed a one-dimensional criteria based on the scalar product of vectors. Not accounting group sources of distortion, which belong to the subject of the energy market, to determine their total factual contribution as the residual of the factual contribution between all sources of distortion. Originality. Simulation mode power supply system was carried out in the phase components space, taking into account the distributed characteristics of distortion sources. Practical value. The results of research can be used to develop methods and tools for distributed measurement and analytical systems assessment of the power quality.

Substructuring in the implicit simulation of single point incremental sheet forming

NARCIS (Netherlands)

Hadoush, A.; van den Boogaard, Antonius H.

2009-01-01

This paper presents a direct substructuring method to reduce the computing time of implicit simulations of single point incremental forming (SPIF). Substructuring is used to divide the finite element (FE) mesh into several non-overlapping parts. Based on the hypothesis that plastic deformation is

Polynomial approach method to solve the neutron point kinetics equations with use of the analytic continuation

Energy Technology Data Exchange (ETDEWEB)

Tumelero, Fernanda; Petersen, Claudio Zen; Goncalves, Glenio Aguiar [Universidade Federal de Pelotas, Capao do Leao, RS (Brazil). Programa de Pos Graduacao em Modelagem Matematica; Schramm, Marcelo [Universidade Federal do Rio Grande do Sul, Porto Alegre, RS (Brazil). Programa de Pos-Graduacao em Engenharia Mecanica

2016-12-15

In this work, we report a solution to solve the Neutron Point Kinetics Equations applying the Polynomial Approach Method. The main idea is to expand the neutron density and delayed neutron precursors as a power series considering the reactivity as an arbitrary function of the time in a relatively short time interval around an ordinary point. In the first interval one applies the initial conditions and the analytical continuation is used to determine the solutions of the next intervals. A genuine error control is developed based on an analogy with the Rest Theorem. For illustration, we also report simulations for different approaches types (linear, quadratic and cubic). The results obtained by numerical simulations for linear approximation are compared with results in the literature.

Field-scale operation of methane biofiltration systems to mitigate point source methane emissions

International Nuclear Information System (INIS)

Hettiarachchi, Vijayamala C.; Hettiaratchi, Patrick J.; Mehrotra, Anil K.; Kumar, Sunil

2011-01-01

Methane biofiltration (MBF) is a novel low-cost technique for reducing low volume point source emissions of methane (CH 4 ). MBF uses a granular medium, such as soil or compost, to support the growth of methanotrophic bacteria responsible for converting CH 4 to carbon dioxide (CO 2 ) and water (H 2 O). A field research program was undertaken to evaluate the potential to treat low volume point source engineered CH 4 emissions using an MBF at a natural gas monitoring station. A new comprehensive three-dimensional numerical model was developed incorporating advection-diffusive flow of gas, biological reactions and heat and moisture flow. The one-dimensional version of this model was used as a guiding tool for designing and operating the MBF. The long-term monitoring results of the field MBF are also presented. The field MBF operated with no control of precipitation, evaporation, and temperature, provided more than 80% of CH 4 oxidation throughout spring, summer, and fall seasons. The numerical model was able to predict the CH 4 oxidation behavior of the field MBF with high accuracy. The numerical model simulations are presented for estimating CH 4 oxidation efficiencies under various operating conditions, including different filter bed depths and CH 4 flux rates. The field observations as well as numerical model simulations indicated that the long-term performance of MBFs is strongly dependent on environmental factors, such as ambient temperature and precipitation. - Highlights: → One-dimensional version of the model was used as a guiding tool for designing and operating the MBF. → Mathematical model predicted CH 4 oxidation behaviors of the field MBF with high accuracy i.e. (> 80 %). → Performance of MBF is dependent on ambient temperature and precipitation. - The developed numerical model simulations and field observations for estimating CH 4 oxidation efficiencies under various operating conditions indicate that the long-term performance of MBFs is strongly

Point Cloud Based Relative Pose Estimation of a Satellite in Close Range

Directory of Open Access Journals (Sweden)

Lujiang Liu

2016-06-01

Full Text Available Determination of the relative pose of satellites is essential in space rendezvous operations and on-orbit servicing missions. The key problems are the adoption of suitable sensor on board of a chaser and efficient techniques for pose estimation. This paper aims to estimate the pose of a target satellite in close range on the basis of its known model by using point cloud data generated by a flash LIDAR sensor. A novel model based pose estimation method is proposed; it includes a fast and reliable pose initial acquisition method based on global optimal searching by processing the dense point cloud data directly, and a pose tracking method based on Iterative Closest Point algorithm. Also, a simulation system is presented in this paper in order to evaluate the performance of the sensor and generate simulated sensor point cloud data. It also provides truth pose of the test target so that the pose estimation error can be quantified. To investigate the effectiveness of the proposed approach and achievable pose accuracy, numerical simulation experiments are performed; results demonstrate algorithm capability of operating with point cloud directly and large pose variations. Also, a field testing experiment is conducted and results show that the proposed method is effective.

Parametric study of two-body floating-point wave absorber

Science.gov (United States)

Amiri, Atena; Panahi, Roozbeh; Radfar, Soheil

2016-03-01

In this paper, we present a comprehensive numerical simulation of a point wave absorber in deep water. Analyses are performed in both the frequency and time domains. The converter is a two-body floating-point absorber (FPA) with one degree of freedom in the heave direction. Its two parts are connected by a linear mass-spring-damper system. The commercial ANSYS-AQWA software used in this study performs well in considering validations. The velocity potential is obtained by assuming incompressible and irrotational flow. As such, we investigated the effects of wave characteristics on energy conversion and device efficiency, including wave height and wave period, as well as the device diameter, draft, geometry, and damping coefficient. To validate the model, we compared our numerical results with those from similar experiments. Our study results can clearly help to maximize the converter's efficiency when considering specific conditions.

Word Length Selection Method for Controller Implementation on FPGAs Using the VHDL-2008 Fixed-Point and Floating-Point Packages

Directory of Open Access Journals (Sweden)

Urriza I

2010-01-01

Full Text Available Abstract This paper presents a word length selection method for the implementation of digital controllers in both fixed-point and floating-point hardware on FPGAs. This method uses the new types defined in the VHDL-2008 fixed-point and floating-point packages. These packages allow customizing the word length of fixed and floating point representations and shorten the design cycle simplifying the design of arithmetic operations. The method performs bit-true simulations in order to determine the word length to represent the constant coefficients and the internal signals of the digital controller while maintaining the control system specifications. A mixed-signal simulation tool is used to simulate the closed loop system as a whole in order to analyze the impact of the quantization effects and loop delays on the control system performance. The method is applied to implement a digital controller for a switching power converter. The digital circuit is implemented on an FPGA, and the simulations are experimentally verified.

Numerical simulation of large deformation polycrystalline plasticity

International Nuclear Information System (INIS)

Inal, K.; Neale, K.W.; Wu, P.D.; MacEwen, S.R.

2000-01-01

A finite element model based on crystal plasticity has been developed to simulate the stress-strain response of sheet metal specimens in uniaxial tension. Each material point in the sheet is considered to be a polycrystalline aggregate of FCC grains. The Taylor theory of crystal plasticity is assumed. The numerical analysis incorporates parallel computing features enabling simulations of realistic models with large number of grains. Simulations have been carried out for the AA3004-H19 aluminium alloy and the results are compared with experimental data. (author)

Use of Green's functions in the numerical solution of two-point boundary value problems

Science.gov (United States)

Gallaher, L. J.; Perlin, I. E.

1974-01-01

This study investigates the use of Green's functions in the numerical solution of the two-point boundary value problem. The first part deals with the role of the Green's function in solving both linear and nonlinear second order ordinary differential equations with boundary conditions and systems of such equations. The second part describes procedures for numerical construction of Green's functions and considers briefly the conditions for their existence. Finally, there is a description of some numerical experiments using nonlinear problems for which the known existence, uniqueness or convergence theorems do not apply. Examples here include some problems in finding rendezvous orbits of the restricted three body system.

Simulation of ultrasonic surface waves with multi-Gaussian and point source beam models

International Nuclear Information System (INIS)

Zhao, Xinyu; Schmerr, Lester W. Jr.; Li, Xiongbing; Sedov, Alexander

2014-01-01

In the past decade, multi-Gaussian beam models have been developed to solve many complicated bulk wave propagation problems. However, to date those models have not been extended to simulate the generation of Rayleigh waves. Here we will combine Gaussian beams with an explicit high frequency expression for the Rayleigh wave Green function to produce a three-dimensional multi-Gaussian beam model for the fields radiated from an angle beam transducer mounted on a solid wedge. Simulation results obtained with this model are compared to those of a point source model. It is shown that the multi-Gaussian surface wave beam model agrees well with the point source model while being computationally much more efficient

Complex saddle points and the sign problem in complex Langevin simulation

International Nuclear Information System (INIS)

Hayata, Tomoya; Hidaka, Yoshimasa; Tanizaki, Yuya

2016-01-01

We show that complex Langevin simulation converges to a wrong result within the semiclassical analysis, by relating it to the Lefschetz-thimble path integral, when the path-integral weight has different phases among dominant complex saddle points. Equilibrium solution of the complex Langevin equation forms local distributions around complex saddle points. Its ensemble average approximately becomes a direct sum of the average in each local distribution, where relative phases among them are dropped. We propose that by taking these phases into account through reweighting, we can solve the wrong convergence problem. However, this prescription may lead to a recurrence of the sign problem in the complex Langevin method for quantum many-body systems.

Point defects and irradiation in oxides: simulations at the atomic scale

International Nuclear Information System (INIS)

Crocombette, J.P.

2005-12-01

The studies done by Jean-Paul Crocombette between 1996 and 2005 in the Service de Recherches de Metallurgie Physique of the Direction de l'Energie Nucleaire in Saclay are presented in this Habilitation thesis. These works were part of the material science researches on the ageing, especially under irradiation, of oxides of interest for the nuclear industry. In this context simulation studies at the atomic scale were performed on two elementary components of ageing under irradiation : point defects and displacement cascades ; using two complementary simulation techniques : ab initio electronic structure calculations and empirical potential molecular dynamics. The first part deals with point defects : self defects (vacancies or interstitials) or hetero-atomic dopants. One first recalls the energetics of such defects in oxides, the specific features of defects calculations and the expected accuracy of these calculations. Then one presents the results obtained on uranium dioxide, oxygen in silver and amorphous silica. The second part tackles the modelling of disintegration recoil nuclei in various?displacement cascades created by crystalline matrices for actinide waste disposal. Cascade calculations give access to the amorphization mechanisms under irradiation of these materials. One thus predicts that the amorphization in zircon takes place directly in the tracks whereas in lanthanum zirconate, the amorphization proceeds through the accumulation of point defects. Finally the prospects of these studies are discussed. (author)

Statistical representation of a spray as a point process

International Nuclear Information System (INIS)

Subramaniam, S.

2000-01-01

The statistical representation of a spray as a finite point process is investigated. One objective is to develop a better understanding of how single-point statistical information contained in descriptions such as the droplet distribution function (ddf), relates to the probability density functions (pdfs) associated with the droplets themselves. Single-point statistical information contained in the droplet distribution function (ddf) is shown to be related to a sequence of single surrogate-droplet pdfs, which are in general different from the physical single-droplet pdfs. It is shown that the ddf contains less information than the fundamental single-point statistical representation of the spray, which is also described. The analysis shows which events associated with the ensemble of spray droplets can be characterized by the ddf, and which cannot. The implications of these findings for the ddf approach to spray modeling are discussed. The results of this study also have important consequences for the initialization and evolution of direct numerical simulations (DNS) of multiphase flows, which are usually initialized on the basis of single-point statistics such as the droplet number density in physical space. If multiphase DNS are initialized in this way, this implies that even the initial representation contains certain implicit assumptions concerning the complete ensemble of realizations, which are invalid for general multiphase flows. Also the evolution of a DNS initialized in this manner is shown to be valid only if an as yet unproven commutation hypothesis holds true. Therefore, it is questionable to what extent DNS that are initialized in this manner constitute a direct simulation of the physical droplets. Implications of these findings for large eddy simulations of multiphase flows are also discussed. (c) 2000 American Institute of Physics

Pseudo-dynamic source modelling with 1-point and 2-point statistics of earthquake source parameters

KAUST Repository

Song, S. G.

2013-12-24

Ground motion prediction is an essential element in seismic hazard and risk analysis. Empirical ground motion prediction approaches have been widely used in the community, but efficient simulation-based ground motion prediction methods are needed to complement empirical approaches, especially in the regions with limited data constraints. Recently, dynamic rupture modelling has been successfully adopted in physics-based source and ground motion modelling, but it is still computationally demanding and many input parameters are not well constrained by observational data. Pseudo-dynamic source modelling keeps the form of kinematic modelling with its computational efficiency, but also tries to emulate the physics of source process. In this paper, we develop a statistical framework that governs the finite-fault rupture process with 1-point and 2-point statistics of source parameters in order to quantify the variability of finite source models for future scenario events. We test this method by extracting 1-point and 2-point statistics from dynamically derived source models and simulating a number of rupture scenarios, given target 1-point and 2-point statistics. We propose a new rupture model generator for stochastic source modelling with the covariance matrix constructed from target 2-point statistics, that is, auto- and cross-correlations. Our sensitivity analysis of near-source ground motions to 1-point and 2-point statistics of source parameters provides insights into relations between statistical rupture properties and ground motions. We observe that larger standard deviation and stronger correlation produce stronger peak ground motions in general. The proposed new source modelling approach will contribute to understanding the effect of earthquake source on near-source ground motion characteristics in a more quantitative and systematic way.

Calculation Method for Equilibrium Points in Dynamical Systems Based on Adaptive Sinchronization

Directory of Open Access Journals (Sweden)

Manuel Prian Rodríguez

2017-12-01

Full Text Available In this work, a control system is proposed as an equivalent numerical procedure whose aim is to obtain the natural equilibrium points of a dynamical system. These equilibrium points may be employed later as setpoint signal for different control techniques. The proposed procedure is based on the adaptive synchronization between an oscillator and a reference model driven by the oscillator state variables. A stability analysis is carried out and a simplified algorithm is proposed. Finally, satisfactory simulation results are shown.

The flux-coordinate independent approach applied to X-point geometries

International Nuclear Information System (INIS)

Hariri, F.; Hill, P.; Ottaviani, M.; Sarazin, Y.

2014-01-01

A Flux-Coordinate Independent (FCI) approach for anisotropic systems, not based on magnetic flux coordinates, has been introduced in Hariri and Ottaviani [Comput. Phys. Commun. 184, 2419 (2013)]. In this paper, we show that the approach can tackle magnetic configurations including X-points. Using the code FENICIA, an equilibrium with a magnetic island has been used to show the robustness of the FCI approach to cases in which a magnetic separatrix is present in the system, either by design or as a consequence of instabilities. Numerical results are in good agreement with the analytic solutions of the sound-wave propagation problem. Conservation properties are verified. Finally, the critical gain of the FCI approach in situations including the magnetic separatrix with an X-point is demonstrated by a fast convergence of the code with the numerical resolution in the direction of symmetry. The results highlighted in this paper show that the FCI approach can efficiently deal with X-point geometries

Identification of discrete chaotic maps with singular points

Directory of Open Access Journals (Sweden)

P. G. Akishin

2001-01-01

Full Text Available We investigate the ability of artificial neural networks to reconstruct discrete chaotic maps with singular points. We use as a simple test model the Cusp map. We compare the traditional Multilayer Perceptron, the Chebyshev Neural Network and the Wavelet Neural Network. The numerical scheme for the accurate determination of a singular point is also developed. We show that combining a neural network with the numerical algorithm for the determination of the singular point we are able to accurately approximate discrete chaotic maps with singularities.

Radar Image Simulation: Validation of the Point Scattering Method. Volume 2

Science.gov (United States)

1977-09-01

the Engineer Topographic Labor - atory (ETL), Fort Belvoir, Virginia. This Radar Simulation Study was performed to validate the point tcattering radar...e.n For radar, the number of Independent samples in a given re.-olution cell is given by 5 ,: N L 2w (16) L Acoso where: 0 Radar incidence angle; w

POD for Real-Time Simulation of Hyperelastic Soft Biological Tissue Using the Point Collocation Method of Finite Spheres

Directory of Open Access Journals (Sweden)

Suleiman Banihani

2013-01-01

Full Text Available The point collocation method of finite spheres (PCMFS is used to model the hyperelastic response of soft biological tissue in real time within the framework of virtual surgery simulation. The proper orthogonal decomposition (POD model order reduction (MOR technique was used to achieve reduced-order model of the problem, minimizing computational cost. The PCMFS is a physics-based meshfree numerical technique for real-time simulation of surgical procedures where the approximation functions are applied directly on the strong form of the boundary value problem without the need for integration, increasing computational efficiency. Since computational speed has a significant role in simulation of surgical procedures, the proposed technique was able to model realistic nonlinear behavior of organs in real time. Numerical results are shown to demonstrate the effectiveness of the new methodology through a comparison between full and reduced analyses for several nonlinear problems. It is shown that the proposed technique was able to achieve good agreement with the full model; moreover, the computational and data storage costs were significantly reduced.

Point-of-care echocardiography in simulation-based education and assessment.

Science.gov (United States)

Amini, Richard; Stolz, Lori A; Javedani, Parisa P; Gaskin, Kevin; Baker, Nicola; Ng, Vivienne; Adhikari, Srikar

2016-01-01

Emergency medicine milestones released by the Accreditation Council for Graduate Medical Education require residents to demonstrate competency in bedside ultrasound (US). The acquisition of these skills necessitates a combination of exposure to clinical pathology, hands-on US training, and feedback. We describe a novel simulation-based educational and assessment tool designed to evaluate emergency medicine residents' competency in point-of-care echocardiography for evaluation of a hypotensive patient with chest pain using bedside US. This was a cross-sectional study conducted at an academic medical center. A simulation-based module was developed to teach and assess the use of point-of-care echocardiography in the evaluation of the hypotensive patient. The focus of this module was sonographic imaging of cardiac pathology, and this focus was incorporated in all components of the session: asynchronous learning, didactic lecture, case-based learning, and hands-on stations. A total of 52 residents with varying US experience participated in this study. Questions focused on knowledge assessment demonstrated improvement across the postgraduate year (PGY) of training. Objective standardized clinical examination evaluation demonstrated improvement between PGY I and PGY III; however, it was noted that there was a small dip in hands-on scanning skills during the PGY II. Clinical diagnosis and management skills also demonstrated incremental improvement across the PGY of training. The 1-day, simulation-based US workshop was an effective educational and assessment tool at our institution.

Theory of Single Point Incremental Forming

DEFF Research Database (Denmark)

Martins, P.A.F.; Bay, Niels; Skjødt, Martin

2008-01-01

This paper presents a closed-form theoretical analysis modelling the fundamentals of single point incremental forming and explaining the experimental and numerical results available in the literature for the past couple of years. The model is based on membrane analysis with bi-directional in-plan......-plane contact friction and is focused on the extreme modes of deformation that are likely to be found in single point incremental forming processes. The overall investigation is supported by experimental work performed by the authors and data retrieved from the literature.......This paper presents a closed-form theoretical analysis modelling the fundamentals of single point incremental forming and explaining the experimental and numerical results available in the literature for the past couple of years. The model is based on membrane analysis with bi-directional in...

Big Rock Point: 35 years of electrical generation

International Nuclear Information System (INIS)

Petrosky, T.D.

1998-01-01

On September 27, 1962, the 75 MWe boiling water reactor, designed and built by General Electric, of the Big Rock Point Nuclear Power Station went critical for the first time. The US Atomic Energy Commission (AEC) and the plant operator, Consumers Power, had designed the plant also as a research reactor. The first studies were devoted to fuel behavior, higher burnup, and materials research. The reactor was also used for medical technology: Co-60 radiation sources were produced for the treatment of more than 120,000 cancer patients. After the accident at the Three Mile Island-2 nuclear generating unit in 1979, Big Rock Point went through an extensive backfitting phase. Personnel from numerous other American nuclear power plants were trained at the simulator of Big Rock Point. The plant was decommissioned permanently on August 29, 1997 after more than 35 years of operation and a cumulated electric power production of 13,291 GWh. A period of five to seven years is estimated for decommissioning and demolition work up to the 'green field' stage. (orig.) [de

Octopuses use a human-like strategy to control precise point-to-point arm movements.

Science.gov (United States)

Sumbre, Germán; Fiorito, Graziano; Flash, Tamar; Hochner, Binyamin

2006-04-18

One of the key problems in motor control is mastering or reducing the number of degrees of freedom (DOFs) through coordination. This problem is especially prominent with hyper-redundant limbs such as the extremely flexible arm of the octopus. Several strategies for simplifying these control problems have been suggested for human point-to-point arm movements. Despite the evolutionary gap and morphological differences, humans and octopuses evolved similar strategies when fetching food to the mouth. To achieve this precise point-to-point-task, octopus arms generate a quasi-articulated structure based on three dynamic joints. A rotational movement around these joints brings the object to the mouth . Here, we describe a peripheral neural mechanism-two waves of muscle activation propagate toward each other, and their collision point sets the medial-joint location. This is a remarkably simple mechanism for adjusting the length of the segments according to where the object is grasped. Furthermore, similar to certain human arm movements, kinematic invariants were observed at the joint level rather than at the end-effector level, suggesting intrinsic control coordination. The evolutionary convergence to similar geometrical and kinematic features suggests that a kinematically constrained articulated limb controlled at the level of joint space is the optimal solution for precise point-to-point movements.

Application of point kinetic model in the study of fluidized bed reactor dynamic

International Nuclear Information System (INIS)

Borges, Volnei; Vilhena, Marco Tullio de; Streck, Elaine E.

1995-01-01

In this work the dynamical behavior of the fluidized bed nuclear reactor is analysed. The main goal consist to study the effect of the acceleration term in the point kinetic equations. Numerical simulations are reported considering constant acceleration. (author). 7 refs, 4 figs

Solving the Richardson equations close to the critical points

Energy Technology Data Exchange (ETDEWEB)

DomInguez, F [Departamento de Matematicas, Universidad de Alcala, 28871 Alcala de Henares (Spain); Esebbag, C [Departamento de Matematicas, Universidad de Alcala, 28871 Alcala de Henares (Spain); Dukelsky, J [Instituto de Estructura de la Materia, CSIC, Serrano 123, 28006 Madrid (Spain)

2006-09-15

We study the Richardson equations close to the critical values of the pairing strength g{sub c}, where the occurrence of divergences precludes numerical solutions. We derive a set of equations for determining the critical g values and the non-collapsing pair energies. Studying the behaviour of the solutions close to the critical points, we develop a procedure to solve numerically the Richardson equations for arbitrary coupling strength.

Evolution and experience with the ATLAS Simulation at Point1 Project

CERN Document Server

AUTHOR|(INSPIRE)INSPIRE-00389536; The ATLAS collaboration; Brasolin, Franco; Kouba, Tomas; Schovancova, Jaroslava; Fazio, Daniel; Di Girolamo, Alessandro; Scannicchio, Diana; Twomey, Matthew Shaun; Wang, Fuquan; Zaytsev, Alexander; Lee, Christopher

2017-01-01

The Simulation at Point1 project is successfully running standard ATLAS simulation jobs on the TDAQ HLT resources. The pool of available resources changes dynamically, therefore we need to be very effective in exploiting the available computing cycles. We present our experience with using the Event Service that provides the event-level granularity of computations. We show the design decisions and overhead time related to the usage of the Event Service. The improved utilization of the resources is also presented with the recent development in monitoring, automatic alerting, deployment and GUI.

Evolution and experience with the ATLAS simulation at Point1 project

CERN Document Server

Ballestrero, Sergio; The ATLAS collaboration; Fazio, Daniel; Di Girolamo, Alessandro; Kouba, Tomas; Lee, Christopher; Scannicchio, Diana; Schovancova, Jaroslava; Twomey, Matthew Shaun; Wang, Fuquan; Zaytsev, Alexander

2016-01-01

The Simulation at Point1 project is successfully running traditional ATLAS simulation jobs on the TDAQ HLT resources. The pool of available resources changes dynamically, therefore we need to be very effective in exploiting the available computing cycles. We will present our experience with using the Event Service that provides the event-level granularity of computations. We will show the design decisions and overhead time related to the usage of the Event Service. The improved utilization of the resources will also be presented with the recent development in monitoring, automatic alerting, deployment and GUI.

A Fixed-point Scheme for the Numerical Construction of Magnetohydrostatic Atmospheres in Three Dimensions

Science.gov (United States)

Gilchrist, S. A.; Braun, D. C.; Barnes, G.

2016-12-01

Magnetohydrostatic models of the solar atmosphere are often based on idealized analytic solutions because the underlying equations are too difficult to solve in full generality. Numerical approaches, too, are often limited in scope and have tended to focus on the two-dimensional problem. In this article we develop a numerical method for solving the nonlinear magnetohydrostatic equations in three dimensions. Our method is a fixed-point iteration scheme that extends the method of Grad and Rubin ( Proc. 2nd Int. Conf. on Peaceful Uses of Atomic Energy 31, 190, 1958) to include a finite gravity force. We apply the method to a test case to demonstrate the method in general and our implementation in code in particular.

Determination of the boiling-point distribution by simulated distillation from n-pentane through n-tetratetracontane in 70 to 80 seconds.

Science.gov (United States)

Lubkowitz, Joaquin A; Meneghini, Roberto I

2002-01-01

This work presents the carrying out of boiling-point distributions by simulated distillation with direct-column heating rather than oven-column heating. Column-heating rates of 300 degrees C/min are obtained yielding retention times of 73 s for n-tetratetracontane. The calibration curves of the retention time versus the boiling point, in the range of n-pentane to n-tetratetracontane, are identical to those obtained by slower oven-heating rates. The boiling-point distribution of the reference gas oil is compared with that obtained with column oven heating at rates of 15 to 40 degrees C/min. The results show boiling-point distribution values nearly the same (1-2 degrees F) as those obtained with oven column heating from the initial boiling point to 80% distilled off. Slightly higher differences are obtained (3-4 degrees F) for the 80% distillation to final boiling-point interval. Nonetheless, allowed consensus differences are never exceeded. Precision of the boiling-point distributions (expressed as standard deviations) are 0.1-0.3% for the data obtained in the direct column-heating mode.

Point-of-care ultrasound education: the increasing role of simulation and multimedia resources.

Science.gov (United States)

Lewiss, Resa E; Hoffmann, Beatrice; Beaulieu, Yanick; Phelan, Mary Beth

2014-01-01

This article reviews the current technology, literature, teaching models, and methods associated with simulation-based point-of-care ultrasound training. Patient simulation appears particularly well suited for learning point-of-care ultrasound, which is a required core competency for emergency medicine and other specialties. Work hour limitations have reduced the opportunities for clinical practice, and simulation enables practicing a skill multiple times before it may be used on patients. Ultrasound simulators can be categorized into 2 groups: low and high fidelity. Low-fidelity simulators are usually static simulators, meaning that they have nonchanging anatomic examples for sonographic practice. Advantages are that the model may be reused over time, and some simulators can be homemade. High-fidelity simulators are usually high-tech and frequently consist of many computer-generated cases of virtual sonographic anatomy that can be scanned with a mock probe. This type of equipment is produced commercially and is more expensive. High-fidelity simulators provide students with an active and safe learning environment and make a reproducible standardized assessment of many different ultrasound cases possible. The advantages and disadvantages of using low- versus high-fidelity simulators are reviewed. An additional concept used in simulation-based ultrasound training is blended learning. Blended learning may include face-to-face or online learning often in combination with a learning management system. Increasingly, with simulation and Web-based learning technologies, tools are now available to medical educators for the standardization of both ultrasound skills training and competency assessment.

Numerical methods used in simulation

International Nuclear Information System (INIS)

Caseau, Paul; Perrin, Michel; Planchard, Jacques

1978-01-01

The fundamental numerical problem posed by simulation problems is the stability of the resolution diagram. The system of the most used equations is defined, since there is a family of models of increasing complexity with 3, 4 or 5 equations although only models with 3 and 4 equations have been used extensively. After defining what is meant by explicit or implicit, the best established stability results is given for one-dimension problems and then for two-dimension problems. It is shown that two types of discretisation may be defined: four and eight point diagrams (in one or two dimensions) and six and ten point diagrams (in one or two dimensions). To end, some results are given on problems that are not usually treated very much, i.e. non-asymptotic stability and the stability of diagrams based on finite elements [fr

Practical considerations in developing numerical simulators for thermal recovery

Energy Technology Data Exchange (ETDEWEB)

Abou-Kassem, J.H. [Chemical and Petroleum Engineering Department, UAE University, Al-Ain (United Arab Emirates)

1996-08-15

Numerical simulation of steam injection and in-situ combustion-based oil recovery processes is of great importance in project design. Development of such numerical simulators is an on-going process, with improvements made as the process description becomes more complete, and also as better methods are devised to resolve certain numerical difficulties. This paper addresses some of the latter, and based on the author`s experience gives useful guidelines for developing more efficient numerical simulators of steam injection and in-situ combustion. The paper takes up a series of questions related to simulating thermal processes. Included are: the elimination of constraint equations at the matrix level, phase change, steam injection rate, alternative treatments of heat loss, relative permeabilities and importance of hysteresis effects, improved solutions to the grid orientation problem and other simulation problems such as potential inversion, grid block size, time-step size control and induced fractures. The points discussed in the paper should be of use to both simulator developers and users alike, and will lead to a better understanding of simulation results

From Contrapuntal Music to Polyphonic Novel: Aldous Huxley’s Point Counter Point

Directory of Open Access Journals (Sweden)

Mevlüde ZENGİN

2015-06-01

Full Text Available Taken at face value Point Counter Point (1928 written by Aldous Huxley seems to be a novel including many stories of various and sundry people and reflecting their points of view about the world in which they live and about the life they have been leading. However, it is this very quality of the novel that provides grounds for the study of the novel as a polyphonic one. The novel presents to its reader an aggregate of strikingly different characters and thus a broad spectrum of contemporary society. The characters in the novel are all characterized by and individualized with easily recognizable physical, intellectual, emotional, psychological and moral qualities. Each of them is well-contrived through their differences in social status, political views, wealth, etc. Thus, many different viewpoints, conflicting voices, contrasting insights and ideas are heard and seen synchronically in Point Counter Point, which makes it polyphonic. Polyphony is a musical motif referring to different notes and chords played at the same time to create a rhythm. It was first adopted by M. M. Bakhtin to analyze F. M. Dostoyevsky’s fiction. The aim of this study is firstly to elucidate, in Bakhtinian thought, polyphony and then dialogism and heteroglossia closely related to his concept of polyphony; and then to put the polyphonic qualities in Point Counter Point forth, studying the novel’s dialogism and heteroglot qualities

A new dipolar potential for numerical simulations of polar fluids on the 4D hypersphere

International Nuclear Information System (INIS)

Caillol, Jean-Michel; Trulsson, Martin

2014-01-01

We present a new method for Monte Carlo or Molecular Dynamics numerical simulations of three-dimensional polar fluids. The simulation cell is defined to be the surface of the northern hemisphere of a four-dimensional (hyper)sphere. The point dipoles are constrained to remain tangent to the sphere and their interactions are derived from the basic laws of electrostatics in this geometry. The dipole-dipole potential has two singularities which correspond to the following boundary conditions: when a dipole leaves the northern hemisphere at some point of the equator, it reappears at the antipodal point bearing the same dipole moment. We derive all the formal expressions needed to obtain the thermodynamic and structural properties of a polar liquid at thermal equilibrium in actual numerical simulation. We notably establish the expression of the static dielectric constant of the fluid as well as the behavior of the pair correlation at large distances. We report and discuss the results of extensive numerical Monte Carlo simulations for two reference states of a fluid of dipolar hard spheres and compare these results with previous methods with a special emphasis on finite size effects

Digital microwave communication engineering point-to-point microwave systems

CERN Document Server

Kizer, George

2013-01-01

The first book to cover all engineering aspects of microwave communication path design for the digital age Fixed point-to-point microwave systems provide moderate-capacity digital transmission between well-defined locations. Most popular in situations where fiber optics or satellite communication is impractical, it is commonly used for cellular or PCS site interconnectivity where digital connectivity is needed but not economically available from other sources, and in private networks where reliability is most important. Until now, no book has adequately treated all en

Comparison of Dose When Prescribed to Point A and Point H for Brachytherapy in Cervical Cancer

Energy Technology Data Exchange (ETDEWEB)

Gang, Ji Hyeong; Gim, Il Hwan; Hwang, Seon Boong; Kim, Woong; Im, Hyeong Seo; Gang, Jin Mook; Gim, Gi Hwan; Lee, Ah Ram [Dept. of Radiation Oncology, Korea Institute of Radiological and Medical Sciences, Seou (Korea, Republic of)

2012-09-15

The purpose of this study is to compare plans prescribed to point A with these prescribed to point H recommended by ABS (American Brachytherapy Society) in high dose rate intracavitary brachytherapy for cervical carcinoma. This study selected 103 patients who received HDR (High Dose Rate) brachytherapy using tandem and ovoids from March 2010 to January 2012. Point A, bladder point, and rectal point conform with Manchester System. Point H conforms with ABS recommendation. Also Sigmoid colon point, and vagina point were established arbitrarily. We examined distance between point A and point H. The percent dose at point A was calculated when 100% dose was prescribed to point H. Additionally, the percent dose at each reference points when dose is prescribed to point H and point A were calculated. The relative dose at point A was lower when point H was located inferior to point A. The relative doses at bladder, rectal, sigmoid colon, and vagina points were higher when point H was located superior to point A, and lower when point H was located inferior to point A. This study found out that as point H got located much superior to point A, the absorbed dose of surrounding normal organs became higher, and as point H got located much inferior to point A, the absorbed dose of surrounding normal organs became lower. This differences dose not seem to affect the treatment. However, we suggest this new point is worth being considered for the treatment of HDR if dose distribution and absorbed dose at normal organs have large differences between prescribed to point A and H.

Present status on numerical algorithms and benchmark tests for point kinetics and quasi-static approximate kinetics

International Nuclear Information System (INIS)

Ise, Takeharu

1976-12-01

Review studies have been made on algorithms of numerical analysis and benchmark tests on point kinetics and quasistatic approximate kinetics computer codes to perform efficiently benchmark tests on space-dependent neutron kinetics codes. Point kinetics methods have now been improved since they can be directly applied to the factorization procedures. Methods based on Pade rational function give numerically stable solutions and methods on matrix-splitting are interested in the fact that they are applicable to the direct integration methods. An improved quasistatic (IQ) approximation is the best and the most practical method; it is numerically shown that the IQ method has a high stability and precision and the computation time which is about one tenth of that of the direct method. IQ method is applicable to thermal reactors as well as fast reactors and especially fitted for fast reactors to which many time steps are necessary. Two-dimensional diffusion kinetics codes are most practicable though there exist also three-dimensional diffusion kinetics code as well as two-dimensional transport kinetics code. On developing a space-dependent kinetics code, in any case, it is desirable to improve the method so as to have a high computing speed for solving static diffusion and transport equations. (auth.)

Variogram based and Multiple - Point Statistical simulation of shallow aquifer structures in the Upper Salzach valley, Austria

Science.gov (United States)

Jandrisevits, Carmen; Marschallinger, Robert

2014-05-01

Quarternary sediments in overdeepened alpine valleys and basins in the Eastern Alps bear substantial groundwater resources. The associated aquifer systems are generally geometrically complex with highly variable hydraulic properties. 3D geological models provide predictions of both geometry and properties of the subsurface required for subsequent modelling of groundwater flow and transport. In hydrology, geostatistical Kriging and Kriging based conditional simulations are widely used to predict the spatial distribution of hydrofacies. In the course of investigating the shallow aquifer structures in the Zell basin in the Upper Salzach valley (Salzburg, Austria), a benchmark of available geostatistical modelling and simulation methods was performed: traditional variogram based geostatistical methods, i.e. Indicator Kriging, Sequential Indicator Simulation and Sequential Indicator Co - Simulation were used as well as Multiple Point Statistics. The ~ 6 km2 investigation area is sampled by 56 drillings with depths of 5 to 50 m; in addition, there are 2 geophysical sections with lengths of 2 km and depths of 50 m. Due to clustered drilling sites, indicator Kriging models failed to consistently model the spatial variability of hydrofacies. Using classical variogram based geostatistical simulation (SIS), equally probable realizations were generated with differences among the realizations providing an uncertainty measure. The yielded models are unstructured from a geological point - they do not portray the shapes and lateral extensions of associated sedimentary units. Since variograms consider only two - point spatial correlations, they are unable to capture the spatial variability of complex geological structures. The Multiple Point Statistics approach overcomes these limitations of two point statistics as it uses a Training image instead of variograms. The 3D Training Image can be seen as a reference facies model where geological knowledge about depositional

Numerical solution of multi groups point kinetic equations by simulink toolbox of Matlab software

International Nuclear Information System (INIS)

Hadad, K.; Mohamadi, A.; Sabet, H.; Ayobian, N.; Khani, M.

2004-01-01

The simulink toolbox of Matlab Software was employed to solve the point kinetics equation with six group delayed neutrons. The method of Adams-Bash ford showed a good convergence in solving the system of simultaneous equations and the obtained results showed good agreements with other numerical schemes. The flexibility of the package in changing the system parameters and the user friendly interface makes this approach a reliable educational package in revealing the affects of reactivity changes on power incursions

Derivation Method for the Foundation Boundaries of Hydraulic Numerical Simulation Models Based on the Elastic Boussinesq Solution

Directory of Open Access Journals (Sweden)

Jintao Song

2015-01-01

Full Text Available The foundation boundaries of numerical simulation models of hydraulic structures dominated by a vertical load are investigated. The method used is based on the stress formula for fundamental solutions to semi-infinite space body elastic mechanics under a vertical concentrated force. The limit method is introduced into the original formula, which is then partitioned and analyzed according to the direction of the depth extension of the foundation. The point load will be changed to a linear load with a length of 2a. Inverse proportion function assumptions are proposed at parameter a and depth l of the calculation points to solve the singularity questions of elastic stress in a semi-infinite space near the ground. Compared with the original formula, changing the point load to a linear load with a length of 2a is more reasonable. Finally, the boundary depth criterion of a hydraulic numerical simulation model is derived and applied to determine the depth boundary formula for gravity dam numerical simulations.

Throughput analysis of point-to-multi-point hybric FSO/RF network

KAUST Repository

Rakia, Tamer; Gebali, Fayez; Yang, Hong-Chuan; Alouini, Mohamed-Slim

2017-01-01

This paper presents and analyzes a point-to-multi-point (P2MP) network that uses a number of free-space optical (FSO) links for data transmission from the central node to the different remote nodes. A common backup radio-frequency (RF) link is used

Bayesian Estimation Of Shift Point In Poisson Model Under Asymmetric Loss Functions

Directory of Open Access Journals (Sweden)

uma srivastava

2012-01-01

Full Text Available The paper deals with estimating  shift point which occurs in any sequence of independent observations  of Poisson model in statistical process control. This shift point occurs in the sequence when  i.e. m  life data are observed. The Bayes estimator on shift point 'm' and before and after shift process means are derived for symmetric and asymmetric loss functions under informative and non informative priors. The sensitivity analysis of Bayes estimators are carried out by simulation and numerical comparisons with  R-programming. The results shows the effectiveness of shift in sequence of Poisson disribution .

Application of fixed point theory to chaotic attractors of forced oscillators

International Nuclear Information System (INIS)

Stewart, H.B.

1990-11-01

A review of the structure of chaotic attractors of periodically forced second order nonlinear oscillators suggests that the theory of fixed points of transformations gives information about the fundamental topological structure of attractors. First a simple extension of the Levinson index formula is proved. Then numerical evidence is used to formulate plausible conjectures about absorbing regions containing chaotic attractors in forced oscillators. Applying the Levinson formula suggests a fundamental relation between the number of fixed points or periodic points in a section of the chaotic attractor on the one hand, and a topological invariant of an absorbing region on the other hand. (author)

Numerical computation of inventory policies, based on the EOQ/sigma-x value for order-point systems

DEFF Research Database (Denmark)

Alstrøm, Poul

2001-01-01

This paper examines the numerical computation of two control parameters, order size and order point in the well-known inventory control model, an (s,Q)system with a beta safety strategy. The aim of the paper is to show that the EOQ/sigma-x value is both sufficient for controlling the system and e...

Numerical computation of inventory policies, based on the EOQ/sigma-x value for order-point systems

DEFF Research Database (Denmark)

Alstrøm, Poul

2000-01-01

This paper examines the numerical computation of two control parameters, order size and order point in the well-known inventory control model, an (s,Q)system with a beta safety strategy. The aim of the paper is to show that the EOQ/sigma-x value is both sufficient for controlling the system and e...

Critical points in magnetic systems

International Nuclear Information System (INIS)

Bongaarts, A.L.M.

1975-01-01

The magnetical phase transitions of CsCoCl 3 .2H 2 O and CsCoCl 3 .2D 2 O are investigated by neutron diffraction techniques with special attention to the critical points in the phase diagrams. CsCoCl 3 .2H 2 O turned out to be a one-dimentional magnetic antiferromagnet with ferromagnetic and antiferromagnetic interactions. In the vicinity of the Neel point, the critical behavior in zero magnetic field could be described as a three-dimentional long range ordering, while the fluctuations in the system are one-dimensional. In the presence of a magnetic field, the behavior of the system in the critical region of the magnetic phase diagram between the Neel temperature at zero field (3.3degK) and 1.85degK, was in good agreement with the theory. Below 1.85degK, the phase transition in a magnetic field changes into a line of triple points whose end point could be identified as a tricritical point, i.e., an intersection of three critical lines. The parameters derived from observations in the neighborhood of this tricritical point obey the scaling laws but are not in numerical agreement with theoretical predictions

Drift-based scrape-off particle width in X-point geometry

Science.gov (United States)

Reiser, D.; Eich, T.

2017-04-01

The Goldston heuristic estimate of the scrape-off layer width (Goldston 2012 Nucl. Fusion 52 013009) is reconsidered using a fluid description for the plasma dynamics. The basic ingredient is the inclusion of a compressible diamagnetic drift for the particle cross field transport. Instead of testing the heuristic model in a sophisticated numerical simulation including several physical mechanisms working together, the purpose of this work is to point out basic consequences for a drift-dominated cross field transport using a reduced fluid model. To evaluate the model equations and prepare them for subsequent numerical solution a specific analytical model for 2D magnetic field configurations with X-points is employed. In a first step parameter scans in high-resolution grids for isothermal plasmas are done to assess the basic formulas of the heuristic model with respect to the functional dependence of the scrape-off width on the poloidal magnetic field and plasma temperature. Particular features in the 2D-fluid calculations—especially the appearance of supersonic parallel flows and shock wave like bifurcational jumps—are discussed and can be understood partly in the framework of a reduced 1D model. The resulting semi-analytical findings might give hints for experimental proof and implementation in more elaborated fluid simulations.

Material-Point Analysis of Large-Strain Problems

DEFF Research Database (Denmark)

Andersen, Søren

The aim of this thesis is to apply and improve the material-point method for modelling of geotechnical problems. One of the geotechnical phenomena that is a subject of active research is the study of landslides. A large amount of research is focused on determining when slopes become unstable. Hence......, it is possible to predict if a certain slope is stable using commercial finite element or finite difference software such as PLAXIS, ABAQUS or FLAC. However, the dynamics during a landslide are less explored. The material-point method (MPM) is a novel numerical method aimed at analysing problems involving...... materials subjected to large strains in a dynamical time–space domain. This thesis explores the material-point method with the specific aim of improving the performance for geotechnical problems. Large-strain geotechnical problems such as landslides pose a major challenge to model numerically. Employing...

The validity of multiphase DNS initialized on the basis of single--point statistics

Science.gov (United States)

Subramaniam, Shankar

1999-11-01

A study of the point--process statistical representation of a spray reveals that single--point statistical information contained in the droplet distribution function (ddf) is related to a sequence of single surrogate--droplet pdf's, which are in general different from the physical single--droplet pdf's. The results of this study have important consequences for the initialization and evolution of direct numerical simulations (DNS) of multiphase flows, which are usually initialized on the basis of single--point statistics such as the average number density in physical space. If multiphase DNS are initialized in this way, this implies that even the initial representation contains certain implicit assumptions concerning the complete ensemble of realizations, which are invalid for general multiphase flows. Also the evolution of a DNS initialized in this manner is shown to be valid only if an as yet unproven commutation hypothesis holds true. Therefore, it is questionable to what extent DNS that are initialized in this manner constitute a direct simulation of the physical droplets.

Dynamics of Multibody Systems Near Lagrangian Points

Science.gov (United States)

Wong, Brian

This thesis examines the dynamics of a physically connected multi-spacecraft system in the vicinity of the Lagrangian points of a Circular Restricted Three-Body System. The spacecraft system is arranged in a wheel-spoke configuration with smaller and less massive satellites connected to a central hub using truss/beams or tether connectors. The kinematics of the system is first defined, and the kinetic, gravitational potential energy and elastic potential energy of the system are derived. The Assumed Modes Method is used to discretize the continuous variables of the system, and a general set of ordinary differential equations describing the dynamics of the connectors and the central hub are obtained using the Lagrangian method. The flexible body dynamics of the tethered and truss connected systems are examined using numerical simulations. The results show that these systems experienced only small elastic deflections when they are naturally librating or rotating at moderate angular velocities, and these deflections have relatively small effect on the attitude dynamics of the systems. Based on these results, it is determined that the connectors can be modeled as rigid when only the attitude dynamics of the system is of interest. The equations of motion of rigid satellites stationed at the Lagrangian points are linearized, and the stability conditions of the satellite are obtained from the linear equations. The required conditions are shown to be similar to those of geocentric satellites. Study of the linear equations also revealed the resonant conditions of rigid Lagrangian point satellites, when a librational natural frequency of the satellite matches the frequency of its station-keeping orbit leading to large attitude motions. For tethered satellites, the linear analysis shows that the tethers are in stable equilibrium when they lie along a line joining the two primary celestial bodies of the Three-Body System. Numerical simulations are used to study the long term

Orbital Maneuvers for Spacecrafts Travelling to/from the Lagrangian Points

Science.gov (United States)

Bertachini, A.

The well-known Lagrangian points that appear in the planar restricted three-body problem (Szebehely, 1967) are very important for astronautical applications. They are five points of equilibrium in the equations of motion, what means that a particle located at one of those points with zero velocity will remain there indefinitely. The collinear points (L1, L2 and L3) are always unstable and the triangular points (L4 and L5) are stable in the present case studied (Sun-Earth system). They are all very good points to locate a space-station, since they require a small amount of V (and fuel), the control to be used for station-keeping. The triangular points are specially good for this purpose, since they are stable equilibrium points. In this paper, the planar restricted three-body problem is regularized (using Lemaître regularization) and combined with numerical integration and gradient methods to solve the two point boundary value problem (the Lambert's three-body problem). This combination is applied to the search of families of transfer orbits between the Lagrangian points and the Earth, in the Sun-Earth system, with the minimum possible cost of the control used. So, the final goal of this paper is to find the magnitude and direction of the two impulses to be applied in the spacecraft to complete the transfer: the first one when leaving/arriving at the Lagrangian point and the second one when arriving/living at the Earth. This paper is a continuation of two previous papers that studied transfers in the Earth-Moon system: Broucke (1979), that studied transfer orbits between the Lagrangian points and the Moon and Prado (1996), that studied transfer orbits between the Lagrangian points and the Earth. So, the equations of motion are: whereis the pseudo-potential given by: To solve the TPBVP in the regularized variables the following steps are used: i) Guess a initial velocity Vi, so together with the initial prescribed position ri the complete initial state is known; ii

Improving multiple-point-based a priori models for inverse problems by combining Sequential Simulation with the Frequency Matching Method

DEFF Research Database (Denmark)

Cordua, Knud Skou; Hansen, Thomas Mejer; Lange, Katrine

In order to move beyond simplified covariance based a priori models, which are typically used for inverse problems, more complex multiple-point-based a priori models have to be considered. By means of marginal probability distributions ‘learned’ from a training image, sequential simulation has...... proven to be an efficient way of obtaining multiple realizations that honor the same multiple-point statistics as the training image. The frequency matching method provides an alternative way of formulating multiple-point-based a priori models. In this strategy the pattern frequency distributions (i.......e. marginals) of the training image and a subsurface model are matched in order to obtain a solution with the same multiple-point statistics as the training image. Sequential Gibbs sampling is a simulation strategy that provides an efficient way of applying sequential simulation based algorithms as a priori...

Direct numerical simulation of electrokinetic instability and transition to chaotic motion

Science.gov (United States)

Demekhin, E. A.; Nikitin, N. V.; Shelistov, V. S.

2013-12-01

A new type of instability—electrokinetic instability—and an unusual transition to chaotic motion near a charge-selective surface (semiselective electric membrane, electrode, or system of micro-/nanochannels) was studied by the numerical integration of the Nernst-Planck-Poisson-Stokes system and a weakly nonlinear analysis near the threshold of instability. A special finite-difference method was used for the space discretization along with a semi-implicit 31/3-step Runge-Kutta scheme for the integration in time. Two kinds of initial conditions were considered: (a) white-noise initial conditions to mimic "room disturbances" and subsequent natural evolution of the solution, and (b) an artificial monochromatic ion distribution with a fixed wave number to simulate regular wave patterns. The results were studied from the viewpoint of hydrodynamic stability and bifurcation theory. The threshold of electroconvective movement was found by the linear spectral stability theory, the results of which were confirmed by numerical simulation of the entire system. Our weakly nonlinear analysis and numerical integration of the entire system predict possibility of both kinds of bifurcations at the critical point, supercritical and subcritical, depending on the system parameters. The following regimes, which replace each other as the potential drop between the selective surfaces increases, were obtained: one-dimensional steady solution, two-dimensional steady electroconvective vortices (stationary point in a proper phase space), unsteady vortices aperiodically changing their parameters (homoclinic contour), periodic motion (limit cycle), and chaotic motion. The transition to chaotic motion does not include Hopf bifurcation. The numerical resolution of the thin concentration polarization layer showed spike-like charge profiles along the surface, which could be, depending on the regime, either steady or aperiodically coalescent. The numerical investigation confirmed the

Direct numerical simulation of electrokinetic instability and transition to chaotic motion

International Nuclear Information System (INIS)

Demekhin, E. A.; Nikitin, N. V.; Shelistov, V. S.

2013-01-01

A new type of instability—electrokinetic instability—and an unusual transition to chaotic motion near a charge-selective surface (semiselective electric membrane, electrode, or system of micro-/nanochannels) was studied by the numerical integration of the Nernst-Planck-Poisson-Stokes system and a weakly nonlinear analysis near the threshold of instability. A special finite-difference method was used for the space discretization along with a semi-implicit 31/3 -step Runge-Kutta scheme for the integration in time. Two kinds of initial conditions were considered: (a) white-noise initial conditions to mimic “room disturbances” and subsequent natural evolution of the solution, and (b) an artificial monochromatic ion distribution with a fixed wave number to simulate regular wave patterns. The results were studied from the viewpoint of hydrodynamic stability and bifurcation theory. The threshold of electroconvective movement was found by the linear spectral stability theory, the results of which were confirmed by numerical simulation of the entire system. Our weakly nonlinear analysis and numerical integration of the entire system predict possibility of both kinds of bifurcations at the critical point, supercritical and subcritical, depending on the system parameters. The following regimes, which replace each other as the potential drop between the selective surfaces increases, were obtained: one-dimensional steady solution, two-dimensional steady electroconvective vortices (stationary point in a proper phase space), unsteady vortices aperiodically changing their parameters (homoclinic contour), periodic motion (limit cycle), and chaotic motion. The transition to chaotic motion does not include Hopf bifurcation. The numerical resolution of the thin concentration polarization layer showed spike-like charge profiles along the surface, which could be, depending on the regime, either steady or aperiodically coalescent. The numerical investigation confirmed the

Direct numerical simulation of electrokinetic instability and transition to chaotic motion

Energy Technology Data Exchange (ETDEWEB)

Demekhin, E. A., E-mail: edemekhi@gmail.com [Laboratory of Micro- and Nanofluidics, Moscow State University, Moscow 119192 (Russian Federation); Department of Computation Mathematics and Computer Science, Kuban State University, Krasnodar 350040 (Russian Federation); Institute of Mechanics, Moscow State University, Moscow 117192 (Russian Federation); Nikitin, N. V. [Institute of Mechanics, Moscow State University, Moscow 117192 (Russian Federation); Shelistov, V. S. [Institute of Mechanics, Moscow State University, Moscow 117192 (Russian Federation); Scientific Research Department, Kuban State University, Krasnodar 350040 (Russian Federation)

2013-12-15

A new type of instability—electrokinetic instability—and an unusual transition to chaotic motion near a charge-selective surface (semiselective electric membrane, electrode, or system of micro-/nanochannels) was studied by the numerical integration of the Nernst-Planck-Poisson-Stokes system and a weakly nonlinear analysis near the threshold of instability. A special finite-difference method was used for the space discretization along with a semi-implicit 31/3 -step Runge-Kutta scheme for the integration in time. Two kinds of initial conditions were considered: (a) white-noise initial conditions to mimic “room disturbances” and subsequent natural evolution of the solution, and (b) an artificial monochromatic ion distribution with a fixed wave number to simulate regular wave patterns. The results were studied from the viewpoint of hydrodynamic stability and bifurcation theory. The threshold of electroconvective movement was found by the linear spectral stability theory, the results of which were confirmed by numerical simulation of the entire system. Our weakly nonlinear analysis and numerical integration of the entire system predict possibility of both kinds of bifurcations at the critical point, supercritical and subcritical, depending on the system parameters. The following regimes, which replace each other as the potential drop between the selective surfaces increases, were obtained: one-dimensional steady solution, two-dimensional steady electroconvective vortices (stationary point in a proper phase space), unsteady vortices aperiodically changing their parameters (homoclinic contour), periodic motion (limit cycle), and chaotic motion. The transition to chaotic motion does not include Hopf bifurcation. The numerical resolution of the thin concentration polarization layer showed spike-like charge profiles along the surface, which could be, depending on the regime, either steady or aperiodically coalescent. The numerical investigation confirmed the

Modeling hard clinical end-point data in economic analyses.

Science.gov (United States)

Kansal, Anuraag R; Zheng, Ying; Palencia, Roberto; Ruffolo, Antonio; Hass, Bastian; Sorensen, Sonja V

2013-11-01

The availability of hard clinical end-point data, such as that on cardiovascular (CV) events among patients with type 2 diabetes mellitus, is increasing, and as a result there is growing interest in using hard end-point data of this type in economic analyses. This study investigated published approaches for modeling hard end-points from clinical trials and evaluated their applicability in health economic models with different disease features. A review of cost-effectiveness models of interventions in clinically significant therapeutic areas (CV diseases, cancer, and chronic lower respiratory diseases) was conducted in PubMed and Embase using a defined search strategy. Only studies integrating hard end-point data from randomized clinical trials were considered. For each study included, clinical input characteristics and modeling approach were summarized and evaluated. A total of 33 articles (23 CV, eight cancer, two respiratory) were accepted for detailed analysis. Decision trees, Markov models, discrete event simulations, and hybrids were used. Event rates were incorporated either as constant rates, time-dependent risks, or risk equations based on patient characteristics. Risks dependent on time and/or patient characteristics were used where major event rates were >1%/year in models with fewer health states (Models of infrequent events or with numerous health states generally preferred constant event rates. The detailed modeling information and terminology varied, sometimes requiring interpretation. Key considerations for cost-effectiveness models incorporating hard end-point data include the frequency and characteristics of the relevant clinical events and how the trial data is reported. When event risk is low, simplification of both the model structure and event rate modeling is recommended. When event risk is common, such as in high risk populations, more detailed modeling approaches, including individual simulations or explicitly time-dependent event rates, are

Dynamic magnetic x-points

International Nuclear Information System (INIS)

Leboeuf, J.N.; Tajima, T.; Dawson, J.M.

1981-03-01

Two-and-one-half dimensional magnetostatic and electromagnetic particle simulations of time-varying magnetic x-points and the associated plasma response are reported. The stability and topology depend on the crossing angle of the field lines at the x-point, irrespective of the plasma β. The electrostatic field and finite Larmor radius effects play an important role in current penetration and shaping of the plasma flow. The snapping of the field lines, and dragging of the plasma into, and confinement of the plasma at, an o-point (magnetic island) is observed. Magnetic island coalescence with explosive growth of the coalescence mode occurs and is accompanied by a large increase of kinetic energy and temperature as well as the formation of hot tails on the distribution functions

Pricing and simulation for real estate index options: Radial basis point interpolation

Science.gov (United States)

Gong, Pu; Zou, Dong; Wang, Jiayue

2018-06-01

This study employs the meshfree radial basis point interpolation (RBPI) for pricing real estate derivatives contingent on real estate index. This method combines radial and polynomial basis functions, which can guarantee the interpolation scheme with Kronecker property and effectively improve accuracy. An exponential change of variables, a mesh refinement algorithm and the Richardson extrapolation are employed in this study to implement the RBPI. Numerical results are presented to examine the computational efficiency and accuracy of our method.

Characterization results and Markov chain Monte Carlo algorithms including exact simulation for some spatial point processes

DEFF Research Database (Denmark)

Häggström, Olle; Lieshout, Marie-Colette van; Møller, Jesper

1999-01-01

The area-interaction process and the continuum random-cluster model are characterized in terms of certain functional forms of their respective conditional intensities. In certain cases, these two point process models can be derived from a bivariate point process model which in many respects...... is simpler to analyse and simulate. Using this correspondence we devise a two-component Gibbs sampler, which can be used for fast and exact simulation by extending the recent ideas of Propp and Wilson. We further introduce a Swendsen-Wang type algorithm. The relevance of the results within spatial statistics...

Dynamics and mission design near libration points

CERN Document Server

Gómez, G; Simo, C; Martínez, R

2001-01-01

In this book the problem of station keeping is studied for orbits near libration points in the solar system. The main focus is on orbits near halo ones in the (Earth+Moon)-Sun system. Taking as starting point the restricted three-body problem, the motion in the full solar system is considered as a perturbation of this simplified model. All the study is done with enough generality to allow easy application to other primary-secondary systems as a simple extension of the analytical and numerical computations. Contents: Bibliographical Survey; Halo Orbits. Analytical and Numerical Study; The Neigh

The universal cardinal ordering of fixed points

International Nuclear Information System (INIS)

San Martin, Jesus; Moscoso, Ma Jose; Gonzalez Gomez, A.

2009-01-01

We present the theorem which determines, by a permutation, the cardinal ordering of fixed points for any orbit of a period doubling cascade. The inverse permutation generates the orbit and the symbolic sequence of the orbit is obtained as a corollary. Interestingly enough, it is important to point that this theorem needs no previous information about any other orbit; also the cardinal ordering is achieved automatically with no need to compare numerical values associated with every point of the orbit (as would be the case if kneading theory were used).

Origin of chaos near critical points of quantum flow.

Science.gov (United States)

Efthymiopoulos, C; Kalapotharakos, C; Contopoulos, G

2009-03-01

The general theory of motion in the vicinity of a moving quantum nodal point (vortex) is studied in the framework of the de Broglie-Bohm trajectory method of quantum mechanics. Using an adiabatic approximation, we find that near any nodal point of an arbitrary wave function psi there is an unstable point (called the X point) in a frame of reference moving with the nodal point. The local phase portrait forms always a characteristic pattern called the "nodal-point- X -point complex." We find general formulas for this complex as well as necessary and sufficient conditions of validity of the adiabatic approximation. We demonstrate that chaos emerges from the consecutive scattering events of the orbits with nodal-point- X -point complexes. The scattering events are of two types (called type I and type II). A theoretical model is constructed yielding the local value of the Lyapunov characteristic numbers in scattering events of both types. The local Lyapunov characteristic number scales as an inverse power of the speed of the nodal point in the rest frame, implying that it scales proportionally to the size of the nodal-point- X -point complex. It is also an inverse power of the distance of a trajectory from the X point's stable manifold far from the complex. This distance plays the role of an effective "impact parameter." The results of detailed numerical experiments with different wave functions, possessing one, two, or three moving nodal points, are reported. Examples are given of regular and chaotic trajectories, and the statistics of the Lyapunov characteristic numbers of the orbits are found and compared to the number of encounter events of each orbit with the nodal-point- X -point complexes. The numerical results are in agreement with the theory, and various phenomena appearing at first as counterintuitive find a straightforward explanation.

Research on numerical simulation technology about regional important pollutant diffusion of haze

Science.gov (United States)

Du, Boying; Ma, Yunfeng; Li, Qiangqiang; Wang, Qi; Hu, Qiongqiong; Bian, Yushan

2018-02-01

In order to analyze the formation of haze in Shenyang and the factors that affect the diffusion of pollutants, the simulation experiment adopted in this paper is based on the numerical model of WRF/CALPUFF coupling. Simulation experiment was conducted to select PM10 of Shenyang City in the period from March 1 to 8, and the PM10 in the regional important haze was simulated. The survey was conducted with more than 120 enterprises section the point of the emission source of this experiment. The contrastive data were analyzed with 11 air quality monitoring points, and the simulation results were compared. Analyze the contribution rate of each typical enterprise to the air quality, verify the correctness of the simulation results, and then use the model to establish the prediction model.

Limits in point to point resolution of MOS based pixels detector arrays

Science.gov (United States)

Fourches, N.; Desforge, D.; Kebbiri, M.; Kumar, V.; Serruys, Y.; Gutierrez, G.; Leprêtre, F.; Jomard, F.

2018-01-01

In high energy physics point-to-point resolution is a key prerequisite for particle detector pixel arrays. Current and future experiments require the development of inner-detectors able to resolve the tracks of particles down to the micron range. Present-day technologies, although not fully implemented in actual detectors, can reach a 5-μm limit, this limit being based on statistical measurements, with a pixel-pitch in the 10 μm range. This paper is devoted to the evaluation of the building blocks for use in pixel arrays enabling accurate tracking of charged particles. Basing us on simulations we will make here a quantitative evaluation of the physical and technological limits in pixel size. Attempts to design small pixels based on SOI technology will be briefly recalled here. A design based on CMOS compatible technologies that allow a reduction of the pixel size below the micrometer is introduced here. Its physical principle relies on a buried carrier-localizing collecting gate. The fabrication process needed by this pixel design can be based on existing process steps used in silicon microelectronics. The pixel characteristics will be discussed as well as the design of pixel arrays. The existing bottlenecks and how to overcome them will be discussed in the light of recent ion implantation and material characterization experiments.

Dynamical Response near Quantum Critical Points.

Science.gov (United States)

Lucas, Andrew; Gazit, Snir; Podolsky, Daniel; Witczak-Krempa, William

2017-02-03

We study high-frequency response functions, notably the optical conductivity, in the vicinity of quantum critical points (QCPs) by allowing for both detuning from the critical coupling and finite temperature. We consider general dimensions and dynamical exponents. This leads to a unified understanding of sum rules. In systems with emergent Lorentz invariance, powerful methods from quantum field theory allow us to fix the high-frequency response in terms of universal coefficients. We test our predictions analytically in the large-N O(N) model and using the gauge-gravity duality and numerically via quantum Monte Carlo simulations on a lattice model hosting the interacting superfluid-insulator QCP. In superfluid phases, interacting Goldstone bosons qualitatively change the high-frequency optical conductivity and the corresponding sum rule.

Algorithms for solving common fixed point problems

CERN Document Server

Zaslavski, Alexander J

2018-01-01

This book details approximate solutions to common fixed point problems and convex feasibility problems in the presence of perturbations. Convex feasibility problems search for a common point of a finite collection of subsets in a Hilbert space; common fixed point problems pursue a common fixed point of a finite collection of self-mappings in a Hilbert space. A variety of algorithms are considered in this book for solving both types of problems, the study of which has fueled a rapidly growing area of research. This monograph is timely and highlights the numerous applications to engineering, computed tomography, and radiation therapy planning. Totaling eight chapters, this book begins with an introduction to foundational material and moves on to examine iterative methods in metric spaces. The dynamic string-averaging methods for common fixed point problems in normed space are analyzed in Chapter 3. Dynamic string methods, for common fixed point problems in a metric space are introduced and discussed in Chapter ...

Confidence in Numerical Simulations

Energy Technology Data Exchange (ETDEWEB)

Hemez, Francois M. [Los Alamos National Lab. (LANL), Los Alamos, NM (United States)

2015-02-23

This PowerPoint presentation offers a high-level discussion of uncertainty, confidence and credibility in scientific Modeling and Simulation (M&S). It begins by briefly evoking M&S trends in computational physics and engineering. The first thrust of the discussion is to emphasize that the role of M&S in decision-making is either to support reasoning by similarity or to “forecast,” that is, make predictions about the future or extrapolate to settings or environments that cannot be tested experimentally. The second thrust is to explain that M&S-aided decision-making is an exercise in uncertainty management. The three broad classes of uncertainty in computational physics and engineering are variability and randomness, numerical uncertainty and model-form uncertainty. The last part of the discussion addresses how scientists “think.” This thought process parallels the scientific method where by a hypothesis is formulated, often accompanied by simplifying assumptions, then, physical experiments and numerical simulations are performed to confirm or reject the hypothesis. “Confidence” derives, not just from the levels of training and experience of analysts, but also from the rigor with which these assessments are performed, documented and peer-reviewed.

A New Approach for the Approximations of Solutions to a Common Fixed Point Problem in Metric Fixed Point Theory

Directory of Open Access Journals (Sweden)

Ishak Altun

2016-01-01

Full Text Available We provide sufficient conditions for the existence of a unique common fixed point for a pair of mappings T,S:X→X, where X is a nonempty set endowed with a certain metric. Moreover, a numerical algorithm is presented in order to approximate such solution. Our approach is different to the usual used methods in the literature.

New gaussian points for the solution of first order ordinary ...

African Journals Online (AJOL)

Numerical experiments carried out using the new Gaussian points revealed there efficiency on stiff differential equations. The results also reveal that methods using the new Gaussian points are more accurate than those using the standard Gaussian points on non-stiff initial value problems. Keywords: Gaussian points ...

Circuit Simulation for Solar Power Maximum Power Point Tracking with Different Buck-Boost Converter Topologies

Directory of Open Access Journals (Sweden)

Jaw-Kuen Shiau

2014-08-01

Full Text Available The power converter is one of the essential elements for effective use of renewable power sources. This paper focuses on the development of a circuit simulation model for maximum power point tracking (MPPT evaluation of solar power that involves using different buck-boost power converter topologies; including SEPIC, Zeta, and four-switch type buck-boost DC/DC converters. The circuit simulation model mainly includes three subsystems: a PV model; a buck-boost converter-based MPPT system; and a fuzzy logic MPPT controller. Dynamic analyses of the current-fed buck-boost converter systems are conducted and results are presented in the paper. The maximum power point tracking function is achieved through appropriate control of the power switches of the power converter. A fuzzy logic controller is developed to perform the MPPT function for obtaining maximum power from the PV panel. The MATLAB-based Simulink piecewise linear electric circuit simulation tool is used to verify the complete circuit simulation model.

A Generic Simulation Approach for the Fast and Accurate Estimation of the Outage Probability of Single Hop and Multihop FSO Links Subject to Generalized Pointing Errors

KAUST Repository

Ben Issaid, Chaouki; Park, Kihong; Alouini, Mohamed-Slim

2017-01-01

When assessing the performance of the free space optical (FSO) communication systems, the outage probability encountered is generally very small, and thereby the use of nave Monte Carlo simulations becomes prohibitively expensive. To estimate these rare event probabilities, we propose in this work an importance sampling approach which is based on the exponential twisting technique to offer fast and accurate results. In fact, we consider a variety of turbulence regimes, and we investigate the outage probability of FSO communication systems, under a generalized pointing error model based on the Beckmann distribution, for both single and multihop scenarios. Selected numerical simulations are presented to show the accuracy and the efficiency of our approach compared to naive Monte Carlo.

A Generic Simulation Approach for the Fast and Accurate Estimation of the Outage Probability of Single Hop and Multihop FSO Links Subject to Generalized Pointing Errors

KAUST Repository

Ben Issaid, Chaouki

2017-07-28

When assessing the performance of the free space optical (FSO) communication systems, the outage probability encountered is generally very small, and thereby the use of nave Monte Carlo simulations becomes prohibitively expensive. To estimate these rare event probabilities, we propose in this work an importance sampling approach which is based on the exponential twisting technique to offer fast and accurate results. In fact, we consider a variety of turbulence regimes, and we investigate the outage probability of FSO communication systems, under a generalized pointing error model based on the Beckmann distribution, for both single and multihop scenarios. Selected numerical simulations are presented to show the accuracy and the efficiency of our approach compared to naive Monte Carlo.

Analysis of point source size on measurement accuracy of lateral point-spread function of confocal Raman microscopy

Science.gov (United States)

Fu, Shihang; Zhang, Li; Hu, Yao; Ding, Xiang

2018-01-01

Confocal Raman Microscopy (CRM) has matured to become one of the most powerful instruments in analytical science because of its molecular sensitivity and high spatial resolution. Compared with conventional Raman Microscopy, CRM can perform three dimensions mapping of tiny samples and has the advantage of high spatial resolution thanking to the unique pinhole. With the wide application of the instrument, there is a growing requirement for the evaluation of the imaging performance of the system. Point-spread function (PSF) is an important approach to the evaluation of imaging capability of an optical instrument. Among a variety of measurement methods of PSF, the point source method has been widely used because it is easy to operate and the measurement results are approximate to the true PSF. In the point source method, the point source size has a significant impact on the final measurement accuracy. In this paper, the influence of the point source sizes on the measurement accuracy of PSF is analyzed and verified experimentally. A theoretical model of the lateral PSF for CRM is established and the effect of point source size on full-width at half maximum of lateral PSF is simulated. For long-term preservation and measurement convenience, PSF measurement phantom using polydimethylsiloxane resin, doped with different sizes of polystyrene microspheres is designed. The PSF of CRM with different sizes of microspheres are measured and the results are compared with the simulation results. The results provide a guide for measuring the PSF of the CRM.

Solution of the reactor point kinetics equations by MATLAB computing

Directory of Open Access Journals (Sweden)

Singh Sudhansu S.

2015-01-01

Full Text Available The numerical solution of the point kinetics equations in the presence of Newtonian temperature feedback has been a challenging issue for analyzing the reactor transients. Reactor point kinetics equations are a system of stiff ordinary differential equations which need special numerical treatments. Although a plethora of numerical intricacies have been introduced to solve the point kinetics equations over the years, some of the simple and straightforward methods still work very efficiently with extraordinary accuracy. As an example, it has been shown recently that the fundamental backward Euler finite difference algorithm with its simplicity has proven to be one of the most effective legacy methods. Complementing the back-ward Euler finite difference scheme, the present work demonstrates the application of ordinary differential equation suite available in the MATLAB software package to solve the stiff reactor point kinetics equations with Newtonian temperature feedback effects very effectively by analyzing various classic benchmark cases. Fair accuracy of the results implies the efficient application of MATLAB ordinary differential equation suite for solving the reactor point kinetics equations as an alternate method for future applications.

On the solution of two-point linear differential eigenvalue problems. [numerical technique with application to Orr-Sommerfeld equation

Science.gov (United States)

Antar, B. N.

1976-01-01

A numerical technique is presented for locating the eigenvalues of two point linear differential eigenvalue problems. The technique is designed to search for complex eigenvalues belonging to complex operators. With this method, any domain of the complex eigenvalue plane could be scanned and the eigenvalues within it, if any, located. For an application of the method, the eigenvalues of the Orr-Sommerfeld equation of the plane Poiseuille flow are determined within a specified portion of the c-plane. The eigenvalues for alpha = 1 and R = 10,000 are tabulated and compared for accuracy with existing solutions.

A Model Stitching Architecture for Continuous Full Flight-Envelope Simulation of Fixed-Wing Aircraft and Rotorcraft from Discrete Point Linear Models

Science.gov (United States)

2016-04-01

AND ROTORCRAFT FROM DISCRETE -POINT LINEAR MODELS Eric L. Tobias and Mark B. Tischler Aviation Development Directorate Aviation and Missile...Stitching Architecture for Continuous Full Flight-Envelope Simulation of Fixed-Wing Aircraft and Rotorcraft from Discrete -Point Linear Models 5...of discrete -point linear models and trim data. The model stitching simulation architecture is applicable to any aircraft configuration readily

Beginning SharePoint 2010 Administration Windows SharePoint Foundation 2010 and Microsoft SharePoint Server 2010

CERN Document Server

Husman, Göran

2010-01-01

Complete coverage on the latest advances in SharePoint 2010 administration. SharePoint 2010 comprises an abundance of new features, and this book shows you how to take advantage of all SharePoint 2010's many improvements. Written by a four-time SharePoint MVP, Beginning SharePoint 2010 Administration begins with a comparison of SharePoint 2010 compared to the previous version and then examines the differences between WSS 4.0 and MSS 2010. Packed with step-by-step instructions, tips and tricks, and real-world examples, this book dives into the basics of how to install, manage, and administrate

Error Mitigation of Point-to-Point Communication for Fault-Tolerant Computing

Science.gov (United States)

Akamine, Robert L.; Hodson, Robert F.; LaMeres, Brock J.; Ray, Robert E.

2011-01-01

Fault tolerant systems require the ability to detect and recover from physical damage caused by the hardware s environment, faulty connectors, and system degradation over time. This ability applies to military, space, and industrial computing applications. The integrity of Point-to-Point (P2P) communication, between two microcontrollers for example, is an essential part of fault tolerant computing systems. In this paper, different methods of fault detection and recovery are presented and analyzed.

Monte Carlo simulation of scatter in non-uniform symmetrical attenuating media for point and distributed sources

International Nuclear Information System (INIS)

Henry, L.J.; Rosenthal, M.S.

1992-01-01

We report results of scatter simulations for both point and distributed sources of 99m Tc in symmetrical non-uniform attenuating media. The simulations utilized Monte Carlo techniques and were tested against experimental phantoms. Both point and ring sources were used inside a 10.5 cm radius acrylic phantom. Attenuating media consisted of combinations of water, ground beef (to simulate muscle mass), air and bone meal (to simulate bone mass). We estimated/measured energy spectra, detector efficiencies and peak height ratios for all cases. In all cases, the simulated spectra agree with the experimentally measured spectra within 2 SD. Detector efficiencies and peak height ratios also are in agreement. The Monte Carlo code is able to properly model the non-uniform attenuating media used in this project. With verification of the simulations, it is possible to perform initial evaluation studies of scatter correction algorithms by evaluating the mechanisms of action of the correction algorithm on the simulated spectra where the magnitude and sources of scatter are known. (author)

Detecting Change-Point via Saddlepoint Approximations

Institute of Scientific and Technical Information of China (English)

Zhaoyuan LI; Maozai TIAN

2017-01-01

It's well-known that change-point problem is an important part of model statistical analysis.Most of the existing methods are not robust to criteria of the evaluation of change-point problem.In this article,we consider "mean-shift" problem in change-point studies.A quantile test of single quantile is proposed based on saddlepoint approximation method.In order to utilize the information at different quantile of the sequence,we further construct a "composite quantile test" to calculate the probability of every location of the sequence to be a change-point.The location of change-point can be pinpointed rather than estimated within a interval.The proposed tests make no assumptions about the functional forms of the sequence distribution and work sensitively on both large and small size samples,the case of change-point in the tails,and multiple change-points situation.The good performances of the tests are confirmed by simulations and real data analysis.The saddlepoint approximation based distribution of the test statistic that is developed in the paper is of independent interest and appealing.This finding may be of independent interest to the readers in this research area.

On Motion Planning for Point-to-Point Maneuvers for a Class of Sailing Vehicles

DEFF Research Database (Denmark)

Xiao, Lin; Jouffroy, Jerome

2011-01-01

Despite their interesting dynamic and controllability properties, sailing vehicles have not been much studied in the control community. In this paper, we investigate motion planning of such vehicles. Starting from a simple dynamic model of sailing vessels in one dimension, this paper first...... considers their associated controllability issues, with the so-called no-sailing zone as a starting point, and it links them with a motion planning strategy using two-point boundary value problems as the main mathematical tool. This perspective is then expanded to do point-to-point maneuvers of sailing...

Fermion-induced quantum critical points.

Science.gov (United States)

Li, Zi-Xiang; Jiang, Yi-Fan; Jian, Shao-Kai; Yao, Hong

2017-08-22

A unified theory of quantum critical points beyond the conventional Landau-Ginzburg-Wilson paradigm remains unknown. According to Landau cubic criterion, phase transitions should be first-order when cubic terms of order parameters are allowed by symmetry in the Landau-Ginzburg free energy. Here, from renormalization group analysis, we show that second-order quantum phase transitions can occur at such putatively first-order transitions in interacting two-dimensional Dirac semimetals. As such type of Landau-forbidden quantum critical points are induced by gapless fermions, we call them fermion-induced quantum critical points. We further introduce a microscopic model of SU(N) fermions on the honeycomb lattice featuring a transition between Dirac semimetals and Kekule valence bond solids. Remarkably, our large-scale sign-problem-free Majorana quantum Monte Carlo simulations show convincing evidences of a fermion-induced quantum critical points for N = 2, 3, 4, 5 and 6, consistent with the renormalization group analysis. We finally discuss possible experimental realizations of the fermion-induced quantum critical points in graphene and graphene-like materials.Quantum phase transitions are governed by Landau-Ginzburg theory and the exceptions are rare. Here, Li et al. propose a type of Landau-forbidden quantum critical points induced by gapless fermions in two-dimensional Dirac semimetals.

On a closed form solution of the point kinetics equations with reactivity feedback of temperature

International Nuclear Information System (INIS)

Silva, Jeronimo J.A.; Vilhena, Marco T.M.B.; Petersen, Claudio Z.; Bodmann, Bardo E.J.; Alvim, Antonio C.M.

2011-01-01

An analytical solution of the point kinetics equations to calculate reactivity as a function of time by the Decomposition method has recently appeared in the literature. In this paper, we go one step forward, by considering the neutron point kinetics equations together with temperature feedback effects. To accomplish that, we extended the point kinetics by a temperature perturbation, obtaining a second order nonlinear ordinary differential equation. This equation is then solved by the Decomposition Method, that is, by expanding the neutron density in a series and the nonlinear terms into Adomian Polynomials. Substituting these expansions into the nonlinear ordinary equation, we construct a recursive set of linear problems that can be solved by the methodology previously mentioned for the point kinetics equation. We also report on numerical simulations and comparisons against literature results. (author)

Numerical investigation of the energy performance of a guideless irregular heat and mass exchanger with corrugated heat transfer surface for dew point cooling

International Nuclear Information System (INIS)

Xu, Peng; Ma, Xiaoli; Diallo, Thierno M.O.; Zhao, Xudong; Fancey, Kevin; Li, Deying; Chen, Hongbing

2016-01-01

The paper presents an investigation into the energy performance of a novel irregular heat and mass exchanger for dew point cooling which, compared to the existing flat-plate heat exchangers, removed the use of the channel supporting guides and implemented the corrugated heat transfer surface, thus expecting to achieve the reduced air flow resistance, increased heat transfer area, and improved energy efficiency (i.e. Coefficient of Performance (COP)) of the air cooling process. CFD simulation was carried out to determine the flow resistance (K) factors of various elements within the dry and wet channels of the exchanger, while the ‘finite-element’ based ‘Newton-iteration’ numerical simulation was undertaken to investigate its cooling capacity, cooling effectiveness and COP at various geometrical and operational conditions. Compared to the existing flat-plate heat and mass exchangers with the same geometrical dimensions and operational conditions, the new irregular exchanger could achieve 32.9%–37% higher cooling capacity, dew-point and wet-bulb effectiveness, 29.7%–33.3% higher COP, and 55.8%–56.2% lower pressure drop. While undertaking dew point air cooling, the irregular heat and mass exchanger had the optimum air velocity of 1 m/s within the flow channels and working-to-intake air ratio of 0.3, which allowed the highest cooling capacity and COP to be achieved. In terms of the exchanger dimensions, the optimum height of the channel was 5 mm while its length was in the range 1–2 m. Overall, the proposed irregular heat and mass exchanger could lead to significant enhanced energy performance compared to the existing flat-plate dew point cooling heat exchanger of the same geometrical dimensions. To achieve the same amount cooling output, the irregular heat and mass exchanger had the reduced size and cost against the flat-plate ones. - Highlights: • Numerical investigation of an irregular heat and mass exchanger was undertaken. • A �

Flat Knitting Loop Deformation Simulation Based on Interlacing Point Model

Directory of Open Access Journals (Sweden)

Jiang Gaoming

2017-12-01

Full Text Available In order to create realistic loop primitives suitable for the faster CAD of the flat-knitted fabric, we have performed research on the model of the loop as well as the variation of the loop surface. This paper proposes an interlacing point-based model for the loop center curve, and uses the cubic Bezier curve to fit the central curve of the regular loop, elongated loop, transfer loop, and irregular deformed loop. In this way, a general model for the central curve of the deformed loop is obtained. The obtained model is then utilized to perform texture mapping, texture interpolation, and brightness processing, simulating a clearly structured and lifelike deformed loop. The computer program LOOP is developed by using the algorithm. The deformed loop is simulated with different yarns, and the deformed loop is applied to design of a cable stitch, demonstrating feasibility of the proposed algorithm. This paper provides a loop primitive simulation method characterized by lifelikeness, yarn material variability, and deformation flexibility, and facilitates the loop-based fast computer-aided design (CAD of the knitted fabric.

Asymmetrical floating point array processors, their application to exploration and exploitation

Energy Technology Data Exchange (ETDEWEB)

Geriepy, B L

1983-01-01

An asymmetrical floating point array processor is a special-purpose scientific computer which operates under asymmetrical control of a host computer. Although an array processor can receive fixed point input and produce fixed point output, its primary mode of operation is floating point. The first generation of array processors was oriented towards time series information. The next generation of array processors has proved much more versatile and their applicability ranges from petroleum reservoir simulation to speech syntheses. Array processors are becoming commonplace in mining, the primary usage being construction of grids-by usual methods or by kriging. The Australian mining community is among the world's leaders in regard to computer-assisted exploration and exploitation systems. Part of this leadership role must be providing guidance to computer vendors in regard to current and future requirements.

The Tipping Points of Technology Development

Directory of Open Access Journals (Sweden)

Tauno Kekäle

2014-07-01

Full Text Available The tipping point, the decisive point in time in the competition between old and new, is an interesting phenomenon in physics of today. This aspect in technology acceptance is connected to many business decisions such as technology investments, product releases, resource allocation, sales forecasts and, ultimately, affects the profitability and even survival of a company. The tipping point itself is based on many stochastic and dynamic variables, and the process may at least partly be described as path-dependent. This paper analyses the tipping point from three aspects: (1 product performance, (2 features of the market and infrastructure (including related technologies and human network externalities, and (3 actions of the incumbents (including customer lock-in, systems lock-in, and sustaining innovation. The paper is based on the Bass s-curve idea and the technology trajectory concept proposed by Dosi. Three illustrative cases are presented to make the point of the multiple factors affecting technology acceptance and, thus, the tipping point. The paper also suggests outlines for further research in field of computer simulation.

The robustness of subcarrier-index modulation in 16-QAM CO-OFDM system with 1024-point FFT.

Science.gov (United States)

Jan, Omar H A; Sandel, David; Puntsri, Kidsanapong; Al-Bermani, Ali; El-Darawy, Mohamed; Noé, Reinhold

2012-12-17

We present in numerical simulations the robustness of subcarrier index modulation (SIM) OFDM to combat laser phase noise. The ability of using DFB lasers with SIM-OFDM in 16-QAM CO-OFDM system with 1024-point FFT has been verified. Although SIM-OFDM has lower spectral efficiency compared to the conventional CO-OFDM system, it is a good candidate for 16-QAM CO-OFDM system with 1024-point FFT which uses a DFB laser of 1 MHz linewidth. In addition, we show the tolerance of SIM-OFDM for mitigation of fiber nonlinearities in long-haul CO-OFDM system. The simulation results show a significant penalty reduction, essentially that due to SPM.

Diffuse scattering from crystals with point defects

International Nuclear Information System (INIS)

Andrushevsky, N.M.; Shchedrin, B.M.; Simonov, V.I.; Malakhova, L.F.

2002-01-01

The analytical expressions for calculating the intensities of X-ray diffuse scattering from a crystal of finite dimensions and monatomic substitutional, interstitial, or vacancy-type point defects have been derived. The method for the determination of the three-dimensional structure by experimental diffuse-scattering data from crystals with point defects having various concentrations is discussed and corresponding numerical algorithms are suggested

The computation of fixed points and applications

CERN Document Server

Todd, Michael J

1976-01-01

Fixed-point algorithms have diverse applications in economics, optimization, game theory and the numerical solution of boundary-value problems. Since Scarf's pioneering work [56,57] on obtaining approximate fixed points of continuous mappings, a great deal of research has been done in extending the applicability and improving the efficiency of fixed-point methods. Much of this work is available only in research papers, although Scarf's book [58] gives a remarkably clear exposition of the power of fixed-point methods. However, the algorithms described by Scarf have been super~eded by the more sophisticated restart and homotopy techniques of Merrill [~8,~9] and Eaves and Saigal [1~,16]. To understand the more efficient algorithms one must become familiar with the notions of triangulation and simplicial approxi- tion, whereas Scarf stresses the concept of primitive set. These notes are intended to introduce to a wider audience the most recent fixed-point methods and their applications. Our approach is therefore ...

Optimal Power Flow by Interior Point and Non Interior Point Modern Optimization Algorithms

Directory of Open Access Journals (Sweden)

Marcin Połomski

2013-03-01

Full Text Available The idea of optimal power flow (OPF is to determine the optimal settings for control variables while respecting various constraints, and in general it is related to power system operational and planning optimization problems. A vast number of optimization methods have been applied to solve the OPF problem, but their performance is highly dependent on the size of a power system being optimized. The development of the OPF recently has tracked significant progress both in numerical optimization techniques and computer techniques application. In recent years, application of interior point methods to solve OPF problem has been paid great attention. This is due to the fact that IP methods are among the fastest algorithms, well suited to solve large-scale nonlinear optimization problems. This paper presents the primal-dual interior point method based optimal power flow algorithm and new variant of the non interior point method algorithm with application to optimal power flow problem. Described algorithms were implemented in custom software. The experiments show the usefulness of computational software and implemented algorithms for solving the optimal power flow problem, including the system model sizes comparable to the size of the National Power System.

Cognitive load predicts point-of-care ultrasound simulator performance.

Science.gov (United States)

Aldekhyl, Sara; Cavalcanti, Rodrigo B; Naismith, Laura M

2018-02-01

The ability to maintain good performance with low cognitive load is an important marker of expertise. Incorporating cognitive load measurements in the context of simulation training may help to inform judgements of competence. This exploratory study investigated relationships between demographic markers of expertise, cognitive load measures, and simulator performance in the context of point-of-care ultrasonography. Twenty-nine medical trainees and clinicians at the University of Toronto with a range of clinical ultrasound experience were recruited. Participants answered a demographic questionnaire then used an ultrasound simulator to perform targeted scanning tasks based on clinical vignettes. Participants were scored on their ability to both acquire and interpret ultrasound images. Cognitive load measures included participant self-report, eye-based physiological indices, and behavioural measures. Data were analyzed using a multilevel linear modelling approach, wherein observations were clustered by participants. Experienced participants outperformed novice participants on ultrasound image acquisition. Ultrasound image interpretation was comparable between the two groups. Ultrasound image acquisition performance was predicted by level of training, prior ultrasound training, and cognitive load. There was significant convergence between cognitive load measurement techniques. A marginal model of ultrasound image acquisition performance including prior ultrasound training and cognitive load as fixed effects provided the best overall fit for the observed data. In this proof-of-principle study, the combination of demographic and cognitive load measures provided more sensitive metrics to predict ultrasound simulator performance. Performance assessments which include cognitive load can help differentiate between levels of expertise in simulation environments, and may serve as better predictors of skill transfer to clinical practice.

Dose point kernels for beta-emitting radioisotopes

International Nuclear Information System (INIS)

Prestwich, W.V.; Chan, L.B.; Kwok, C.S.; Wilson, B.

1986-01-01

Knowledge of the dose point kernel corresponding to a specific radionuclide is required to calculate the spatial dose distribution produced in a homogeneous medium by a distributed source. Dose point kernels for commonly used radionuclides have been calculated previously using as a basis monoenergetic dose point kernels derived by numerical integration of a model transport equation. The treatment neglects fluctuations in energy deposition, an effect which has been later incorporated in dose point kernels calculated using Monte Carlo methods. This work describes new calculations of dose point kernels using the Monte Carlo results as a basis. An analytic representation of the monoenergetic dose point kernels has been developed. This provides a convenient method both for calculating the dose point kernel associated with a given beta spectrum and for incorporating the effect of internal conversion. An algebraic expression for allowed beta spectra has been accomplished through an extension of the Bethe-Bacher approximation, and tested against the exact expression. Simplified expression for first-forbidden shape factors have also been developed. A comparison of the calculated dose point kernel for 32 P with experimental data indicates good agreement with a significant improvement over the earlier results in this respect. An analytic representation of the dose point kernel associated with the spectrum of a single beta group has been formulated. 9 references, 16 figures, 3 tables

Seeing your error alters my pointing: observing systematic pointing errors induces sensori-motor after-effects.

Directory of Open Access Journals (Sweden)

Roberta Ronchi

Full Text Available During the procedure of prism adaptation, subjects execute pointing movements to visual targets under a lateral optical displacement: as consequence of the discrepancy between visual and proprioceptive inputs, their visuo-motor activity is characterized by pointing errors. The perception of such final errors triggers error-correction processes that eventually result into sensori-motor compensation, opposite to the prismatic displacement (i.e., after-effects. Here we tested whether the mere observation of erroneous pointing movements, similar to those executed during prism adaptation, is sufficient to produce adaptation-like after-effects. Neurotypical participants observed, from a first-person perspective, the examiner's arm making incorrect pointing movements that systematically overshot visual targets location to the right, thus simulating a rightward optical deviation. Three classical after-effect measures (proprioceptive, visual and visual-proprioceptive shift were recorded before and after first-person's perspective observation of pointing errors. Results showed that mere visual exposure to an arm that systematically points on the right-side of a target (i.e., without error correction produces a leftward after-effect, which mostly affects the observer's proprioceptive estimation of her body midline. In addition, being exposed to such a constant visual error induced in the observer the illusion "to feel" the seen movement. These findings indicate that it is possible to elicit sensori-motor after-effects by mere observation of movement errors.

Wrox SharePoint 2010 SharePoint911 three-pack

CERN Document Server

Klindt, Todd; Mason, Jennifer; Rogers, Laura; Drisgill, Randy; Ross, John; Riemann, Larry; Perran, Amanda; Perran, Shane; Sanford, Jacob J; Stubbs, Paul; Caravajal, Steve

2012-01-01

The Wrox SharePoint 2010 SharePoint911 Three-Pack combines the contents of three full e-books written by the experts from SharePoint911.  That's over 1800 pages of hands-on advice from Todd Klindt, Shane Young, Laura Rogers, Randy Drisgill, Jennifer Mason, John Ross, and Larry Riemann, among others. In Beginning SharePoint 2010: Building Business Solutions with SharePoint (ISBN 978-0-470-61789-2) by Amanda Perran, Shane Perran, Jennifer Mason, and Laura Rogers, readers learn the core concepts, terminology, and features of SharePoint 2010. In Professiona

The timing of control signals underlying fast point-to-point arm movements.

Science.gov (United States)

Ghafouri, M; Feldman, A G

2001-04-01

It is known that proprioceptive feedback induces muscle activation when the facilitation of appropriate motoneurons exceeds their threshold. In the suprathreshold range, the muscle-reflex system produces torques depending on the position and velocity of the joint segment(s) that the muscle spans. The static component of the torque-position relationship is referred to as the invariant characteristic (IC). According to the equilibrium-point (EP) hypothesis, control systems produce movements by changing the activation thresholds and thus shifting the IC of the appropriate muscles in joint space. This control process upsets the balance between muscle and external torques at the initial limb configuration and, to regain the balance, the limb is forced to establish a new configuration or, if the movement is prevented, a new level of static torques. Taken together, the joint angles and the muscle torques generated at an equilibrium configuration define a single variable called the EP. Thus by shifting the IC, control systems reset the EP. Muscle activation and movement emerge following the EP resetting because of the natural physical tendency of the system to reach equilibrium. Empirical and simulation studies support the notion that the control IC shifts and the resulting EP shifts underlying fast point-to-point arm movements are gradual rather than step-like. However, controversies exist about the duration of these shifts. Some studies suggest that the IC shifts cease with the movement offset. Other studies propose that the IC shifts end early in comparison to the movement duration (approximately, at peak velocity). The purpose of this study was to evaluate the duration of the IC shifts underlying fast point-to-point arm movements. Subjects made fast (hand peak velocity about 1.3 m/s) planar arm movements toward different targets while grasping a handle. Hand forces applied to the handle and shoulder/elbow torques were, respectively, measured from a force sensor placed

Beginning SharePoint 2010 Building Business Solutions with SharePoint

CERN Document Server

Perran, Amanda; Mason, Jennifer; Rogers, Laura

2010-01-01

Two SharePoint MVPs provide the ultimate introduction to SharePoint 2010Beginning SharePoint 2010: Building Team Solutions with SharePoint provides information workers and site managers with extensive knowledge and expert advice, empowering them to become SharePoint champions within their organizations.Provides expansive coverage of SharePoint topics, as well as specialty areas such as forms, excel services, records management, and web content managementDetails realistic usage scenarios, and includes practice examples that highlight best practices for configuration and customizationIncludes de

The algorithm to generate color point-cloud with the registration between panoramic image and laser point-cloud

International Nuclear Information System (INIS)

Zeng, Fanyang; Zhong, Ruofei

2014-01-01

Laser point cloud contains only intensity information and it is necessary for visual interpretation to obtain color information from other sensor. Cameras can provide texture, color, and other information of the corresponding object. Points with color information of corresponding pixels in digital images can be used to generate color point-cloud and is conducive to the visualization, classification and modeling of point-cloud. Different types of digital cameras are used in different Mobile Measurement Systems (MMS).the principles and processes for generating color point-cloud in different systems are not the same. The most prominent feature of the panoramic images is the field of 360 degrees view angle in the horizontal direction, to obtain the image information around the camera as much as possible. In this paper, we introduce a method to generate color point-cloud with panoramic image and laser point-cloud, and deduce the equation of the correspondence between points in panoramic images and laser point-clouds. The fusion of panoramic image and laser point-cloud is according to the collinear principle of three points (the center of the omnidirectional multi-camera system, the image point on the sphere, the object point). The experimental results show that the proposed algorithm and formulae in this paper are correct

Effect of Antenna Pointing Errors on SAR Imaging Considering the Change of the Point Target Location

Science.gov (United States)

Zhang, Xin; Liu, Shijie; Yu, Haifeng; Tong, Xiaohua; Huang, Guoman

2018-04-01

Towards spaceborne spotlight SAR, the antenna is regulated by the SAR system with specific regularity, so the shaking of the internal mechanism is inevitable. Moreover, external environment also has an effect on the stability of SAR platform. Both of them will cause the jitter of the SAR platform attitude. The platform attitude instability will introduce antenna pointing error on both the azimuth and range directions, and influence the acquisition of SAR original data and ultimate imaging quality. In this paper, the relations between the antenna pointing errors and the three-axis attitude errors are deduced, then the relations between spaceborne spotlight SAR imaging of the point target and antenna pointing errors are analysed based on the paired echo theory, meanwhile, the change of the azimuth antenna gain is considered as the spotlight SAR platform moves ahead. The simulation experiments manifest the effects on spotlight SAR imaging caused by antenna pointing errors are related to the target location, that is, the pointing errors of the antenna beam will severely influence the area far away from the scene centre of azimuth direction in the illuminated scene.

Modelling The Effect of Changing Point Systems to Teams’ Competition Standing in A Malaysian Soccer Super League

Science.gov (United States)

Mat Yusof, Muhammad; Khalid, Ruzelan; Hamid, Mohamad Shukri Abdul; Mansor, Rosnalini; Sulaiman, Tajularipin

2018-05-01

In a sports league such as in a soccer league, the teams’ competition standing is based on a cumulative point system. Typically, the standard point system is given to every single match for win, draw and lose teams is the 3-1-0 point system. In this paper, we explore the effect of changing point systems to teams’ competition standing by changing the weightage values for win, draw and lose teams. Three types of point systems are explored in our soccer simulation model; firstly the 3-1-0, secondly the 2-1-0 and thirdly the 4-1-0 point system. Based on the teams participating in a Malaysian soccer Super League, our simulation result shows that there are small changes in term of teams’ competition standing when we compared the actual rank and the simulation rank position. However, the 4-1-0 point system recorded the highest Pearson correlation value which is 0.97, followed by the 2-1-0 point system (0.95) and thirdly the 3-1-0 point system (0.94).

Confidence in Numerical Simulations

International Nuclear Information System (INIS)

Hemez, Francois M.

2015-01-01

This PowerPoint presentation offers a high-level discussion of uncertainty, confidence and credibility in scientific Modeling and Simulation (M&S). It begins by briefly evoking M&S trends in computational physics and engineering. The first thrust of the discussion is to emphasize that the role of M&S in decision-making is either to support reasoning by similarity or to ''forecast,'' that is, make predictions about the future or extrapolate to settings or environments that cannot be tested experimentally. The second thrust is to explain that M&S-aided decision-making is an exercise in uncertainty management. The three broad classes of uncertainty in computational physics and engineering are variability and randomness, numerical uncertainty and model-form uncertainty. The last part of the discussion addresses how scientists ''think.'' This thought process parallels the scientific method where by a hypothesis is formulated, often accompanied by simplifying assumptions, then, physical experiments and numerical simulations are performed to confirm or reject the hypothesis. ''Confidence'' derives, not just from the levels of training and experience of analysts, but also from the rigor with which these assessments are performed, documented and peer-reviewed.

Thinning spatial point processes into Poisson processes

DEFF Research Database (Denmark)

Møller, Jesper; Schoenberg, Frederic Paik

, and where one simulates backwards and forwards in order to obtain the thinned process. In the case of a Cox process, a simple independent thinning technique is proposed. In both cases, the thinning results in a Poisson process if and only if the true Papangelou conditional intensity is used, and thus can......This paper describes methods for randomly thinning certain classes of spatial point processes. In the case of a Markov point process, the proposed method involves a dependent thinning of a spatial birth-and-death process, where clans of ancestors associated with the original points are identified...... be used as a diagnostic for assessing the goodness-of-fit of a spatial point process model. Several examples, including clustered and inhibitive point processes, are considered....

Accelerating simulation for the multiple-point statistics algorithm using vector quantization

Science.gov (United States)

Zuo, Chen; Pan, Zhibin; Liang, Hao

2018-03-01

Multiple-point statistics (MPS) is a prominent algorithm to simulate categorical variables based on a sequential simulation procedure. Assuming training images (TIs) as prior conceptual models, MPS extracts patterns from TIs using a template and records their occurrences in a database. However, complex patterns increase the size of the database and require considerable time to retrieve the desired elements. In order to speed up simulation and improve simulation quality over state-of-the-art MPS methods, we propose an accelerating simulation for MPS using vector quantization (VQ), called VQ-MPS. First, a variable representation is presented to make categorical variables applicable for vector quantization. Second, we adopt a tree-structured VQ to compress the database so that stationary simulations are realized. Finally, a transformed template and classified VQ are used to address nonstationarity. A two-dimensional (2D) stationary channelized reservoir image is used to validate the proposed VQ-MPS. In comparison with several existing MPS programs, our method exhibits significantly better performance in terms of computational time, pattern reproductions, and spatial uncertainty. Further demonstrations consist of a 2D four facies simulation, two 2D nonstationary channel simulations, and a three-dimensional (3D) rock simulation. The results reveal that our proposed method is also capable of solving multifacies, nonstationarity, and 3D simulations based on 2D TIs.

Metal Carbon Eutectics to Extend the Use of the Fixed-Point Technique in Precision IR Thermometry

Science.gov (United States)

Battuello, M.; Girard, F.; Florio, M.

2008-06-01

The high-temperature extension of the fixed-point technique for primary calibration of precision infrared (IR) thermometers was investigated both through mathematical simulations and laboratory investigations. Simulations were performed with Co C (1,324°C) and Pd C (1, 492°C) eutectic fixed points, and a precision IR thermometer was calibrated from the In point (156.5985°C) up to the Co C point. Mathematical simulations suggested the possibility of directly deriving the transition temperature of the Co C and Pd C points by extrapolating the calibration derived from fixed-point measurements from In to the Cu point. Both temperatures, as a result of the low uncertainty associated with the In Cu calibration and the high number of fixed points involved in the calibration process, can be derived with an uncertainty of 0.11°C for Co C and 0.18°C for Pd C. A transition temperature of 1,324.3°C for Co C was determined from the experimental verification, a value higher than, but compatible with, the one proposed by the thermometry community for inclusion as a secondary reference point for ITS-90 dissemination, i.e., 1,324.0°C.

Network Competition - the Coexistence of Hub-and-Spoke and Point-to-Point Systems

NARCIS (Netherlands)

Alderighi, M.; Cento, A.; Nijkamp, P.; Rietveld, P.

2005-01-01

The paper identifies conditions under which asymmetric equilibria may exist when carriers compete in designing their network configurations in a game-theoretical framework. Two carriers are assumed here, which are allowed to play three different strategies: Point-to-point (PP), hub-and-spoke (HS) or

Thermal Analysis of a Cracked Half-plane under Moving Point Heat Source

Directory of Open Access Journals (Sweden)

He Kuanfang

2017-09-01

Full Text Available The heat conduction in half-plane with an insulated crack subjected to moving point heat source is investigated. The analytical solution and the numerical means are combined to analyze the transient temperature distribution of a cracked half-plane under moving point heat source. The transient temperature distribution of the half plane structure under moving point heat source is obtained by the moving coordinate method firstly, then the heat conduction equation with thermal boundary of an insulated crack face is changed to singular integral equation by applying Fourier transforms and solved by the numerical method. The numerical examples of the temperature distribution on the cracked half-plane structure under moving point heat source are presented and discussed in detail.

SIMULATION OF NEW SIMPLE FUZZY LOGIC MAXIMUM POWER POINT TRACKER FOR PHOTOVOLTAIC ARRAY

Directory of Open Access Journals (Sweden)

H. Serhoud

2015-08-01

Full Text Available A new simple fuzzy method used for tracking the maximum power point tracker (MPPT for photovoltaic systems is proposed. The input parameters   and duty cycle D are used to generate the optimal MPPT under different operating conditions, The photovoltaic system simulated and constructed by photovoltaic arrays, a DC/DC boost converter, a fuzzy MPPT control and a resistive load, The Fuzzy control law designed and the results in a simulation platform will be presented and compare to Perturbation and observation (P&O controller.

A highly accurate algorithm for the solution of the point kinetics equations

International Nuclear Information System (INIS)

Ganapol, B.D.

2013-01-01

Highlights: • Point kinetics equations for nuclear reactor transient analysis are numerically solved to extreme accuracy. • Results for classic benchmarks found in the literature are given to 9-digit accuracy. • Recent results of claimed accuracy are shown to be less accurate than claimed. • Arguably brings a chapter of numerical evaluation of the PKEs to a close. - Abstract: Attempts to resolve the point kinetics equations (PKEs) describing nuclear reactor transients have been the subject of numerous articles and texts over the past 50 years. Some very innovative methods, such as the RTS (Reactor Transient Simulation) and CAC (Continuous Analytical Continuation) methods of G.R. Keepin and J. Vigil respectively, have been shown to be exceptionally useful. Recently however, several authors have developed methods they consider accurate without a clear basis for their assertion. In response, this presentation will establish a definitive set of benchmarks to enable those developing PKE methods to truthfully assess the degree of accuracy of their methods. Then, with these benchmarks, two recently published methods, found in this journal will be shown to be less accurate than claimed and a legacy method from 1984 will be confirmed

Determine point-to-point networking interactions using regular expressions

Directory of Open Access Journals (Sweden)

Konstantin S. Deev

2015-06-01

Full Text Available As Internet growth and becoming more popular, the number of concurrent data flows start to increasing, which makes sense in bandwidth requested. Providers and corporate customers need ability to identify point-to-point interactions. The best is to use special software and hardware implementations that distribute the load in the internals of the complex, using the principles and approaches, in particular, described in this paper. This paper represent the principles of building system, which searches for a regular expression match using computing on graphics adapter in server station. A significant computing power and capability to parallel execution on modern graphic processor allows inspection of large amounts of data through sets of rules. Using the specified characteristics can lead to increased computing power in 30…40 times compared to the same setups on the central processing unit. The potential increase in bandwidth capacity could be used in systems that provide packet analysis, firewalls and network anomaly detectors.

Point spread function due to multiple scattering of light in the atmosphere

International Nuclear Information System (INIS)

Pękala, J.; Wilczyński, H.

2013-01-01

The atmospheric scattering of light has a significant influence on the results of optical observations of air showers. It causes attenuation of direct light from the shower, but also contributes a delayed signal to the observed light. The scattering of light therefore should be accounted for, both in simulations of air shower detection and reconstruction of observed events. In this work a Monte Carlo simulation of multiple scattering of light has been used to determine the contribution of the scattered light in observations of a point source of light. Results of the simulations and a parameterization of the angular distribution of the scattered light contribution to the observed signal (the point spread function) are presented. -- Author-Highlights: •Analysis of atmospheric scattering of light from an isotropic point source. •Different geometries and atmospheric conditions were investigated. •A parameterization of scattered light distribution has been developed. •The parameterization allows one to easily account for the light scattering in air. •The results will be useful in analyses of observations of extensive air shower

Maximum entropy state of the quasi-geostrophic bi-disperse point vortex system: bifurcation phenomena under periodic boundary conditions

Energy Technology Data Exchange (ETDEWEB)

Funakoshi, Satoshi; Sato, Tomoyoshi; Miyazaki, Takeshi, E-mail: funakosi@miyazaki.mce.uec.ac.jp, E-mail: miyazaki@mce.uec.ac.jp [Department of Mechanical Engineering and Intelligent Systems, University of Electro-Communications, 1-5-1, Chofugaoka, Chofu, Tokyo 182-8585 (Japan)

2012-06-01

We investigate the statistical mechanics of quasi-geostrophic point vortices of mixed sign (bi-disperse system) numerically and theoretically. Direct numerical simulations under periodic boundary conditions are performed using a fast special-purpose computer for molecular dynamics (GRAPE-DR). Clustering of point vortices of like sign is observed and two-dimensional (2D) equilibrium states are formed. It is shown that they are the solutions of the 2D mean-field equation, i.e. the sinh-Poisson equation. The sinh-Poisson equation is generalized to study the 3D nature of the equilibrium states, and a new mean-field equation with the 3D Laplace operator is derived based on the maximum entropy theory. 3D solutions are obtained at very low energy level. These solution branches, however, cannot be traced up to the higher energy level at which the direct numerical simulations are performed, and transitions to 2D solution branches take place when the energy is increased. (paper)

Synthesis of Numerical Methods for Modeling Wave Energy Converter-Point Absorbers: Preprint

Energy Technology Data Exchange (ETDEWEB)

Li, Y.; Yu, Y. H.

2012-05-01

During the past few decades, wave energy has received significant attention among all ocean energy formats. Industry has proposed hundreds of prototypes such as an oscillating water column, a point absorber, an overtopping system, and a bottom-hinged system. In particular, many researchers have focused on modeling the floating-point absorber as the technology to extract wave energy. Several modeling methods have been used such as the analytical method, the boundary-integral equation method, the Navier-Stokes equations method, and the empirical method. However, no standardized method has been decided. To assist the development of wave energy conversion technologies, this report reviews the methods for modeling the floating-point absorber.

Thinning spatial point processes into Poisson processes

DEFF Research Database (Denmark)

Møller, Jesper; Schoenberg, Frederic Paik

2010-01-01

are identified, and where we simulate backwards and forwards in order to obtain the thinned process. In the case of a Cox process, a simple independent thinning technique is proposed. In both cases, the thinning results in a Poisson process if and only if the true Papangelou conditional intensity is used, and......In this paper we describe methods for randomly thinning certain classes of spatial point processes. In the case of a Markov point process, the proposed method involves a dependent thinning of a spatial birth-and-death process, where clans of ancestors associated with the original points......, thus, can be used as a graphical exploratory tool for inspecting the goodness-of-fit of a spatial point process model. Several examples, including clustered and inhibitive point processes, are considered....

Two-point versus multiple-point geostatistics: the ability of geostatistical methods to capture complex geobodies and their facies associations—an application to a channelized carbonate reservoir, southwest Iran

International Nuclear Information System (INIS)

Hashemi, Seyyedhossein; Javaherian, Abdolrahim; Ataee-pour, Majid; Khoshdel, Hossein

2014-01-01

Facies models try to explain facies architectures which have a primary control on the subsurface heterogeneities and the fluid flow characteristics of a given reservoir. In the process of facies modeling, geostatistical methods are implemented to integrate different sources of data into a consistent model. The facies models should describe facies interactions; the shape and geometry of the geobodies as they occur in reality. Two distinct categories of geostatistical techniques are two-point and multiple-point (geo) statistics (MPS). In this study, both of the aforementioned categories were applied to generate facies models. A sequential indicator simulation (SIS) and a truncated Gaussian simulation (TGS) represented two-point geostatistical methods, and a single normal equation simulation (SNESIM) selected as an MPS simulation representative. The dataset from an extremely channelized carbonate reservoir located in southwest Iran was applied to these algorithms to analyze their performance in reproducing complex curvilinear geobodies. The SNESIM algorithm needs consistent training images (TI) in which all possible facies architectures that are present in the area are included. The TI model was founded on the data acquired from modern occurrences. These analogies delivered vital information about the possible channel geometries and facies classes that are typically present in those similar environments. The MPS results were conditioned to both soft and hard data. Soft facies probabilities were acquired from a neural network workflow. In this workflow, seismic-derived attributes were implemented as the input data. Furthermore, MPS realizations were conditioned to hard data to guarantee the exact positioning and continuity of the channel bodies. A geobody extraction workflow was implemented to extract the most certain parts of the channel bodies from the seismic data. These extracted parts of the channel bodies were applied to the simulation workflow as hard data

Calculus on Surfaces with General Closest Point Functions

KAUST Repository

Mä rz, Thomas; Macdonald, Colin B.

2012-01-01

The closest point method for solving partial differential equations (PDEs) posed on surfaces was recently introduced by Ruuth and Merriman [J. Comput. Phys., 227 (2008), pp. 1943- 1961] and successfully applied to a variety of surface PDEs. In this paper we study the theoretical foundations of this method. The main idea is that surface differentials of a surface function can be replaced with Cartesian differentials of its closest point extension, i.e., its composition with a closest point function. We introduce a general class of these closest point functions (a subset of differentiable retractions), show that these are exactly the functions necessary to satisfy the above idea, and give a geometric characterization of this class. Finally, we construct some closest point functions and demonstrate their effectiveness numerically on surface PDEs. © 2012 Society for Industrial and Applied Mathematics.

An approximate factorization procedure for solving nine-point elliptic difference equations. Application for a fast 2-D relativistic Fokker-Planck solver

Energy Technology Data Exchange (ETDEWEB)

Peysson, Y. [Association Euratom-CEA, CEA Grenoble, 38 (France). Dept. de Recherches sur la Fusion Controlee; Choucri, M. [Centre Canadien de Fusion Magnetique, Varennes, PQ (Canada)

1997-09-01

A full implicit numerical procedure based on the use of a nine-point difference operator is presented to solve the two dimensional (2{sub D}) relativistic Fokker-Planck equation for the current drive problem and synergetic effects between the lower hybrid and the electron cyclotron waves in tokamaks. As compared to the standard approach based on the use of a five-point difference operator [M. Shoucri, I. Shkarofsky, Comput. Phys. Comm. 82 (1994) 287], the convergence rate towards the steady state solution may be significantly enhanced with no loss of accuracy on the distribution function. Moreover, it is shown that the numerical stability may be strongly improved without a large degradation of the CPU time consumption as in the five-point scheme, making this approach very attractive for a fast solution of the 2-D Fokker-Planck equation on a fine grid in conjunction with other numerical codes for realistic plasma simulations. This new algorithm, based on an approximate matrix factorization technique, may be applied to all numerical problems with large sets of equations which involve nine-point difference operators. (author) 21 refs.

An approximate factorization procedure for solving nine-point elliptic difference equations. Application for a fast 2-D relativistic Fokker-Planck solver

International Nuclear Information System (INIS)

Peysson, Y.

1997-09-01

A full implicit numerical procedure based on the use of a nine-point difference operator is presented to solve the two dimensional (2 D ) relativistic Fokker-Planck equation for the current drive problem and synergetic effects between the lower hybrid and the electron cyclotron waves in tokamaks. As compared to the standard approach based on the use of a five-point difference operator [M. Shoucri, I. Shkarofsky, Comput. Phys. Comm. 82 (1994) 287], the convergence rate towards the steady state solution may be significantly enhanced with no loss of accuracy on the distribution function. Moreover, it is shown that the numerical stability may be strongly improved without a large degradation of the CPU time consumption as in the five-point scheme, making this approach very attractive for a fast solution of the 2-D Fokker-Planck equation on a fine grid in conjunction with other numerical codes for realistic plasma simulations. This new algorithm, based on an approximate matrix factorization technique, may be applied to all numerical problems with large sets of equations which involve nine-point difference operators. (author)

Effects of the second X-point on hot VDE in HL-2M

Science.gov (United States)

Xue, L.; Duan, X. R.; Zheng, G. Y.; Liu, Y. Q.; Dokuka, V. N.; Lukash, V. E.; Khayrutdinov, R. R.

2017-05-01

Study of the hot-plasma vertical displacement event (VDE) in advanced divertor configurations is of significant importance for ITER and for future fusion reactors. The newly designed, medium-sized copper-conductor machine HL-2M has the capability of generating the second X-point for various advanced divertor configurations. In this paper, effects of the second X-point on the hot VDE in HL-2M are numerically investigated by utilizing the non-linear time-dependent DINA code. The simulation results show that the existence of the second X-point at certain special locations appears to have a better stability in the vertical direction, compared to the standard configuration with the same main plasma parameters. Meanwhile, the peak halo current during the current quench tends to increase as the second X-point changes in the horizontal direction. The same quantity decreases as the second X-point changes in the vertical direction away from the dominant X-point. From the view point of minimizing the halo current, the tripod is better than the standard configuration, followed by the snowflake-plus and the exact snowflake (SF) configuration. The SF-minus is the worst scenario. On the other hand, the tripod configuration, as well as the SF minus configurations, results in relatively higher peak electromagnetic force acting on the vacuum vessel, when compared to other aforementioned configurations.

Analysis and research on Maximum Power Point Tracking of Photovoltaic Array with Fuzzy Logic Control and Three-point Weight Comparison Method

Institute of Scientific and Technical Information of China (English)

LIN; Kuang-Jang; LIN; Chii-Ruey

2010-01-01

The Photovoltaic Array has a best optimal operating point where the array operating can obtain the maximum power.However, the optimal operating point can be compromised by the strength of solar radiation,angle,and by the change of environment and load.Due to the constant changes in these conditions,it has become very difficult to locate the optimal operating point by following a mathematical model.Therefore,this study will focus mostly on the application of Fuzzy Logic Control theory and Three-point Weight Comparison Method in effort to locate the optimal operating point of solar panel and achieve maximum efficiency in power generation. The Three-point Weight Comparison Method is the comparison between the characteristic curves of the voltage of photovoltaic array and output power;it is a rather simple way to track the maximum power.The Fuzzy Logic Control,on the other hand,can be used to solve problems that cannot be effectively dealt with by calculation rules,such as concepts,contemplation, deductive reasoning,and identification.Therefore,this paper uses these two kinds of methods to make simulation successively. The simulation results show that,the Three-point Comparison Method is more effective under the environment with more frequent change of solar radiation;however,the Fuzzy Logic Control has better tacking efficiency under the environment with violent change of solar radiation.

Simulation of surface crack initiation induced by slip localization and point defects kinetics

International Nuclear Information System (INIS)

Sauzay, Maxime; Liu, Jia; Rachdi, Fatima

2014-01-01

Crack initiation along surface persistent slip bands (PSBs) has been widely observed and modelled. Nevertheless, from our knowledge, no physically-based fracture modelling has been proposed and validated with respect to the numerous recent experimental data showing the strong relationship between extrusion and microcrack initiation. The whole FE modelling accounts for: - localized plastic slip in PSBs; - production and annihilation of vacancies induced by cyclic slip. If temperature is high enough, point defects may diffuse in the surrounding matrix due to large concentration gradients, allowing continuous extrusion growth in agreement with Polak's model. At each cycle, the additional atoms diffusing from the matrix are taken into account by imposing an incremental free dilatation; - brittle fracture at the interfaces between PSBs and their surrounding matrix which is simulated using cohesive zone modelling. Any inverse fitting of parameter is avoided. Only experimental single crystal data are used such as hysteresis loops and resistivity values. Two fracture parameters are required: the {111} surface energy which depends on environment and the cleavage stress which is predicted by the universal binding energy relationship. The predicted extrusion growth curves agree rather well with the experimental data published for copper and the 316L steel. A linear dependence with respect to PSB length, thickness and slip plane angle is predicted in agreement with recent AFM measurement results. Crack initiation simulations predict fairly well the effects of PSB length and environment for copper single and poly-crystals. (authors)

Adding-point strategy for reduced-order hypersonic aerothermodynamics modeling based on fuzzy clustering

Science.gov (United States)

Chen, Xin; Liu, Li; Zhou, Sida; Yue, Zhenjiang

2016-09-01

Reduced order models(ROMs) based on the snapshots on the CFD high-fidelity simulations have been paid great attention recently due to their capability of capturing the features of the complex geometries and flow configurations. To improve the efficiency and precision of the ROMs, it is indispensable to add extra sampling points to the initial snapshots, since the number of sampling points to achieve an adequately accurate ROM is generally unknown in prior, but a large number of initial sampling points reduces the parsimony of the ROMs. A fuzzy-clustering-based adding-point strategy is proposed and the fuzzy clustering acts an indicator of the region in which the precision of ROMs is relatively low. The proposed method is applied to construct the ROMs for the benchmark mathematical examples and a numerical example of hypersonic aerothermodynamics prediction for a typical control surface. The proposed method can achieve a 34.5% improvement on the efficiency than the estimated mean squared error prediction algorithm and shows same-level prediction accuracy.

Sequential function approximation on arbitrarily distributed point sets

Science.gov (United States)

Wu, Kailiang; Xiu, Dongbin

2018-02-01

We present a randomized iterative method for approximating unknown function sequentially on arbitrary point set. The method is based on a recently developed sequential approximation (SA) method, which approximates a target function using one data point at each step and avoids matrix operations. The focus of this paper is on data sets with highly irregular distribution of the points. We present a nearest neighbor replacement (NNR) algorithm, which allows one to sample the irregular data sets in a near optimal manner. We provide mathematical justification and error estimates for the NNR algorithm. Extensive numerical examples are also presented to demonstrate that the NNR algorithm can deliver satisfactory convergence for the SA method on data sets with high irregularity in their point distributions.

Statistical aspects of determinantal point processes

DEFF Research Database (Denmark)

Lavancier, Frédéric; Møller, Jesper; Rubak, Ege

The statistical aspects of determinantal point processes (DPPs) seem largely unexplored. We review the appealing properties of DDPs, demonstrate that they are useful models for repulsiveness, detail a simulation procedure, and provide freely available software for simulation and statistical infer...

Anderson Acceleration for Fixed-Point Iterations

Energy Technology Data Exchange (ETDEWEB)

Walker, Homer F. [Worcester Polytechnic Institute, MA (United States)

2015-08-31

The purpose of this grant was to support research on acceleration methods for fixed-point iterations, with applications to computational frameworks and simulation problems that are of interest to DOE.

Address Points - COUNTY_ADDRESS_POINTS_IDHS_IN: Address Points Maintained by County Agencies in Indiana (Indiana Department of Homeland Security, Point feature class)

Data.gov (United States)

NSGIC State | GIS Inventory — COUNTY_ADDRESS_POINTS_IDHS_IN is an ESRI Geodatabase point feature class that contains address points maintained by county agencies in Indiana, provided by personnel...

Analysis of Stress Updates in the Material-point Method

DEFF Research Database (Denmark)

Andersen, Søren; Andersen, Lars

2009-01-01

The material-point method (MPM) is a new numerical method for analysis of large strain engineering problems. The MPM applies a dual formulation, where the state of the problem (mass, stress, strain, velocity etc.) is tracked using a finite set of material points while the governing equations...... are solved on a background computational grid. Several references state, that one of the main advantages of the material-point method is the easy application of complicated material behaviour as the constitutive response is updated individually for each material point. However, as discussed here, the MPM way...

QCD and the chiral critical point

International Nuclear Information System (INIS)

Gavin, S.; Gocksch, A.; Pisarski, R.D.

1994-01-01

As an extension of QCD, consider a theory with ''2+1'' flavors, where the current quark masses are held in a fixed ratio as the overall scale of the quark masses is varied. At nonzero temperature and baryon density it is expected that in the chiral limit the chiral phase transition is of first order. Increasing the quark mass from zero, the chiral transition becomes more weakly first order, and can end in a chiral critical point. We show that the only massless field at the chiral critical point is a σ meson, with the universality class that of the Ising model. Present day lattice simulations indicate that QCD is (relatively) near to the chiral critical point

On the Required Number of Antennas in a Point-to-Point Large-but-Finite MIMO System: Outage-Limited Scenario

KAUST Repository

Makki, Behrooz

2016-03-22

This paper investigates the performance of the point-To-point multiple-input-multiple-output (MIMO) systems in the presence of a large but finite numbers of antennas at the transmitters and/or receivers. Considering the cases with and without hybrid automatic repeat request (HARQ) feedback, we determine the minimum numbers of the transmit/receive antennas, which are required to satisfy different outage probability constraints. Our results are obtained for different fading conditions and the effect of the power amplifiers efficiency/feedback error probability on the performance of the MIMO-HARQ systems is analyzed. Then, we use some recent results on the achievable rates of finite block-length codes, to analyze the effect of the codewords lengths on the system performance. Moreover, we derive closed-form expressions for the asymptotic performance of the MIMO-HARQ systems when the number of antennas increases. Our analytical and numerical results show that different outage requirements can be satisfied with relatively few transmit/receive antennas. © 1972-2012 IEEE.

Recursive representation of the torus 1-point conformal block

Science.gov (United States)

Hadasz, Leszek; Jaskólski, Zbigniew; Suchanek, Paulina

2010-01-01

The recursive relation for the 1-point conformal block on a torus is derived and used to prove the identities between conformal blocks recently conjectured by Poghossian in [1]. As an illustration of the efficiency of the recurrence method the modular invariance of the 1-point Liouville correlation function is numerically analyzed.

Hardware-accelerated Point Generation and Rendering of Point-based Impostors

DEFF Research Database (Denmark)

Bærentzen, Jakob Andreas

2005-01-01

This paper presents a novel scheme for generating points from triangle models. The method is fast and lends itself well to implementation using graphics hardware. The triangle to point conversion is done by rendering the models, and the rendering may be performed procedurally or by a black box API....... I describe the technique in detail and discuss how the generated point sets can easily be used as impostors for the original triangle models used to create the points. Since the points reside solely in GPU memory, these impostors are fairly efficient. Source code is available online....

Active point out-of-plane ultrasound calibration

Science.gov (United States)

Cheng, Alexis; Guo, Xiaoyu; Zhang, Haichong K.; Kang, Hyunjae; Etienne-Cummings, Ralph; Boctor, Emad M.

2015-03-01

Image-guided surgery systems are often used to provide surgeons with informational support. Due to several unique advantages such as ease of use, real-time image acquisition, and no ionizing radiation, ultrasound is a common intraoperative medical imaging modality used in image-guided surgery systems. To perform advanced forms of guidance with ultrasound, such as virtual image overlays or automated robotic actuation, an ultrasound calibration process must be performed. This process recovers the rigid body transformation between a tracked marker attached to the transducer and the ultrasound image. Point-based phantoms are considered to be accurate, but their calibration framework assumes that the point is in the image plane. In this work, we present the use of an active point phantom and a calibration framework that accounts for the elevational uncertainty of the point. Given the lateral and axial position of the point in the ultrasound image, we approximate a circle in the axial-elevational plane with a radius equal to the axial position. The standard approach transforms all of the imaged points to be a single physical point. In our approach, we minimize the distances between the circular subsets of each image, with them ideally intersecting at a single point. We simulated in noiseless and noisy cases, presenting results on out-of-plane estimation errors, calibration estimation errors, and point reconstruction precision. We also performed an experiment using a robot arm as the tracker, resulting in a point reconstruction precision of 0.64mm.

December 2002 Lidar Point Data of Southern California Coastline: Dana Point to Point La Jolla

Data.gov (United States)

National Oceanic and Atmospheric Administration, Department of Commerce — This data set contains lidar point data (latitude and longitude) from a strip of Southern California coastline (including water, beach, cliffs, and top of cliffs)...

April 2004 Lidar Point Data of Southern California Coastline: Dana Point to Point La Jolla

Data.gov (United States)

National Oceanic and Atmospheric Administration, Department of Commerce — This data set contains lidar point data (Geodetic Coordinates) from a strip of Southern California coastline (including water, beach, cliffs, and top of cliffs) from...

March 2003 Lidar Point Data of Southern California Coastline: Dana Point to Point La Jolla

Data.gov (United States)

National Oceanic and Atmospheric Administration, Department of Commerce — This data set contains lidar point data (Geodetic Coordinates) from a strip of Southern California coastline (including water, beach, cliffs, and top of cliffs) from...

Four points function fitted and first derivative procedure for determining the end points in potentiometric titration curves: statistical analysis and method comparison.

Science.gov (United States)

Kholeif, S A

2001-06-01

A new method that belongs to the differential category for determining the end points from potentiometric titration curves is presented. It uses a preprocess to find first derivative values by fitting four data points in and around the region of inflection to a non-linear function, and then locate the end point, usually as a maximum or minimum, using an inverse parabolic interpolation procedure that has an analytical solution. The behavior and accuracy of the sigmoid and cumulative non-linear functions used are investigated against three factors. A statistical evaluation of the new method using linear least-squares method validation and multifactor data analysis are covered. The new method is generally applied to symmetrical and unsymmetrical potentiometric titration curves, and the end point is calculated using numerical procedures only. It outperforms the "parent" regular differential method in almost all factors levels and gives accurate results comparable to the true or estimated true end points. Calculated end points from selected experimental titration curves compatible with the equivalence point category of methods, such as Gran or Fortuin, are also compared with the new method.

Acceleration of saddle-point searches with machine learning

Energy Technology Data Exchange (ETDEWEB)

Peterson, Andrew A., E-mail: andrew-peterson@brown.edu [School of Engineering, Brown University, Providence, Rhode Island 02912 (United States)

2016-08-21

In atomistic simulations, the location of the saddle point on the potential-energy surface (PES) gives important information on transitions between local minima, for example, via transition-state theory. However, the search for saddle points often involves hundreds or thousands of ab initio force calls, which are typically all done at full accuracy. This results in the vast majority of the computational effort being spent calculating the electronic structure of states not important to the researcher, and very little time performing the calculation of the saddle point state itself. In this work, we describe how machine learning (ML) can reduce the number of intermediate ab initio calculations needed to locate saddle points. Since machine-learning models can learn from, and thus mimic, atomistic simulations, the saddle-point search can be conducted rapidly in the machine-learning representation. The saddle-point prediction can then be verified by an ab initio calculation; if it is incorrect, this strategically has identified regions of the PES where the machine-learning representation has insufficient training data. When these training data are used to improve the machine-learning model, the estimates greatly improve. This approach can be systematized, and in two simple example problems we demonstrate a dramatic reduction in the number of ab initio force calls. We expect that this approach and future refinements will greatly accelerate searches for saddle points, as well as other searches on the potential energy surface, as machine-learning methods see greater adoption by the atomistics community.

Acceleration of saddle-point searches with machine learning

International Nuclear Information System (INIS)

Peterson, Andrew A.

2016-01-01

In atomistic simulations, the location of the saddle point on the potential-energy surface (PES) gives important information on transitions between local minima, for example, via transition-state theory. However, the search for saddle points often involves hundreds or thousands of ab initio force calls, which are typically all done at full accuracy. This results in the vast majority of the computational effort being spent calculating the electronic structure of states not important to the researcher, and very little time performing the calculation of the saddle point state itself. In this work, we describe how machine learning (ML) can reduce the number of intermediate ab initio calculations needed to locate saddle points. Since machine-learning models can learn from, and thus mimic, atomistic simulations, the saddle-point search can be conducted rapidly in the machine-learning representation. The saddle-point prediction can then be verified by an ab initio calculation; if it is incorrect, this strategically has identified regions of the PES where the machine-learning representation has insufficient training data. When these training data are used to improve the machine-learning model, the estimates greatly improve. This approach can be systematized, and in two simple example problems we demonstrate a dramatic reduction in the number of ab initio force calls. We expect that this approach and future refinements will greatly accelerate searches for saddle points, as well as other searches on the potential energy surface, as machine-learning methods see greater adoption by the atomistics community.

Acceleration of saddle-point searches with machine learning.

Science.gov (United States)

Peterson, Andrew A

2016-08-21

In atomistic simulations, the location of the saddle point on the potential-energy surface (PES) gives important information on transitions between local minima, for example, via transition-state theory. However, the search for saddle points often involves hundreds or thousands of ab initio force calls, which are typically all done at full accuracy. This results in the vast majority of the computational effort being spent calculating the electronic structure of states not important to the researcher, and very little time performing the calculation of the saddle point state itself. In this work, we describe how machine learning (ML) can reduce the number of intermediate ab initio calculations needed to locate saddle points. Since machine-learning models can learn from, and thus mimic, atomistic simulations, the saddle-point search can be conducted rapidly in the machine-learning representation. The saddle-point prediction can then be verified by an ab initio calculation; if it is incorrect, this strategically has identified regions of the PES where the machine-learning representation has insufficient training data. When these training data are used to improve the machine-learning model, the estimates greatly improve. This approach can be systematized, and in two simple example problems we demonstrate a dramatic reduction in the number of ab initio force calls. We expect that this approach and future refinements will greatly accelerate searches for saddle points, as well as other searches on the potential energy surface, as machine-learning methods see greater adoption by the atomistics community.

Numerical and experimental analysis of a point focus solar collector using high concentration imaging PMMA Fresnel lens

International Nuclear Information System (INIS)

Xie, W.T.; Dai, Y.J.; Wang, R.Z.

2011-01-01

Research highlights: → We studied a point focus Fresnel solar collector using different cavity receivers. → The collector heat removal factors are derived to find the optimal cavity shape. → Numerical and experimental analysis shows that the conical cavity is optimum. -- Abstract: A high concentration imaging Fresnel solar collector provided with different cavity receivers was developed and its behavior was investigated. Round copper pipes winded into different spring shapes were used as receiver by placing in the cylindrical cavity to absorb concentrated solar energy and transfer it to a heat transfer fluid (HTF). The collector efficiency factor and collector heat removal factor were derived for the cavity receivers to find out heat transfer mechanism and to propose an effective way for evaluating the performance of Fresnel solar collector and determining the optimal cavity structure. The problem of Fresnel solar collector with synthetic heat transfer oil flow was simulated and analyzed to investigate heat loss from different cavity receivers. Solar irradiation as well as convection and heat transfer in the circulating fluid and between the internal surfaces of the cavity and the environment are considered in the model. The temperature distribution over its area as well as the collector thermal efficiency at nominal flow rate was used in order to validate the simulation results. It was found that the simulated temperature distribution during operation and the average collector efficiency are in good agreement with the experimental data. Finally, the optimal shape of solar cavity receiver, as well as its thermal performance, are deeply analyzed and discussed.

Fractional neutron point kinetics equations for nuclear reactor dynamics

International Nuclear Information System (INIS)

Espinosa-Paredes, Gilberto; Polo-Labarrios, Marco-A.; Espinosa-Martinez, Erick-G.; Valle-Gallegos, Edmundo del

2011-01-01

The fractional point-neutron kinetics model for the dynamic behavior in a nuclear reactor is derived and analyzed in this paper. The fractional model retains the main dynamic characteristics of the neutron motion in which the relaxation time associated with a rapid variation in the neutron flux contains a fractional order, acting as exponent of the relaxation time, to obtain the best representation of a nuclear reactor dynamics. The physical interpretation of the fractional order is related with non-Fickian effects from the neutron diffusion equation point of view. The numerical approximation to the solution of the fractional neutron point kinetics model, which can be represented as a multi-term high-order linear fractional differential equation, is calculated by reducing the problem to a system of ordinary and fractional differential equations. The numerical stability of the fractional scheme is investigated in this work. Results for neutron dynamic behavior for both positive and negative reactivity and for different values of fractional order are shown and compared with the classic neutron point kinetic equations. Additionally, a related review with the neutron point kinetics equations is presented, which encompasses papers written in English about this research topic (as well as some books and technical reports) published since 1940 up to 2010.

Graphene Dirac point tuned by ferroelectric polarization field

Science.gov (United States)

Wang, Xudong; Chen, Yan; Wu, Guangjian; Wang, Jianlu; Tian, Bobo; Sun, Shuo; Shen, Hong; Lin, Tie; Hu, Weida; Kang, Tingting; Tang, Minghua; Xiao, Yongguang; Sun, Jinglan; Meng, Xiangjian; Chu, Junhao

2018-04-01

Graphene has received numerous attention for future nanoelectronics and optoelectronics. The Dirac point is a key parameter of graphene that provides information about its carrier properties. There are lots of methods to tune the Dirac point of graphene, such as chemical doping, impurities, defects, and disorder. In this study, we report a different approach to tune the Dirac point of graphene using a ferroelectric polarization field. The Dirac point can be adjusted to near the ferroelectric coercive voltage regardless its original position. We have ensured this phenomenon by temperature-dependent experiments, and analyzed its mechanism with the theory of impurity correlation in graphene. Additionally, with the modulation of ferroelectric polymer, the current on/off ratio and mobility of graphene transistor both have been improved. This work provides an effective method to tune the Dirac point of graphene, which can be readily used to configure functional devices such as p-n junctions and inverters.

A model for the two-point velocity correlation function in turbulent channel flow

International Nuclear Information System (INIS)

Sahay, A.; Sreenivasan, K.R.

1996-01-01

A relatively simple analytical expression is presented to approximate the equal-time, two-point, double-velocity correlation function in turbulent channel flow. To assess the accuracy of the model, we perform the spectral decomposition of the integral operator having the model correlation function as its kernel. Comparisons of the empirical eigenvalues and eigenfunctions with those constructed from direct numerical simulations data show good agreement. copyright 1996 American Institute of Physics

Simulation Investigation on Combustion Characteristics in a Four-Point Lean Direct Injection Combustor with Hydrogen/Air

Directory of Open Access Journals (Sweden)

Jianzhong Li

2017-06-01

Full Text Available To investigate the combustion characteristics in multi-point lean direct injection (LDI combustors with hydrogen/air, two swirl–venturi 2 × 2 array four-point LDI combustors were designed. The four-point LDI combustor consists of injector assembly, swirl–venturi array and combustion chamber. The injector, swirler and venturi together govern the rapid mixing of hydrogen and air to form the mixture for combustion. Using clockwise swirlers and anticlockwise swirlers, the co-swirling and count-swirling swirler arrays LDI combustors were achieved. Using Reynolds-Averaged Navier–Stokes (RANS code for steady-state reacting flow computations, the four-point LDI combustors with hydrogen/air were simulated with an 11 species and 23 lumped reaction steps H2/Air reaction mechanism. The axial velocity, turbulence kinetic energy, total pressure drop coefficient, outlet temperature, mass fraction of OH and emission of pollutant NO of four-point LDI combustors, with different equivalence ratios, are here presented and discussed. As the equivalence ratios increased, the total pressure drop coefficient became higher because of increasing heat loss. Increasing equivalence ratios also corresponded with the rise in outlet temperature of the four-point LDI combustors, as well as an increase in the emission index of NO EINO in the four-point LDI combustors. Along the axial distance, the EINO always increased and was at maximum at the exit of the dump. Along the chamber, the EINO gradually increased, maximizing at the exit of chamber. The total temperature of four-point LDI combustors with different equivalence ratios was identical to the theoretical equilibrium temperature. The EINO was an exponential function of the equivalence ratio.

Focusing: coming to the point in metamaterials

Science.gov (United States)

Guenneau, S.; Diatta, A.; McPhedran, R. C.

2010-04-01

This paper reviews some properties of lenses in curved and folded optical spaces. The point of the paper is to show some limitations of geometrical optics in the analysis of subwavelength focusing. We first provide a comprehensive derivation for the equation of geodesics in curved optical spaces, which is a tool of choice to design metamaterials in transformation optics. We then analyse the resolution of the image of a line source radiating in the Maxwell fisheye and the Veselago-Pendry slab lens. The former optical medium is deduced from the stereographic projection of a virtual sphere and displays a heterogeneous refractive index n(r) which is proportional to the inverse of 1 + r 2. The latter is described by a homogeneous, but negative, refractive index. It has been suggested that the fisheye makes a perfect lens without negative refraction [Leonhardt, Philbin arxiv:0805.4778v2]. However, we point out that the definition of super-resolution in such a heterogeneous medium should be computed with respect to the wavelength in a homogenised medium, and it is perhaps more adequate to talk about a conjugate image rather than a perfect image (the former does not necessarily contain the evanescent components of the source). We numerically find that both the Maxwell fisheye and a thick silver slab lens lead to a resolution close to λ/3 in transverse magnetic polarisation (electric field pointing orthogonal to the plane). We note a shift of the image plane in the latter lens. We also observe that two sources lead to multiple secondary images in the former lens, as confirmed from light rays travelling along geodesics of the virtual sphere. We further observe resolutions ranging from λ/2 to nearly λ/4 for magnetic dipoles of varying orientations of dipole moments within the fisheye in transverse electric polarisation (magnetic field pointing orthogonal to the plane). Finally, we analyse the Eaton lens for which the source and its image are either located within a unit

A Bayesian MCMC method for point process models with intractable normalising constants

DEFF Research Database (Denmark)

Berthelsen, Kasper Klitgaard; Møller, Jesper

2004-01-01

to simulate from the "unknown distribution", perfect simulation algorithms become useful. We illustrate the method in cases whre the likelihood is given by a Markov point process model. Particularly, we consider semi-parametric Bayesian inference in connection to both inhomogeneous Markov point process models...... and pairwise interaction point processes....

Point Information Gain and Multidimensional Data Analysis

Directory of Open Access Journals (Sweden)

Renata Rychtáriková

2016-10-01

Full Text Available We generalize the point information gain (PIG and derived quantities, i.e., point information gain entropy (PIE and point information gain entropy density (PIED, for the case of the Rényi entropy and simulate the behavior of PIG for typical distributions. We also use these methods for the analysis of multidimensional datasets. We demonstrate the main properties of PIE/PIED spectra for the real data with the examples of several images and discuss further possible utilizations in other fields of data processing.

Tipping Point

Medline Plus

Full Text Available ... en español Blog About OnSafety CPSC Stands for Safety The Tipping Point Home > 60 Seconds of Safety (Videos) > The Tipping Point The Tipping Point by ... danger death electrical fall furniture head injury product safety television tipover tv Watch the video in Adobe ...

Point charges optimally placed to represent the multipole expansion of charge distributions.

Directory of Open Access Journals (Sweden)

Ramu Anandakrishnan

Full Text Available We propose an approach for approximating electrostatic charge distributions with a small number of point charges to optimally represent the original charge distribution. By construction, the proposed optimal point charge approximation (OPCA retains many of the useful properties of point multipole expansion, including the same far-field asymptotic behavior of the approximate potential. A general framework for numerically computing OPCA, for any given number of approximating charges, is described. We then derive a 2-charge practical point charge approximation, PPCA, which approximates the 2-charge OPCA via closed form analytical expressions, and test the PPCA on a set of charge distributions relevant to biomolecular modeling. We measure the accuracy of the new approximations as the RMS error in the electrostatic potential relative to that produced by the original charge distribution, at a distance 2x the extent of the charge distribution--the mid-field. The error for the 2-charge PPCA is found to be on average 23% smaller than that of optimally placed point dipole approximation, and comparable to that of the point quadrupole approximation. The standard deviation in RMS error for the 2-charge PPCA is 53% lower than that of the optimal point dipole approximation, and comparable to that of the point quadrupole approximation. We also calculate the 3-charge OPCA for representing the gas phase quantum mechanical charge distribution of a water molecule. The electrostatic potential calculated by the 3-charge OPCA for water, in the mid-field (2.8 Å from the oxygen atom, is on average 33.3% more accurate than the potential due to the point multipole expansion up to the octupole order. Compared to a 3 point charge approximation in which the charges are placed on the atom centers, the 3-charge OPCA is seven times more accurate, by RMS error. The maximum error at the oxygen-Na distance (2.23 Å is half that of the point multipole expansion up to the octupole

Quench dynamics across quantum critical points

International Nuclear Information System (INIS)

Sengupta, K.; Powell, Stephen; Sachdev, Subir

2004-01-01

We study the quantum dynamics of a number of model systems as their coupling constants are changed rapidly across a quantum critical point. The primary motivation is provided by the recent experiments of Greiner et al. [Nature (London) 415, 39 (2002)] who studied the response of a Mott insulator of ultracold atoms in an optical lattice to a strong potential gradient. In a previous work, it had been argued that the resonant response observed at a critical potential gradient could be understood by proximity to an Ising quantum critical point describing the onset of density wave order. Here we obtain numerical results on the evolution of the density wave order as the potential gradient is scanned across the quantum critical point. This is supplemented by studies of the integrable quantum Ising spin chain in a transverse field, where we obtain exact results for the evolution of the Ising order correlations under a time-dependent transverse field. We also study the evolution of transverse superfluid order in the three-dimensional case. In all cases, the order parameter is best enhanced in the vicinity of the quantum critical point

Critical point of Nf=3 QCD from lattice simulations in the canonical ensemble

International Nuclear Information System (INIS)

Li Anyi; Alexandru, Andrei; Liu, Keh-Fei

2011-01-01

A canonical ensemble algorithm is employed to study the phase diagram of N f =3 QCD using lattice simulations. We lock in the desired quark number sector using an exact Fourier transform of the fermion determinant. We scan the phase space below T c and look for an S-shape structure in the chemical potential, which signals the coexistence phase of a first order phase transition in finite volume. Applying Maxwell construction, we determine the boundaries of the coexistence phase at three temperatures and extrapolate them to locate the critical point. Using an improved gauge action and improved Wilson fermions on lattices with a spatial extent of 1.8 fm and quark masses close to that of the strange, we find the critical point at T E =0.925(5)T c and baryon chemical potential μ B E =2.60(8)T c .

A simple model for correcting the zero point energy problem in classical trajectory simulations of polyatomic molecules

International Nuclear Information System (INIS)

Miller, W.H.; Hase, W.L.; Darling, C.L.

1989-01-01

A simple model is proposed for correcting problems with zero point energy in classical trajectory simulations of dynamical processes in polyatomic molecules. The ''problems'' referred to are that classical mechanics allows the vibrational energy in a mode to decrease below its quantum zero point value, and since the total energy is conserved classically this can allow too much energy to pool in other modes. The proposed model introduces hard sphere-like terms in action--angle variables that prevent the vibrational energy in any mode from falling below its zero point value. The algorithm which results is quite simple in terms of the cartesian normal modes of the system: if the energy in a mode k, say, decreases below its zero point value at time t, then at this time the momentum P k for that mode has its sign changed, and the trajectory continues. This is essentially a time reversal for mode k (only exclamation point), and it conserves the total energy of the system. One can think of the model as supplying impulsive ''quantum kicks'' to a mode whose energy attempts to fall below its zero point value, a kind of ''Planck demon'' analogous to a Brownian-like random force. The model is illustrated by application to a model of CH overtone relaxation

Interval Mathematics Applied to Critical Point Transitions

Directory of Open Access Journals (Sweden)

Benito A. Stradi

2012-03-01

Full Text Available The determination of critical points of mixtures is important for both practical and theoretical reasons in the modeling of phase behavior, especially at high pressure. The equations that describe the behavior of complex mixtures near critical points are highly nonlinear and with multiplicity of solutions to the critical point equations. Interval arithmetic can be used to reliably locate all the critical points of a given mixture. The method also verifies the nonexistence of a critical point if a mixture of a given composition does not have one. This study uses an interval Newton/Generalized Bisection algorithm that provides a mathematical and computational guarantee that all mixture critical points are located. The technique is illustrated using several example problems. These problems involve cubic equation of state models; however, the technique is general purpose and can be applied in connection with other nonlinear problems.

Performance of Closely Spaced Point Absorbers with Constrained Floater Motion

DEFF Research Database (Denmark)

De Backer, G.; Vantorre, M.; Beels, C.

2009-01-01

The performance of a wave energy converter array of twelve heaving point absorbers has been assessed numerically in a frequency domain model. Each point absorber is assumed to have its own linear power take-off. The impact of slamming, stroke and force restrictions on the power absorption...

Application of the nudged elastic band method to the point-to-point radio wave ray tracing in IRI modeled ionosphere

Science.gov (United States)

Nosikov, I. A.; Klimenko, M. V.; Bessarab, P. F.; Zhbankov, G. A.

2017-07-01

Point-to-point ray tracing is an important problem in many fields of science. While direct variational methods where some trajectory is transformed to an optimal one are routinely used in calculations of pathways of seismic waves, chemical reactions, diffusion processes, etc., this approach is not widely known in ionospheric point-to-point ray tracing. We apply the Nudged Elastic Band (NEB) method to a radio wave propagation problem. In the NEB method, a chain of points which gives a discrete representation of the radio wave ray is adjusted iteratively to an optimal configuration satisfying the Fermat's principle, while the endpoints of the trajectory are kept fixed according to the boundary conditions. Transverse displacements define the radio ray trajectory, while springs between the points control their distribution along the ray. The method is applied to a study of point-to-point ionospheric ray tracing, where the propagation medium is obtained with the International Reference Ionosphere model taking into account traveling ionospheric disturbances. A 2-dimensional representation of the optical path functional is developed and used to gain insight into the fundamental difference between high and low rays. We conclude that high and low rays are minima and saddle points of the optical path functional, respectively.

Close pairs of relative equilibria for identical point vortices

DEFF Research Database (Denmark)

Dirksen, Tobias; Aref, Hassan

2011-01-01

Numerical solution of the classical problem of relative equilibria for identical point vortices on the unbounded plane reveals configurations that are very close to the analytically known, centered, symmetrically arranged, nested equilateral triangles. New numerical solutions of this kind are fou...

Developing a Tool Point Control Scheme for a Hydraulic Crane Using Interactive Real-time Dynamic Simulation

DEFF Research Database (Denmark)

Pedersen, Mikkel Melters; Hansen, Michael Rygaard; Ballebye, Morten

2010-01-01

This paper describes the implementation of an interactive real-time dynamic simulation model of a hydraulic crane. The user input to the model is given continuously via joystick and output is presented continuously in a 3D animation. Using this simulation model, a tool point control scheme...... is developed for the specific crane, considering the saturation phenomena of the system and practical implementation....

Control of the electrical operating point of an electric discharge by mean of its power supply; Controle du point de fonctionnement des decharges electriques par l'intermediaire de leur alimentation

Energy Technology Data Exchange (ETDEWEB)

Salanne, J.Ph.

2005-11-15

The operating points obtained by the coupling of the power supply with the electric discharge system can be unstable because of the dynamical behaviour of the discharge or because of a change in its length. In this work, the different possible couplings existing between the discharge and the characteristics of its power supply are analyzed in order to optimize the design and control of the power supply and to control the operating point. Analytical and numerical modeling of the system are proposed which allow to simulate the couplings between the power supply and the discharge. This approach is completed by experimental investigations allowing to consider the cases of peak/peak discharges, dielectric barrier discharges (DBD), and gliding arcs. (J.S.)

Track length estimation applied to point detectors

International Nuclear Information System (INIS)

Rief, H.; Dubi, A.; Elperin, T.

1984-01-01

The concept of the track length estimator is applied to the uncollided point flux estimator (UCF) leading to a new algorithm of calculating fluxes at a point. It consists essentially of a line integral of the UCF, and although its variance is unbounded, the convergence rate is that of a bounded variance estimator. In certain applications, involving detector points in the vicinity of collimated beam sources, it has a lower variance than the once-more-collided point flux estimator, and its application is more straightforward

OFF-Stagnation point testing in plasma facility

Science.gov (United States)

Viladegut, A.; Chazot, O.

2015-06-01

Reentry space vehicles face extreme conditions of heat flux when interacting with the atmosphere at hypersonic velocities. Stagnation point heat flux is normally used as a reference for Thermal Protection Material (TPS) design; however, many critical phenomena also occur at off-stagnation point. This paper adresses the implementation of an offstagnation point methodology able to duplicate in ground facility the hypersonic boundary layer over a flat plate model. The first analysis using two-dimensional (2D) computational fluid dynamics (CFD) simulations is carried out to understand the limitations of this methodology when applying it in plasma wind tunnel. The results from the testing campaign at VKI Plasmatron are also presented.

Recent advances in lossy compression of scientific floating-point data

Science.gov (United States)

Lindstrom, P.

2017-12-01

With a continuing exponential trend in supercomputer performance, ever larger data sets are being generated through numerical simulation. Bandwidth and storage capacity are, however, not keeping pace with this increase in data size, causing significant data movement bottlenecks in simulation codes and substantial monetary costs associated with archiving vast volumes of data. Worse yet, ever smaller fractions of data generated can be stored for further analysis, where scientists frequently rely on decimating or averaging large data sets in time and/or space. One way to mitigate these problems is to employ data compression to reduce data volumes. However, lossless compression of floating-point data can achieve only very modest size reductions on the order of 10-50%. We present ZFP and FPZIP, two state-of-the-art lossy compressors for structured floating-point data that routinely achieve one to two orders of magnitude reduction with little to no impact on the accuracy of visualization and quantitative data analysis. We provide examples of the use of such lossy compressors in climate and seismic modeling applications to effectively accelerate I/O and reduce storage requirements. We further discuss how the design decisions behind these and other compressors impact error distributions and other statistical and differential properties, including derived quantities of interest relevant to each science application.

Point to point multispectral light projection applied to cultural heritage

Science.gov (United States)

Vázquez, D.; Alvarez, A.; Canabal, H.; Garcia, A.; Mayorga, S.; Muro, C.; Galan, T.

2017-09-01

Use of new of light sources based on LED technology should allow the develop of systems that combine conservation and exhibition requirements and allow to make these art goods available to the next generations according to sustainability principles. The goal of this work is to develop light systems and sources with an optimized spectral distribution for each specific point of the art piece. This optimization process implies to maximize the color fidelity reproduction and the same time to minimize the photochemical damage. Perceived color under these sources will be similar (metameric) to technical requirements given by the restoration team uncharged of the conservation and exhibition of the goods of art. Depending of the fragility of the exposed art objects (i.e. spectral responsivity of the material) the irradiance must be kept under a critical level. Therefore, it is necessary to develop a mathematical model that simulates with enough accuracy both the visual effect of the illumination and the photochemical impact of the radiation. Spectral reflectance of a reference painting The mathematical model is based on a merit function that optimized the individual intensity of the LED-light sources taking into account the damage function of the material and color space coordinates. Moreover the algorithm used weights for damage and color fidelity in order to adapt the model to a specific museal application. In this work we show a sample of this technology applied to a picture of Sorolla (1863-1923) an important Spanish painter title "woman walking at the beach".

Determination of point of maximum likelihood in failure domain using genetic algorithms

International Nuclear Information System (INIS)

Obadage, A.S.; Harnpornchai, N.

2006-01-01

The point of maximum likelihood in a failure domain yields the highest value of the probability density function in the failure domain. The maximum-likelihood point thus represents the worst combination of random variables that contribute in the failure event. In this work Genetic Algorithms (GAs) with an adaptive penalty scheme have been proposed as a tool for the determination of the maximum likelihood point. The utilization of only numerical values in the GAs operation makes the algorithms applicable to cases of non-linear and implicit single and multiple limit state function(s). The algorithmic simplicity readily extends its application to higher dimensional problems. When combined with Monte Carlo Simulation, the proposed methodology will reduce the computational complexity and at the same time will enhance the possibility in rare-event analysis under limited computational resources. Since, there is no approximation done in the procedure, the solution obtained is considered accurate. Consequently, GAs can be used as a tool for increasing the computational efficiency in the element and system reliability analyses

Modeling of Landslides with the Material Point Method

DEFF Research Database (Denmark)

Andersen, Søren Mikkel; Andersen, Lars

2008-01-01

A numerical model for studying the dynamic evolution of landslides is presented. The numerical model is based on the Generalized Interpolation Material Point Method. A simplified slope with a house placed on top is analysed. An elasto-plastic material model based on the Mohr-Coulomb yield criterion...

Modelling of Landslides with the Material-point Method

DEFF Research Database (Denmark)

Andersen, Søren; Andersen, Lars

2009-01-01

A numerical model for studying the dynamic evolution of landslides is presented. The numerical model is based on the Generalized Interpolation Material Point Method. A simplified slope with a house placed on top is analysed. An elasto-plastic material model based on the Mohr-Coulomb yield criterion...

When should we recommend use of dual time-point and delayed time-point imaging techniques in FDG PET?

Energy Technology Data Exchange (ETDEWEB)

Cheng, Gang [Philadelphia VA Medical Center, Department of Radiology, Philadelphia, PA (United States); Hospital of the University of Pennsylvania, Department of Radiology, Philadelphia, PA (United States); Torigian, Drew A.; Alavi, Abass [Hospital of the University of Pennsylvania, Department of Radiology, Philadelphia, PA (United States); Zhuang, Hongming [Children' s Hospital of Philadelphia, Department of Radiology, Philadelphia, PA (United States)

2013-05-15

FDG PET and PET/CT are now widely used in oncological imaging for tumor characterization, staging, restaging, and response evaluation. However, numerous benign etiologies may cause increased FDG uptake indistinguishable from that of malignancy. Multiple studies have shown that dual time-point imaging (DTPI) of FDG PET may be helpful in differentiating malignancy from benign processes. However, exceptions exist, and some studies have demonstrated significant overlap of FDG uptake patterns between benign and malignant lesions on delayed time-point images. In this review, we summarize our experience and opinions on the value of DTPI and delayed time-point imaging in oncology, with a review of the relevant literature. We believe that the major value of DTPI and delayed time-point imaging is the increased sensitivity due to continued clearance of background activity and continued FDG accumulation in malignant lesions, if the same diagnostic criteria (as in the initial standard single time-point imaging) are used. The specificity of DTPI and delayed time-point imaging depends on multiple factors, including the prevalence of malignancies, the patient population, and the cut-off values (either SUV or retention index) used to define a malignancy. Thus, DTPI and delayed time-point imaging would be more useful if performed for evaluation of lesions in regions with significant background activity clearance over time (such as the liver, the spleen, the mediastinum), and if used in the evaluation of the extent of tumor involvement rather than in the characterization of the nature of any specific lesion. Acute infectious and non-infectious inflammatory lesions remain as the major culprit for diminished diagnostic performance of these approaches (especially in tuberculosis-endemic regions). Tumor heterogeneity may also contribute to inconsistent performance of DTPI. The authors believe that selective use of DTPI and delayed time-point imaging will improve diagnostic accuracy and

Conversion of fracture toughness testing values from small scale three point bending test specimens to small scale yielding state (SSY) by elastic-plastic stress analysis

International Nuclear Information System (INIS)

Ikonen, K.

1993-07-01

The report describes the work performed for achieving readiness to calculate fracture toughness dependence on dimension effects and loading conditions in fracture test specimens and real structures. In the report two- and three-dimensional computer codes developed and calculational methods applied are described. One of the main goals is to converse fracture toughness from small scale three point bending test specimens to case of a depth crack in plane strain i.e. to small scale yielding state (SSY) by numerical elastic-plastic stress analysis. Thickness effect of a test specimens and effect of a crack depth are separately investigated. Tests of three point bending specimens with and without sidegrooves and curved crack front are numerically simulated and experimental and computed results are compared. J-integral is calculated along crack front and also from force-deflection dependence of the beam. For the analyses the computing system was thoroughly automatized. Measuring capacity of three point bending test specimens was tried to evaluate. (orig.) (7 refs., 54 figs.)

Modelling occupants’ heating set-point prefferences

DEFF Research Database (Denmark)

Andersen, Rune Vinther; Olesen, Bjarne W.; Toftum, Jørn

2011-01-01

consumption. Simultaneous measurement of the set-point of thermostatic radiator valves (trv), and indoor and outdoor environment characteristics was carried out in 15 dwellings in Denmark in 2008. Linear regression was used to infer a model of occupants’ interactions with trvs. This model could easily...... be implemented in most simulation software packages to increase the validity of the simulation outcomes....

Numerical simulation of collision-free plasma using Vlasov hybrid simulation

International Nuclear Information System (INIS)

Nunn, D.

1990-01-01

A novel scheme for the numerical simulation of wave particle interactions in space plasmas has been developed. The method, termed VHS or Vlasov Hybrid Simulation, is applicable to hot collision free plasmas in which the unperturbed distribution functions is smooth and free of delta function singularities. The particle population is described as a continuous Vlasov fluid in phase space-granularity and collisional effects being ignored. In traditional PIC/CIC codes the charge/current due to each simulation particle is assigned to a fixed spatial grid. In the VHS method the simulation particles sample the Vlasov fluid and provide information about the value of distribution function (F(r,v) at random points in phase space. Values of F are interpolated from the simulation particles onto a fixed grid in velocity/position or phase space. With distribution function defined on a phase space grid the plasma charge/current field is quickly calculated. The simulation particles serve only to provide information, and thus the particle population may be dynamic. Particles no longer resonant with the wavefield may be discarded from the simulation, and new particles may be inserted into the Vlasov fluid where required

Numerical analysis of sandwich beam with corrugated core under three-point bending

Energy Technology Data Exchange (ETDEWEB)

Wittenbeck, Leszek [Poznan University of Technology, Institute of Mathematics Piotrowo Street No. 5, 60-965 Poznan (Poland); Grygorowicz, Magdalena; Paczos, Piotr [Poznan University of Technology, Institute of Applied Mechanics Jana Pawla IIStreet No. 24, 60-965 Poznan (Poland)

2015-03-10

The strength problem of sandwich beam with corrugated core under three-point bending is presented.The beam are made of steel and formed by three mutually orthogonal corrugated layers. The finite element analysis (FEA) of the sandwich beam is performed with the use of the FEM system - ABAQUS. The relationship between the applied load and deflection in three-point bending is considered.

Comparison of solar photospheric bright points between Sunrise observations and MHD simulations

Science.gov (United States)

Riethmüller, T. L.; Solanki, S. K.; Berdyugina, S. V.; Schüssler, M.; Martínez Pillet, V.; Feller, A.; Gandorfer, A.; Hirzberger, J.

2014-08-01

Bright points (BPs) in the solar photosphere are thought to be the radiative signatures (small-scale brightness enhancements) of magnetic elements described by slender flux tubes or sheets located in the darker intergranular lanes in the solar photosphere. They contribute to the ultraviolet (UV) flux variations over the solar cycle and hence may play a role in influencing the Earth's climate. Here we aim to obtain a better insight into their properties by combining high-resolution UV and spectro-polarimetric observations of BPs by the Sunrise Observatory with 3D compressible radiation magnetohydrodynamical (MHD) simulations. To this end, full spectral line syntheses are performed with the MHD data and a careful degradation is applied to take into account all relevant instrumental effects of the observations. In a first step it is demonstrated that the selected MHD simulations reproduce the measured distributions of intensity at multiple wavelengths, line-of-sight velocity, spectral line width, and polarization degree rather well. The simulated line width also displays the correct mean, but a scatter that is too small. In the second step, the properties of observed BPs are compared with synthetic ones. Again, these are found to match relatively well, except that the observations display a tail of large BPs with strong polarization signals (most likely network elements) not found in the simulations, possibly due to the small size of the simulation box. The higher spatial resolution of the simulations has a significant effect, leading to smaller and more numerous BPs. The observation that most BPs are weakly polarized is explained mainly by the spatial degradation, the stray light contamination, and the temperature sensitivity of the Fe i line at 5250.2 Å. Finally, given that the MHD simulations are highly consistent with the observations, we used the simulations to explore the properties of BPs further. The Stokes V asymmetries increase with the distance to the

Numerical simulation of particle settling and cohesion in liquid

Energy Technology Data Exchange (ETDEWEB)

Johno, Y; Nakashima, K; Shigematsu, T; Ono, B [SASEBO National College of Technology, 1-1 Okishin, Sasebo, Nagasaki, 857-1193 (Japan); Satomi, M, E-mail: yjohno@post.cc.sasebo.ac.j [Sony Semiconductor Kyushu Corporation, Kikuchigun, Kumamoto (Japan)

2009-02-01

In this study, the motions of particles and particle clusters in liquid were numerically simulated. The particles of two sizes (Dp=40mum and 20mum) settle while repeating cohesion and dispersion, and finally the sediment of particles are formed at the bottom of a hexahedron container which is filled up with pure water. The flow field was solved with the Navier-Stokes equations and the particle motions were solved with the Lagrangian-type motion equations, where the interaction between fluid and particles due to drag forces were taken into account. The collision among particles was calculated using Distinct Element Method (DEM), and the effects of cohesive forces by van der Waals force acting on particle contact points were taken into account. Numerical simulations were performed under conditions in still flow and in shear flow. It was found that the simulation results enable us to know the state of the particle settling and the particle condensation.

A method for computing the stationary points of a function subject to linear equality constraints

International Nuclear Information System (INIS)

Uko, U.L.

1989-09-01

We give a new method for the numerical calculation of stationary points of a function when it is subject to equality constraints. An application to the solution of linear equations is given, together with a numerical example. (author). 5 refs

Controllable resonant tunnelling through single-point potentials: A point triode

International Nuclear Information System (INIS)

Zolotaryuk, A.V.; Zolotaryuk, Yaroslav

2015-01-01

A zero-thickness limit of three-layer heterostructures under two bias voltages applied externally, where one of which is supposed to be a gate parameter, is studied. As a result, an effect of controllable resonant tunnelling of electrons through single-point potentials is shown to exist. Therefore the limiting structure may be termed a “point triode” and considered in the theory of point interactions as a new object. The simple limiting analytical expressions adequately describe the resonant behaviour in the transistor with realistic parameter values and thus one can conclude that the zero-range limit of multi-layer structures may be used in fabricating nanodevices. The difference between the resonant tunnelling across single-point potentials and the Fabry–Pérot interference effect is also emphasized. - Highlights: • The zero-thickness limit of three-layer heterostructures is described in terms of point interactions. • The effect of resonant tunnelling through these single-point potentials is established. • The resonant tunnelling is shown to be controlled by a gate voltage

Interesting Interest Points

DEFF Research Database (Denmark)

Aanæs, Henrik; Dahl, Anders Lindbjerg; Pedersen, Kim Steenstrup

2012-01-01

on spatial invariance of interest points under changing acquisition parameters by measuring the spatial recall rate. The scope of this paper is to investigate the performance of a number of existing well-established interest point detection methods. Automatic performance evaluation of interest points is hard......Not all interest points are equally interesting. The most valuable interest points lead to optimal performance of the computer vision method in which they are employed. But a measure of this kind will be dependent on the chosen vision application. We propose a more general performance measure based...... position. The LED illumination provides the option for artificially relighting the scene from a range of light directions. This data set has given us the ability to systematically evaluate the performance of a number of interest point detectors. The highlights of the conclusions are that the fixed scale...

Numerical study of unsteady MHD oblique stagnation point flow and heat transfer due to an oscillating stream

Science.gov (United States)

Javed, T.; Ghaffari, A.; Ahmad, H.

2016-05-01

The unsteady stagnation point flow impinging obliquely on a flat plate in presence of a uniform applied magnetic field due to an oscillating stream has been studied. The governing partial differential equations are transformed into dimensionless form and the stream function is expressed in terms of Hiemenz and tangential components. The dimensionless partial differential equations are solved numerically by using well-known implicit finite difference scheme named as Keller-box method. The obtained results are compared with those available in the literature. It is observed that the results are in excellent agreement with the previous studies. The effects of pertinent parameters involved in the problem namely magnetic parameter, Prandtl number and impinging angle on flow and heat transfer characteristics are illustrated through graphs. It is observed that the influence of magnetic field strength increases the fluid velocity and by the increase of obliqueness parameter, the skin friction increases.

Three points of view in transport theory

International Nuclear Information System (INIS)

Ruben, Panta Pazos; Tilio de Vilhena, M.

2001-01-01

A lot of efforts in Transport Theory is used to develop numerical methods or hybrid numerical-analytical techniques. We present in this work three points of view about transport problems. First the C0 semigroup approach, in which the free transport operator ψ → μ ∇ generates an strongly continuous semigroup. The operators operator ψ → σt and operator ψ → ∫ ∇ k(x,μ,μ' ψ(x,μ') dμ' are bounded operators, and by perturbation the transport operator ψ → μ ∇ ψ + σt ψ - K ψ also generates an strongly continuous semigroup. To prove the convergence of the approximations of a numerical methods to the exact solution we use the approximation theorem of C0 semi-groups in canonical form. In other way, the discrete schemes theory is employed in searching the rate of convergence of numerical techniques in transport theory. For 1D dependent of time transport problem and two-dimensional steady state problem we summarize some estimates, incorporating different boundary conditions. Finally we give a survey about the dynamical behavior of the SN approximations. In order to give a unified approach, some results illustrates the equivalence of the three points of views for the case of the steady-state transport problem for slab geometry. (author)

Characterizing fixed points

Directory of Open Access Journals (Sweden)

Sanjo Zlobec

2017-04-01

Full Text Available A set of sufficient conditions which guarantee the existence of a point x⋆ such that f(x⋆ = x⋆ is called a "fixed point theorem". Many such theorems are named after well-known mathematicians and economists. Fixed point theorems are among most useful ones in applied mathematics, especially in economics and game theory. Particularly important theorem in these areas is Kakutani's fixed point theorem which ensures existence of fixed point for point-to-set mappings, e.g., [2, 3, 4]. John Nash developed and applied Kakutani's ideas to prove the existence of (what became known as "Nash equilibrium" for finite games with mixed strategies for any number of players. This work earned him a Nobel Prize in Economics that he shared with two mathematicians. Nash's life was dramatized in the movie "Beautiful Mind" in 2001. In this paper, we approach the system f(x = x differently. Instead of studying existence of its solutions our objective is to determine conditions which are both necessary and sufficient that an arbitrary point x⋆ is a fixed point, i.e., that it satisfies f(x⋆ = x⋆. The existence of solutions for continuous function f of the single variable is easy to establish using the Intermediate Value Theorem of Calculus. However, characterizing fixed points x⋆, i.e., providing answers to the question of finding both necessary and sufficient conditions for an arbitrary given x⋆ to satisfy f(x⋆ = x⋆, is not simple even for functions of the single variable. It is possible that constructive answers do not exist. Our objective is to find them. Our work may require some less familiar tools. One of these might be the "quadratic envelope characterization of zero-derivative point" recalled in the next section. The results are taken from the author's current research project "Studying the Essence of Fixed Points". They are believed to be original. The author has received several feedbacks on the preliminary report and on parts of the project

Tipping Point

Medline Plus

Full Text Available ... OnSafety CPSC Stands for Safety The Tipping Point Home > 60 Seconds of Safety (Videos) > The Tipping Point ... 24 hours a day. For young children whose home is a playground, it’s the best way to ...

Experiment and simulation study on unidirectional carbon fiber composite component under dynamic 3 point bending loading

Energy Technology Data Exchange (ETDEWEB)

Zhou, Guowei; Sun, Qingping; Zeng, Danielle; Li, Dayong; Su, Xuming

2018-04-10

In current work, unidirectional (UD) carbon fiber composite hatsection component with two different layups are studied under dynamic 3 point bending loading. The experiments are performed at various impact velocities, and the effects of impactor velocity and layup on acceleration histories are compared. A macro model is established with LS-Dyna for more detailed study. The simulation results show that the delamination plays an important role during dynamic 3 point bending test. Based on the analysis with high speed camera, the sidewall of hatsection shows significant buckling rather than failure. Without considering the delamination, current material model cannot capture the post failure phenomenon correctly. The sidewall delamination is modeled by assumption of larger failure strain together with slim parameters, and the simulation results of different impact velocities and layups match the experimental results reasonable well.

Point-point and point-line moving-window correlation spectroscopy and its applications

Science.gov (United States)

Zhou, Qun; Sun, Suqin; Zhan, Daqi; Yu, Zhiwu

2008-07-01

In this paper, we present a new extension of generalized two-dimensional (2D) correlation spectroscopy. Two new algorithms, namely point-point (P-P) correlation and point-line (P-L) correlation, have been introduced to do the moving-window 2D correlation (MW2D) analysis. The new method has been applied to a spectral model consisting of two different processes. The results indicate that P-P correlation spectroscopy can unveil the details and re-constitute the entire process, whilst the P-L can provide general feature of the concerned processes. Phase transition behavior of dimyristoylphosphotidylethanolamine (DMPE) has been studied using MW2D correlation spectroscopy. The newly proposed method verifies that the phase transition temperature is 56 °C, same as the result got from a differential scanning calorimeter. To illustrate the new method further, a lysine and lactose mixture has been studied under thermo perturbation. Using the P-P MW2D, the Maillard reaction of the mixture was clearly monitored, which has been very difficult using conventional display of FTIR spectra.

Numerical simulation of the knotted nylon netting panel

Directory of Open Access Journals (Sweden)

Li Yuwei

2016-01-01

Full Text Available A piece of netting, consists of the 8 8 meshes, fixed on a square frame, was simulated and the tensions and their distribution, the positions of knots and netting shape were calculated by means of MATLAB in computer. The dynamic mathematic model was developed based on lumped mass method, the netting was treated as spring-mass system, the Runge-Kutta fifth-order and sixth-order method was used to solve the differential equations for every step, then the displacement and tension of each mass point were obtained. For verify this model, the tests have been carried out in a flume tank. The results of the numerical simulation fully agreed with the experiments.

Zero-point oscillations, zero-point fluctuations, and fluctuations of zero-point oscillations

International Nuclear Information System (INIS)

Khalili, Farit Ya

2003-01-01

Several physical effects and methodological issues relating to the ground state of an oscillator are considered. Even in the simplest case of an ideal lossless harmonic oscillator, its ground state exhibits properties that are unusual from the classical point of view. In particular, the mean value of the product of two non-negative observables, kinetic and potential energies, is negative in the ground state. It is shown that semiclassical and rigorous quantum approaches yield substantially different results for the ground state energy fluctuations of an oscillator with finite losses. The dependence of zero-point fluctuations on the boundary conditions is considered. Using this dependence, it is possible to transmit information without emitting electromagnetic quanta. Fluctuations of electromagnetic pressure of zero-point oscillations are analyzed, and the corresponding mechanical friction is considered. This friction can be viewed as the most fundamental mechanism limiting the quality factor of mechanical oscillators. Observation of these effects exceeds the possibilities of contemporary experimental physics but almost undoubtedly will be possible in the near future. (methodological notes)

Reduced-complexity adaptive multi-channel assignment for shared access points in over-loaded small-cell networks

KAUST Repository

Radaydeh, Redha Mahmoud

2013-06-01

This paper proposes a reduced-complexity downlink multi-channel assignment scheme when feedback links are capacity-limited. The system model treats the case when multiple access points are allocated to serve scheduled users in over-loaded (i.e. dense) pico/femtocell networks. It assumes that the deployed access points can be shared simultaneously and employ isotropic antenna arrays of arbitrary sizes. Moreover, they transmit their data on a common physical channel and can not coordinate their transmissions. On the other hand, each scheduled user can be served by single transmit channel from each active access point at a time, and it lacks coordination with concurrent active users. The scheme operates according to the occupancy of available transmit channels, wherein extensively occupied access points are avoided adaptively, while reducing the load of processing. The operation is linked to a target performance via controlling the observed aggregate interference from the projected set of serving points. Through the analysis, results for the scheduled user outage performance, and the average number of active access points are presented. Numerical and simulations studies clarify the gains of the proposed scheme for different operating conditions. © 2013 IEEE.

Reduced-complexity adaptive multi-channel assignment for shared access points in over-loaded small-cell networks

KAUST Repository

Radaydeh, Redha Mahmoud; Qaraqe, Khalid A.; Alouini, Mohamed-Slim

2013-01-01

This paper proposes a reduced-complexity downlink multi-channel assignment scheme when feedback links are capacity-limited. The system model treats the case when multiple access points are allocated to serve scheduled users in over-loaded (i.e. dense) pico/femtocell networks. It assumes that the deployed access points can be shared simultaneously and employ isotropic antenna arrays of arbitrary sizes. Moreover, they transmit their data on a common physical channel and can not coordinate their transmissions. On the other hand, each scheduled user can be served by single transmit channel from each active access point at a time, and it lacks coordination with concurrent active users. The scheme operates according to the occupancy of available transmit channels, wherein extensively occupied access points are avoided adaptively, while reducing the load of processing. The operation is linked to a target performance via controlling the observed aggregate interference from the projected set of serving points. Through the analysis, results for the scheduled user outage performance, and the average number of active access points are presented. Numerical and simulations studies clarify the gains of the proposed scheme for different operating conditions. © 2013 IEEE.

Realworld maximum power point tracking simulation of PV system based on Fuzzy Logic control

Science.gov (United States)

Othman, Ahmed M.; El-arini, Mahdi M. M.; Ghitas, Ahmed; Fathy, Ahmed

2012-12-01

In the recent years, the solar energy becomes one of the most important alternative sources of electric energy, so it is important to improve the efficiency and reliability of the photovoltaic (PV) systems. Maximum power point tracking (MPPT) plays an important role in photovoltaic power systems because it maximize the power output from a PV system for a given set of conditions, and therefore maximize their array efficiency. This paper presents a maximum power point tracker (MPPT) using Fuzzy Logic theory for a PV system. The work is focused on the well known Perturb and Observe (P&O) algorithm and is compared to a designed fuzzy logic controller (FLC). The simulation work dealing with MPPT controller; a DC/DC Ćuk converter feeding a load is achieved. The results showed that the proposed Fuzzy Logic MPPT in the PV system is valid.

Numerical flow simulation methods and additive manufacturing methods for the development of a flow optimised design of a novel point-of-care diagnostic device

Directory of Open Access Journals (Sweden)

Dethloff Manuel

2017-09-01

Full Text Available For the development of a novel, user-friendly and low cost point-of-care diagnostic device for the detection of disease specific biomarker a flow optimised design of the test system has to be investigated. The resulting test system is characterised by a reduced execution period, a reduction of execution steps and an integrated waste management. Based on previous results, the current study focused on the design implementation of the fluidic requirements, e. g. tightness, inside the test device. With the help of fluid flow simulations (CFD – computational fluid dynamics the flow behaviour inside the test device was analysed for different designs and arrangements. Prototypes generated from additive manufacturing technologies (PolyJet modeling are used for validating the simulation results and further experimental tests.

Fixed point of the parabolic renormalization operator

CERN Document Server

Lanford III, Oscar E

2014-01-01

This monograph grew out of the authors' efforts to provide a natural geometric description for the class of maps invariant under parabolic renormalization and for the Inou-Shishikura fixed point itself as well as to carry out a computer-assisted study of the parabolic renormalization operator. It introduces a renormalization-invariant class of analytic maps with a maximal domain of analyticity and rigid covering properties and presents a numerical scheme for computing parabolic renormalization of a germ, which is used to compute the Inou-Shishikura renormalization fixed point.   Inside, readers will find a detailed introduction into the theory of parabolic bifurcation,  Fatou coordinates, Écalle-Voronin conjugacy invariants of parabolic germs, and the definition and basic properties of parabolic renormalization.   The systematic view of parabolic renormalization developed in the book and the numerical approach to its study will be interesting to both experts in the field as well as graduate students wishi...

Temperature distribution model for the semiconductor dew point detector

Science.gov (United States)

Weremczuk, Jerzy; Gniazdowski, Z.; Jachowicz, Ryszard; Lysko, Jan M.

2001-08-01

The simulation results of temperature distribution in the new type silicon dew point detector are presented in this paper. Calculations were done with use of the SMACEF simulation program. Fabricated structures, apart from the impedance detector used to the dew point detection, contained the resistive four terminal thermometer and two heaters. Two detector structures, the first one located on the silicon membrane and the second one placed on the bulk materials were compared in this paper.

Bottom Slamming on Heaving Point Absorber Wave Energy Devices

DEFF Research Database (Denmark)

De Backer, Griet; Vantorre, Marc; Frigaard, Peter

2010-01-01

shapes are considered: a hemisphere and two conical shapes with deadrise angles of 30 and 45, with a waterline diameter of 5 m. The simulations indicate that the risk of rising out of the water is largely dependent on the buoy draft and sea state. Although associated with power losses, emergence......Oscillating point absorber buoys may rise out of the water and be subjected to bottom slamming upon re-entering the water. Numerical simulations are performed to estimate the power absorption, the impact velocities and the corresponding slamming forces for various slamming constraints. Three buoy...... occurrence probabilities can be significantly reduced by adapting the control parameters. The magnitude of the slamming load is severely influenced by the buoy shape. The ratio between the peak impact load on the hemisphere and that on the 45 cone is approximately 2, whereas the power absorption is only 4...

Pairwise contact energy statistical potentials can help to find probability of point mutations.

Science.gov (United States)

Saravanan, K M; Suvaithenamudhan, S; Parthasarathy, S; Selvaraj, S

2017-01-01

To adopt a particular fold, a protein requires several interactions between its amino acid residues. The energetic contribution of these residue-residue interactions can be approximated by extracting statistical potentials from known high resolution structures. Several methods based on statistical potentials extracted from unrelated proteins are found to make a better prediction of probability of point mutations. We postulate that the statistical potentials extracted from known structures of similar folds with varying sequence identity can be a powerful tool to examine probability of point mutation. By keeping this in mind, we have derived pairwise residue and atomic contact energy potentials for the different functional families that adopt the (α/β) 8 TIM-Barrel fold. We carried out computational point mutations at various conserved residue positions in yeast Triose phosphate isomerase enzyme for which experimental results are already reported. We have also performed molecular dynamics simulations on a subset of point mutants to make a comparative study. The difference in pairwise residue and atomic contact energy of wildtype and various point mutations reveals probability of mutations at a particular position. Interestingly, we found that our computational prediction agrees with the experimental studies of Silverman et al. (Proc Natl Acad Sci 2001;98:3092-3097) and perform better prediction than i Mutant and Cologne University Protein Stability Analysis Tool. The present work thus suggests deriving pairwise contact energy potentials and molecular dynamics simulations of functionally important folds could help us to predict probability of point mutations which may ultimately reduce the time and cost of mutation experiments. Proteins 2016; 85:54-64. © 2016 Wiley Periodicals, Inc. © 2016 Wiley Periodicals, Inc.

Fixed Points

Indian Academy of Sciences (India)

Home; Journals; Resonance – Journal of Science Education; Volume 5; Issue 5. Fixed Points - From Russia with Love - A Primer of Fixed Point Theory. A K Vijaykumar. Book Review Volume 5 Issue 5 May 2000 pp 101-102. Fulltext. Click here to view fulltext PDF. Permanent link:

Effect of Finite Particle Size on Convergence of Point Particle Models in Euler-Lagrange Multiphase Dispersed Flow

Science.gov (United States)

Nili, Samaun; Park, Chanyoung; Haftka, Raphael T.; Kim, Nam H.; Balachandar, S.

2017-11-01

Point particle methods are extensively used in simulating Euler-Lagrange multiphase dispersed flow. When particles are much smaller than the Eulerian grid the point particle model is on firm theoretical ground. However, this standard approach of evaluating the gas-particle coupling at the particle center fails to converge as the Eulerian grid is reduced below particle size. We present an approach to model the interaction between particles and fluid for finite size particles that permits convergence. We use the generalized Faxen form to compute the force on a particle and compare the results against traditional point particle method. We apportion the different force components on the particle to fluid cells based on the fraction of particle volume or surface in the cell. The application is to a one-dimensional model of shock propagation through a particle-laden field at moderate volume fraction, where the convergence is achieved for a well-formulated force model and back coupling for finite size particles. Comparison with 3D direct fully resolved numerical simulations will be used to check if the approach also improves accuracy compared to the point particle model. Work supported by the U.S. Department of Energy, National Nuclear Security Administration, Advanced Simulation and Computing Program, as a Cooperative Agreement under the Predictive Science Academic Alliance Program, under Contract No. DE-NA0002378.

Numerical Analysis Of The Resistance To Pullout Test Of Clinched Assemblies Of Thin Metal Sheets

International Nuclear Information System (INIS)

Jomaa, Moez; Billardon, Rene

2007-01-01

This paper presents the finite element analysis of the resistance of a clinch point to pullout test -that follows the numerical analysis of the forming process of the point-. The simulations have been validated by comparison with experimental evidences. The influence on the numerical predictions of various computation and process parameters have been evaluated

Csf Based Non-Ground Points Extraction from LIDAR Data

Science.gov (United States)

Shen, A.; Zhang, W.; Shi, H.

2017-09-01

Region growing is a classical method of point cloud segmentation. Based on the idea of collecting the pixels with similar properties to form regions, region growing is widely used in many fields such as medicine, forestry and remote sensing. In this algorithm, there are two core problems. One is the selection of seed points, the other is the setting of the growth constraints, in which the selection of the seed points is the foundation. In this paper, we propose a CSF (Cloth Simulation Filtering) based method to extract the non-ground seed points effectively. The experiments have shown that this method can obtain a group of seed spots compared with the traditional methods. It is a new attempt to extract seed points

Demerit points systems.

NARCIS (Netherlands)

2006-01-01

In 2012, 21 of the 27 EU Member States had some form of demerit points system. In theory, demerit points systems contribute to road safety through three mechanisms: 1) prevention of unsafe behaviour through the risk of receiving penalty points, 2) selection and suspension of the most frequent

Pointing and the Evolution of Language: An Applied Evolutionary Epistemological Approach

Directory of Open Access Journals (Sweden)

Nathalie Gontier

2013-07-01

Full Text Available Numerous evolutionary linguists have indicated that human pointing behaviour might be associated with the evolution of language. At an ontogenetic level, and in normal individuals, pointing develops spontaneously and the onset of human pointing precedes as well as facilitates phases in speech and language development. Phylogenetically, pointing behaviour might have preceded and facilitated the evolutionary origin of both gestural and vocal language. Contrary to wild non-human primates, captive and human-reared nonhuman primates also demonstrate pointing behaviour. In this article, we analyse the debates on pointing and its role it might have played in language evolution from a meta-level. From within an Applied Evolutionary Epistemological approach, we examine how exactly we can determine whether pointing has been a unit, a level or a mechanism in language evolution.

High-Temperature Gas-Cooled Test Reactor Point Design

Energy Technology Data Exchange (ETDEWEB)

Sterbentz, James William [Idaho National Laboratory; Bayless, Paul David [Idaho National Laboratory; Nelson, Lee Orville [Idaho National Laboratory; Gougar, Hans David [Idaho National Laboratory; Kinsey, James Carl [Idaho National Laboratory; Strydom, Gerhard [Idaho National Laboratory; Kumar, Akansha [Idaho National Laboratory

2016-04-01

A point design has been developed for a 200 MW high-temperature gas-cooled test reactor. The point design concept uses standard prismatic blocks and 15.5% enriched UCO fuel. Reactor physics and thermal-hydraulics simulations have been performed to characterize the capabilities of the design. In addition to the technical data, overviews are provided on the technological readiness level, licensing approach and costs.

Analytical and numerical construction of vertical periodic orbits about triangular libration points based on polynomial expansion relations among directions

Science.gov (United States)

Qian, Ying-Jing; Yang, Xiao-Dong; Zhai, Guan-Qiao; Zhang, Wei

2017-08-01

Innovated by the nonlinear modes concept in the vibrational dynamics, the vertical periodic orbits around the triangular libration points are revisited for the Circular Restricted Three-body Problem. The ζ -component motion is treated as the dominant motion and the ξ and η -component motions are treated as the slave motions. The slave motions are in nature related to the dominant motion through the approximate nonlinear polynomial expansions with respect to the ζ -position and ζ -velocity during the one of the periodic orbital motions. By employing the relations among the three directions, the three-dimensional system can be transferred into one-dimensional problem. Then the approximate three-dimensional vertical periodic solution can be analytically obtained by solving the dominant motion only on ζ -direction. To demonstrate the effectiveness of the proposed method, an accuracy study was carried out to validate the polynomial expansion (PE) method. As one of the applications, the invariant nonlinear relations in polynomial expansion form are used as constraints to obtain numerical solutions by differential correction. The nonlinear relations among the directions provide an alternative point of view to explore the overall dynamics of periodic orbits around libration points with general rules.

Magnetic neutral point stretching and coalescence in tearing-generated magnetohydrodynamic structures

International Nuclear Information System (INIS)

Malara, F.; Veltri, P.; Carbone, V.

1991-01-01

The time evolution of the instability of a sheet pinch is numerically studied using a sufficiently high ratio of system length to width in order to allow the simultaneous growth of several unstable wavelengths. This numerical simulation provides new insights into the nonlinear development of the tearing instability. Before the instability saturates, the nonlinear interactions among the unstable modes produce local coalescence phenomena that destroy the weaker current pinches and reduce the number of magnetic islands. In contrast with the usual picture, this coalescence is not due to the attraction between the current maxima, but is due to the stretching of the X-neutral points associated with the most intense current pinches. The global perturbation growth rate remains essentially unchanged in time, being of the order of the resistive instability growth rate

Nonuniform fast Fourier transform method for numerical diffraction simulation on tilted planes.

Science.gov (United States)

Xiao, Yu; Tang, Xiahui; Qin, Yingxiong; Peng, Hao; Wang, Wei; Zhong, Lijing

2016-10-01

The method, based on the rotation of the angular spectrum in the frequency domain, is generally used for the diffraction simulation between the tilted planes. Due to the rotation of the angular spectrum, the interval between the sampling points in the Fourier domain is not even. For the conventional fast Fourier transform (FFT)-based methods, a spectrum interpolation is needed to get the approximate sampling value on the equidistant sampling points. However, due to the numerical error caused by the spectrum interpolation, the calculation accuracy degrades very quickly as the rotation angle increases. Here, the diffraction propagation between the tilted planes is transformed into a problem about the discrete Fourier transform on the uneven sampling points, which can be evaluated effectively and precisely through the nonuniform fast Fourier transform method (NUFFT). The most important advantage of this method is that the conventional spectrum interpolation is avoided and the high calculation accuracy can be guaranteed for different rotation angles, even when the rotation angle is close to π/2. Also, its calculation efficiency is comparable with that of the conventional FFT-based methods. Numerical examples as well as a discussion about the calculation accuracy and the sampling method are presented.

Point spread functions and deconvolution of ultrasonic images.

Science.gov (United States)

Dalitz, Christoph; Pohle-Fröhlich, Regina; Michalk, Thorsten

2015-03-01

This article investigates the restoration of ultrasonic pulse-echo C-scan images by means of deconvolution with a point spread function (PSF). The deconvolution concept from linear system theory (LST) is linked to the wave equation formulation of the imaging process, and an analytic formula for the PSF of planar transducers is derived. For this analytic expression, different numerical and analytic approximation schemes for evaluating the PSF are presented. By comparing simulated images with measured C-scan images, we demonstrate that the assumptions of LST in combination with our formula for the PSF are a good model for the pulse-echo imaging process. To reconstruct the object from a C-scan image, we compare different deconvolution schemes: the Wiener filter, the ForWaRD algorithm, and the Richardson-Lucy algorithm. The best results are obtained with the Richardson-Lucy algorithm with total variation regularization. For distances greater or equal twice the near field distance, our experiments show that the numerically computed PSF can be replaced with a simple closed analytic term based on a far field approximation.

Statistical aspects of determinantal point processes

DEFF Research Database (Denmark)

Lavancier, Frédéric; Møller, Jesper; Rubak, Ege Holger

The statistical aspects of determinantal point processes (DPPs) seem largely unexplored. We review the appealing properties of DDPs, demonstrate that they are useful models for repulsiveness, detail a simulation procedure, and provide freely available software for simulation and statistical...... inference. We pay special attention to stationary DPPs, where we give a simple condition ensuring their existence, construct parametric models, describe how they can be well approximated so that the likelihood can be evaluated and realizations can be simulated, and discuss how statistical inference...

An unsteady point vortex method for coupled fluid-solid problems

Energy Technology Data Exchange (ETDEWEB)

Michelin, Sebastien [Jacobs School of Engineering, UCSD, Department of Mechanical and Aerospace Engineering, La Jolla, CA (United States); Ecole Nationale Superieure des Mines de Paris, Paris (France); Llewellyn Smith, Stefan G. [Jacobs School of Engineering, UCSD, Department of Mechanical and Aerospace Engineering, La Jolla, CA (United States)

2009-06-15

A method is proposed for the study of the two-dimensional coupled motion of a general sharp-edged solid body and a surrounding inviscid flow. The formation of vorticity at the body's edges is accounted for by the shedding at each corner of point vortices whose intensity is adjusted at each time step to satisfy the regularity condition on the flow at the generating corner. The irreversible nature of vortex shedding is included in the model by requiring the vortices' intensity to vary monotonically in time. A conservation of linear momentum argument is provided for the equation of motion of these point vortices (Brown-Michael equation). The forces and torques applied on the solid body are computed as explicit functions of the solid body velocity and the vortices' position and intensity, thereby providing an explicit formulation of the vortex-solid coupled problem as a set of non-linear ordinary differential equations. The example of a falling card in a fluid initially at rest is then studied using this method. The stability of broadside-on fall is analysed and the shedding of vorticity from both plate edges is shown to destabilize this position, consistent with experimental studies and numerical simulations of this problem. The reduced-order representation of the fluid motion in terms of point vortices is used to understand the physical origin of this destabilization. (orig.)

Exploring Foundation Concepts in Introductory Statistics Using Dynamic Data Points

Science.gov (United States)

Ekol, George

2015-01-01

This paper analyses introductory statistics students' verbal and gestural expressions as they interacted with a dynamic sketch (DS) designed using "Sketchpad" software. The DS involved numeric data points built on the number line whose values changed as the points were dragged along the number line. The study is framed on aggregate…

Mechanistic spatio-temporal point process models for marked point processes, with a view to forest stand data

DEFF Research Database (Denmark)

Møller, Jesper; Ghorbani, Mohammad; Rubak, Ege Holger

We show how a spatial point process, where to each point there is associated a random quantitative mark, can be identified with a spatio-temporal point process specified by a conditional intensity function. For instance, the points can be tree locations, the marks can express the size of trees......, and the conditional intensity function can describe the distribution of a tree (i.e., its location and size) conditionally on the larger trees. This enable us to construct parametric statistical models which are easily interpretable and where likelihood-based inference is tractable. In particular, we consider maximum...

Three points of view in transport theory

Energy Technology Data Exchange (ETDEWEB)

Ruben, Panta Pazos [Faculdade de Matematica, PUCRS, Porto Alegre, RS (Brazil); Tilio de Vilhena, M. [Instituto de Matematica, UFRGS, Porto Alegre, RS (Brazil)

2001-07-01

A lot of efforts in Transport Theory is used to develop numerical methods or hybrid numerical-analytical techniques. We present in this work three points of view about transport problems. First the C0 semigroup approach, in which the free transport operator {psi} {yields} {mu} {nabla} generates an strongly continuous semigroup. The operators operator {psi} {yields} {sigma}t and operator {psi} {yields} {integral} {nabla} k(x,{mu},{mu}') {psi}(x,{mu}') d{mu}' are bounded operators, and by perturbation the transport operator {psi} {yields} {mu} {nabla} {psi} + {sigma}t {psi} - K {psi} also generates an strongly continuous semigroup. To prove the convergence of the approximations of a numerical methods to the exact solution we use the approximation theorem of C0 semi-groups in canonical form. In other way, the discrete schemes theory is employed in searching the rate of convergence of numerical techniques in transport theory. For 1D dependent of time transport problem and two-dimensional steady state problem we summarize some estimates, incorporating different boundary conditions. Finally we give a survey about the dynamical behavior of the SN approximations. In order to give a unified approach, some results illustrates the equivalence of the three points of views for the case of the steady-state transport problem for slab geometry. (author)

Realworld maximum power point tracking simulation of PV system based on Fuzzy Logic control

Directory of Open Access Journals (Sweden)

Ahmed M. Othman

2012-12-01

Full Text Available In the recent years, the solar energy becomes one of the most important alternative sources of electric energy, so it is important to improve the efficiency and reliability of the photovoltaic (PV systems. Maximum power point tracking (MPPT plays an important role in photovoltaic power systems because it maximize the power output from a PV system for a given set of conditions, and therefore maximize their array efficiency. This paper presents a maximum power point tracker (MPPT using Fuzzy Logic theory for a PV system. The work is focused on the well known Perturb and Observe (P&O algorithm and is compared to a designed fuzzy logic controller (FLC. The simulation work dealing with MPPT controller; a DC/DC Ćuk converter feeding a load is achieved. The results showed that the proposed Fuzzy Logic MPPT in the PV system is valid.

Ecological association between HIV and concurrency point ...

African Journals Online (AJOL)

Ecological association between HIV and concurrency point-prevalence in South Africa's ethnic groups. Chris Kenyon. Abstract. HIV prevalence between different ethnic groups within South Africa exhibits considerable variation. Numerous authors believe that elevated sexual partner concurrency rates are important in the ...

Function parametrization by using 4-point transforms

International Nuclear Information System (INIS)

Dikusar, N.D.

1996-01-01

A continuous parametrization of the smooth curve f(x)=f(x;R) is suggested on a basis of four-point transformations. Coordinates of three reference points of the curve are chosen as parameters R. This approach allows to derive a number of advantages in function approximation and fitting of empiric data. The transformations have made possible to derive a new class of polynomials (monosplines) having the better approximation quality than monomials {x n }. A behaviour of an error of the approximation has a uniform character. A three-point model of the cubic spline (TPS) is proposed. The model allows to reduce a number of unknown parameters in twice and to obtain an advantage in a computing aspect. The new approach to the function approximation and fitting are shown on a number of examples. The proposed approach gives a new mathematical tool and a new possibility in both practical applications and theoretical research of numerical and computational methods. 13 refs., 13 figs., 2 tabs

Point specificity in acupuncture

Directory of Open Access Journals (Sweden)

Choi Emma M

2012-02-01

Full Text Available Abstract The existence of point specificity in acupuncture is controversial, because many acupuncture studies using this principle to select control points have found that sham acupoints have similar effects to those of verum acupoints. Furthermore, the results of pain-related studies based on visual analogue scales have not supported the concept of point specificity. In contrast, hemodynamic, functional magnetic resonance imaging and neurophysiological studies evaluating the responses to stimulation of multiple points on the body surface have shown that point-specific actions are present. This review article focuses on clinical and laboratory studies supporting the existence of point specificity in acupuncture and also addresses studies that do not support this concept. Further research is needed to elucidate the point-specific actions of acupuncture.

Comprehensive Interpretation of a Three-Point Gauss Quadrature with Variable Sampling Points and Its Application to Integration for Discrete Data

Directory of Open Access Journals (Sweden)

Young-Doo Kwon

2013-01-01

Full Text Available This study examined the characteristics of a variable three-point Gauss quadrature using a variable set of weighting factors and corresponding optimal sampling points. The major findings were as follows. The one-point, two-point, and three-point Gauss quadratures that adopt the Legendre sampling points and the well-known Simpson’s 1/3 rule were found to be special cases of the variable three-point Gauss quadrature. In addition, the three-point Gauss quadrature may have out-of-domain sampling points beyond the domain end points. By applying the quadratically extrapolated integrals and nonlinearity index, the accuracy of the integration could be increased significantly for evenly acquired data, which is popular with modern sophisticated digital data acquisition systems, without using higher-order extrapolation polynomials.

New position of Dirac points in the strained graphene reciprocal lattice

Directory of Open Access Journals (Sweden)

Cui-Lian Li

2014-08-01

Full Text Available In the strained graphene, Fermi velocity shows space-dependent and it changes as the position of Dirac point shifts. In this paper, we apply the tight-binding approach within linear elasticity theory to investigate the shifting of Dirac points in the strained graphene reciprocal lattice space. Based on this, we derive the analytical expression on the new positions of the Dirac points as the strain parameter varies. Comparing the data from our analytical expression, ones from Eq. (20 in Phys. Rev. B 80, 045401 (2009, and those from numerical calculation, we find that our analytical expression raises the effective prediction range of the strain parameter from 3% to 15%. i.e., our analytical expression is practicable until the strain parameter is larger than 15%. This almost includes the whole range where the Dirac points present and the energy gap is zero. Moreover, we further calculate the energy gap by numerical method when the shear strain parameter varies from 0 to 20%, and find that the energy gap can not open until the strain parameter is larger than 16%. After this, the energy gap open and the Dirac points disappear.

Near-point string: Simple method to demonstrate anticipated near point for multifocal and accommodating intraocular lenses.

Science.gov (United States)

George, Monica C; Lazer, Zane P; George, David S

2016-05-01

We present a technique that uses a near-point string to demonstrate the anticipated near point of multifocal and accommodating intraocular lenses (IOLs). Beads are placed on the string at distances corresponding to the near points for diffractive and accommodating IOLs. The string is held up to the patient's eye to demonstrate where each of the IOLs is likely to provide the best near vision. None of the authors has a financial or proprietary interest in any material or method mentioned. Copyright © 2016 ASCRS and ESCRS. Published by Elsevier Inc. All rights reserved.

Publication point indicators

DEFF Research Database (Denmark)

Elleby, Anita; Ingwersen, Peter

2010-01-01

; the Cumulated Publication Point Indicator (CPPI), which graphically illustrates the cumulated gain of obtained vs. ideal points, both seen as vectors; and the normalized Cumulated Publication Point Index (nCPPI) that represents the cumulated gain of publication success as index values, either graphically......The paper presents comparative analyses of two publication point systems, The Norwegian and the in-house system from the interdisciplinary Danish Institute of International Studies (DIIS), used as case in the study for publications published 2006, and compares central citation-based indicators...... with novel publication point indicators (PPIs) that are formalized and exemplified. Two diachronic citation windows are applied: 2006-07 and 2006-08. Web of Science (WoS) as well as Google Scholar (GS) are applied to observe the cite delay and citedness for the different document types published by DIIS...

Point defect relaxation volumes for copper

International Nuclear Information System (INIS)

Miller, K.M.

1979-11-01

The methods used for the determination of point defect relaxation volumes are discussed and it is shown that a previous interatomic potential derived for copper is inaccurate and results obtained using it are invalid. A new interatomic potential has been produced for copper and a computer simulation of point and planar defects carried out. A vacancy relaxation volume of -0.33 atomic volumes has been found with interstitial values in the range 1.7 to 2.0 atomic volumes. It is shown that these values in current theories of irradiation induced swelling lead to an anomalously high value for dislocation bias compared with that determined experimentally. (author)

Magnetic Reconnection at a Three-dimensional Solar Null Point

DEFF Research Database (Denmark)

Frederiksen, Jacob Trier; Baumann, Gisela; Galsgaard, Klaus

2012-01-01

Using a specific solar null point reconnection case studied by Masson et al (2009; ApJ 700, 559) we investigate the dependence of the reconnection rate on boundary driving speed, numerical resolution, type of resistivity (constant or numerical), and assumed stratification (constant density or sol...

Point-to-point radio link variation at E-band and its effect on antenna design

NARCIS (Netherlands)

Al-Rawi, A.; Dubok, A.; Herben, M.H.A.J.; Smolders, A.B.

2015-01-01

Radio propagation will strongly influence the design of the antenna and front-end components of E-band point-to-point communication systems. Based on the ITU rain model, the rain attenuation is estimated in a statistical sense and it is concluded that for backhaul links of 1–10 km, antennas with a

Continuation of connecting orbits in 3d-ODEs' (i) point-to-cycle connections.

NARCIS (Netherlands)

Doedel, E.J.; Kooi, B.W.; van Voorn, G.A.K.; Kuznetzov, Y.A.

2008-01-01

We propose new methods for the numerical continuation of point-to-cycle connecting orbits in three-dimensional autonomous ODE's using projection boundary conditions. In our approach, the projection boundary conditions near the cycle are formulated using an eigenfunction of the associated adjoint

Delay analysis of a point-to-multipoint spectrum sharing network with CSI based power allocation

KAUST Repository

Khan, Fahd Ahmed

2012-10-01

In this paper, we analyse the delay performance of a point-to-multipoint cognitive radio network which is sharing the spectrum with a point-to-multipoint primary network. The channel is assumed to be independent but not identically distributed and has Nakagami-m fading. A constraint on the peak transmit power of the secondary user transmitter (SU-Tx) is also considered in addition to the peak interference power constraint. Based on the constraints, a power allocation scheme which requires knowledge of the instantaneous channel state information (CSI) of the interference links is derived. The SU-Tx is assumed to be equipped with a buffer and is modelled using the M/G/1 queueing model. Closed form expressions for the probability distribution function (PDF) and cumulative distribution function (CDF) of the packet transmission time is derived. Using the PDF, the expressions for the moments of transmission time are obtained. In addition, using the moments, the expressions for the performance measures such as the total average waiting time of packets and the average number of packets waiting in the buffer of the SU-Tx are also obtained. Numerical simulations corroborate the theoretical results. © 2012 IEEE.

Single-point incremental forming and formability-failure diagrams

DEFF Research Database (Denmark)

Silva, M.B.; Skjødt, Martin; Atkins, A.G.

2008-01-01

In a recent work [1], the authors constructed a closed-form analytical model that is capable of dealing with the fundamentals of single point incremental forming and explaining the experimental and numerical results published in the literature over the past couple of years. The model is based...... of deformation that are commonly found in general single point incremental forming processes; and (ii) to investigate the formability limits of SPIF in terms of ductile damage mechanics and the question of whether necking does, or does not, precede fracture. Experimentation by the authors together with data...

Determination of point of incidence for the case of reflection or refraction at spherical surface knowing two points lying on the ray.

Science.gov (United States)

Mikš, Antonín; Novák, Pavel

2017-09-01

The paper is focused on the problem of determination of the point of incidence of a light ray for the case of reflection or refraction at the spherical optical surface, assuming that two fixed points in space that the sought light ray should go through are given. The requirement is that one of these points lies on the incident ray and the other point on the reflected/refracted ray. Although at first glance it seems to be a simple problem, it will be shown that it has no simple analytical solution. The basic idea of the solution is given, and it is shown that the problem leads to a nonlinear equation in one variable. The roots of the resulting nonlinear equation can be found by numerical methods of mathematical optimization. The proposed methods were implemented in MATLAB, and the proper function of these algorithms was verified on several examples.

Tipping Point

Medline Plus

Full Text Available ... 60 Seconds of Safety (Videos) > The Tipping Point The Tipping Point by CPSC Blogger September 22, 2009 appliance child Childproofing CPSC danger death electrical fall furniture head injury product safety television tipover tv Watch the video in Adobe Flash ...

High Temperature Gas-Cooled Test Reactor Point Design: Summary Report

Energy Technology Data Exchange (ETDEWEB)

Sterbentz, James William [Idaho National Lab. (INL), Idaho Falls, ID (United States); Bayless, Paul David [Idaho National Lab. (INL), Idaho Falls, ID (United States); Nelson, Lee Orville [Idaho National Lab. (INL), Idaho Falls, ID (United States); Gougar, Hans David [Idaho National Lab. (INL), Idaho Falls, ID (United States); Kinsey, J. [Idaho National Lab. (INL), Idaho Falls, ID (United States); Strydom, Gerhard [Idaho National Lab. (INL), Idaho Falls, ID (United States)

2016-03-01

A point design has been developed for a 200-MW high-temperature gas-cooled test reactor. The point design concept uses standard prismatic blocks and 15.5% enriched uranium oxycarbide fuel. Reactor physics and thermal-hydraulics simulations have been performed to characterize the capabilities of the design. In addition to the technical data, overviews are provided on the technology readiness level, licensing approach, and costs of the test reactor point design.

High Temperature Gas-Cooled Test Reactor Point Design: Summary Report

Energy Technology Data Exchange (ETDEWEB)

Sterbentz, James William [Idaho National Lab. (INL), Idaho Falls, ID (United States); Bayless, Paul David [Idaho National Lab. (INL), Idaho Falls, ID (United States); Nelson, Lee Orville [Idaho National Lab. (INL), Idaho Falls, ID (United States); Gougar, Hans David [Idaho National Lab. (INL), Idaho Falls, ID (United States); Strydom, Gerhard [Idaho National Lab. (INL), Idaho Falls, ID (United States)

2016-01-01

A point design has been developed for a 200-MW high-temperature gas-cooled test reactor. The point design concept uses standard prismatic blocks and 15.5% enriched uranium oxycarbide fuel. Reactor physics and thermal-hydraulics simulations have been performed to characterize the capabilities of the design. In addition to the technical data, overviews are provided on the technology readiness level, licensing approach, and costs of the test reactor point design.

Point Genetics: A New Concept to Assess Neutron Kinetics

International Nuclear Information System (INIS)

Klein Meulekamp, R.; Kuijper, J.C.; Schikorr, M.

2005-01-01

Point genetic equations are introduced. These equations are similar to the well-known point kinetic equations but characterize and couple individual fission generations in subcritical systems. Point genetic equations are able to describe dynamic behavior of source-driven subcritical systems on shorter timescales than is possible with point kinetic equations. Point genetic parameters can be used as a first-order characterization of the system and can be calculated using standard Monte Carlo techniques; the implementation in other calculational schemes seems straightforward. A Godiva sphere is considered to show the applicability of the point genetic equations in describing a detector response on short timescales. For this system the point genetic parameters are calculated and compared with reference calculations. Typical dynamic source behavior is considered by studying a transient in which the neutron source energy decreases from 20 to 1 MeV. For all cases studied, the point genetic equations are compared to full space-time kinetic solutions, and it is shown that point genetics performs well

An integral constraint for the evolution of the galaxy two-point correlation function

International Nuclear Information System (INIS)

Peebles, P.J.E.; Groth, E.J.

1976-01-01

Under some conditions an integral over the galaxy two-point correlation function, xi(x,t), evolves with the expansion of the universe in a simple manner easily computed from linear perturbation theory.This provides a useful constraint on the possible evolution of xi(x,t) itself. We test the integral constraint with both an analytic model and numerical N-body simulations for the evolution of irregularities in an expanding universe. Some applications are discussed. (orig.) [de

The succinonitrile triple-point standard: a fixed point to improve the accuracy of temperature measurements in the clinical laboratory.

Science.gov (United States)

Mangum, B W

1983-07-01

In an investigation of the melting and freezing behavior of succinonitrile, the triple-point temperature was determined to be 58.0805 degrees C, with an estimated uncertainty of +/- 0.0015 degrees C relative to the International Practical Temperature Scale of 1968 (IPTS-68). The triple-point temperature of this material is evaluated as a temperature-fixed point, and some clinical laboratory applications of this fixed point are proposed. In conjunction with the gallium and ice points, the availability of succinonitrile permits thermistor thermometers to be calibrated accurately and easily on the IPTS-68.

Anomalous transparency in photonic crystals and its application to point-by-point grating inscription in photonic crystal fibers.

Science.gov (United States)

Baghdasaryan, Tigran; Geernaert, Thomas; Chah, Karima; Caucheteur, Christophe; Schuster, Kay; Kobelke, Jens; Thienpont, Hugo; Berghmans, Francis

2018-04-03

It is common belief that photonic crystals behave similarly to isotropic and transparent media only when their feature sizes are much smaller than the wavelength of light. Here, we counter that belief and we report on photonic crystals that are transparent for anomalously high normalized frequencies up to 0.9, where the crystal's feature sizes are comparable with the free space wavelength. Using traditional photonic band theory, we demonstrate that the isofrequency curves can be circular in the region above the first stop band for triangular lattice photonic crystals. In addition, by simulating how efficiently a tightly focused Gaussian beam propagates through the photonic crystal slab, we judge on the photonic crystal's transparency rather than on isotropy only. Using this approach, we identified a wide range of photonic crystal parameters that provide anomalous transparency. Our findings indicate the possibility to scale up the features of photonic crystals and to extend their operational wavelength range for applications including optical cloaking and graded index guiding. We applied our result in the domain of femtosecond laser micromachining, by demonstrating what we believe to be the first point-by-point grating inscribed in a multi-ring photonic crystal fiber.

Design of four-point SENB specimens with stable crack growth

DEFF Research Database (Denmark)

Jørgensen, Jeppe Bjørn; Kildegaard, Casper; Sørensen, Bent F.

2018-01-01

A four-point single-edge-notch-beam (SENB) test specimen loaded in displacement control (fixed grip) is proposed for studying crack deflection at bi-material interfaces. In order to ensure stable crack growth, a novel analytical model of the four-point SENB specimen in fixed grip is derived...... and compared with numerical models. Model results show that the specimen should be short and thick, and the start-crack length should be deep for the crack to propagate stable towards the bi-material interface. Observations from experimental tests of four-point SENB specimens with different start-crack lengths...

An easy-to-build, low-budget point-of-care ultrasound simulator: from Linux to a web-based solution.

Science.gov (United States)

Damjanovic, Domagoj; Goebel, Ulrich; Fischer, Benedikt; Huth, Martin; Breger, Hartmut; Buerkle, Hartmut; Schmutz, Axel

2017-12-01

Hands-on training in point-of-care ultrasound (POC-US) should ideally comprise bedside teaching, as well as simulated clinical scenarios. High-fidelity phantoms and portable ultrasound simulation systems are commercially available, however, at considerable costs. This limits their suitability for medical schools. A Linux-based software for Emergency Department Ultrasound Simulation (edus2TM) was developed by Kulyk and Olszynski in 2011. Its feasibility for POC-US education has been well-documented, and shows good acceptance. An important limitation to an even more widespread use of edus2, however, may be due to the need for a virtual machine for WINDOWS ® systems. Our aim was to adapt the original software toward an HTML-based solution, thus making it affordable and applicable in any simulation setting. We created an HTML browser-based ultrasound simulation application, which reads the input of different sensors, triggering an ultrasound video to be displayed on a respective device. RFID tags, NFC tags, and QR Codes™ have been integrated into training phantoms or were attached to standardized patients. The RFID antenna was hidden in a mock ultrasound probe. The application is independent from the respective device. Our application was used successfully with different trigger/scanner combinations and mounted readily into simulated training scenarios. The application runs independently from operating systems or electronic devices. This low-cost, browser-based ultrasound simulator is easy-to-build, very adaptive, and independent from operating systems. It has the potential to facilitate POC-US training throughout the world, especially in resource-limited areas.

Numerical simulation of aircraft crash on nuclear containment structure

Energy Technology Data Exchange (ETDEWEB)

Iqbal, M.A., E-mail: iqbalfce@iitr.ernet.in [Department of Civil Engineering, Indian Institute of Technology Roorkee, Roorkee 247667 (India); Rai, S.; Sadique, M.R.; Bhargava, P. [Department of Civil Engineering, Indian Institute of Technology Roorkee, Roorkee 247667 (India)

2012-02-15

Highlights: Black-Right-Pointing-Pointer The deformation was more localised at the center of cylindrical portion. Black-Right-Pointing-Pointer The peak deflection at the junction of dome and cylinder was found to be 67 mm. Black-Right-Pointing-Pointer The peak deflection at midpoint of the cylindrical portion was found to be 88.9 mm. Black-Right-Pointing-Pointer The strain rate was found to be an important parameter to effect the deformation. Black-Right-Pointing-Pointer The model without strain rate and 290 s{sup -1} strain rate predicted very high deformations. - Abstract: Numerical simulations were carried with ABAQUS/Explicit finite element code in order to predict the response of BWR Mark III type nuclear containment against Boeing 707-320 aircraft crash. The load of the aircraft was applied using and force history curve. The damaged plasticity model was used to predict the behavior of concrete while the Johnson-Cook elasto-viscoplastic material model was used to incorporate the behavior of steel reinforcement. The crash was considered to occur at two different locations i.e., the midpoint of the cylindrical portion and the junction of dome and cylinder. The midpoint of the cylindrical portion experienced more deformation. The strain rate in the material model was varied and found to have a significant effect on the response of containment. The results of the present investigation were compared with those of the studies available in literature and a close agreement with the previous results was found in terms of maximum target deformation.

Critical-point nuclei

International Nuclear Information System (INIS)

Clark, R.M.

2004-01-01

It has been suggested that a change of nuclear shape may be described in terms of a phase transition and that specific nuclei may lie close to the critical point of the transition. Analytical descriptions of such critical-point nuclei have been introduced recently and they are described briefly. The results of extensive searches for possible examples of critical-point behavior are presented. Alternative pictures, such as describing bands in the candidate nuclei using simple ΔK = 0 and ΔK = 2 rotational-coupling models, are discussed, and the limitations of the different approaches highlighted. A possible critical-point description of the transition from a vibrational to rotational pairing phase is suggested

Scattering and absorption of particles emitted by a point source in a cluster of point scatterers

International Nuclear Information System (INIS)

Liljequist, D.

2012-01-01

A theory for the scattering and absorption of particles isotropically emitted by a point source in a cluster of point scatterers is described and related to the theory for the scattering of an incident particle beam. The quantum mechanical probability of escape from the cluster in different directions is calculated, as well as the spatial distribution of absorption events within the cluster. A source strength renormalization procedure is required. The average quantum scattering in clusters with randomly shifting scatterer positions is compared to trajectory simulation with the aim of studying the validity of the trajectory method. Differences between the results of the quantum and trajectory methods are found primarily for wavelengths larger than the average distance between nearest neighbour scatterers. The average quantum results include, for example, a local minimum in the number of absorption events at the location of the point source and interference patterns in the angle-dependent escape probability as well as in the distribution of absorption events. The relative error of the trajectory method is in general, though not generally, of similar magnitude as that obtained for beam scattering.

Computer simulations of channel meandering and the formation of point bars: Linking channel dynamics to the preserved stratigraphy

Science.gov (United States)

Sun, T.; Covault, J. A.; Pyrcz, M.; Sullivan, M.

2012-12-01

Meandering rivers are probably one of the most recognizable geomorphic features on earth. As they meander across alluvial and delta plains, channels migrate laterally and develop point bars, splays, levees and other geomorphic and sedimentary features that compose substantial portions of the fill within many sedimentary basins. These basins can include hydrocarbon producing fields. Therefore, a good understanding of the processes of meandering channels and their associated deposits is critical for exploiting these reservoirs in the subsurface. In the past couple of decades, significant progress has been made in our understanding of the morphodynamics of channel meandering. Basic fluid dynamics and sediment transport (Ikeda and Parker, 1981; Howard, 1992) has shown that many characteristic features of meandering rivers, such as the meandering wavelength, growth rate and downstream migration rate, can be predicted quantitatively. As a result, a number of variations and improvement of the theory have emerged (e.g., Blondeaux and Seminara, 1985; Parker and Andrews, 1985, 1986; and Sun et al., 2001a, b).The main improvements include the recognition of so called "bar-bend" interactions, where the development of bars on the channel bed and their interactions with the channel bend is recognized as a primary cause for meandering channels to develop greater complexity than the classic goose-neck meander bend shapes, such as compound bend. Recently, Sun and others have shown that the spatial patterns of width variations in meandering channels can be explained by an extrinsic periodic flow variations coupled with the intrinsic bend instability dynamics. In contrast to the significant improvement of our understanding of channel meandering, little work has been done to link the geomorphic features of meandering channels to the geometry and heterogeneity of the deposits they form and ultimately preserves. A computer simulation model based on the work of Sun and others (1996, 2001

Non-parametric Bayesian inference for inhomogeneous Markov point processes

DEFF Research Database (Denmark)

Berthelsen, Kasper Klitgaard; Møller, Jesper; Johansen, Per Michael

is a shot noise process, and the interaction function for a pair of points depends only on the distance between the two points and is a piecewise linear function modelled by a marked Poisson process. Simulation of the resulting posterior using a Metropolis-Hastings algorithm in the "conventional" way...

IMAGE TO POINT CLOUD METHOD OF 3D-MODELING

Directory of Open Access Journals (Sweden)

A. G. Chibunichev

2012-07-01

Full Text Available This article describes the method of constructing 3D models of objects (buildings, monuments based on digital images and a point cloud obtained by terrestrial laser scanner. The first step is the automated determination of exterior orientation parameters of digital image. We have to find the corresponding points of the image and point cloud to provide this operation. Before the corresponding points searching quasi image of point cloud is generated. After that SIFT algorithm is applied to quasi image and real image. SIFT algorithm allows to find corresponding points. Exterior orientation parameters of image are calculated from corresponding points. The second step is construction of the vector object model. Vectorization is performed by operator of PC in an interactive mode using single image. Spatial coordinates of the model are calculated automatically by cloud points. In addition, there is automatic edge detection with interactive editing available. Edge detection is performed on point cloud and on image with subsequent identification of correct edges. Experimental studies of the method have demonstrated its efficiency in case of building facade modeling.

Genetic Algorithm Based Economic Dispatch with Valve Point Effect

Energy Technology Data Exchange (ETDEWEB)

Park, Jong Nam; Park, Kyung Won; Kim, Ji Hong; Kim, Jin O [Hanyang University (Korea, Republic of)

1999-03-01

This paper presents a new approach on genetic algorithm to economic dispatch problem for valve point discontinuities. Proposed approach in this paper on genetic algorithms improves the performance to solve economic dispatch problem for valve point discontinuities through improved death penalty method, generation-apart elitism, atavism and sexual selection with sexual distinction. Numerical results on a test system consisting of 13 thermal units show that the proposed approach is faster, more robust and powerful than conventional genetic algorithms. (author). 8 refs., 10 figs.

Multiview point clouds denoising based on interference elimination

Science.gov (United States)

Hu, Yang; Wu, Qian; Wang, Le; Jiang, Huanyu

2018-03-01

Newly emerging low-cost depth sensors offer huge potentials for three-dimensional (3-D) modeling, but existing high noise restricts these sensors from obtaining accurate results. Thus, we proposed a method for denoising registered multiview point clouds with high noise to solve that problem. The proposed method is aimed at fully using redundant information to eliminate the interferences among point clouds of different views based on an iterative procedure. In each iteration, noisy points are either deleted or moved to their weighted average targets in accordance with two cases. Simulated data and practical data captured by a Kinect v2 sensor were tested in experiments qualitatively and quantitatively. Results showed that the proposed method can effectively reduce noise and recover local features from highly noisy multiview point clouds with good robustness, compared to truncated signed distance function and moving least squares (MLS). Moreover, the resulting low-noise point clouds can be further smoothed by the MLS to achieve improved results. This study provides the feasibility of obtaining fine 3-D models with high-noise devices, especially for depth sensors, such as Kinect.

The registration of non-cooperative moving targets laser point cloud in different view point

Science.gov (United States)

Wang, Shuai; Sun, Huayan; Guo, Huichao

2018-01-01

Non-cooperative moving target multi-view cloud registration is the key technology of 3D reconstruction of laser threedimension imaging. The main problem is that the density changes greatly and noise exists under different acquisition conditions of point cloud. In this paper, firstly, the feature descriptor is used to find the most similar point cloud, and then based on the registration algorithm of region segmentation, the geometric structure of the point is extracted by the geometric similarity between point and point, The point cloud is divided into regions based on spectral clustering, feature descriptors are created for each region, searching to find the most similar regions in the most similar point of view cloud, and then aligning the pair of point clouds by aligning their minimum bounding boxes. Repeat the above steps again until registration of all point clouds is completed. Experiments show that this method is insensitive to the density of point clouds and performs well on the noise of laser three-dimension imaging.

[Dancing with Pointe Shoes: Characteristics and Assessment Criteria for Pointe Readiness].

Science.gov (United States)

Wanke, Eileen M; Exner-Grave, Elisabeth

2017-12-01

Training with pointe shoes is an integral part of professional dance education and ambitious hobby dancing. Pointe shoes - developed more than hundred years ago and almost unaltered since then - are highly specific and strike a balance between aesthetics, function, protection, and health care. Therefore, pointe readiness should be tested prior to all dance training or career training. Medical specialists are often confronted with this issue. Specific anatomical dance technique-orientated general conditional and coordinative preconditions as well as dance-technical prerequisites must be met by pointe readiness tests in order to keep traumatic injuries or long-term damage at a minimum. In addition to a (training) history, medical counselling sessions have come to include various tests that enable a reliable decision for or against pointe work. This article suggests adequate testing procedures (STT TEST), taking account of professional dancing as well as hobby dancing. © Georg Thieme Verlag KG Stuttgart · New York.

Numerical simulation of a short RFQ resonator using the MAFIA codes

International Nuclear Information System (INIS)

Wang, H.; Ben-Zvi, I.; Jain, A.; Paul, P.; Lombardi, A.

1991-01-01

The electrical characteristics of a short (2βλ=0.4 m) resonator with large modulation (m=4) have been studied using the three dimensional codes, MAFIA. The complete resonator, including the modulated electrodes and a complex support structure, has been simulated using ∼ 350,000 mesh points. Important characteristics studied include the resonant frequency, electric and magnetic fields distributions, quality factor and stored energy. The results of the numerical simulations are compared with the measurements of an actual resonator and analytical approximations. 7 refs., 3 figs., 1 tab

Performances improvement of maximum power point tracking perturb and observe method

Energy Technology Data Exchange (ETDEWEB)

Egiziano, L.; Femia, N.; Granozio, D.; Petrone, G.; Spagnuolo, G. [Salermo Univ., Salermo (Italy); Vitelli, M. [Seconda Univ. di Napoli, Napoli (Italy)

2006-07-01

Perturb and observe best operation conditions were investigated in order to identify edge efficiency performance capabilities of a maximum power point (MPP) tracking technique for photovoltaic (PV) applications. The strategy was developed to ensure a 3-points behavior across the MPP under a fixed irradiation level with a central point blocked on the MPP and 2 operating points operating at voltage values that guaranteed the same power levels. The system was also devised to quickly detect the MPP movement in the presence of varying atmospheric conditions by increasing the perturbation so that the MPP was guaranteed within a few sampling periods. A perturbation equation was selected where amplitude was represented as a function of the actual power drawn from the PV field together with the adoption of a parabolic interpolation of the sequence of the final 3 acquired voltage power couples corresponding to as many operating points. The technique was developed to ensure that the power difference between 2 consecutive operating points was higher than the power quantization error. Simulations were conducted to demonstrate that the proposed technique arranged operating points symmetrically around the MPP. The average power of the 3-points set was achieved by means of the parabolic prediction. Experiments conducted to validate the simulation showed a reduced power oscillation below the MPP and a real power gain. 2 refs., 8 figs.

THE GROWTH POINTS OF STATISTICAL METHODS

OpenAIRE

Orlov A. I.

2014-01-01

On the basis of a new paradigm of applied mathematical statistics, data analysis and economic-mathematical methods are identified; we have also discussed five topical areas in which modern applied statistics is developing as well as the other statistical methods, i.e. five "growth points" – nonparametric statistics, robustness, computer-statistical methods, statistics of interval data, statistics of non-numeric data

Ferrimagnetism and compensation points in a decorated 3D Ising model

International Nuclear Information System (INIS)

Oitmaa, J.; Zheng, W.

2003-01-01

Full text: Ferrimagnets are materials where ions on different sublattices have opposing magnetic moments which do not exactly cancel even at zero temperature. An intriguing possibility then is the existence of a compensation point, below the Curie temperature, where the net moment changes sign. This has obvious technological significance. Most theoretical studies of such systems have used mean-field approaches, making it difficult to distinguish real properties of the model from artefacts of the approximation. For this reason a number of simpler models have been proposed, where treatments beyond mean-field theory are possible. Of particular interest are decorated systems, which can be mapped exactly onto simpler models and, in this way, either solved exactly or to a high degree of numerical precision. We use this approach to study a ferrimagnetic Ising system with spins 1/2 at the sites of a simple cubic lattice and spins S=1 or 3/2 located on the bonds. Our results, which are exact to high numerical precision, show a number of surprising and interesting features: for S=1 the possibility of zero, one or two compensation points, re-entrant behaviour, and up to three critical points; for S=3/2 always a simple critical point and zero or one compensation point

Point Clouds to Indoor/outdoor Accessibility Diagnosis

Science.gov (United States)

Balado, J.; Díaz-Vilariño, L.; Arias, P.; Garrido, I.

2017-09-01

This work presents an approach to automatically detect structural floor elements such as steps or ramps in the immediate environment of buildings, elements that may affect the accessibility to buildings. The methodology is based on Mobile Laser Scanner (MLS) point cloud and trajectory information. First, the street is segmented in stretches along the trajectory of the MLS to work in regular spaces. Next, the lower region of each stretch (the ground zone) is selected as the ROI and normal, curvature and tilt are calculated for each point. With this information, points in the ROI are classified in horizontal, inclined or vertical. Points are refined and grouped in structural elements using raster process and connected components in different phases for each type of previously classified points. At last, the trajectory data is used to distinguish between road and sidewalks. Adjacency information is used to classify structural elements in steps, ramps, curbs and curb-ramps. The methodology is tested in a real case study, consisting of 100 m of an urban street. Ground elements are correctly classified in an acceptable computation time. Steps and ramps also are exported to GIS software to enrich building models from Open Street Map with information about accessible/inaccessible entrances and their locations.

Hyperchaos Numerical Simulation and Control in a 4D Hyperchaotic System

Directory of Open Access Journals (Sweden)

Junhai Ma

2013-01-01

Full Text Available A hyperchaotic system is introduced, and the complex dynamical behaviors of such system are investigated by means of numerical simulations. The bifurcation diagrams, Lyapunov exponents, hyperchaotic attractors, the power spectrums, and time charts are mapped out through the theory analysis and dynamic simulations. The chaotic and hyper-chaotic attractors exist and alter over a wide range of parameters according to the variety of Lyapunov exponents and bifurcation diagrams. Furthermore, linear feedback controllers are designed for stabilizing the hyperchaos to the unstable equilibrium points; thus, we achieve the goal of a second control which is more useful in application.

Comparison of methods for accurate end-point detection of potentiometric titrations

Science.gov (United States)

Villela, R. L. A.; Borges, P. P.; Vyskočil, L.

2015-01-01

Detection of the end point in potentiometric titrations has wide application on experiments that demand very low measurement uncertainties mainly for certifying reference materials. Simulations of experimental coulometric titration data and consequential error analysis of the end-point values were conducted using a programming code. These simulations revealed that the Levenberg-Marquardt method is in general more accurate than the traditional second derivative technique used currently as end-point detection for potentiometric titrations. Performance of the methods will be compared and presented in this paper.

Comparison of methods for accurate end-point detection of potentiometric titrations

International Nuclear Information System (INIS)

Villela, R L A; Borges, P P; Vyskočil, L

2015-01-01

Detection of the end point in potentiometric titrations has wide application on experiments that demand very low measurement uncertainties mainly for certifying reference materials. Simulations of experimental coulometric titration data and consequential error analysis of the end-point values were conducted using a programming code. These simulations revealed that the Levenberg-Marquardt method is in general more accurate than the traditional second derivative technique used currently as end-point detection for potentiometric titrations. Performance of the methods will be compared and presented in this paper

Efficient 3D Volume Reconstruction from a Point Cloud Using a Phase-Field Method

Directory of Open Access Journals (Sweden)

Darae Jeong

2018-01-01

Full Text Available We propose an explicit hybrid numerical method for the efficient 3D volume reconstruction from unorganized point clouds using a phase-field method. The proposed three-dimensional volume reconstruction algorithm is based on the 3D binary image segmentation method. First, we define a narrow band domain embedding the unorganized point cloud and an edge indicating function. Second, we define a good initial phase-field function which speeds up the computation significantly. Third, we use a recently developed explicit hybrid numerical method for solving the three-dimensional image segmentation model to obtain efficient volume reconstruction from point cloud data. In order to demonstrate the practical applicability of the proposed method, we perform various numerical experiments.

Estimating Function Approaches for Spatial Point Processes

Science.gov (United States)

Deng, Chong

Spatial point pattern data consist of locations of events that are often of interest in biological and ecological studies. Such data are commonly viewed as a realization from a stochastic process called spatial point process. To fit a parametric spatial point process model to such data, likelihood-based methods have been widely studied. However, while maximum likelihood estimation is often too computationally intensive for Cox and cluster processes, pairwise likelihood methods such as composite likelihood, Palm likelihood usually suffer from the loss of information due to the ignorance of correlation among pairs. For many types of correlated data other than spatial point processes, when likelihood-based approaches are not desirable, estimating functions have been widely used for model fitting. In this dissertation, we explore the estimating function approaches for fitting spatial point process models. These approaches, which are based on the asymptotic optimal estimating function theories, can be used to incorporate the correlation among data and yield more efficient estimators. We conducted a series of studies to demonstrate that these estmating function approaches are good alternatives to balance the trade-off between computation complexity and estimating efficiency. First, we propose a new estimating procedure that improves the efficiency of pairwise composite likelihood method in estimating clustering parameters. Our approach combines estimating functions derived from pairwise composite likeli-hood estimation and estimating functions that account for correlations among the pairwise contributions. Our method can be used to fit a variety of parametric spatial point process models and can yield more efficient estimators for the clustering parameters than pairwise composite likelihood estimation. We demonstrate its efficacy through a simulation study and an application to the longleaf pine data. Second, we further explore the quasi-likelihood approach on fitting

PowerPoint 2010 Bible

CERN Document Server

Wempen, Faithe

2010-01-01

Master PowerPoint and improve your presentation skills-with one book!. It's no longer enough to have slide after slide of text, bullets, and charts. It's not even enough to have good speaking skills if your PowerPoint slides bore your audience. Get the very most out of all that PowerPoint 2010 has to offer while also learning priceless tips and techniques for making good presentations in this new PowerPoint 2010 Bible. Well-known PowerPoint expert and author Faithe Wempen provides formatting tips; shows you how to work with drawings, tables, and SmartArt; introduces new collaboration tools; wa

Numerical Simulation of Tsunami Hazard Mitigation by Mangrove Forest in North Coast Bali, Indonesia

Directory of Open Access Journals (Sweden)

Putu Harry Gunawan

2015-06-01

Full Text Available Mangrove forest or known as bakau forest is important forest as a natural wave barrier or tsunami wave mitigation. Some advantages of mangrove forest to reduce the water waves are already studied. Mangrove forest in north coast of Bali’s island, Buleleng regency, Indonesia is in damaged condition. The aim of this paper is to present the importance of mangrove forest as the water wave mitigation in numerical simulation point of view. Moreover, the results also show the effect of tsunami propagation to the coastal area with and without mangrove resistance. Here, the nonlinear shallow water equations are used to govern the model of numerical simulation.

Numerical Simulation of Tsunami Hazard Mitigation by Mangrove Forest in North Coast Bali, Indonesia

Directory of Open Access Journals (Sweden)

Putu Harry Gunawan