Extreme trajectory approach to the problem of determining the degree of stability of a laser surface scanning system

NASA Astrophysics Data System (ADS)

Nikitin, V. N.; Chemodanov, V. B.

2018-02-01

The degree of stability of a laser system for surface scanning with nonlinear multiplicative crosstalks is discussed. To determine its stability, the action functional is introduced, which is defined on the set of virtual (achievable) trajectories. The action functional is a measure of external action, which should be applied to a system to move it along a predetermined trial trajectory in the state space.The degree of stability of the system depends on the minimum value of the action functional which is reached on the extreme trajectory transferring the laser scanning system from equilibrium to the limit of the normal operation range. Numerical methods are proposed for calculating the degree of stability.

Weight Multispectral Reconstruction Strategy for Enhanced Reconstruction Accuracy and Stability With Cerenkov Luminescence Tomography.

PubMed

Hongbo Guo; Xiaowei He; Muhan Liu; Zeyu Zhang; Zhenhua Hu; Jie Tian

2017-06-01

Cerenkov luminescence tomography (CLT) provides a novel technique for 3-D noninvasive detection of radiopharmaceuticals in living subjects. However, because of the severe scattering of Cerenkov light, the reconstruction accuracy and stability of CLT is still unsatisfied. In this paper, a modified weight multispectral CLT (wmCLT) reconstruction strategy was developed which split the Cerenkov radiation spectrum into several sub-spectral bands and weighted the sub-spectral results to obtain the final result. To better evaluate the property of the wmCLT reconstruction strategy in terms of accuracy, stability and practicability, several numerical simulation experiments and in vivo experiments were conducted and the results obtained were compared with the traditional multispectral CLT (mCLT) and hybrid-spectral CLT (hCLT) reconstruction strategies. The numerical simulation results indicated that wmCLT strategy significantly improved the accuracy of Cerenkov source localization and intensity quantitation and exhibited good stability in suppressing noise in numerical simulation experiments. And the comparison of the results achieved from different in vivo experiments further indicated significant improvement of the wmCLT strategy in terms of the shape recovery of the bladder and the spatial resolution of imaging xenograft tumors. Overall the strategy reported here will facilitate the development of nuclear and optical molecular tomography in theoretical study.

Development of Theoretical and Numerical Techniques for Achieving Stability in Gyrotron Traveling-Wave Amplifiers.

DTIC Science & Technology

1989-02-01

analysis methods diverge significantly. The electron current density found in Eq. 2.106 may be evaluated" as I J ...S..Y.v Yvt r t) (2.107) 0 ZO where 10...will be specified by the geometry and mode under consider- ation. It was noted earlier that the point of divergence between the two principle...techniques lies in the methods used to calculate the current density. Actually, the divergence is present only in theory. Theoreti- cally and numerically, Eq

Precise attitude rate estimation using star images obtained by mission telescope for satellite missions

NASA Astrophysics Data System (ADS)

Inamori, Takaya; Hosonuma, Takayuki; Ikari, Satoshi; Saisutjarit, Phongsatorn; Sako, Nobutada; Nakasuka, Shinichi

2015-02-01

Recently, small satellites have been employed in various satellite missions such as astronomical observation and remote sensing. During these missions, the attitudes of small satellites should be stabilized to a higher accuracy to obtain accurate science data and images. To achieve precise attitude stabilization, these small satellites should estimate their attitude rate under the strict constraints of mass, space, and cost. This research presents a new method for small satellites to precisely estimate angular rate using star blurred images by employing a mission telescope to achieve precise attitude stabilization. In this method, the angular velocity is estimated by assessing the quality of a star image, based on how blurred it appears to be. Because the proposed method utilizes existing mission devices, a satellite does not require additional precise rate sensors, which makes it easier to achieve precise stabilization given the strict constraints possessed by small satellites. The research studied the relationship between estimation accuracy and parameters used to achieve an attitude rate estimation, which has a precision greater than 1 × 10-6 rad/s. The method can be applied to all attitude sensors, which use optics systems such as sun sensors and star trackers (STTs). Finally, the method is applied to the nano astrometry satellite Nano-JASMINE, and we investigate the problems that are expected to arise with real small satellites by performing numerical simulations.

influx_s: increasing numerical stability and precision for metabolic flux analysis in isotope labelling experiments.

PubMed

Sokol, Serguei; Millard, Pierre; Portais, Jean-Charles

2012-03-01

The problem of stationary metabolic flux analysis based on isotope labelling experiments first appeared in the early 1950s and was basically solved in early 2000s. Several algorithms and software packages are available for this problem. However, the generic stochastic algorithms (simulated annealing or evolution algorithms) currently used in these software require a lot of time to achieve acceptable precision. For deterministic algorithms, a common drawback is the lack of convergence stability for ill-conditioned systems or when started from a random point. In this article, we present a new deterministic algorithm with significantly increased numerical stability and accuracy of flux estimation compared with commonly used algorithms. It requires relatively short CPU time (from several seconds to several minutes with a standard PC architecture) to estimate fluxes in the central carbon metabolism network of Escherichia coli. The software package influx_s implementing this algorithm is distributed under an OpenSource licence at http://metasys.insa-toulouse.fr/software/influx/. Supplementary data are available at Bioinformatics online.

Wake Numerical Simulation Based on the Park-Gauss Model and Considering Atmospheric Stability

NASA Astrophysics Data System (ADS)

Yang, Xiangsheng; Zhao, Ning; Tian, Linlin; Zhu, Jun

2016-06-01

In this paper, a new Park-Gauss model based on the assumption of the Park model and the Eddy-viscosity model is investigated to conduct the wake numerical simulation for solving a single wind turbine problem. The initial wake radius has been modified to improve the model’s numerical accuracy. Then the impact of the atmospheric stability based on the Park-Gauss model has been studied in the wake region. By the comparisons and the analyses of the test results, it turns out that the new Park-Gauss model could achieve better effects of the wind velocity simulation in the wake region. The wind velocity in the wake region recovers quickly under the unstable atmospheric condition provided the wind velocity is closest to the test result, and recovers slowly under stable atmospheric condition in case of the wind velocity is lower than the test result. Meanwhile, the wind velocity recovery falls in between the unstable and stable neutral atmospheric conditions.

Investigation of the burning configuration of a coaxial injector in a combustion chamber

NASA Technical Reports Server (NTRS)

Ohara, J.

1978-01-01

An analytical investigation was made into the stability of the burning configuration of a single coaxial injector surrounded by similar injectors. The stability criteria was based on an average pressure difference along the boundaries of the adjacent stream tubes as calculated using Spaulding's numerical method. The results indicate qualitatively that there is a tendency for the injectors to have different burning configurations. It is believed that the configuration achieved is random, however once the burning configuration is established, it is believed to persist.

Finite Differences and Collocation Methods for the Solution of the Two Dimensional Heat Equation

NASA Technical Reports Server (NTRS)

Kouatchou, Jules

1999-01-01

In this paper we combine finite difference approximations (for spatial derivatives) and collocation techniques (for the time component) to numerically solve the two dimensional heat equation. We employ respectively a second-order and a fourth-order schemes for the spatial derivatives and the discretization method gives rise to a linear system of equations. We show that the matrix of the system is non-singular. Numerical experiments carried out on serial computers, show the unconditional stability of the proposed method and the high accuracy achieved by the fourth-order scheme.

A 3-D enlarged cell technique (ECT) for elastic wave modelling of a curved free surface

NASA Astrophysics Data System (ADS)

Wei, Songlin; Zhou, Jianyang; Zhuang, Mingwei; Liu, Qing Huo

2016-09-01

The conventional finite-difference time-domain (FDTD) method for elastic waves suffers from the staircasing error when applied to model a curved free surface because of its structured grid. In this work, an improved, stable and accurate 3-D FDTD method for elastic wave modelling on a curved free surface is developed based on the finite volume method and enlarged cell technique (ECT). To achieve a sufficiently accurate implementation, a finite volume scheme is applied to the curved free surface to remove the staircasing error; in the mean time, to achieve the same stability as the FDTD method without reducing the time step increment, the ECT is introduced to preserve the solution stability by enlarging small irregular cells into adjacent cells under the condition of conservation of force. This method is verified by several 3-D numerical examples. Results show that the method is stable at the Courant stability limit for a regular FDTD grid, and has much higher accuracy than the conventional FDTD method.

Scale-Limited Lagrange Stability and Finite-Time Synchronization for Memristive Recurrent Neural Networks on Time Scales.

PubMed

Xiao, Qiang; Zeng, Zhigang

2017-10-01

The existed results of Lagrange stability and finite-time synchronization for memristive recurrent neural networks (MRNNs) are scale-free on time evolvement, and some restrictions appear naturally. In this paper, two novel scale-limited comparison principles are established by means of inequality techniques and induction principle on time scales. Then the results concerning Lagrange stability and global finite-time synchronization of MRNNs on time scales are obtained. Scaled-limited Lagrange stability criteria are derived, in detail, via nonsmooth analysis and theory of time scales. Moreover, novel criteria for achieving the global finite-time synchronization are acquired. In addition, the derived method can also be used to study global finite-time stabilization. The proposed results extend or improve the existed ones in the literatures. Two numerical examples are chosen to show the effectiveness of the obtained results.

Design of Rock Slope Reinforcement: An Himalayan Case Study

NASA Astrophysics Data System (ADS)

Tiwari, Gaurav; Latha, Gali Madhavi

2016-06-01

The stability analysis of the two abutment slopes of a railway bridge proposed at about 359 m above the ground level, crossing a river and connecting two hill faces in the Himalayas, India, is presented. The bridge is located in a zone of high seismic activity. The rock slopes are composed of a heavily jointed rock mass and the spacing, dip and dip direction of joint sets are varying at different locations. Geological mapping was carried out to characterize all discontinuities present along the slopes. Laboratory and field investigations were conducted to assess the geotechnical properties of the intact rock, rock mass and joint infill. Stability analyses of these rock slopes were carried out using numerical programmes. Loads from the foundations resting on the slopes and seismic accelerations estimated from site-specific ground response analysis were considered. The proposed slope profile with several berms between successive foundations was simulated in the numerical model. An equivalent continuum approach with Hoek and Brown failure criterion was initially used in a finite element model to assess the global stability of the slope abutments. In the second stage, finite element analysis of rock slopes with all joint sets with their orientations, spacing and properties explicitly incorporated into the numerical model was taken up using continuum with joints approach. It was observed that the continuum with joints approach was able to capture the local failures in some of the slope sections, which were verified using wedge failure analysis and stereographic projections. Based on the slope deformations and failure patterns observed from the numerical analyses, rock anchors were designed to achieve the target factors of safety against failure while keeping the deformations within the permissible limits. Detailed design of rock anchors and comparison of the stability of slopes with and without reinforcement are presented.

Discrete-time systems with random switches: From systems stability to networks synchronization.

PubMed

Guo, Yao; Lin, Wei; Ho, Daniel W C

2016-03-01

In this article, we develop some approaches, which enable us to more accurately and analytically identify the essential patterns that guarantee the almost sure stability of discrete-time systems with random switches. We allow for the case that the elements in the switching connection matrix even obey some unbounded and continuous-valued distributions. In addition to the almost sure stability, we further investigate the almost sure synchronization in complex dynamical networks consisting of randomly connected nodes. Numerical examples illustrate that a chaotic dynamics in the synchronization manifold is preserved when statistical parameters enter some almost sure synchronization region established by the developed approach. Moreover, some delicate configurations are considered on probability space for ensuring synchronization in networks whose nodes are described by nonlinear maps. Both theoretical and numerical results on synchronization are presented by setting only a few random connections in each switch duration. More interestingly, we analytically find it possible to achieve almost sure synchronization in the randomly switching complex networks even with very large population sizes, which cannot be easily realized in non-switching but deterministically connected networks.

Discrete-time systems with random switches: From systems stability to networks synchronization

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

Guo, Yao; Lin, Wei, E-mail: wlin@fudan.edu.cn; Shanghai Key Laboratory of Contemporary Applied Mathematics, LMNS, and Shanghai Center for Mathematical Sciences, Shanghai 200433

2016-03-15

In this article, we develop some approaches, which enable us to more accurately and analytically identify the essential patterns that guarantee the almost sure stability of discrete-time systems with random switches. We allow for the case that the elements in the switching connection matrix even obey some unbounded and continuous-valued distributions. In addition to the almost sure stability, we further investigate the almost sure synchronization in complex dynamical networks consisting of randomly connected nodes. Numerical examples illustrate that a chaotic dynamics in the synchronization manifold is preserved when statistical parameters enter some almost sure synchronization region established by the developedmore » approach. Moreover, some delicate configurations are considered on probability space for ensuring synchronization in networks whose nodes are described by nonlinear maps. Both theoretical and numerical results on synchronization are presented by setting only a few random connections in each switch duration. More interestingly, we analytically find it possible to achieve almost sure synchronization in the randomly switching complex networks even with very large population sizes, which cannot be easily realized in non-switching but deterministically connected networks.« less

Tilt angle measurement with a Gaussian-shaped laser beam tracking

NASA Astrophysics Data System (ADS)

Šarbort, Martin; Řeřucha, Šimon; Jedlička, Petr; Lazar, Josef; Číp, Ondrej

2014-05-01

We have addressed the challenge to carry out the angular tilt stabilization of a laser guiding mirror which is intended to route a laser beam with a high energy density. Such an application requires good angular accuracy as well as large operating range, long term stability and absolute positioning. We have designed an instrument for such a high precision angular tilt measurement based on a triangulation method where a laser beam with Gaussian profile is reflected off the stabilized mirror and detected by an image sensor. As the angular deflection of the mirror causes a change of the beam spot position, the principal task is to measure the position on the image chip surface. We have employed a numerical analysis of the Gaussian intensity pattern which uses the nonlinear regression algorithm. The feasibility and performance of the method were tested by numeric modeling as well as experimentally. The experimental results indicate that the assembled instrument achieves a measurement error of 0.13 microradian in the range +/-0.65 degrees over the period of one hour. This corresponds to the dynamic range of 1:170 000.

A second-order cell-centered Lagrangian ADER-MOOD finite volume scheme on multidimensional unstructured meshes for hydrodynamics

NASA Astrophysics Data System (ADS)

Boscheri, Walter; Dumbser, Michael; Loubère, Raphaël; Maire, Pierre-Henri

2018-04-01

In this paper we develop a conservative cell-centered Lagrangian finite volume scheme for the solution of the hydrodynamics equations on unstructured multidimensional grids. The method is derived from the Eucclhyd scheme discussed in [47,43,45]. It is second-order accurate in space and is combined with the a posteriori Multidimensional Optimal Order Detection (MOOD) limiting strategy to ensure robustness and stability at shock waves. Second-order of accuracy in time is achieved via the ADER (Arbitrary high order schemes using DERivatives) approach. A large set of numerical test cases is proposed to assess the ability of the method to achieve effective second order of accuracy on smooth flows, maintaining an essentially non-oscillatory behavior on discontinuous profiles, general robustness ensuring physical admissibility of the numerical solution, and precision where appropriate.

New exponential stability criteria for stochastic BAM neural networks with impulses

NASA Astrophysics Data System (ADS)

Sakthivel, R.; Samidurai, R.; Anthoni, S. M.

2010-10-01

In this paper, we study the global exponential stability of time-delayed stochastic bidirectional associative memory neural networks with impulses and Markovian jumping parameters. A generalized activation function is considered, and traditional assumptions on the boundedness, monotony and differentiability of activation functions are removed. We obtain a new set of sufficient conditions in terms of linear matrix inequalities, which ensures the global exponential stability of the unique equilibrium point for stochastic BAM neural networks with impulses. The Lyapunov function method with the Itô differential rule is employed for achieving the required result. Moreover, a numerical example is provided to show that the proposed result improves the allowable upper bound of delays over some existing results in the literature.

DYNAMIC STABILITY OF THE SOLAR SYSTEM: STATISTICALLY INCONCLUSIVE RESULTS FROM ENSEMBLE INTEGRATIONS

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

Zeebe, Richard E., E-mail: zeebe@soest.hawaii.edu

Due to the chaotic nature of the solar system, the question of its long-term stability can only be answered in a statistical sense, for instance, based on numerical ensemble integrations of nearby orbits. Destabilization of the inner planets, leading to close encounters and/or collisions can be initiated through a large increase in Mercury's eccentricity, with a currently assumed likelihood of ∼1%. However, little is known at present about the robustness of this number. Here I report ensemble integrations of the full equations of motion of the eight planets and Pluto over 5 Gyr, including contributions from general relativity. The resultsmore » show that different numerical algorithms lead to statistically different results for the evolution of Mercury's eccentricity (e{sub M}). For instance, starting at present initial conditions (e{sub M}≃0.21), Mercury's maximum eccentricity achieved over 5 Gyr is, on average, significantly higher in symplectic ensemble integrations using heliocentric rather than Jacobi coordinates and stricter error control. In contrast, starting at a possible future configuration (e{sub M}≃0.53), Mercury's maximum eccentricity achieved over the subsequent 500 Myr is, on average, significantly lower using heliocentric rather than Jacobi coordinates. For example, the probability for e{sub M} to increase beyond 0.53 over 500 Myr is >90% (Jacobi) versus only 40%-55% (heliocentric). This poses a dilemma because the physical evolution of the real system—and its probabilistic behavior—cannot depend on the coordinate system or the numerical algorithm chosen to describe it. Some tests of the numerical algorithms suggest that symplectic integrators using heliocentric coordinates underestimate the odds for destabilization of Mercury's orbit at high initial e{sub M}.« less

High-Order Implicit-Explicit Multi-Block Time-stepping Method for Hyperbolic PDEs

NASA Technical Reports Server (NTRS)

Nielsen, Tanner B.; Carpenter, Mark H.; Fisher, Travis C.; Frankel, Steven H.

2014-01-01

This work seeks to explore and improve the current time-stepping schemes used in computational fluid dynamics (CFD) in order to reduce overall computational time. A high-order scheme has been developed using a combination of implicit and explicit (IMEX) time-stepping Runge-Kutta (RK) schemes which increases numerical stability with respect to the time step size, resulting in decreased computational time. The IMEX scheme alone does not yield the desired increase in numerical stability, but when used in conjunction with an overlapping partitioned (multi-block) domain significant increase in stability is observed. To show this, the Overlapping-Partition IMEX (OP IMEX) scheme is applied to both one-dimensional (1D) and two-dimensional (2D) problems, the nonlinear viscous Burger's equation and 2D advection equation, respectively. The method uses two different summation by parts (SBP) derivative approximations, second-order and fourth-order accurate. The Dirichlet boundary conditions are imposed using the Simultaneous Approximation Term (SAT) penalty method. The 6-stage additive Runge-Kutta IMEX time integration schemes are fourth-order accurate in time. An increase in numerical stability 65 times greater than the fully explicit scheme is demonstrated to be achievable with the OP IMEX method applied to 1D Burger's equation. Results from the 2D, purely convective, advection equation show stability increases on the order of 10 times the explicit scheme using the OP IMEX method. Also, the domain partitioning method in this work shows potential for breaking the computational domain into manageable sizes such that implicit solutions for full three-dimensional CFD simulations can be computed using direct solving methods rather than the standard iterative methods currently used.

On the stability analysis of a pair of van der Pol oscillators with delayed self-connection, position and velocity couplings

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

Hu, Kun; Department of Mathematics, City University of Hong Kong, 83 Tat Chee Avenue, Kowloon; Chung, Kwok-wai, E-mail: makchung@cityu.edu.hk

2013-11-15

In this paper, we perform a stability analysis of a pair of van der Pol oscillators with delayed self-connection, position and velocity couplings. Bifurcation diagram of the damping, position and velocity coupling strengths is constructed, which gives insight into how stability boundary curves come into existence and how these curves evolve from small closed loops into open-ended curves. The van der Pol oscillator has been considered by many researchers as the nodes for various networks. It is inherently unstable at the zero equilibrium. Stability control of a network is always an important problem. Currently, the stabilization of the zero equilibriummore » of a pair of van der Pol oscillators can be achieved only for small damping strength by using delayed velocity coupling. An interesting question arises naturally: can the zero equilibrium be stabilized for an arbitrarily large value of the damping strength? We prove that it can be. In addition, a simple condition is given on how to choose the feedback parameters to achieve such goal. We further investigate how the in-phase mode or the out-of-phase mode of a periodic solution is related to the stability boundary curve that it emerges from a Hopf bifurcation. Analytical expression of a periodic solution is derived using an integration method. Some illustrative examples show that the theoretical prediction and numerical simulation are in good agreement.« less

Analysis of torque transmitting behavior and wheel slip prevention control during regenerative braking for high speed EMU trains

NASA Astrophysics Data System (ADS)

Xu, Kun; Xu, Guo-Qing; Zheng, Chun-Hua

2016-04-01

The wheel-rail adhesion control for regenerative braking systems of high speed electric multiple unit trains is crucial to maintaining the stability, improving the adhesion utilization, and achieving deep energy recovery. There remain technical challenges mainly because of the nonlinear, uncertain, and varying features of wheel-rail contact conditions. This research analyzes the torque transmitting behavior during regenerative braking, and proposes a novel methodology to detect the wheel-rail adhesion stability. Then, applications to the wheel slip prevention during braking are investigated, and the optimal slip ratio control scheme is proposed, which is based on a novel optimal reference generation of the slip ratio and a robust sliding mode control. The proposed methodology achieves the optimal braking performance without the wheel-rail contact information. Numerical simulation results for uncertain slippery rails verify the effectiveness of the proposed methodology.

A novel double-convection chaotic attractor, its adaptive control and circuit simulation

NASA Astrophysics Data System (ADS)

Mamat, M.; Vaidyanathan, S.; Sambas, A.; Mujiarto; Sanjaya, W. S. M.; Subiyanto

2018-03-01

A 3-D novel double-convection chaotic system with three nonlinearities is proposed in this research work. The dynamical properties of the new chaotic system are described in terms of phase portraits, Lyapunov exponents, Kaplan-Yorke dimension, dissipativity, stability analysis of equilibria, etc. Adaptive control and synchronization of the new chaotic system with unknown parameters are achieved via nonlinear controllers and the results are established using Lyapunov stability theory. Furthermore, an electronic circuit realization of the new 3-D novel chaotic system is presented in detail. Finally, the circuit experimental results of the 3-D novel chaotic attractor show agreement with the numerical simulations.

Streakline-based closed-loop control of a bluff body flow

NASA Astrophysics Data System (ADS)

Roca, Pablo; Cammilleri, Ada; Duriez, Thomas; Mathelin, Lionel; Artana, Guillermo

2014-04-01

A novel closed-loop control methodology is introduced to stabilize a cylinder wake flow based on images of streaklines. Passive scalar tracers are injected upstream the cylinder and their concentration is monitored downstream at certain image sectors of the wake. An AutoRegressive with eXogenous inputs mathematical model is built from these images and a Generalized Predictive Controller algorithm is used to compute the actuation required to stabilize the wake by adding momentum tangentially to the cylinder wall through plasma actuators. The methodology is new and has real-world applications. It is demonstrated on a numerical simulation and the provided results show that good performances are achieved.

An Energy Decaying Scheme for Nonlinear Dynamics of Shells

NASA Technical Reports Server (NTRS)

Bottasso, Carlo L.; Bauchau, Olivier A.; Choi, Jou-Young; Bushnell, Dennis M. (Technical Monitor)

2000-01-01

A novel integration scheme for nonlinear dynamics of geometrically exact shells is developed based on the inextensible director assumption. The new algorithm is designed so as to imply the strict decay of the system total mechanical energy at each time step, and consequently unconditional stability is achieved in the nonlinear regime. Furthermore, the scheme features tunable high frequency numerical damping and it is therefore stiffly accurate. The method is tested for a finite element spatial formulation of shells based on mixed interpolations of strain tensorial components and on a two-parameter representation of director rotations. The robustness of the, scheme is illustrated with the help of numerical examples.

Current Status of Astrometry Satellite missions in Japan: JASMINE project series

NASA Astrophysics Data System (ADS)

Yano, T.; Gouda, N.; Kobayashi, Y.; Tsujimoto, T.; Hatsutori, Y.; Murooka, J.; Niwa, Y.; Yamada, Y.

Astrometry satellites have common technological issues. (A) Astrometry satellites are required to measure the positions of stars with high accuracy from the huge amount of data during the observational period. (B) The high stabilization of the thermal environment in the telescope is required. (C) The attitude-pointing stability of these satellites with sub-pixel accuracy is also required. Measurement of the positions of stars from a huge amount of data is the essence of astrometry. It is needed to exclude the systematic errors adequately for each image of stars in order to obtain the accurate positions. We have carried out a centroiding experiment for determining the positions of stars from about 10 000 image data. The following two points are important issues for the mission system of JASMINE in order to achieve our aim. For the small-JASMINE, we require the thermal stabilization of the telescope in order to obtain high astrometric accuracy of about 10 micro-arcsec. In order to accomplish a measurement of positions of stars with high accuracy, we must make a model of the distortion of the image on the focal plane with the accuracy of less than 0.1 nm. We have investigated numerically that the above requirement is achieved if the thermal variation is within about 1 K / 0.75 h. We also require the accuracy of the attitude-pointing stability of about 200 mas / 7 s. The utilization of the Tip-tilt mirror will make it possible to achieve such a stable pointing.

A mimetic finite difference method for the Stokes problem with elected edge bubbles

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

Lipnikov, K; Berirao, L

2009-01-01

A new mimetic finite difference method for the Stokes problem is proposed and analyzed. The unstable P{sub 1}-P{sub 0} discretization is stabilized by adding a small number of bubble functions to selected mesh edges. A simple strategy for selecting such edges is proposed and verified with numerical experiments. The discretizations schemes for Stokes and Navier-Stokes equations must satisfy the celebrated inf-sup (or the LBB) stability condition. The stability condition implies a balance between discrete spaces for velocity and pressure. In finite elements, this balance is frequently achieved by adding bubble functions to the velocity space. The goal of this articlemore » is to show that the stabilizing edge bubble functions can be added only to a small set of mesh edges. This results in a smaller algebraic system and potentially in a faster calculations. We employ the mimetic finite difference (MFD) discretization technique that works for general polyhedral meshes and can accomodate non-uniform distribution of stabilizing bubbles.« less

A study of numerical methods for hyperbolic conservation laws with stiff source terms

NASA Technical Reports Server (NTRS)

Leveque, R. J.; Yee, H. C.

1988-01-01

The proper modeling of nonequilibrium gas dynamics is required in certain regimes of hypersonic flow. For inviscid flow this gives a system of conservation laws coupled with source terms representing the chemistry. Often a wide range of time scales is present in the problem, leading to numerical difficulties as in stiff systems of ordinary differential equations. Stability can be achieved by using implicit methods, but other numerical difficulties are observed. The behavior of typical numerical methods on a simple advection equation with a parameter-dependent source term was studied. Two approaches to incorporate the source term were utilized: MacCormack type predictor-corrector methods with flux limiters, and splitting methods in which the fluid dynamics and chemistry are handled in separate steps. Various comparisons over a wide range of parameter values were made. In the stiff case where the solution contains discontinuities, incorrect numerical propagation speeds are observed with all of the methods considered. This phenomenon is studied and explained.

An approach toward the numerical evaluation of multi-loop Feynman diagrams

NASA Astrophysics Data System (ADS)

Passarino, Giampiero

2001-12-01

A scheme for systematically achieving accurate numerical evaluation of multi-loop Feynman diagrams is developed. This shows the feasibility of a project aimed to produce a complete calculation for two-loop predictions in the Standard Model. As a first step an algorithm, proposed by F.V. Tkachov and based on the so-called generalized Bernstein functional relation, is applied to one-loop multi-leg diagrams with particular emphasis to the presence of infrared singularities, to the problem of tensorial reduction and to the classification of all singularities of a given diagram. Successively, the extension of the algorithm to two-loop diagrams is examined. The proposed solution consists in applying the functional relation to the one-loop sub-diagram which has the largest number of internal lines. In this way the integrand can be made smooth, a part from a factor which is a polynomial in xS, the vector of Feynman parameters needed for the complementary sub-diagram with the smallest number of internal lines. Since the procedure does not introduce new singularities one can distort the xS-integration hyper-contour into the complex hyper-plane, thus achieving numerical stability. The algorithm is then modified to deal with numerical evaluation around normal thresholds. Concise and practical formulas are assembled and presented, numerical results and comparisons with the available literature are shown and discussed for the so-called sunset topology.

Impact of aerobic stabilization on the characteristics of treatment sludge in the leather tanning industry.

PubMed

Cokgor, Emine Ubay; Aydinli, Ebru; Tas, Didem Okutman; Zengin, Gulsum Emel; Orhon, Derin

2014-01-01

The efficiency of aerobic stabilization on the treatment sludge generated from the leather industry was investigated to meet the expected characteristics and conditions of sludge prior to landfill. The sludge types subjected to aerobic stabilization were chemical treatment sludge, biological excess sludge, and the mixture of both chemical and biological sludges. At the end of 23 days of stabilization, suspended solids, volatile suspended solids and total organic carbon removal efficiencies were determined as 17%, 19% and 23% for biological sludge 31%, 35% and 54% for chemical sludge, and 32%, 34% and 63% for the mixture of both chemical and biological sludges, respectively. Model simulations of the respirometric oxygen uptake rate measurements showed that the ratio of active biomass remained the same at the end of the stabilization for all the sludge samples. Although mixing the chemical and biological sludges resulted in a relatively effective organic carbon and solids removal, the level of stabilization achieved remained clearly below the required level of organic carbon content for landfill. These findings indicate the potential risk of setting numerical restrictions without referring to proper scientific support.

Polyhedral meshing in numerical analysis of conjugate heat transfer

NASA Astrophysics Data System (ADS)

Sosnowski, Marcin; Krzywanski, Jaroslaw; Grabowska, Karolina; Gnatowska, Renata

2018-06-01

Computational methods have been widely applied in conjugate heat transfer analysis. The very first and crucial step in such research is the meshing process which consists in dividing the analysed geometry into numerous small control volumes (cells). In Computational Fluid Dynamics (CFD) applications it is desirable to use the hexahedral cells as the resulting mesh is characterized by low numerical diffusion. Unfortunately generating such mesh can be a very time-consuming task and in case of complicated geometry - it may not be possible to generate cells of good quality. Therefore tetrahedral cells have been implemented into commercial pre-processors. Their advantage is the ease of its generation even in case of very complex geometry. On the other hand tetrahedrons cannot be stretched excessively without decreasing the mesh quality factor, so significantly larger number of cells has to be used in comparison to hexahedral mesh in order to achieve a reasonable accuracy. Moreover the numerical diffusion of tetrahedral elements is significantly higher. Therefore the polyhedral cells are proposed within the paper in order to combine the advantages of hexahedrons (low numerical diffusion resulting in accurate solution) and tetrahedrons (rapid semi-automatic generation) as well as to overcome the disadvantages of both the above mentioned mesh types. The major benefit of polyhedral mesh is that each individual cell has many neighbours, so gradients can be well approximated. Polyhedrons are also less sensitive to stretching than tetrahedrons which results in better mesh quality leading to improved numerical stability of the model. In addition, numerical diffusion is reduced due to mass exchange over numerous faces. This leads to a more accurate solution achieved with a lower cell count. Therefore detailed comparison of numerical modelling results concerning conjugate heat transfer using tetrahedral and polyhedral meshes is presented in the paper.

Fault stability under conditions of variable normal stress

USGS Publications Warehouse

Dieterich, J.H.; Linker, M.F.

1992-01-01

The stability of fault slip under conditions of varying normal stress is modelled as a spring and slider system with rate- and state-dependent friction. Coupling of normal stress to shear stress is achieved by inclining the spring at an angle, ??, to the sliding surface. Linear analysis yields two conditions for unstable slip. The first, of a type previously identified for constant normal stress systems, results in instability if stiffness is below a critical value. Critical stiffness depends on normal stress, constitutive parameters, characteristic sliding distance and the spring angle. Instability of the first type is possible only for velocity-weakening friction. The second condition yields instability if spring angle ?? <-cot-1??ss, where ??ss is steady-state sliding friction. The second condition can arise under conditions of velocity strengthening or weakening. Stability fields for finite perturbations are investigated by numerical simulation. -Authors

Observer enhanced control for spin-stabilized tethered formation in earth orbit

NASA Astrophysics Data System (ADS)

Guang, Zhai; Yuyang, Li; Liang, Bin

2018-04-01

This paper addresses the issues relevant to control of spin-stabilized tethered formation in circular orbit. Due to the dynamic complexities and nonlinear perturbations, it is challenging to promote the control precision for the formation deployment and maintenance. In this work, the formation dynamics are derived with considering the spinning rate of the central body, then major attention is dedicated to develop the nonlinear disturbance observer. To achieve better control performance, the observer-enhanced controller is designed by incorporating the disturbance observer into the control loop, benefits from the disturbance compensation are demonstrated, and also, the dependences of the disturbance observer performance on some important parameters are theoretically and numerically analyzed.

A novel single thruster control strategy for spacecraft attitude stabilization

NASA Astrophysics Data System (ADS)

Godard; Kumar, Krishna Dev; Zou, An-Min

2013-05-01

Feasibility of achieving three axis attitude stabilization using a single thruster is explored in this paper. Torques are generated using a thruster orientation mechanism with which the thrust vector can be tilted on a two axis gimbal. A robust nonlinear control scheme is developed based on the nonlinear kinematic and dynamic equations of motion of a rigid body spacecraft in the presence of gravity gradient torque and external disturbances. The spacecraft, controlled using the proposed concept, constitutes an underactuated system (a system with fewer independent control inputs than degrees of freedom) with nonlinear dynamics. Moreover, using thruster gimbal angles as control inputs make the system non-affine (control terms appear nonlinearly in the state equation). This necessitates the control algorithms to be developed based on nonlinear control theory since linear control methods are not directly applicable. The stability conditions for the spacecraft attitude motion for robustness against uncertainties and disturbances are derived to establish the regions of asymptotic 3-axis attitude stabilization. Several numerical simulations are presented to demonstrate the efficacy of the proposed controller and validate the theoretical results. The control algorithm is shown to compensate for time-varying external disturbances including solar radiation pressure, aerodynamic forces, and magnetic disturbances; and uncertainties in the spacecraft inertia parameters. The numerical results also establish the robustness of the proposed control scheme to negate disturbances caused by orbit eccentricity.

Effect of Swirl on Turbulent Structures in Supersonic Jets

NASA Technical Reports Server (NTRS)

Rao, Ram Mohan; Lundgren, Thomas S.

1998-01-01

Direct Numerical Simulation (DNS) is used to study the mechanism of generation and evolution of turbulence structures in a temporally evolving supersonic swirling round jet and also to examine the resulting acoustic radiations. Fourier spectral expansions are used in the streamwise and azimuthal directions and a 1-D b-spline Galerkin representation is used in the radial direction. Spectral-like accuracy is achieved using this numerical scheme. Direct numerical simulations, using the b-spline spectral method, are carried out starting from mean flow initial conditions which are perturbed by the most unstable linear stability eigenfunctions. It is observed that the initial helical instability waves evolve into helical vortices which eventually breakdown into smaller scales of turbulence. 'Rib' structures similar to those seen in incompressible mixing layer flow of Rogers and Moserl are observed. The jet core breakdown stage exhibits increased acoustic radiations.

Numerical investigations of passive scalar transport in Taylor-Couette flows: Counter-rotation effect

NASA Astrophysics Data System (ADS)

Ouazib, Nabila; Salhi, Yacine; Si-Ahmed, El-Khider; Legrand, Jack; Degrez, G.

2017-07-01

Numerical methods for solving convection-diffusion-reaction (CDR) scalar transport equation in three-dimensional flow are used in the present investigation. The flow is confined between two concentric cylinders both the inner cylinder and the outer one are allowed to rotate. Direct numerical simulations (DNS) have been achieved to study the effects of the gravitational and the centrifugal potentials on the stability of incompressible Taylor-Couette flow. The Navier-Stokes equations and the uncoupled convection-diffusion-reaction equation are solved using a spectral development in one direction combined together with a finite element discretization in the two remaining directions. The complexity of the patterns is highlighted. Since, it increases as the rotation rates of the cylinders increase. In addition, the effect of the counter-rotation of the cylinders on the mass transfer is pointed out.

Tempest - Efficient Computation of Atmospheric Flows Using High-Order Local Discretization Methods

NASA Astrophysics Data System (ADS)

Ullrich, P. A.; Guerra, J. E.

2014-12-01

The Tempest Framework composes several compact numerical methods to easily facilitate intercomparison of atmospheric flow calculations on the sphere and in rectangular domains. This framework includes the implementations of Spectral Elements, Discontinuous Galerkin, Flux Reconstruction, and Hybrid Finite Element methods with the goal of achieving optimal accuracy in the solution of atmospheric problems. Several advantages of this approach are discussed such as: improved pressure gradient calculation, numerical stability by vertical/horizontal splitting, arbitrary order of accuracy, etc. The local numerical discretization allows for high performance parallel computation and efficient inclusion of parameterizations. These techniques are used in conjunction with a non-conformal, locally refined, cubed-sphere grid for global simulations and standard Cartesian grids for simulations at the mesoscale. A complete implementation of the methods described is demonstrated in a non-hydrostatic setting.

Coordination of fractional-order nonlinear multi-agent systems via distributed impulsive control

NASA Astrophysics Data System (ADS)

Ma, Tiedong; Li, Teng; Cui, Bing

2018-01-01

The coordination of fractional-order nonlinear multi-agent systems via distributed impulsive control method is studied in this paper. Based on the theory of impulsive differential equations, algebraic graph theory, Lyapunov stability theory and Mittag-Leffler function, two novel sufficient conditions for achieving the cooperative control of a class of fractional-order nonlinear multi-agent systems are derived. Finally, two numerical simulations are verified to illustrate the effectiveness and feasibility of the proposed method.

Influence of thermal agitation on the electric field induced precessional magnetization reversal with perpendicular easy axis

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

Cheng, Hongguang, E-mail: chenghg7932@gmail.com; Deng, Ning

2013-12-15

We investigated the influence of thermal agitation on the electric field induced precessional magnetization switching probability with perpendicular easy axis by solving the Fokker-Planck equation numerically with finite difference method. The calculated results show that the thermal agitation during the reversal process crucially influences the switching probability. The switching probability can be achieved is only determined by the thermal stability factor Δ of the free layer, it is independent on the device dimension, which is important for the high density device application. Ultra-low error rate down to the order of 10{sup −9} can be achieved for the device of thermalmore » stability factor Δ of 40. Low damping factor α material should be used for the free layer for high reliability device applications. These results exhibit potential of electric field induced precessional magnetization switching with perpendicular easy axis for ultra-low power, high speed and high density magnetic random access memory (MRAM) applications.« less

Pitching stability analysis of half-rotating wing air vehicle

NASA Astrophysics Data System (ADS)

Wang, Xiaoyi; Wu, Yang; Li, Qian; Li, Congmin; Qiu, Zhizhen

2017-06-01

Half-Rotating Wing (HRW) is a new power wing which had been developed by our work team using rotating-type flapping instead of oscillating-type flapping. Half-Rotating Wing Air Vehicle (HRWAV) is similar as Bionic Flapping Wing Air Vehicle (BFWAV). It is necessary to guarantee pitching stability of HRWAV to maintain flight stability. The working principle of HRW was firstly introduced in this paper. The rule of motion indicated that the fuselage of HRWAV without empennage would overturn forward as it generated increased pitching movement. Therefore, the empennage was added on the tail of HRWAV to balance the additional moment generated by aerodynamic force during flight. The stability analysis further shows that empennage could weaken rapidly the pitching disturbance on HRWAV and a new balance of fuselage could be achieved in a short time. Case study using numerical analysis verified correctness and validity of research results mentioned above, which could provide theoretical guidance to design and control HRWAV.

Simulation of Nonlinear Instabilities in an Attachment-Line Boundary Layer

NASA Technical Reports Server (NTRS)

Joslin, Ronald D.

1996-01-01

The linear and the nonlinear stability of disturbances that propagate along the attachment line of a three-dimensional boundary layer is considered. The spatially evolving disturbances in the boundary layer are computed by direct numerical simulation (DNS) of the unsteady, incompressible Navier-Stokes equations. Disturbances are introduced either by forcing at the in ow or by applying suction and blowing at the wall. Quasi-parallel linear stability theory and a nonparallel theory yield notably different stability characteristics for disturbances near the critical Reynolds number; the DNS results con rm the latter theory. Previously, a weakly nonlinear theory and computations revealed a high wave-number region of subcritical disturbance growth. More recent computations have failed to achieve this subcritical growth. The present computational results indicate the presence of subcritically growing disturbances; the results support the weakly nonlinear theory. Furthermore, an explanation is provided for the previous theoretical and computational discrepancy. In addition, the present results demonstrate that steady suction can be used to stabilize disturbances that otherwise grow subcritically along the attachment line.

A locally conservative stabilized continuous Galerkin finite element method for two-phase flow in poroelastic subsurfaces

NASA Astrophysics Data System (ADS)

Deng, Q.; Ginting, V.; McCaskill, B.; Torsu, P.

2017-10-01

We study the application of a stabilized continuous Galerkin finite element method (CGFEM) in the simulation of multiphase flow in poroelastic subsurfaces. The system involves a nonlinear coupling between the fluid pressure, subsurface's deformation, and the fluid phase saturation, and as such, we represent this coupling through an iterative procedure. Spatial discretization of the poroelastic system employs the standard linear finite element in combination with a numerical diffusion term to maintain stability of the algebraic system. Furthermore, direct calculation of the normal velocities from pressure and deformation does not entail a locally conservative field. To alleviate this drawback, we propose an element based post-processing technique through which local conservation can be established. The performance of the method is validated through several examples illustrating the convergence of the method, the effectivity of the stabilization term, and the ability to achieve locally conservative normal velocities. Finally, the efficacy of the method is demonstrated through simulations of realistic multiphase flow in poroelastic subsurfaces.

Olive-oil nanocapsules stabilized by HSA: influence of processing variables on particle properties

NASA Astrophysics Data System (ADS)

Molina-Bolívar, J. A.; Galisteo-González, F.

2015-10-01

Liquid lipid nanocapsules (LLN) are considered to be promising drug carriers in the medical field. The size and the surface charge of these nanocarriers are of major importance, affecting their bioavailability and the in vivo behaviour after intravenous injection. This research provides a comprehensive study on the preparation of olive-oil nanocapsules stabilized with a human serum albumin shell (HSA). LLN were prepared by modified solvent-displacement method. Numerous experimental variables were examined in order to characterize their impact on LLN size, distribution, and electrophoretic mobility. Physicochemical parameters of LLN were controlled by adjusting the nanodroplet stabilizing shell of adsorbed protein molecules, which was affected by the oil:HSA ratio, pH, and ionic strength of aqueous medium. The stronger the repulsion between adsorbed HSA molecules, the smaller and more monodisperse the particles proved. Other process parameters, including the ethanol:acetone ratio, organic:aqueous phase ratio, speed of organic-phase injection, and stirring rate were examined to achieve optimum preparation conditions. LLN produced by our standardized formulation were in the range of 170-175 nm with low polydispersity index (<0.1). Long-term colloidal stability of samples was evaluated after 6 months of storage. Efficient incorporation of curcumin, a model for a water-insoluble drug, into olive-oil nanocapsules was achieved (90 %). Encapsulation of curcumin into LLN had a stabilizing effect with respect to drug photodecomposition compared to that of the free molecule in solution.

Biomechanical Concepts for Fracture Fixation

PubMed Central

Bottlang, Michael; Schemitsch, Christine E.; Nauth, Aaron; Routt, Milton; Egol, Kenneth; Cook, Gillian E.; Schemitsch, Emil H.

2015-01-01

Application of the correct fixation construct is critical for fracture healing and long-term stability; however, it is a complex issue with numerous significant factors. This review describes a number of common fracture types, and evaluates their currently available fracture fixation constructs. In the setting of complex elbow instability, stable fixation or radial head replacement with an appropriately sized implant in conjunction with ligamentous repair is required to restore stability. For unstable sacral fractures, “standard” iliosacral screw fixation is not sufficient for fractures with vertical or multiplanar instabilities. Periprosthetic femur fractures, in particular Vancouver B1 fractures, have increased stability when using 90/90 fixation versus a single locking plate. Far Cortical Locking combines the concept of dynamization with locked plating in order to achieve superior healing of a distal femur fracture. Finally, there is no ideal construct for syndesmotic fracture stabilization; however, these fractures should be fixed using a device that allows for sufficient motion in the syndesmosis. In general, orthopaedic surgeons should select a fracture fixation construct that restores stability and promotes healing at the fracture site, while reducing the potential for fixation failure. PMID:26584263

Entropy Stable Wall Boundary Conditions for the Compressible Navier-Stokes Equations

NASA Technical Reports Server (NTRS)

Parsani, Matteo; Carpenter, Mark H.; Nielsen, Eric J.

2014-01-01

Non-linear entropy stability and a summation-by-parts framework are used to derive entropy stable wall boundary conditions for the compressible Navier-Stokes equations. A semi-discrete entropy estimate for the entire domain is achieved when the new boundary conditions are coupled with an entropy stable discrete interior operator. The data at the boundary are weakly imposed using a penalty flux approach and a simultaneous-approximation-term penalty technique. Although discontinuous spectral collocation operators are used herein for the purpose of demonstrating their robustness and efficacy, the new boundary conditions are compatible with any diagonal norm summation-by-parts spatial operator, including finite element, finite volume, finite difference, discontinuous Galerkin, and flux reconstruction schemes. The proposed boundary treatment is tested for three-dimensional subsonic and supersonic flows. The numerical computations corroborate the non-linear stability (entropy stability) and accuracy of the boundary conditions.

Entropy Stable Wall Boundary Conditions for the Three-Dimensional Compressible Navier-Stokes Equations

NASA Technical Reports Server (NTRS)

Parsani, Matteo; Carpenter, Mark H.; Nielsen, Eric J.

2015-01-01

Non-linear entropy stability and a summation-by-parts framework are used to derive entropy stable wall boundary conditions for the three-dimensional compressible Navier-Stokes equations. A semi-discrete entropy estimate for the entire domain is achieved when the new boundary conditions are coupled with an entropy stable discrete interior operator. The data at the boundary are weakly imposed using a penalty flux approach and a simultaneous-approximation-term penalty technique. Although discontinuous spectral collocation operators on unstructured grids are used herein for the purpose of demonstrating their robustness and efficacy, the new boundary conditions are compatible with any diagonal norm summation-by-parts spatial operator, including finite element, finite difference, finite volume, discontinuous Galerkin, and flux reconstruction/correction procedure via reconstruction schemes. The proposed boundary treatment is tested for three-dimensional subsonic and supersonic flows. The numerical computations corroborate the non-linear stability (entropy stability) and accuracy of the boundary conditions.

Enhanced ethylene separation and plasticization resistance in polymer membranes incorporating metal-organic framework nanocrystals.

PubMed

Bachman, Jonathan E; Smith, Zachary P; Li, Tao; Xu, Ting; Long, Jeffrey R

2016-08-01

The implementation of membrane-based separations in the petrochemical industry has the potential to reduce energy consumption significantly relative to conventional separation processes. Achieving this goal, however, requires the development of new membrane materials with greater selectivity, permeability and stability than available at present. Here, we report composite materials consisting of nanocrystals of metal-organic frameworks dispersed within a high-performance polyimide, which can exhibit enhanced selectivity for ethylene over ethane, greater ethylene permeability and improved membrane stability. Our results suggest that framework-polymer interactions reduce chain mobility of the polymer while simultaneously boosting membrane separation performance. The increased stability, or plasticization resistance, is expected to improve membrane utility under real process conditions for petrochemical separations and natural gas purification. Furthermore, this approach can be broadly applied to numerous polymers that encounter aggressive environments, potentially making gas separations possible that were previously inaccessible to membranes.

Lag and anticipating synchronization without time-delay coupling.

PubMed

Corron, Ned J; Blakely, Jonathan N; Pethel, Shawn D

2005-06-01

We describe a new method for achieving approximate lag and anticipating synchronization in unidirectionally coupled chaotic oscillators. The method uses a specific parameter mismatch between the drive and response that is a first-order approximation to true time-delay coupling. As a result, an adjustable lag or anticipation effect can be achieved without the need for a variable delay line, making the method simpler and more economical to implement in many physical systems. We present a stability analysis, demonstrate the method numerically, and report experimental observation of the effect in radio-frequency electronic oscillators. In the circuit experiments, both lag and anticipation are controlled by tuning a single capacitor in the response oscillator.

Finite-time fault tolerant attitude stabilization control for rigid spacecraft.

PubMed

Huo, Xing; Hu, Qinglei; Xiao, Bing

2014-03-01

A sliding mode based finite-time control scheme is presented to address the problem of attitude stabilization for rigid spacecraft in the presence of actuator fault and external disturbances. More specifically, a nonlinear observer is first proposed to reconstruct the amplitude of actuator faults and external disturbances. It is proved that precise reconstruction with zero observer error is achieved in finite time. Then, together with the system states, the reconstructed information is used to synthesize a nonsingular terminal sliding mode attitude controller. The attitude and the angular velocity are asymptotically governed to zero with finite-time convergence. A numerical example is presented to demonstrate the effectiveness of the proposed scheme. © 2013 Published by ISA on behalf of ISA.

Front and pulse solutions for the complex Ginzburg-Landau equation with higher-order terms.

PubMed

Tian, Huiping; Li, Zhonghao; Tian, Jinping; Zhou, Guosheng

2002-12-01

We investigate one-dimensional complex Ginzburg-Landau equation with higher-order terms and discuss their influences on the multiplicity of solutions. An exact analytic front solution is presented. By stability analysis for the original partial differential equation, we derive its necessary stability condition for amplitude perturbations. This condition together with the exact front solution determine the region of parameter space where the uniformly translating front solution can exist. In addition, stable pulses, chaotic pulses, and attenuation pulses appear generally if the parameters are out of the range. Finally, applying these analysis into the optical transmission system numerically we find that the stable transmission of optical pulses can be achieved if the parameters are appropriately chosen.

Numerical study on the Welander oscillatory natural circulation problem using high-order numerical methods

DOE PAGES

Zou, Ling; Zhao, Haihua; Kim, Seung Jun

2016-11-16

In this study, the classical Welander’s oscillatory natural circulation problem is investigated using high-order numerical methods. As originally studied by Welander, the fluid motion in a differentially heated fluid loop can exhibit stable, weakly instable, and strongly instable modes. A theoretical stability map has also been originally derived from the stability analysis. Numerical results obtained in this paper show very good agreement with Welander’s theoretical derivations. For stable cases, numerical results from both the high-order and low-order numerical methods agree well with the non-dimensional flow rate analytically derived. The high-order numerical methods give much less numerical errors compared to themore » low-order methods. For stability analysis, the high-order numerical methods could perfectly predict the stability map, while the low-order numerical methods failed to do so. For all theoretically unstable cases, the low-order methods predicted them to be stable. The result obtained in this paper is a strong evidence to show the benefits of using high-order numerical methods over the low-order ones, when they are applied to simulate natural circulation phenomenon that has already gain increasing interests in many future nuclear reactor designs.« less

Elastic Differential Cross Sections

NASA Technical Reports Server (NTRS)

Werneth, Charles M.; Maung, Khin M.; Ford, William P.; Norbury, John W.; Vera, Michael D.

2014-01-01

The eikonal, partial wave (PW) Lippmann-Schwinger, and three-dimensional Lippmann-Schwinger (LS3D) methods are compared for nuclear reactions that are relevant for space radiation applications. Numerical convergence of the eikonal method is readily achieved when exact formulas of the optical potential are used for light nuclei (A less than or equal to 16) and the momentum-space optical potential is used for heavier nuclei. The PW solution method is known to be numerically unstable for systems that require a large number of partial waves, and, as a result, the LS3D method is employed. The effect of relativistic kinematics is studied with the PW and LS3D methods and is compared to eikonal results. It is recommended that the LS3D method be used for high energy nucleon- nucleus reactions and nucleus-nucleus reactions at all energies because of its rapid numerical convergence and stability.

Nuclear Cross Sections for Space Radiation Applications

NASA Technical Reports Server (NTRS)

Werneth, C. M.; Maung, K. M.; Ford, W. P.; Norbury, J. W.; Vera, M. D.

2015-01-01

The eikonal, partial wave (PW) Lippmann-Schwinger, and three-dimensional Lippmann-Schwinger (LS3D) methods are compared for nuclear reactions that are relevant for space radiation applications. Numerical convergence of the eikonal method is readily achieved when exact formulas of the optical potential are used for light nuclei (A = 16) and the momentum-space optical potential is used for heavier nuclei. The PW solution method is known to be numerically unstable for systems that require a large number of partial waves, and, as a result, the LS3D method is employed. The effect of relativistic kinematics is studied with the PW and LS3D methods and is compared to eikonal results. It is recommended that the LS3D method be used for high energy nucleon-nucleus reactions and nucleus-nucleus reactions at all energies because of its rapid numerical convergence and stability for both non-relativistic and relativistic kinematics.

The DUV Stability of Superlattice-Doped CMOS Detector Arrays

NASA Technical Reports Server (NTRS)

Hoenk, M. E.; Carver, A.; Jones, T.; Dickie, M.; Cheng, P.; Greer, H. F.; Nikzad, S.; Sgro, J.

2013-01-01

In this paper, we present experimental results and band structure calculations that illuminate the unique properties of superlattice-doped detectors. Numerical band structure calculations are presented to analyze the dependencies of surface passivation on dopant profiles and interface trap densities (Figure 3). Experiments and calculations show that quantum-engineered surfaces, grown at JPL by low temperature molecular beam epitaxy, achieve a qualitative as well as quantitative uniqueness in their near-immunity to high densities of surface and interface traps.

Cluster-modified function projective synchronisation of complex networks with asymmetric coupling

NASA Astrophysics Data System (ADS)

Wang, Shuguo

2018-02-01

This paper investigates the cluster-modified function projective synchronisation (CMFPS) of a generalised linearly coupled network with asymmetric coupling and nonidentical dynamical nodes. A novel synchronisation scheme is proposed to achieve CMFPS in community networks. We use adaptive control method to derive CMFPS criteria based on Lyapunov stability theory. Each cluster of networks is synchronised with target system by state transformation with scaling function matrix. Numerical simulation results are presented finally to illustrate the effectiveness of this method.

Numerical investigation of the pseudopotential lattice Boltzmann modeling of liquid-vapor for multi-phase flows

NASA Astrophysics Data System (ADS)

Nemati, Maedeh; Shateri Najaf Abady, Ali Reza; Toghraie, Davood; Karimipour, Arash

2018-01-01

The incorporation of different equations of state into single-component multiphase lattice Boltzmann model is considered in this paper. The original pseudopotential model is first detailed, and several cubic equations of state, the Redlich-Kwong, Redlich-Kwong-Soave, and Peng-Robinson are then incorporated into the lattice Boltzmann model. A comparison of the numerical simulation achievements on the basis of density ratios and spurious currents is used for presentation of the details of phase separation in these non-ideal single-component systems. The paper demonstrates that the scheme for the inter-particle interaction force term as well as the force term incorporation method matters to achieve more accurate and stable results. The velocity shifting method is demonstrated as the force term incorporation method, among many, with accuracy and stability results. Kupershtokh scheme also makes it possible to achieve large density ratio (up to 104) and to reproduce the coexistence curve with high accuracy. Significant reduction of the spurious currents at vapor-liquid interface is another observation. High-density ratio and spurious current reduction resulted from the Redlich-Kwong-Soave and Peng-Robinson EOSs, in higher accordance with the Maxwell construction results.

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

Zou, Ling; Zhao, Haihua; Kim, Seung Jun

In this study, the classical Welander’s oscillatory natural circulation problem is investigated using high-order numerical methods. As originally studied by Welander, the fluid motion in a differentially heated fluid loop can exhibit stable, weakly instable, and strongly instable modes. A theoretical stability map has also been originally derived from the stability analysis. Numerical results obtained in this paper show very good agreement with Welander’s theoretical derivations. For stable cases, numerical results from both the high-order and low-order numerical methods agree well with the non-dimensional flow rate analytically derived. The high-order numerical methods give much less numerical errors compared to themore » low-order methods. For stability analysis, the high-order numerical methods could perfectly predict the stability map, while the low-order numerical methods failed to do so. For all theoretically unstable cases, the low-order methods predicted them to be stable. The result obtained in this paper is a strong evidence to show the benefits of using high-order numerical methods over the low-order ones, when they are applied to simulate natural circulation phenomenon that has already gain increasing interests in many future nuclear reactor designs.« less

Lens-mount stability trade-off: a survey exemplified for DUV wafer inspection objectives

NASA Astrophysics Data System (ADS)

Bouazzam, Achmed; Erbe, Torsten; Fahr, Stephan; Werschnik, Jan

2015-09-01

The position stability of optical elements is an essential part of the tolerance budget of an optical system because its compensation would require an alignment step after the lens has left the factory. In order to achieve a given built performance the stability error contribution needs to be known and accounted for. Given a high-end lens touching the edge of technology not knowing, under- or overestimating this contribution becomes a serious cost and risk factor. If overestimated the remaining parts of the budget need to be tighter. If underestimated the total project might fail. For many mounting principles the stability benchmark is based on previous systems or information gathered by elaborated testing of complete optical systems. This renders the development of a new system into a risky endeavour, because these experiences are not sufficiently precise and tend to be not transferable when scaling of the optical elements is intended. This contribution discusses the influences of different optical mounting concepts on the position stability using the example of high numerical aperture (HNA) inspection lenses working in the deep ultraviolet (DUV) spectrum. A method to investigate the positional stability is presented for selected mounting examples typical for inspection lenses.

Exact and Approximate Stability of Solutions to Traveling Salesman Problems.

PubMed

Niendorf, Moritz; Girard, Anouck R

2018-02-01

This paper presents the stability analysis of an optimal tour for the symmetric traveling salesman problem (TSP) by obtaining stability regions. The stability region of an optimal tour is the set of all cost changes for which that solution remains optimal and can be understood as the margin of optimality for a solution with respect to perturbations in the problem data. It is known that it is not possible to test in polynomial time whether an optimal tour remains optimal after the cost of an arbitrary set of edges changes. Therefore, this paper develops tractable methods to obtain under and over approximations of stability regions based on neighborhoods and relaxations. The application of the results to the two-neighborhood and the minimum 1 tree (M1T) relaxation are discussed in detail. For Euclidean TSPs, stability regions with respect to vertex location perturbations and the notion of safe radii and location criticalities are introduced. Benefits of this paper include insight into robustness properties of tours, minimum spanning trees, M1Ts, and fast methods to evaluate optimality after perturbations occur. Numerical examples are given to demonstrate the methods and achievable approximation quality.

Centralized Multi-Sensor Square Root Cubature Joint Probabilistic Data Association

PubMed Central

Liu, Jun; Li, Gang; Qi, Lin; Li, Yaowen; He, You

2017-01-01

This paper focuses on the tracking problem of multiple targets with multiple sensors in a nonlinear cluttered environment. To avoid Jacobian matrix computation and scaling parameter adjustment, improve numerical stability, and acquire more accurate estimated results for centralized nonlinear tracking, a novel centralized multi-sensor square root cubature joint probabilistic data association algorithm (CMSCJPDA) is proposed. Firstly, the multi-sensor tracking problem is decomposed into several single-sensor multi-target tracking problems, which are sequentially processed during the estimation. Then, in each sensor, the assignment of its measurements to target tracks is accomplished on the basis of joint probabilistic data association (JPDA), and a weighted probability fusion method with square root version of a cubature Kalman filter (SRCKF) is utilized to estimate the targets’ state. With the measurements in all sensors processed CMSCJPDA is derived and the global estimated state is achieved. Experimental results show that CMSCJPDA is superior to the state-of-the-art algorithms in the aspects of tracking accuracy, numerical stability, and computational cost, which provides a new idea to solve multi-sensor tracking problems. PMID:29113085

Centralized Multi-Sensor Square Root Cubature Joint Probabilistic Data Association.

PubMed

Liu, Yu; Liu, Jun; Li, Gang; Qi, Lin; Li, Yaowen; He, You

2017-11-05

This paper focuses on the tracking problem of multiple targets with multiple sensors in a nonlinear cluttered environment. To avoid Jacobian matrix computation and scaling parameter adjustment, improve numerical stability, and acquire more accurate estimated results for centralized nonlinear tracking, a novel centralized multi-sensor square root cubature joint probabilistic data association algorithm (CMSCJPDA) is proposed. Firstly, the multi-sensor tracking problem is decomposed into several single-sensor multi-target tracking problems, which are sequentially processed during the estimation. Then, in each sensor, the assignment of its measurements to target tracks is accomplished on the basis of joint probabilistic data association (JPDA), and a weighted probability fusion method with square root version of a cubature Kalman filter (SRCKF) is utilized to estimate the targets' state. With the measurements in all sensors processed CMSCJPDA is derived and the global estimated state is achieved. Experimental results show that CMSCJPDA is superior to the state-of-the-art algorithms in the aspects of tracking accuracy, numerical stability, and computational cost, which provides a new idea to solve multi-sensor tracking problems.

Stabilizing canonical-ensemble calculations in the auxiliary-field Monte Carlo method

NASA Astrophysics Data System (ADS)

Gilbreth, C. N.; Alhassid, Y.

2015-03-01

Quantum Monte Carlo methods are powerful techniques for studying strongly interacting Fermi systems. However, implementing these methods on computers with finite-precision arithmetic requires careful attention to numerical stability. In the auxiliary-field Monte Carlo (AFMC) method, low-temperature or large-model-space calculations require numerically stabilized matrix multiplication. When adapting methods used in the grand-canonical ensemble to the canonical ensemble of fixed particle number, the numerical stabilization increases the number of required floating-point operations for computing observables by a factor of the size of the single-particle model space, and thus can greatly limit the systems that can be studied. We describe an improved method for stabilizing canonical-ensemble calculations in AFMC that exhibits better scaling, and present numerical tests that demonstrate the accuracy and improved performance of the method.

Surface mediated assembly of small, metastable gold nanoclusters

NASA Astrophysics Data System (ADS)

Pettibone, John M.; Osborn, William A.; Rykaczewski, Konrad; Talin, A. Alec; Bonevich, John E.; Hudgens, Jeffrey W.; Allendorf, Mark D.

2013-06-01

The unique properties of metallic nanoclusters are attractive for numerous commercial and industrial applications but are generally less stable than nanocrystals. Thus, developing methodologies for stabilizing nanoclusters and retaining their enhanced functionality is of great interest. We report the assembly of PPh3-protected Au9 clusters from a heterogeneous mixture into films consisting of sub 3 nm nanocluster assemblies. The depositing nanoclusters are metastable in solution, but the resulting nanocluster assemblies are stabilized indefinitely in air or fresh solvent. The films exhibit distinct structure from Au nanoparticles observed by X-ray diffraction, and film dissolution data support the preservation of small nanoclusters. UV-Vis spectroscopy, electrospray ionization mass spectrometry, X-ray photoelectron spectroscopy and electron microscopy are used to elucidate information regarding the nanocluster formation and assembly mechanism. Preferential deposition of nanocluster assemblies can be achieved on multiple substrates, including polymer, Cr, Si, SiO2, SiNx, and metal-organic frameworks (MOFs). Unlike other vapor phase coating processes, nanocluster assembly on the MIL-68(In) MOF crystal is capable of preferentially coating the external surface and stabilizing the crystal structure in hydrothermal conditions, which should enhance their storage, separation and delivery capabilities.The unique properties of metallic nanoclusters are attractive for numerous commercial and industrial applications but are generally less stable than nanocrystals. Thus, developing methodologies for stabilizing nanoclusters and retaining their enhanced functionality is of great interest. We report the assembly of PPh3-protected Au9 clusters from a heterogeneous mixture into films consisting of sub 3 nm nanocluster assemblies. The depositing nanoclusters are metastable in solution, but the resulting nanocluster assemblies are stabilized indefinitely in air or fresh solvent. The films exhibit distinct structure from Au nanoparticles observed by X-ray diffraction, and film dissolution data support the preservation of small nanoclusters. UV-Vis spectroscopy, electrospray ionization mass spectrometry, X-ray photoelectron spectroscopy and electron microscopy are used to elucidate information regarding the nanocluster formation and assembly mechanism. Preferential deposition of nanocluster assemblies can be achieved on multiple substrates, including polymer, Cr, Si, SiO2, SiNx, and metal-organic frameworks (MOFs). Unlike other vapor phase coating processes, nanocluster assembly on the MIL-68(In) MOF crystal is capable of preferentially coating the external surface and stabilizing the crystal structure in hydrothermal conditions, which should enhance their storage, separation and delivery capabilities. Electronic supplementary information (ESI) available: Further details on stored plating solution preparation, film characterization, solution processing, MOF crystal FIB reconstruction and stability are available. See DOI: 10.1039/c3nr01708g

Transonic airfoil analysis and design in nonuniform flow

NASA Technical Reports Server (NTRS)

Chang, J. F.; Lan, C. E.

1986-01-01

A nonuniform transonic airfoil code is developed for applications in analysis, inverse design and direct optimization involving an airfoil immersed in propfan slipstream. Problems concerning the numerical stability, convergence, divergence and solution oscillations are discussed. The code is validated by comparing with some known results in incompressible flow. A parametric investigation indicates that the airfoil lift-drag ratio can be increased by decreasing the thickness ratio. A better performance can be achieved if the airfoil is located below the slipstream center. Airfoil characteristics designed by the inverse method and a direct optimization are compared. The airfoil designed with the method of direct optimization exhibits better characteristics and achieves a gain of 22 percent in lift-drag ratio with a reduction of 4 percent in thickness.

Stabilizing photoassociated Cs2 molecules by optimal control

NASA Astrophysics Data System (ADS)

Zhang, Wei; Xie, Ting; Huang, Yin; Wang, Gao-Ren; Cong, Shu-Lin

2013-01-01

We demonstrate theoretically that photoassociated molecules can be stabilized to deeply bound states. This process is achieved by transferring the population from the outer well to the inner well using the optimal control theory, the Cs2 molecule is taken as an example. Numerical calculations show that weakly bound molecules formed in the outer well by a pump pulse can be compressed to the inner well via a vibrational level of the ground electronic state as an intermediary by an additionally optimized laser pulse. The positively chirped pulse can enhance the population of the target state. With a transform-limited dump pulse, nearly all the photoassociated molecules in the inner well of the excited electronic state can be transferred to the deeply vibrational level of the ground electronic state.

Aperiodic Robust Model Predictive Control for Constrained Continuous-Time Nonlinear Systems: An Event-Triggered Approach.

PubMed

Liu, Changxin; Gao, Jian; Li, Huiping; Xu, Demin

2018-05-01

The event-triggered control is a promising solution to cyber-physical systems, such as networked control systems, multiagent systems, and large-scale intelligent systems. In this paper, we propose an event-triggered model predictive control (MPC) scheme for constrained continuous-time nonlinear systems with bounded disturbances. First, a time-varying tightened state constraint is computed to achieve robust constraint satisfaction, and an event-triggered scheduling strategy is designed in the framework of dual-mode MPC. Second, the sufficient conditions for ensuring feasibility and closed-loop robust stability are developed, respectively. We show that robust stability can be ensured and communication load can be reduced with the proposed MPC algorithm. Finally, numerical simulations and comparison studies are performed to verify the theoretical results.

Disturbance observer based active and adaptive synchronization of energy resource chaotic system.

PubMed

Wei, Wei; Wang, Meng; Li, Donghai; Zuo, Min; Wang, Xiaoyi

2016-11-01

In this paper, synchronization of a three-dimensional energy resource chaotic system is considered. For the sake of achieving the synchronization between the drive and response systems, two different nonlinear control approaches, i.e. active control with known parameters and adaptive control with unknown parameters, have been designed. In order to guarantee the transient performance, finite-time boundedness (FTB) and finite-time stability (FTS) are introduced in the design of active control and adaptive control, respectively. Simultaneously, in view of the existence of disturbances, a new disturbance observer is proposed to estimate the disturbance. The conditions of the asymptotic stability for the closed-loop system are obtained. Numerical simulations are provided to illustrate the proposed approaches. Copyright © 2016 ISA. Published by Elsevier Ltd. All rights reserved.

Boattail juncture shaping for spin-stabilized rounds in supersonic flight

NASA Astrophysics Data System (ADS)

Jiajan, W.; Chue, R. S. M.; Nguyen, T.; Yu, S. C. M.

2015-03-01

In this paper, the effects of boattail junction shaping on aerodynamic drag and stability of supersonic spin-stabilized rounds are investigated using computational fluid dynamics. For a generic round body comprising of a secant-ogive nose, a cylindrical body and a conical boattail, the shaping technique was achieved by adding a convex surface of varying degrees of radius of curvature to the junction between the cylindrical body and the boattail. It was shown through numerical simulations that this shaping technique can provide a reduction in aerodynamic drag of up to 5.4 % without destabilizing the round bodies when the radius of curvature is less than 8.8 times the diameter of the cylindrical body. The more gradual change of the flow characteristics, e.g., the pressure over the convex surface, was identified as the main reason for the drag reduction. A unique aspect of the current work is that stability is treated as an integral part of the performance assessment. It was also found that the dynamic instability encountered at large radii of curvature is due to the Magnus effects.

Stability analysis of multigrid acceleration methods for the solution of partial differential equations

NASA Technical Reports Server (NTRS)

Fay, John F.

1990-01-01

A calculation is made of the stability of various relaxation schemes for the numerical solution of partial differential equations. A multigrid acceleration method is introduced, and its effects on stability are explored. A detailed stability analysis of a simple case is carried out and verified by numerical experiment. It is shown that the use of multigrids can speed convergence by several orders of magnitude without adversely affecting stability.

Dynamical behavior and Jacobi stability analysis of wound strings

NASA Astrophysics Data System (ADS)

Lake, Matthew J.; Harko, Tiberiu

2016-06-01

We numerically solve the equations of motion (EOM) for two models of circular cosmic string loops with windings in a simply connected internal space. Since the windings cannot be topologically stabilized, stability must be achieved (if at all) dynamically. As toy models for realistic compactifications, we consider windings on a small section of mathbb {R}^2, which is valid as an approximation to any simply connected internal manifold if the winding radius is sufficiently small, and windings on an S^2 of constant radius mathcal {R}. We then use Kosambi-Cartan-Chern (KCC) theory to analyze the Jacobi stability of the string equations and determine bounds on the physical parameters that ensure dynamical stability of the windings. We find that, for the same initial conditions, the curvature and topology of the internal space have nontrivial effects on the microscopic behavior of the string in the higher dimensions, but that the macroscopic behavior is remarkably insensitive to the details of the motion in the compact space. This suggests that higher-dimensional signatures may be extremely difficult to detect in the effective (3+1)-dimensional dynamics of strings compactified on an internal space, even if configurations with nontrivial windings persist over long time periods.

Effect of Rotational Speed on the Stability of Two Rotating Side-by-side Circular Cylinders at Low Reynolds Number

NASA Astrophysics Data System (ADS)

Dou, Huashu; Zhang, Shuo; Yang, Hui; Setoguchi, Toshiaki; Kinoue, Yoichi

2018-04-01

Flow around two rotating side-by-side circular cylinders of equal diameter D is numerically studied at the Reynolds number 40≤ Re ≤200 and various rotation rate θ i . The incoming flow is assumed to be two-dimensional laminar flow. The governing equations are the incompressible Navier-Stokes equations and solved by the finite volume method (FVM). The ratio of the center-to-center spacing to the cylinder diameter is T/D=2. The objective of the present work is to investigate the effect of rotational speed and Reynolds number on the stability of the flow. The simulation results are compared with the experimental data and a good agreement is achieved. The stability of the flow is analyzed by using the energy gradient theory, which produces the energy gradient function K to identify the region where the flow is the most prone to be destabilized and the degree of the destabilization. Numerical results reveal that K is the most significant at the separated shear layers of the cylinder pair. With Re increases, the length of the wake is shorter and the vortex shedding generally exhibits a symmetrical distribution for θ i < θ crit . It is also shown that the unsteady vortex shedding can be suppressed by rotating the cylinders in the counter-rotating mode.

Tuning the synchronization of a network of weakly coupled self-oscillating gels via capacitors.

PubMed

Fang, Yan; Yashin, Victor V; Dickerson, Samuel J; Balazs, Anna C

2018-05-01

We consider a network of coupled oscillating units, where each unit comprises a self-oscillating polymer gel undergoing the Belousov-Zhabotinsky (BZ) reaction and an overlaying piezoelectric (PZ) cantilever. Through chemo-mechano-electrical coupling, the oscillations of the networked BZ-PZ units achieve in-phase or anti-phase synchronization, enabling, for example, the storage of information within the system. Herein, we develop numerical and computational models to show that the introduction of capacitors into the BZ-PZ system enhances the dynamical behavior of the oscillating network by yielding additional stable synchronization modes. We specifically show that the capacitors lead to a redistribution of charge in the system and alteration of the force that the PZ cantilevers apply to the underlying gel. Hence, the capacitors modify the strength of the coupling between the oscillators in the network. We utilize a linear stability analysis to determine the phase behavior of BZ-PZ networks encompassing different capacitances, force polarities, and number of units and then verify our findings with numerical simulations. Thus, through analytical calculations and numerical simulations, we determine the impact of the capacitors on the existence of the synchronization modes, their stability, and the rate of synchronization within these complex dynamical systems. The findings from our study can be used to design robotic materials that harness the materials' intrinsic, responsive properties to perform such functions as sensing, actuation, and information storage.

Tuning the synchronization of a network of weakly coupled self-oscillating gels via capacitors

NASA Astrophysics Data System (ADS)

Fang, Yan; Yashin, Victor V.; Dickerson, Samuel J.; Balazs, Anna C.

2018-05-01

We consider a network of coupled oscillating units, where each unit comprises a self-oscillating polymer gel undergoing the Belousov-Zhabotinsky (BZ) reaction and an overlaying piezoelectric (PZ) cantilever. Through chemo-mechano-electrical coupling, the oscillations of the networked BZ-PZ units achieve in-phase or anti-phase synchronization, enabling, for example, the storage of information within the system. Herein, we develop numerical and computational models to show that the introduction of capacitors into the BZ-PZ system enhances the dynamical behavior of the oscillating network by yielding additional stable synchronization modes. We specifically show that the capacitors lead to a redistribution of charge in the system and alteration of the force that the PZ cantilevers apply to the underlying gel. Hence, the capacitors modify the strength of the coupling between the oscillators in the network. We utilize a linear stability analysis to determine the phase behavior of BZ-PZ networks encompassing different capacitances, force polarities, and number of units and then verify our findings with numerical simulations. Thus, through analytical calculations and numerical simulations, we determine the impact of the capacitors on the existence of the synchronization modes, their stability, and the rate of synchronization within these complex dynamical systems. The findings from our study can be used to design robotic materials that harness the materials' intrinsic, responsive properties to perform such functions as sensing, actuation, and information storage.

LISA on Table: an optical simulator for LISA

NASA Astrophysics Data System (ADS)

Halloin, H.; Jeannin, O.; Argence, B.; Bourrier, V.; de Vismes, E.; Prat, P.

2017-11-01

LISA, the first space project for detecting gravitational waves, relies on two main technical challenges: the free falling masses and an outstanding precision on phase shift measurements (a few pm on 5 Mkm in the LISA band). The technology of the free falling masses, i.e. their isolation to forces other than gravity and the capability for the spacecraft to precisely follow the test masses, will soon be tested with the technological LISA Pathfinder mission. The performance of the phase measurement will be achieved by at least two stabilization stages: a pre-stabilisation of the laser frequency at a level of 10-13 (relative frequency stability) will be further improved by using numerical algorithms, such as Time Delay Interferometry, which have been theoretically and numerically demonstrated to reach the required performance level (10-21). Nevertheless, these algorithms, though already tested with numerical model of LISA, require experimental validation, including `realistic' hardware elements. Such an experiment would allow to evaluate the expected noise level and the possible interactions between subsystems. To this end, the APC is currently developing an optical benchtop experiment, called LISA On Table (LOT), which is representative of the three LISA spacecraft. A first module of the LOT experiment has been mounted and is being characterized. After completion this facility may be used by the LISA community to test hardware (photodiodes, phasemeters) or software (reconstruction algorithms) components.

H∞ control of combustion in diesel engines using a discrete dynamics model

NASA Astrophysics Data System (ADS)

Hirata, Mitsuo; Ishizuki, Sota; Suzuki, Masayasu

2016-09-01

This paper proposes a control method for combustion in diesel engines using a discrete dynamics model. The proposed two-degree-of-freedom control scheme achieves not only good feedback properties such as disturbance suppression and robust stability but also a good transient response. The method includes a feedforward controller constructed from the inverse model of the plant, and a feedback controller designed by an Hcontrol method, which reduces the effect of the turbocharger lag. The effectiveness of the proposed method is evaluated via numerical simulations.

Finite-time hybrid projective synchronization of the drive-response complex networks with distributed-delay via adaptive intermittent control

NASA Astrophysics Data System (ADS)

Cheng, Lin; Yang, Yongqing; Li, Li; Sui, Xin

2018-06-01

This paper studies the finite-time hybrid projective synchronization of the drive-response complex networks. In the model, general transmission delays and distributed delays are also considered. By designing the adaptive intermittent controllers, the response network can achieve hybrid projective synchronization with the drive system in finite time. Based on finite-time stability theory and several differential inequalities, some simple finite-time hybrid projective synchronization criteria are derived. Two numerical examples are given to illustrate the effectiveness of the proposed method.

Mobility of discrete multibreathers in the exciton dynamics of the Davydov model with saturable nonlinearities.

PubMed

Tchinang Tchameu, J D; Togueu Motcheyo, A B; Tchawoua, C

2014-10-01

We show that the state of amide-I excitations in proteins is modeled by the discrete nonlinear Schrödinger equation with saturable nonlinearities. This is done by extending the Davydov model to take into account the competition between local compression and local dilatation of the lattice, thus leading to the interplay between self-focusing and defocusing saturable nonlinearities. Site-centered (sc) mode and/or bond-centered mode like discrete multihump soliton (DMHS) solutions are found numerically and their stability is analyzed. As a result, we obtained the existence and stability diagrams for all observed types of sc DMHS solutions. We also note that the stability of sc DMHS solutions depends not only on the value of the interpeak separation but also on the number of peaks, while their counterpart having at least one intersite soliton is instable. A study of mobility is achieved and it appears that, depending on the higher-order saturable nonlinearity, DMHS-like mechanism for vibrational energy transport along the protein chain is possible.

Optimal subinterval selection approach for power system transient stability simulation

DOE PAGES

Kim, Soobae; Overbye, Thomas J.

2015-10-21

Power system transient stability analysis requires an appropriate integration time step to avoid numerical instability as well as to reduce computational demands. For fast system dynamics, which vary more rapidly than what the time step covers, a fraction of the time step, called a subinterval, is used. However, the optimal value of this subinterval is not easily determined because the analysis of the system dynamics might be required. This selection is usually made from engineering experiences, and perhaps trial and error. This paper proposes an optimal subinterval selection approach for power system transient stability analysis, which is based on modalmore » analysis using a single machine infinite bus (SMIB) system. Fast system dynamics are identified with the modal analysis and the SMIB system is used focusing on fast local modes. An appropriate subinterval time step from the proposed approach can reduce computational burden and achieve accurate simulation responses as well. As a result, the performance of the proposed method is demonstrated with the GSO 37-bus system.« less

Higher-order compositional modeling of three-phase flow in 3D fractured porous media based on cross-flow equilibrium

NASA Astrophysics Data System (ADS)

Moortgat, Joachim; Firoozabadi, Abbas

2013-10-01

Numerical simulation of multiphase compositional flow in fractured porous media, when all the species can transfer between the phases, is a real challenge. Despite the broad applications in hydrocarbon reservoir engineering and hydrology, a compositional numerical simulator for three-phase flow in fractured media has not appeared in the literature, to the best of our knowledge. In this work, we present a three-phase fully compositional simulator for fractured media, based on higher-order finite element methods. To achieve computational efficiency, we invoke the cross-flow equilibrium (CFE) concept between discrete fractures and a small neighborhood in the matrix blocks. We adopt the mixed hybrid finite element (MHFE) method to approximate convective Darcy fluxes and the pressure equation. This approach is the most natural choice for flow in fractured media. The mass balance equations are discretized by the discontinuous Galerkin (DG) method, which is perhaps the most efficient approach to capture physical discontinuities in phase properties at the matrix-fracture interfaces and at phase boundaries. In this work, we account for gravity and Fickian diffusion. The modeling of capillary effects is discussed in a separate paper. We present the mathematical framework, using the implicit-pressure-explicit-composition (IMPEC) scheme, which facilitates rigorous thermodynamic stability analyses and the computation of phase behavior effects to account for transfer of species between the phases. A deceptively simple CFL condition is implemented to improve numerical stability and accuracy. We provide six numerical examples at both small and larger scales and in two and three dimensions, to demonstrate powerful features of the formulation.

Damping in Space Constructions

NASA Astrophysics Data System (ADS)

de Vreugd, Jan; de Lange, Dorus; Winters, Jasper; Human, Jet; Kamphues, Fred; Tabak, Erik

2014-06-01

Monolithic structures are often used in optomechanical designs for space applications to achieve high dimensional stability and to prevent possible backlash and friction phenomena. The capacity of monolithic structures to dissipate mechanical energy is however limited due to the high Q-factor, which might result in high stresses during dynamic launch loads like random vibration, sine sweeps and shock. To reduce the Q-factor in space applications, the effect of constrained layer damping (CLD) is investigated in this work. To predict the damping increase, the CLD effect is implemented locally at the supporting struts in an existing FE model of an optical instrument. Numerical simulations show that the effect of local damping treatment in this instrument could reduce the vibrational stresses with 30-50%. Validation experiments on a simple structure showed good agreement between measured and predicted damping properties. This paper presents material characterization, material modeling, numerical implementation of damping models in finite element code, numerical results on space hardware and the results of validation experiments.

A numerical comparison of discrete Kalman filtering algorithms: An orbit determination case study

NASA Technical Reports Server (NTRS)

Thornton, C. L.; Bierman, G. J.

1976-01-01

The numerical stability and accuracy of various Kalman filter algorithms are thoroughly studied. Numerical results and conclusions are based on a realistic planetary approach orbit determination study. The case study results of this report highlight the numerical instability of the conventional and stabilized Kalman algorithms. Numerical errors associated with these algorithms can be so large as to obscure important mismodeling effects and thus give misleading estimates of filter accuracy. The positive result of this study is that the Bierman-Thornton U-D covariance factorization algorithm is computationally efficient, with CPU costs that differ negligibly from the conventional Kalman costs. In addition, accuracy of the U-D filter using single-precision arithmetic consistently matches the double-precision reference results. Numerical stability of the U-D filter is further demonstrated by its insensitivity of variations in the a priori statistics.

Numerical stability of the error diffusion concept

NASA Astrophysics Data System (ADS)

Weissbach, Severin; Wyrowski, Frank

1992-10-01

The error diffusion algorithm is an easy implementable mean to handle nonlinearities in signal processing, e.g. in picture binarization and coding of diffractive elements. The numerical stability of the algorithm depends on the choice of the diffusion weights. A criterion for the stability of the algorithm is presented and evaluated for some examples.

Experimental and theoretical developments in the Mochi project

NASA Astrophysics Data System (ADS)

You, Setthivoine; von der Linden, Jens; Vereen, Keon; Lavine, Eric Sander; Carroll, Evan; Card, Alexander; Azuara-Rosales, Manuel; Quinley, Morgan; Yun, Gunsu

2015-11-01

The Mochi project investigates the interaction between magnetic fields and plasma flows in cylindrical and toroidal geometries. The configuration is designed to tailor the radial electric field profile with three annular electrodes and allow for shear helical flows in magnetized plasma jets or merging spheromaks. First plasma has been achieved and characterization is in progress with images, magnetic probes, an energy analyzer, an interferometer, a fast ion gauge, and optical and RF spectroscopy. Vector tomography of ion Doppler spectroscopy is progressing with the design of the custom fiber bundle and implementation of the numerical code. The first experiments are investigating the coupling of sausage and kink instabilities, comparing measurements to a new stability criterion and a numerical stability code. A new canonical field theory has been developed to help interpret the dynamics of plasma self-organization. The theory augments the Lagrangian of general dynamical systems to rigourously demonstrate that canonical helicity transport is valid across single particle, kinetic and fluid regimes, that dynamical equations can be re-formulated as a form of Maxwell's equations, and that helicity is conserved only when density gradients are shallow. This work is supported by US DOE Grant DE-SC0010340.

Ares-I Bending Filter Design using a Constrained Optimization Approach

NASA Technical Reports Server (NTRS)

Hall, Charles; Jang, Jiann-Woei; Hall, Robert; Bedrossian, Nazareth

2008-01-01

The Ares-I launch vehicle represents a challenging flex-body structural environment for control system design. Software filtering of the inertial sensor output is required to ensure adequate stable response to guidance commands while minimizing trajectory deviations. This paper presents a design methodology employing numerical optimization to develop the Ares-I bending filters. The design objectives include attitude tracking accuracy and robust stability with respect to rigid body dynamics, propellant slosh, and flex. Under the assumption that the Ares-I time-varying dynamics and control system can be frozen over a short period of time, the bending filters are designed to stabilize all the selected frozen-time launch control systems in the presence of parameter uncertainty. To ensure adequate response to guidance command, step response specifications are introduced as constraints in the optimization problem. Imposing these constrains minimizes performance degradation caused by the addition of the bending filters. The first stage bending filter design achieves stability by adding lag to the first structural frequency to phase stabilize the first flex mode while gain stabilizing the higher modes. The upper stage bending filter design gain stabilizes all the flex bending modes. The bending filter designs provided here have been demonstrated to provide stable first and second stage control systems in both Draper Ares Stability Analysis Tool (ASAT) and the MSFC MAVERIC 6DOF nonlinear time domain simulation.

Ares I Flight Control System Design

NASA Technical Reports Server (NTRS)

Jang, Jiann-Woei; Alaniz, Abran; Hall, Robert; Bedrossian, Nazareth; Hall, Charles; Ryan, Stephen; Jackson, Mark

2010-01-01

The Ares I launch vehicle represents a challenging flex-body structural environment for flight control system design. This paper presents a design methodology for employing numerical optimization to develop the Ares I flight control system. The design objectives include attitude tracking accuracy and robust stability with respect to rigid body dynamics, propellant slosh, and flex. Under the assumption that the Ares I time-varying dynamics and control system can be frozen over a short period of time, the flight controllers are designed to stabilize all selected frozen-time launch control systems in the presence of parametric uncertainty. Flex filters in the flight control system are designed to minimize the flex components in the error signals before they are sent to the attitude controller. To ensure adequate response to guidance command, step response specifications are introduced as constraints in the optimization problem. Imposing these constraints minimizes performance degradation caused by the addition of the flex filters. The first stage bending filter design achieves stability by adding lag to the first structural frequency to phase stabilize the first flex mode while gain stabilizing the higher modes. The upper stage bending filter design gain stabilizes all the flex bending modes. The flight control system designs provided here have been demonstrated to provide stable first and second stage control systems in both Draper Ares Stability Analysis Tool (ASAT) and the MSFC 6DOF nonlinear time domain simulation.

Does surface roughness influence the primary stability of acetabular cups? A numerical and experimental biomechanical evaluation.

PubMed

Le Cann, Sophie; Galland, Alexandre; Rosa, Benoît; Le Corroller, Thomas; Pithioux, Martine; Argenson, Jean-Noël; Chabrand, Patrick; Parratte, Sébastien

2014-09-01

Most acetabular cups implanted today are press-fit impacted cementless. Anchorage begins with the primary stability given by insertion of a slightly oversized cup. This primary stability is key to obtaining bone ingrowth and secondary stability. We tested the hypothesis that primary stability of the cup is related to surface roughness of the implant, using both an experimental and a numerical models to analyze how three levels of surface roughness (micro, macro and combined) affect the primary stability of the cup. We also investigated the effect of differences in diameter between the cup and its substrate, and of insertion force, on the cups' primary stability. The results of our study show that primary stability depends on the surface roughness of the cup. The presence of macro-roughness on the peripheral ring is found to decrease primary stability; there was excessive abrasion of the substrate, damaging it and leading to poor primary stability. Numerical modeling indicates that oversizing the cup compared to its substrate has an impact on primary stability, as has insertion force. Copyright © 2014 IPEM. Published by Elsevier Ltd. All rights reserved.

An extended continuum model considering optimal velocity change with memory and numerical tests

NASA Astrophysics Data System (ADS)

Qingtao, Zhai; Hongxia, Ge; Rongjun, Cheng

2018-01-01

In this paper, an extended continuum model of traffic flow is proposed with the consideration of optimal velocity changes with memory. The new model's stability condition and KdV-Burgers equation considering the optimal velocities change with memory are deduced through linear stability theory and nonlinear analysis, respectively. Numerical simulation is carried out to study the extended continuum model, which explores how optimal velocity changes with memory affected velocity, density and energy consumption. Numerical results show that when considering the effects of optimal velocity changes with memory, the traffic jams can be suppressed efficiently. Both the memory step and sensitivity parameters of optimal velocity changes with memory will enhance the stability of traffic flow efficiently. Furthermore, numerical results demonstrates that the effect of optimal velocity changes with memory can avoid the disadvantage of historical information, which increases the stability of traffic flow on road, and so it improve the traffic flow stability and minimize cars' energy consumptions.

Composite Robust $$H_\\infty$$ Control for Uncertain Stochastic Nonlinear Systems with State Delay via Disturbance Observer

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

Liu, Yunlong; Wang, Hong; Guo, Lei

Here in this note, the robust stochastic stabilization and robust H_infinity control problems are investigated for uncertain stochastic time-delay systems with nonlinearity and multiple disturbances. By estimating the disturbance, which can be described by an exogenous model, a composite hierarchical control scheme is proposed that integrates the output of the disturbance observer with state feedback control law. Sufficient conditions for the existence of the disturbance observer and composite hierarchical controller are established in terms of linear matrix inequalities, which ensure the mean-square asymptotic stability of the resulting closed-loop system and the disturbance attenuation. It has been shown that the disturbancemore » rejection performance can also be achieved. A numerical example is provided to show the potential of the proposed techniques and encouraging results have been obtained.« less

Robust adaptive cruise control of high speed trains.

PubMed

Faieghi, Mohammadreza; Jalali, Aliakbar; Mashhadi, Seyed Kamal-e-ddin Mousavi

2014-03-01

The cruise control problem of high speed trains in the presence of unknown parameters and external disturbances is considered. In particular a Lyapunov-based robust adaptive controller is presented to achieve asymptotic tracking and disturbance rejection. The system under consideration is nonlinear, MIMO and non-minimum phase. To deal with the limitations arising from the unstable zero-dynamics we do an output redefinition such that the zero-dynamics with respect to new outputs becomes stable. Rigorous stability analyses are presented which establish the boundedness of all the internal states and simultaneously asymptotic stability of the tracking error dynamics. The results are presented for two common configurations of high speed trains, i.e. the DD and PPD designs, based on the multi-body model and are verified by several numerical simulations. Copyright © 2013 ISA. Published by Elsevier Ltd. All rights reserved.

Composite Robust $$H_\\infty$$ Control for Uncertain Stochastic Nonlinear Systems with State Delay via Disturbance Observer

DOE PAGES

Liu, Yunlong; Wang, Hong; Guo, Lei

2018-03-26

Here in this note, the robust stochastic stabilization and robust H_infinity control problems are investigated for uncertain stochastic time-delay systems with nonlinearity and multiple disturbances. By estimating the disturbance, which can be described by an exogenous model, a composite hierarchical control scheme is proposed that integrates the output of the disturbance observer with state feedback control law. Sufficient conditions for the existence of the disturbance observer and composite hierarchical controller are established in terms of linear matrix inequalities, which ensure the mean-square asymptotic stability of the resulting closed-loop system and the disturbance attenuation. It has been shown that the disturbancemore » rejection performance can also be achieved. A numerical example is provided to show the potential of the proposed techniques and encouraging results have been obtained.« less

Stability Limits and Dynamics of Nonaxisymmetric Liquid Bridges

NASA Technical Reports Server (NTRS)

Alexander, J. Iwan D.

1996-01-01

Theoretical and experimental investigation of the stability of nonaxisymmetric and axisymmetric bridges contained between equal and unequal radii disks as a function of Bond and Weber number with emphasis on the transition from unstable axisymmetric to stable nonaxisymmetric shapes. Numerical analysis of the stability of nonaxisymmetric bridges between unequal disks for various orientations of the gravity vector Experimental and numerical investigation of bridge stability (nonaxisymmetric and axisymmetric), large amplitude (nonaxisymmetric) oscillations and breaking.

Reliability-Based Stability Analysis of Rock Slopes Using Numerical Analysis and Response Surface Method

NASA Astrophysics Data System (ADS)

Dadashzadeh, N.; Duzgun, H. S. B.; Yesiloglu-Gultekin, N.

2017-08-01

While advanced numerical techniques in slope stability analysis are successfully used in deterministic studies, they have so far found limited use in probabilistic analyses due to their high computation cost. The first-order reliability method (FORM) is one of the most efficient probabilistic techniques to perform probabilistic stability analysis by considering the associated uncertainties in the analysis parameters. However, it is not possible to directly use FORM in numerical slope stability evaluations as it requires definition of a limit state performance function. In this study, an integrated methodology for probabilistic numerical modeling of rock slope stability is proposed. The methodology is based on response surface method, where FORM is used to develop an explicit performance function from the results of numerical simulations. The implementation of the proposed methodology is performed by considering a large potential rock wedge in Sumela Monastery, Turkey. The accuracy of the developed performance function to truly represent the limit state surface is evaluated by monitoring the slope behavior. The calculated probability of failure is compared with Monte Carlo simulation (MCS) method. The proposed methodology is found to be 72% more efficient than MCS, while the accuracy is decreased with an error of 24%.

Robust control design with real parameter uncertainty using absolute stability theory. Ph.D. Thesis

NASA Technical Reports Server (NTRS)

How, Jonathan P.; Hall, Steven R.

1993-01-01

The purpose of this thesis is to investigate an extension of mu theory for robust control design by considering systems with linear and nonlinear real parameter uncertainties. In the process, explicit connections are made between mixed mu and absolute stability theory. In particular, it is shown that the upper bounds for mixed mu are a generalization of results from absolute stability theory. Both state space and frequency domain criteria are developed for several nonlinearities and stability multipliers using the wealth of literature on absolute stability theory and the concepts of supply rates and storage functions. The state space conditions are expressed in terms of Riccati equations and parameter-dependent Lyapunov functions. For controller synthesis, these stability conditions are used to form an overbound of the H2 performance objective. A geometric interpretation of the equivalent frequency domain criteria in terms of off-axis circles clarifies the important role of the multiplier and shows that both the magnitude and phase of the uncertainty are considered. A numerical algorithm is developed to design robust controllers that minimize the bound on an H2 cost functional and satisfy an analysis test based on the Popov stability multiplier. The controller and multiplier coefficients are optimized simultaneously, which avoids the iteration and curve-fitting procedures required by the D-K procedure of mu synthesis. Several benchmark problems and experiments on the Middeck Active Control Experiment at M.I.T. demonstrate that these controllers achieve good robust performance and guaranteed stability bounds.

Analysis of stability for stochastic delay integro-differential equations.

PubMed

Zhang, Yu; Li, Longsuo

2018-01-01

In this paper, we concern stability of numerical methods applied to stochastic delay integro-differential equations. For linear stochastic delay integro-differential equations, it is shown that the mean-square stability is derived by the split-step backward Euler method without any restriction on step-size, while the Euler-Maruyama method could reproduce the mean-square stability under a step-size constraint. We also confirm the mean-square stability of the split-step backward Euler method for nonlinear stochastic delay integro-differential equations. The numerical experiments further verify the theoretical results.

Comparison between numerical and analytical results on the required rf current for stabilizing neoclassical tearing modes

NASA Astrophysics Data System (ADS)

Wang, Xiaojing; Yu, Qingquan; Zhang, Xiaodong; Zhang, Yang; Zhu, Sizheng; Wang, Xiaoguang; Wu, Bin

2018-04-01

Numerical studies on the stabilization of neoclassical tearing modes (NTMs) by electron cyclotron current drive (ECCD) have been carried out based on reduced MHD equations, focusing on the amount of the required driven current for mode stabilization and the comparison with analytical results. The dependence of the minimum driven current required for NTM stabilization on some parameters, including the bootstrap current density, radial width of the driven current, radial deviation of the driven current from the resonant surface, and the island width when applying ECCD, are studied. By fitting the numerical results, simple expressions for these dependences are obtained. Analysis based on the modified Rutherford equation (MRE) has also been carried out, and the corresponding results have the same trend as numerical ones, while a quantitative difference between them exists. This difference becomes smaller when the applied radio frequency (rf) current is smaller.

Truncation effect on Taylor-Aris dispersion in lattice Boltzmann schemes: Accuracy towards stability

NASA Astrophysics Data System (ADS)

Ginzburg, Irina; Roux, Laetitia

2015-10-01

The Taylor dispersion in parabolic velocity field provides a well-known benchmark for advection-diffusion (ADE) schemes and serves as a first step towards accurate modeling of the high-order non-Gaussian effects in heterogeneous flow. While applying the Lattice Boltzmann ADE two-relaxation-times (TRT) scheme for a transport with given Péclet number (Pe) one should select six free-tunable parameters, namely, (i) molecular-diffusion-scale, equilibrium parameter; (ii) three families of equilibrium weights, assigned to the terms of mass, velocity and numerical-diffusion-correction, and (iii) two relaxation rates. We analytically and numerically investigate the respective roles of all these degrees of freedom in the accuracy and stability in the evolution of a Gaussian plume. For this purpose, the third- and fourth-order transient multi-dimensional analysis of the recurrence equations of the TRT ADE scheme is extended for a spatially-variable velocity field. The key point is in the coupling of the truncation and Taylor dispersion analysis which allows us to identify the second-order numerical correction δkT to Taylor dispersivity coefficient kT. The procedure is exemplified for a straight Poiseuille flow where δkT is given in a closed analytical form in equilibrium and relaxation parameter spaces. The predicted longitudinal dispersivity is in excellent agreement with the numerical experiments over a wide parameter range. In relatively small Pe-range, the relative dispersion error increases with Péclet number. This deficiency reduces in the intermediate and high Pe-range where it becomes Pe-independent and velocity-amplitude independent. Eliminating δkT by a proper parameter choice and employing specular reflection for zero flux condition on solid boundaries, the d2Q9 TRT ADE scheme may reproduce the Taylor-Aris result quasi-exactly, from very coarse to fine grids, and from very small to arbitrarily high Péclet numbers. Since free-tunable product of two eigenfunctions also controls stability of the model, the validity of the analytically established von Neumann stability diagram is examined in Poiseuille profile. The simplest coordinate-stencil subclass, which is the d2Q5 TRT bounce-back scheme, demonstrates the best performance and achieves the maximum accuracy for most stable relaxation parameters.

Study report on guidelines and test procedures for investigating stability of nonlinear cardiovascular control system models

NASA Technical Reports Server (NTRS)

Fitzjerrell, D. G.

1974-01-01

A general study of the stability of nonlinear as compared to linear control systems is presented. The analysis is general and, therefore, applies to other types of nonlinear biological control systems as well as the cardiovascular control system models. Both inherent and numerical stability are discussed for corresponding analytical and graphic methods and numerical methods.

Finite element analysis and simulation of rheological properties of bulk molding compound (BMC)

NASA Astrophysics Data System (ADS)

Ergin, M. Fatih; Aydin, Ismail

2013-12-01

Bulk molding compound (BMC) is one of the important composite materials with various engineering applications. BMC is a thermoset plastic resin blend of various inert fillers, fiber reinforcements, catalysts, stabilizers and pigments that form a viscous, molding compound. Depending on the end-use application, bulk molding compounds are formulated to achieve close dimensional control, flame and scratch resistance, electrical insulation, corrosion and stain resistance, superior mechanical properties, low shrink and color stability. Its excellent flow characteristics, dielectric properties, and flame resistance make this thermoset material well-suited to a wide variety of applications requiring precision in detail and dimensions as well as high performance. When a BMC is used for these purposes, the rheological behavior and properties of the BMC is the main concern. In this paper, finite element analysis of rheological properties of bulk molding composite material was studied. For this purpose, standard samples of composite material were obtained by means of uniaxial hot pressing. 3 point flexural tests were then carried out by using a universal testing machine. Finite element analyses were then performed with defined material properties within a specific constitutive material behavior. Experimental and numerical results were then compared. Good correlation between the numerical simulation and the experimental results was obtained. It was expected with this study that effects of various process parameters and boundary conditions on the rheological behavior of bulk molding compounds could be determined by means of numerical analysis without detailed experimental work.

Stable scalable control of soliton propagation in broadband nonlinear optical waveguides

NASA Astrophysics Data System (ADS)

Peleg, Avner; Nguyen, Quan M.; Huynh, Toan T.

2017-02-01

We develop a method for achieving scalable transmission stabilization and switching of N colliding soliton sequences in optical waveguides with broadband delayed Raman response and narrowband nonlinear gain-loss. We show that dynamics of soliton amplitudes in N-sequence transmission is described by a generalized N-dimensional predator-prey model. Stability and bifurcation analysis for the predator-prey model are used to obtain simple conditions on the physical parameters for robust transmission stabilization as well as on-off and off-on switching of M out of N soliton sequences. Numerical simulations for single-waveguide transmission with a system of N coupled nonlinear Schrödinger equations with 2 ≤ N ≤ 4 show excellent agreement with the predator-prey model's predictions and stable propagation over significantly larger distances compared with other broadband nonlinear single-waveguide systems. Moreover, stable on-off and off-on switching of multiple soliton sequences and stable multiple transmission switching events are demonstrated by the simulations. We discuss the reasons for the robustness and scalability of transmission stabilization and switching in waveguides with broadband delayed Raman response and narrowband nonlinear gain-loss, and explain their advantages compared with other broadband nonlinear waveguides.

De novo self-assembling collagen heterotrimers using explicit positive and negative design.

PubMed

Xu, Fei; Zhang, Lei; Koder, Ronald L; Nanda, Vikas

2010-03-23

We sought to computationally design model collagen peptides that specifically associate as heterotrimers. Computational design has been successfully applied to the creation of new protein folds and functions. Despite the high abundance of collagen and its key role in numerous biological processes, fibrous proteins have received little attention as computational design targets. Collagens are composed of three polypeptide chains that wind into triple helices. We developed a discrete computational model to design heterotrimer-forming collagen-like peptides. Stability and specificity of oligomerization were concurrently targeted using a combined positive and negative design approach. The sequences of three 30-residue peptides, A, B, and C, were optimized to favor charge-pair interactions in an ABC heterotrimer, while disfavoring the 26 competing oligomers (i.e., AAA, ABB, BCA). Peptides were synthesized and characterized for thermal stability and triple-helical structure by circular dichroism and NMR. A unique A:B:C-type species was not achieved. Negative design was partially successful, with only A + B and B + C competing mixtures formed. Analysis of computed versus experimental stabilities helps to clarify the role of electrostatics and secondary-structure propensities determining collagen stability and to provide important insight into how subsequent designs can be improved.

A hive of activity: the future of public health.

PubMed

Kisely, S; Jones, J

1995-07-01

In consultation with primary and secondary care, public health has a crucial role in improving the health of the population through health promotion, advice on appropriateness and effectiveness of services and communicable disease control. The speciality of public health medicine can only achieve this goal if it is allowed a degree of stability after the numerous organisational changes of the last two decades. The 'beehive model' of public health may offer an organisational structure that is adaptable to further health service change, reflects the rise in multidisciplinary working, guarantees freedom of speech and is sustainable within existing manpower numbers.

Numerical methods of solving a system of multi-dimensional nonlinear equations of the diffusion type

NASA Technical Reports Server (NTRS)

Agapov, A. V.; Kolosov, B. I.

1979-01-01

The principles of conservation and stability of difference schemes achieved using the iteration control method were examined. For the schemes obtained of the predictor-corrector type, the conversion was proved for the control sequences of approximate solutions to the precise solutions in the Sobolev metrics. Algorithms were developed for reducing the differential problem to integral relationships, whose solution methods are known, were designed. The algorithms for the problem solution are classified depending on the non-linearity of the diffusion coefficients, and practical recommendations for their effective use are given.

Adaptive extended-state observer-based fault tolerant attitude control for spacecraft with reaction wheels

NASA Astrophysics Data System (ADS)

Ran, Dechao; Chen, Xiaoqian; de Ruiter, Anton; Xiao, Bing

2018-04-01

This study presents an adaptive second-order sliding control scheme to solve the attitude fault tolerant control problem of spacecraft subject to system uncertainties, external disturbances and reaction wheel faults. A novel fast terminal sliding mode is preliminarily designed to guarantee that finite-time convergence of the attitude errors can be achieved globally. Based on this novel sliding mode, an adaptive second-order observer is then designed to reconstruct the system uncertainties and the actuator faults. One feature of the proposed observer is that the design of the observer does not necessitate any priori information of the upper bounds of the system uncertainties and the actuator faults. In view of the reconstructed information supplied by the designed observer, a second-order sliding mode controller is developed to accomplish attitude maneuvers with great robustness and precise tracking accuracy. Theoretical stability analysis proves that the designed fault tolerant control scheme can achieve finite-time stability of the closed-loop system, even in the presence of reaction wheel faults and system uncertainties. Numerical simulations are also presented to demonstrate the effectiveness and superiority of the proposed control scheme over existing methodologies.

Simultaneous stability and sensitivity in model cortical networks is achieved through anti-correlations between the in- and out-degree of connectivity

PubMed Central

Vasquez, Juan C.; Houweling, Arthur R.; Tiesinga, Paul

2013-01-01

Neuronal networks in rodent barrel cortex are characterized by stable low baseline firing rates. However, they are sensitive to the action potentials of single neurons as suggested by recent single-cell stimulation experiments that reported quantifiable behavioral responses in response to short spike trains elicited in single neurons. Hence, these networks are stable against internally generated fluctuations in firing rate but at the same time remain sensitive to similarly-sized externally induced perturbations. We investigated stability and sensitivity in a simple recurrent network of stochastic binary neurons and determined numerically the effects of correlation between the number of afferent (“in-degree”) and efferent (“out-degree”) connections in neurons. The key advance reported in this work is that anti-correlation between in-/out-degree distributions increased the stability of the network in comparison to networks with no correlation or positive correlations, while being able to achieve the same level of sensitivity. The experimental characterization of degree distributions is difficult because all pre-synaptic and post-synaptic neurons have to be identified and counted. We explored whether the statistics of network motifs, which requires the characterization of connections between small subsets of neurons, could be used to detect evidence for degree anti-correlations. We find that the sample frequency of the 3-neuron “ring” motif (1→2→3→1), can be used to detect degree anti-correlation for sub-networks of size 30 using about 50 samples, which is of significance because the necessary measurements are achievable experimentally in the near future. Taken together, we hypothesize that barrel cortex networks exhibit degree anti-correlations and specific network motif statistics. PMID:24223550

Study on C-S and P-R EOS in pseudo-potential lattice Boltzmann model for two-phase flows

NASA Astrophysics Data System (ADS)

Peng, Yong; Mao, Yun Fei; Wang, Bo; Xie, Bo

Equations of State (EOS) is crucial in simulating multiphase flows by the pseudo-potential lattice Boltzmann method (LBM). In the present study, the Peng and Robinson (P-R) and Carnahan and Starling (C-S) EOS in the pseudo-potential LBM with Exact Difference Method (EDM) scheme for two-phase flows have been compared. Both of P-R and C-S EOS have been used to study the two-phase separation, surface tension, the maximum two-phase density ratio and spurious currents. The study shows that both of P-R and C-S EOS agree with the analytical solutions although P-R EOS may perform better. The prediction of liquid phase by P-R EOS is more accurate than that of air phase and the contrary is true for C-S EOS. Predictions by both of EOS conform with the Laplace’s law. Besides, adjustment of surface tension is achieved by adjusting T. The P-R EOS can achieve larger maximum density ratio than C-S EOS under the same τ. Besides, no matter the C-S EOS or the P-R EOS, if τ tends to 0.5, the computation is prone to numerical instability. The maximum spurious current for P-R is larger than that of C-S. The multiple-relaxation-time LBM still can improve obviously the numerical stability and can achieve larger maximum density ratio.

Matrix-product-operator approach to the nonequilibrium steady state of driven-dissipative quantum arrays

NASA Astrophysics Data System (ADS)

Mascarenhas, Eduardo; Flayac, Hugo; Savona, Vincenzo

2015-08-01

We develop a numerical procedure to efficiently model the nonequilibrium steady state of one-dimensional arrays of open quantum systems based on a matrix-product operator ansatz for the density matrix. The procedure searches for the null eigenvalue of the Liouvillian superoperator by sweeping along the system while carrying out a partial diagonalization of the single-site stationary problem. It bears full analogy to the density-matrix renormalization-group approach to the ground state of isolated systems, and its numerical complexity scales as a power law with the bond dimension. The method brings considerable advantage when compared to the integration of the time-dependent problem via Trotter decomposition, as it can address arbitrarily long-ranged couplings. Additionally, it ensures numerical stability in the case of weakly dissipative systems thanks to a slow tuning of the dissipation rates along the sweeps. We have tested the method on a driven-dissipative spin chain, under various assumptions for the Hamiltonian, drive, and dissipation parameters, and compared the results to those obtained both by Trotter dynamics and Monte Carlo wave function methods. Accurate and numerically stable convergence was always achieved when applying the method to systems with a gapped Liouvillian and a nondegenerate steady state.

Effects of Cavity Configurations on Flameholding and Performances of Kerosene Fueled Scramjet Combustor

NASA Astrophysics Data System (ADS)

Shi, Deyong; Song, Wenyan; Wang, Yuhang; Wang, Yanhua

2017-08-01

In this work, the effects of cavity flameholder configurations on flameholding and performances of kerosene fueled scramjet combustor were studied experimentally and numerically. For experiments, a directly connected ground facility was used and clean high enthalpy air, with a total temperature of 800 K and a total pressure of 800 Kpa, was provided by an electricity resistance heater. To investigate the effects of cavity configurations on flameholding capacity and reacting-flow characteristics, three different flameholders, one single cavity flameholder and two tandem cavity flameholders, were used in experiments. For the two combustors with tandem cavity flameholders, the location and configurations of its up-stream cavity were same with the single cavity flameholder, and the length-to-depth ratios for down-stream cavities were 9 and 11 respectively. The experimental results showed that stabilize kerosene combustion were achieved for combustor with tandem cavity flameholders mounted, and none for that with single cavity flameholder. The none-reacting and reacting flows of combustor models with tandem cavity flameholders were compared and studied with numerical and experimental results. The results showed that higher combustion efficiencies and pressure recovery ratios were achieved for the combustor with down-stream cavity length-to-depth ratio of 9.

Determination of Global Stability of the Slosh Motion in a Spacecraft via Num Erical Experiment

NASA Astrophysics Data System (ADS)

Kang, Ja-Young

2003-12-01

The global stability of the attitude motion of a spin-stabilized space vehicle is investigated by performing numerical experiment. In the previous study, a stationary solution and a particular resonant condition for a given model were found by using analytical method but failed to represent the system stability over parameter values near and off the stationary points. Accordingly, as an extension of the previous work, this study performs numerical experiment to investigate the stability of the system across the parameter space and determines stable and unstable regions of the design parameters of the system.

Numerical modelling on stabilizing large magnetic island by RF current for disruption avoidance

NASA Astrophysics Data System (ADS)

Wang, Xiaojing; Yu, Qingquan; Zhang, Xiaodong; Zhu, Sizheng; Wang, Xiaoguang; Wu, Bin

2018-01-01

Numerical modelling on tearing mode stabilization by RF current due to electron cyclotron current drive (ECCD) has been carried out for the purposes of disruption avoidance, focusing on stabilizing the magnetic island which can grow to a large width and therefore, might cause plasma disruption. When the island has become large, a threshold in driven current for fully stabilizing the mode is found; below this threshold, the island width only slightly decreases. The island’s O-point shifts radially towards the magnetic axis as the mode grows, as a result, applying ECCD at the minor radius of the island’s O-point has a stronger effect than that at the original equilibrium rational surface for stabilizing a large island. During the island growth, the required driven current for mode stabilization increases with the island’s width, indicating that it is more effective to apply ECCD as early as possible for disruption avoidance, as observed in experiments. The numerical results have been compared with those obtained from the modified Rutherford equation.

Models for small-scale structure on cosmic strings. II. Scaling and its stability

NASA Astrophysics Data System (ADS)

Vieira, J. P. P.; Martins, C. J. A. P.; Shellard, E. P. S.

2016-11-01

We make use of the formalism described in a previous paper [Martins et al., Phys. Rev. D 90, 043518 (2014)] to address general features of wiggly cosmic string evolution. In particular, we highlight the important role played by poorly understood energy loss mechanisms and propose a simple Ansatz which tackles this problem in the context of an extended velocity-dependent one-scale model. We find a general procedure to determine all the scaling solutions admitted by a specific string model and study their stability, enabling a detailed comparison with future numerical simulations. A simpler comparison with previous Goto-Nambu simulations supports earlier evidence that scaling is easier to achieve in the matter era than in the radiation era. In addition, we also find that the requirement that a scaling regime be stable seems to notably constrain the allowed range of energy loss parameters.

Robust fuzzy control subject to state variance and passivity constraints for perturbed nonlinear systems with multiplicative noises.

PubMed

Chang, Wen-Jer; Huang, Bo-Jyun

2014-11-01

The multi-constrained robust fuzzy control problem is investigated in this paper for perturbed continuous-time nonlinear stochastic systems. The nonlinear system considered in this paper is represented by a Takagi-Sugeno fuzzy model with perturbations and state multiplicative noises. The multiple performance constraints considered in this paper include stability, passivity and individual state variance constraints. The Lyapunov stability theory is employed to derive sufficient conditions to achieve the above performance constraints. By solving these sufficient conditions, the contribution of this paper is to develop a parallel distributed compensation based robust fuzzy control approach to satisfy multiple performance constraints for perturbed nonlinear systems with multiplicative noises. At last, a numerical example for the control of perturbed inverted pendulum system is provided to illustrate the applicability and effectiveness of the proposed multi-constrained robust fuzzy control method. Copyright © 2014 ISA. Published by Elsevier Ltd. All rights reserved.

Symmetry breaking in occupation number based slave-particle methods

NASA Astrophysics Data System (ADS)

Georgescu, Alexandru B.; Ismail-Beigi, Sohrab

2017-10-01

We describe a theoretical approach to finding spontaneously symmetry-broken electronic phases due to strong electronic interactions when using recently developed slave-particle (slave-boson) approaches based on occupation numbers. We describe why, to date, spontaneous symmetry breaking has proven difficult to achieve in such approaches. We then provide a total energy based approach for introducing auxiliary symmetry-breaking fields into the solution of the slave-particle problem that leads to lowered total energies for symmetry-broken phases. We point out that not all slave-particle approaches yield energy lowering: the slave-particle model being used must explicitly describe the degrees of freedom that break symmetry. Finally, our total energy approach permits us to greatly simplify the formalism used to achieve a self-consistent solution between spinon and slave modes while increasing the numerical stability and greatly speeding up the calculations.

Non-hydrostatic semi-elastic hybrid-coordinate SISL extension of HIRLAM. Part II: numerical testing

NASA Astrophysics Data System (ADS)

Rõõm, Rein; Männik, Aarne; Luhamaa, Andres; Zirk, Marko

2007-10-01

The semi-implicit semi-Lagrangian (SISL), two-time-level, non-hydrostatic numerical scheme, based on the non-hydrostatic, semi-elastic pressure-coordinate equations, is tested in model experiments with flow over given orography (elliptical hill, mountain ridge, system of successive ridges) in a rectangular domain with emphasis on the numerical accuracy and non-hydrostatic effect presentation capability. Comparison demonstrates good (in strong primary wave generation) to satisfactory (in weak secondary wave reproduction in some cases) consistency of the numerical modelling results with known stationary linear test solutions. Numerical stability of the developed model is investigated with respect to the reference state choice, modelling dynamics of a stationary front. The horizontally area-mean reference temperature proves to be the optimal stability warrant. The numerical scheme with explicit residual in the vertical forcing term becomes unstable for cross-frontal temperature differences exceeding 30 K. Stability is restored, if the vertical forcing is treated implicitly, which enables to use time steps, comparable with the hydrostatic SISL.

Analysis of ELM stability with extended MHD models in JET, JT-60U and future JT-60SA tokamak plasmas

NASA Astrophysics Data System (ADS)

Aiba, N.; Pamela, S.; Honda, M.; Urano, H.; Giroud, C.; Delabie, E.; Frassinetti, L.; Lupelli, I.; Hayashi, N.; Huijsmans, G.; JET Contributors, the; Research Unit, JT-60SA

2018-01-01

The stability with respect to a peeling-ballooning mode (PBM) was investigated numerically with extended MHD simulation codes in JET, JT-60U and future JT-60SA plasmas. The MINERVA-DI code was used to analyze the linear stability, including the effects of rotation and ion diamagnetic drift ({ω }* {{i}}), in JET-ILW and JT-60SA plasmas, and the JOREK code was used to simulate nonlinear dynamics with rotation, viscosity and resistivity in JT-60U plasmas. It was validated quantitatively that the ELM trigger condition in JET-ILW plasmas can be reasonably explained by taking into account both the rotation and {ω }* {{i}} effects in the numerical analysis. When deuterium poloidal rotation is evaluated based on neoclassical theory, an increase in the effective charge of plasma destabilizes the PBM because of an acceleration of rotation and a decrease in {ω }* {{i}}. The difference in the amount of ELM energy loss in JT-60U plasmas rotating in opposite directions was reproduced qualitatively with JOREK. By comparing the ELM affected areas with linear eigenfunctions, it was confirmed that the difference in the linear stability property, due not to the rotation direction but to the plasma density profile, is thought to be responsible for changing the ELM energy loss just after the ELM crash. A predictive study to determine the pedestal profiles in JT-60SA was performed by updating the EPED1 model to include the rotation and {ω }* {{i}} effects in the PBM stability analysis. It was shown that the plasma rotation predicted with the neoclassical toroidal viscosity degrades the pedestal performance by about 10% by destabilizing the PBM, but the pressure pedestal height will be high enough to achieve the target parameters required for the ITER-like shape inductive scenario in JT-60SA.

SToRM: A Model for 2D environmental hydraulics

USGS Publications Warehouse

Simões, Francisco J. M.

2017-01-01

A two-dimensional (depth-averaged) finite volume Godunov-type shallow water model developed for flow over complex topography is presented. The model, SToRM, is based on an unstructured cell-centered finite volume formulation and on nonlinear strong stability preserving Runge-Kutta time stepping schemes. The numerical discretization is founded on the classical and well established shallow water equations in hyperbolic conservative form, but the convective fluxes are calculated using auto-switching Riemann and diffusive numerical fluxes. Computational efficiency is achieved through a parallel implementation based on the OpenMP standard and the Fortran programming language. SToRM’s implementation within a graphical user interface is discussed. Field application of SToRM is illustrated by utilizing it to estimate peak flow discharges in a flooding event of the St. Vrain Creek in Colorado, U.S.A., in 2013, which reached 850 m3/s (~30,000 f3 /s) at the location of this study.

Advanced Computational Modeling of Vapor Deposition in a High-Pressure Reactor

NASA Technical Reports Server (NTRS)

Cardelino, Beatriz H.; Moore, Craig E.; McCall, Sonya D.; Cardelino, Carlos A.; Dietz, Nikolaus; Bachmann, Klaus

2004-01-01

In search of novel approaches to produce new materials for electro-optic technologies, advances have been achieved in the development of computer models for vapor deposition reactors in space. Numerical simulations are invaluable tools for costly and difficult processes, such as those experiments designed for high pressures and microgravity conditions. Indium nitride is a candidate compound for high-speed laser and photo diodes for optical communication system, as well as for semiconductor lasers operating into the blue and ultraviolet regions. But InN and other nitride compounds exhibit large thermal decomposition at its optimum growth temperature. In addition, epitaxy at lower temperatures and subatmospheric pressures incorporates indium droplets into the InN films. However, surface stabilization data indicate that InN could be grown at 900 K in high nitrogen pressures, and microgravity could provide laminar flow conditions. Numerical models for chemical vapor deposition have been developed, coupling complex chemical kinetics with fluid dynamic properties.

Advanced Computational Modeling of Vapor Deposition in a High-pressure Reactor

NASA Technical Reports Server (NTRS)

Cardelino, Beatriz H.; Moore, Craig E.; McCall, Sonya D.; Cardelino, Carlos A.; Dietz, Nikolaus; Bachmann, Klaus

2004-01-01

In search of novel approaches to produce new materials for electro-optic technologies, advances have been achieved in the development of computer models for vapor deposition reactors in space. Numerical simulations are invaluable tools for costly and difficult processes, such as those experiments designed for high pressures and microgravity conditions. Indium nitride is a candidate compound for high-speed laser and photo diodes for optical communication system, as well as for semiconductor lasers operating into the blue and ultraviolet regions. But InN and other nitride compounds exhibit large thermal decomposition at its optimum growth temperature. In addition, epitaxy at lower temperatures and subatmospheric pressures incorporates indium droplets into the InN films. However, surface stabilization data indicate that InN could be grown at 900 K in high nitrogen pressures, and microgravity could provide laminar flow conditions. Numerical models for chemical vapor deposition have been developed, coupling complex chemical kinetics with fluid dynamic properties.

Summary of FY15 results of benchmark modeling activities

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

Arguello, J. Guadalupe

2015-08-01

Sandia is participating in the third phase of an is a contributing partner to a U.S.-German "Joint Project" entitled "Comparison of current constitutive models and simulation procedures on the basis of model calculations of the thermo-mechanical behavior and healing of rock salt." The first goal of the project is to check the ability of numerical modeling tools to correctly describe the relevant deformation phenomena in rock salt under various influences. Achieving this goal will lead to increased confidence in the results of numerical simulations related to the secure storage of radioactive wastes in rock salt, thereby enhancing the acceptance ofmore » the results. These results may ultimately be used to make various assertions regarding both the stability analysis of an underground repository in salt, during the operating phase, and the long-term integrity of the geological barrier against the release of harmful substances into the biosphere, in the post-operating phase.« less

Counterrotating prop-fan simulations which feature a relative-motion multiblock grid decomposition enabling arbitrary time-steps

NASA Technical Reports Server (NTRS)

Janus, J. Mark; Whitfield, David L.

1990-01-01

Improvements are presented of a computer algorithm developed for the time-accurate flow analysis of rotating machines. The flow model is a finite volume method utilizing a high-resolution approximate Riemann solver for interface flux definitions. The numerical scheme is a block LU implicit iterative-refinement method which possesses apparent unconditional stability. Multiblock composite gridding is used to orderly partition the field into a specified arrangement of blocks exhibiting varying degrees of similarity. Block-block relative motion is achieved using local grid distortion to reduce grid skewness and accommodate arbitrary time step selection. A general high-order numerical scheme is applied to satisfy the geometric conservation law. An even-blade-count counterrotating unducted fan configuration is chosen for a computational study comparing solutions resulting from altering parameters such as time step size and iteration count. The solutions are compared with measured data.

Ionic liquid and nanoparticle hybrid systems: Emerging applications.

PubMed

He, Zhiqi; Alexandridis, Paschalis

2017-06-01

Having novel electronic and optical properties that emanate from their nano-scale dimensions, nanoparticles are central to numerous applications. Ionic liquids can confer to nanoparticle chemical protection and physicochemical property enhancement through intermolecular interactions and can consequently improve the stability and reusability of nanoparticle for various operations. With an aim to combine the novel properties of nanoparticles and ionic liquids, different structures have been generated, based on a balance of several intermolecular interactions. Such ionic liquid and nanoparticle hybrids are showing great potential in diverse applications. In this review, we first introduce various types of ionic liquid and nanoparticle hybrids, including nanoparticle colloidal dispersions in ionic liquids, ionic liquid-grafted nanoparticles, and nanoparticle-stabilized ionic liquid-based emulsions. Such hybrid materials exhibit interesting synergisms. We then highlight representative applications of ionic liquid and nanoparticle hybrids in the catalysis, electrochemistry and separations fields. Such hybrids can attain better stability and higher efficiency under a broad range of conditions. Novel and enhanced performance can be achieved in these applications by combining desired properties of ionic liquids and of nanoparticles within an appropriate hybrid nanostructure. Copyright © 2016 Elsevier B.V. All rights reserved.

Craton stability and continental lithosphere dynamics during plume-plate interaction

NASA Astrophysics Data System (ADS)

Wang, H.; Van Hunen, J.; Pearson, D.

2013-12-01

Survival of thick cratonic roots in a vigorously convecting mantle system for billions of years has long been studied by the geodynamical community. A high cratonic root strength is generally considered to be the most important factor. We first perform and discuss new numerical models to investigate craton stability in both Newtonian and non-Newtonian rheology in the stagnant lid regime. The results show that only a modest compositional rheological factor of Δη=10 with non-Newtonian rheology is required for the survival of cratonic roots in a stagnant lid regime. A larger rheological factor (100 or more) is needed to maintain similar craton longevity in a Newtonian rheology environment. Furthermore, chemical buoyancy plays an important role on craton stability and its evolution, but could only work with suitable compositional rheology. During their long lifespan, cratons experienced a suite of dynamic, tectonothermal events, such as nearby subduction and mantle plume activity. Cratonic nuclei are embedded in shorter-lived, more vulnerable continental areas of different thickness, composition and rheology, which would influence the lithosphere dynamic when tectonothermal events happen nearby. South Africa provides a very good example to investigate such dynamic processes as it hosts several cratons and there are many episodic thermal events since the Mesozoic as indicated by a spectrum of magmatic activity. We numerically investigate such an integrated system using the topographic evolution of cratons and surrounding lithosphere as a diagnostic observable. The post-70Ma thinning of pericratonic lithosphere by ~50km around Kaapvaal craton (Mather et al., 2011) is also investigated through our numerical models. The results show that the pericratonic lithosphere cools and grows faster than cratons do, but is also more likely to be effected by episodic thermal events. This leads to surface topography change that is significantly larger around the craton than within the craton itself. Given the considerable debate on the uplift history of southern African plateau (Nyblade and Sleep, 2003), our numerical models that encompass lithospheric heterogeneity within cratons could help to achieve a better understanding of this issue.

Designing Adaptive Low-Dissipative High Order Schemes for Long-Time Integrations. Chapter 1

NASA Technical Reports Server (NTRS)

Yee, Helen C.; Sjoegreen, B.; Mansour, Nagi N. (Technical Monitor)

2001-01-01

A general framework for the design of adaptive low-dissipative high order schemes is presented. It encompasses a rather complete treatment of the numerical approach based on four integrated design criteria: (1) For stability considerations, condition the governing equations before the application of the appropriate numerical scheme whenever it is possible; (2) For consistency, compatible schemes that possess stability properties, including physical and numerical boundary condition treatments, similar to those of the discrete analogue of the continuum are preferred; (3) For the minimization of numerical dissipation contamination, efficient and adaptive numerical dissipation control to further improve nonlinear stability and accuracy should be used; and (4) For practical considerations, the numerical approach should be efficient and applicable to general geometries, and an efficient and reliable dynamic grid adaptation should be used if necessary. These design criteria are, in general, very useful to a wide spectrum of flow simulations. However, the demand on the overall numerical approach for nonlinear stability and accuracy is much more stringent for long-time integration of complex multiscale viscous shock/shear/turbulence/acoustics interactions and numerical combustion. Robust classical numerical methods for less complex flow physics are not suitable or practical for such applications. The present approach is designed expressly to address such flow problems, especially unsteady flows. The minimization of employing very fine grids to overcome the production of spurious numerical solutions and/or instability due to under-resolved grids is also sought. The incremental studies to illustrate the performance of the approach are summarized. Extensive testing and full implementation of the approach is forthcoming. The results shown so far are very encouraging.

High-power hybrid plasma spraying of large yttria-stabilized zirconia powder

NASA Astrophysics Data System (ADS)

Huang, Heji; Eguchi, Keisuke; Yoshida, Toyonobu

2006-03-01

To testify to the advantage of large ceramic powder spraying, numerical simulations and experimental studies on the behavior of large yttria-stabilized zirconia (YSZ) powder in a high-power hybrid plasma spraying process have been carried out. Numeric predictions and experimental results showed that, with the high radio frequency (RF) input power of 100 kW, the most refractory YSZ powder with particle sizes as large as 88 μm could be fully melted and well-flattened splats could be formed. A large degree of flattening (ξ) of 4.7 has been achieved. The improved adhesive strength between the large splat and the substrate was confirmed based on the measurement of the crack density inside of the splats. A thick YSZ coating >300 μm was successfully deposited on a large CoNiCrAlY-coated Inconel substrate (50×50×4 mm in size). The ultradense microstructure without clear boundaries between the splats and the clean and crack-free interface between the top-coat and the bond-coat also indicate the good adhesion. These results showed that highpower hybrid plasma spraying of large ceramic powder is a very promising process for deposition of highquality coatings, especially in the application of thermal barrier coatings (TBCs).

Stability of the mode-locking regime in tapered quantum-dot lasers

NASA Astrophysics Data System (ADS)

Bardella, P.; Drzewietzki, L.; Rossetti, M.; Weber, C.; Breuer, S.

2018-02-01

We study numerically and experimentally the role of the injection current and reverse bias voltage on the pulse stability of tapered, passively mode-locked, Quantum Dot (QD) lasers. By using a multi-section delayed differential equation and introducing in the model the QD inhomogenous broadening, we are able to predict the onset of leading and trailing edge instabilities in the emitted pulse trains and to identify specific trends of stability in dependence on the laser biasing conditions. The numerical results are confirmed experimentally trough amplitude and timing stability analysis of the pulses.

Numerical solution of the general coupled nonlinear Schrödinger equations on unbounded domains.

PubMed

Li, Hongwei; Guo, Yue

2017-12-01

The numerical solution of the general coupled nonlinear Schrödinger equations on unbounded domains is considered by applying the artificial boundary method in this paper. In order to design the local absorbing boundary conditions for the coupled nonlinear Schrödinger equations, we generalize the unified approach previously proposed [J. Zhang et al., Phys. Rev. E 78, 026709 (2008)PLEEE81539-375510.1103/PhysRevE.78.026709]. Based on the methodology underlying the unified approach, the original problem is split into two parts, linear and nonlinear terms, and we then achieve a one-way operator to approximate the linear term to make the wave out-going, and finally we combine the one-way operator with the nonlinear term to derive the local absorbing boundary conditions. Then we reduce the original problem into an initial boundary value problem on the bounded domain, which can be solved by the finite difference method. The stability of the reduced problem is also analyzed by introducing some auxiliary variables. Ample numerical examples are presented to verify the accuracy and effectiveness of our proposed method.

An Laudau-Lifschitz theory based algorithm on calculating post-buckling configuration of a rod buckling in elastic media

NASA Astrophysics Data System (ADS)

Huang, Shicheng; Tan, Likun; Hu, Nan; Grover, Hannah; Chu, Kevin; Chen, Zi

This reserach introduces a new numerical approach of calculating the post-buckling configuration of a thin rod embedded in elastic media. The theoretical base is the governing ODEs describing the balance of forces and moments, the length conservation, and the physics of bending and twisting by Laudau and Lifschitz. The numerical methods applied in the calculation are continuation method and Newton's method of iteration in combination with spectrum method. To the authors' knowledge, it is the first trial of directly applying the L-L theory to numerically studying the phenomenon of rod buckling in elastic medium. This method accounts for nonlinearity of geometry, thus is capable of calculating large deformation. The stability of this method is another advantage achieved by expressing the governing equations in a set of first-order derivative form. The wave length, amplitude, and decay effect all agree with the experiment without any further assumptions. This program can be applied to different occasions with varying stiffness of the elastic medai and rigidity of the rod.

Turing patterns in parabolic systems of conservation laws and numerically observed stability of periodic waves

NASA Astrophysics Data System (ADS)

Barker, Blake; Jung, Soyeun; Zumbrun, Kevin

2018-03-01

Turing patterns on unbounded domains have been widely studied in systems of reaction-diffusion equations. However, up to now, they have not been studied for systems of conservation laws. Here, we (i) derive conditions for Turing instability in conservation laws and (ii) use these conditions to find families of periodic solutions bifurcating from uniform states, numerically continuing these families into the large-amplitude regime. For the examples studied, numerical stability analysis suggests that stable periodic waves can emerge either from supercritical Turing bifurcations or, via secondary bifurcation as amplitude is increased, from subcritical Turing bifurcations. This answers in the affirmative a question of Oh-Zumbrun whether stable periodic solutions of conservation laws can occur. Determination of a full small-amplitude stability diagram - specifically, determination of rigorous Eckhaus-type stability conditions - remains an interesting open problem.

Complete synchronization of uncertain chaotic systems via a single proportional adaptive controller: A comparative study

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

Ahmad, Israr, E-mail: iak-2000plus@yahoo.com; Saaban, Azizan Bin, E-mail: azizan.s@uum.edu.my; Ibrahim, Adyda Binti, E-mail: adyda@uum.edu.my

This paper addresses a comparative computational study on the synchronization quality, cost and converging speed for two pairs of identical chaotic and hyperchaotic systems with unknown time-varying parameters. It is assumed that the unknown time-varying parameters are bounded. Based on the Lyapunov stability theory and using the adaptive control method, a single proportional controller is proposed to achieve the goal of complete synchronizations. Accordingly, appropriate adaptive laws are designed to identify the unknown time-varying parameters. The designed control strategy is easy to implement in practice. Numerical simulations results are provided to verify the effectiveness of the proposed synchronization scheme.

Asynchronous variational integration using continuous assumed gradient elements.

PubMed

Wolff, Sebastian; Bucher, Christian

2013-03-01

Asynchronous variational integration (AVI) is a tool which improves the numerical efficiency of explicit time stepping schemes when applied to finite element meshes with local spatial refinement. This is achieved by associating an individual time step length to each spatial domain. Furthermore, long-term stability is ensured by its variational structure. This article presents AVI in the context of finite elements based on a weakened weak form (W2) Liu (2009) [1], exemplified by continuous assumed gradient elements Wolff and Bucher (2011) [2]. The article presents the main ideas of the modified AVI, gives implementation notes and a recipe for estimating the critical time step.

Structural analysis of an off-grid tiny house

NASA Astrophysics Data System (ADS)

Calluari, Karina Arias; Alonso-Marroquín, Fernando

2017-06-01

The off-grid technologies and tiny house movement have experimented an unprecedented growth in recent years. Putting both sides together, we are trying to achieve an economic and environmental friendly solution to the higher cost of residential properties. This solution is the construction of off-grid tiny houses. This article presents a design for a small modular off-grid house made by pine timber. A numerical analysis of the proposed tiny house was performed to ensure its structural stability. The results were compared with the suggested serviceability limit state criteria, which are contended in the Australia Guidelines Standards making this design reliable for construction.

Generalized Projective Synchronization between Two Complex Networks with Time-Varying Coupling Delay

NASA Astrophysics Data System (ADS)

Sun, Mei; Zeng, Chang-Yan; Tian, Li-Xin

2009-01-01

Generalized projective synchronization (GPS) between two complex networks with time-varying coupling delay is investigated. Based on the Lyapunov stability theory, a nonlinear controller and adaptive updated laws are designed. Feasibility of the proposed scheme is proven in theory. Moreover, two numerical examples are presented, using the energy resource system and Lü's system [Physica A 382 (2007) 672] as the nodes of the networks. GPS between two energy resource complex networks with time-varying coupling delay is achieved. This study can widen the application range of the generalized synchronization methods and will be instructive for the demand-supply of energy resource in some regions of China.

Generalized Hill-stability criteria for hierarchical three-body systems at arbitrary inclinations

NASA Astrophysics Data System (ADS)

Grishin, Evgeni; Perets, Hagai B.; Zenati, Yossef; Michaely, Erez

2017-04-01

A fundamental aspect of the three-body problem is its stability. Most stability studies have focused on the co-planar three-body problem, deriving analytic criteria for the dynamical stability of such pro/retrograde systems. Numerical studies of inclined systems phenomenologically mapped their stability regions, but neither complement it by theoretical framework, nor provided satisfactory fit for their dependence on mutual inclinations. Here we present a novel approach to study the stability of hierarchical three-body systems at arbitrary inclinations, which accounts not only for the instantaneous stability of such systems, but also for the secular stability and evolution through Lidov-Kozai cycles and evection. We generalize the Hill-stability criteria to arbitrarily inclined triple systems, explain the existence of quasi-stable regimes and characterize the inclination dependence of their stability. We complement the analytic treatment with an extensive numerical study, to test our analytic results. We find excellent correspondence up to high inclinations (˜120°), beyond which the agreement is marginal. At such high inclinations, the stability radius is larger, the ratio between the outer and inner periods becomes comparable and our secular averaging approach is no longer strictly valid. We therefore combine our analytic results with polynomial fits to the numerical results to obtain a generalized stability formula for triple systems at arbitrary inclinations. Besides providing a generalized secular-based physical explanation for the stability of non-co-planar systems, our results have direct implications for any triple systems and, in particular, binary planets and moon/satellite systems; we briefly discuss the latter as a test case for our models.

Predictive Relation between Early Numerical Competencies and Mathematics Achievement in First Grade Portuguese Children.

PubMed

Marcelino, Lilia; de Sousa, Óscar; Lopes, António

2017-01-01

Early numerical competencies (ENC) (counting, number relations, and basic arithmetic operations) have a central position in the initial learning of mathematics, and their assessment is useful for predicting later mathematics achievement. Using a regression model, this study aims to analyze the correlational and predictive evidence between ENC and mathematics achievement in first grade Portuguese children ( n = 123). The children's ENC were examined at the point of school entry. Three criterion groups (low, moderate, and high ENC) were formed based on the results of the early numerical brief screener and mathematics achievement measured at the end of first grade. The following hypotheses were tested: children who started first grade with low numerical competencies remained low mathematics achievement at the end of first grade; and children who started with high numerical competencies, finished the first grade with high mathematics achievement. The results showed that ENC contributed to a significant amount of explained variance in mathematics achievement at the end of the first grade. Children with low numerical competencies performed lower than children with moderate and high numerical competencies. Findings suggest that ENC are meaningful for predicting first-grade mathematics difficulties.

Predictive Relation between Early Numerical Competencies and Mathematics Achievement in First Grade Portuguese Children

PubMed Central

Marcelino, Lilia; de Sousa, Óscar; Lopes, António

2017-01-01

Early numerical competencies (ENC) (counting, number relations, and basic arithmetic operations) have a central position in the initial learning of mathematics, and their assessment is useful for predicting later mathematics achievement. Using a regression model, this study aims to analyze the correlational and predictive evidence between ENC and mathematics achievement in first grade Portuguese children (n = 123). The children’s ENC were examined at the point of school entry. Three criterion groups (low, moderate, and high ENC) were formed based on the results of the early numerical brief screener and mathematics achievement measured at the end of first grade. The following hypotheses were tested: children who started first grade with low numerical competencies remained low mathematics achievement at the end of first grade; and children who started with high numerical competencies, finished the first grade with high mathematics achievement. The results showed that ENC contributed to a significant amount of explained variance in mathematics achievement at the end of the first grade. Children with low numerical competencies performed lower than children with moderate and high numerical competencies. Findings suggest that ENC are meaningful for predicting first-grade mathematics difficulties. PMID:28713308

A modified Holly-Preissmann scheme for simulating sharp concentration fronts in streams with steep velocity gradients using RIV1Q

NASA Astrophysics Data System (ADS)

Liu, Zhao-wei; Zhu, De-jun; Chen, Yong-can; Wang, Zhi-gang

2014-12-01

RIV1Q is the stand-alone water quality program of CE-QUAL-RIV1, a hydraulic and water quality model developed by U.S. Army Corps of Engineers Waterways Experiment Station. It utilizes an operator-splitting algorithm and the advection term in governing equation is treated using the explicit two-point, fourth-order accurate, Holly-Preissmann scheme, in order to preserve numerical accuracy for advection of sharp gradients in concentration. In the scheme, the spatial derivative of the transport equation, where the derivative of velocity is included, is introduced to update the first derivative of dependent variable. In the stream with larger cross-sectional variation, steep velocity gradient can be easily found and should be estimated correctly. In the original version of RIV1Q, however, the derivative of velocity is approximated by a finite difference which is first-order accurate. Its leading truncation error leads to the numerical error of concentration which is related with the velocity and concentration gradients and increases with the decreasing Courant number. The simulation may also be unstable when a sharp velocity drop occurs. In the present paper, the derivative of velocity is estimated with a modified second-order accurate scheme and the corresponding numerical error of concentration decreases. Additionally, the stability of the simulation is improved. The modified scheme is verified with a hypothetical channel case and the results demonstrate that satisfactory accuracy and stability can be achieved even when the Courant number is very low. Finally, the applicability of the modified scheme is discussed.

Stabilization of dynamics of oscillatory systems by nonautonomous perturbation.

PubMed

Lucas, Maxime; Newman, Julian; Stefanovska, Aneta

2018-04-01

Synchronization and stability under periodic oscillatory driving are well understood, but little is known about the effects of aperiodic driving, despite its abundance in nature. Here, we consider oscillators subject to driving with slowly varying frequency, and investigate both short-term and long-term stability properties. For a phase oscillator, we find that, counterintuitively, such variation is guaranteed to enlarge the Arnold tongue in parameter space. Using analytical and numerical methods that provide information on time-variable dynamical properties, we find that the growth of the Arnold tongue is specifically due to the growth of a region of intermittent synchronization where trajectories alternate between short-term stability and short-term neutral stability, giving rise to stability on average. We also present examples of higher-dimensional nonlinear oscillators where a similar stabilization phenomenon is numerically observed. Our findings help support the case that in general, deterministic nonautonomous perturbation is a very good candidate for stabilizing complex dynamics.

Stabilization of dynamics of oscillatory systems by nonautonomous perturbation

NASA Astrophysics Data System (ADS)

Lucas, Maxime; Newman, Julian; Stefanovska, Aneta

2018-04-01

Synchronization and stability under periodic oscillatory driving are well understood, but little is known about the effects of aperiodic driving, despite its abundance in nature. Here, we consider oscillators subject to driving with slowly varying frequency, and investigate both short-term and long-term stability properties. For a phase oscillator, we find that, counterintuitively, such variation is guaranteed to enlarge the Arnold tongue in parameter space. Using analytical and numerical methods that provide information on time-variable dynamical properties, we find that the growth of the Arnold tongue is specifically due to the growth of a region of intermittent synchronization where trajectories alternate between short-term stability and short-term neutral stability, giving rise to stability on average. We also present examples of higher-dimensional nonlinear oscillators where a similar stabilization phenomenon is numerically observed. Our findings help support the case that in general, deterministic nonautonomous perturbation is a very good candidate for stabilizing complex dynamics.

A High-Authority/Low-Authority Control Strategy for Coupled Aircraft-Style Bays

NASA Technical Reports Server (NTRS)

Schiller, N. H.; Fuller, C. R.; Cabell, R. H.

2006-01-01

This paper presents a numerical investigation of an active structural acoustic control strategy for coupled aircraft-style bays. While structural coupling can destabilize or limit the performance of some model-based decentralized control systems, fullycoupled centralized control strategies are impractical for typical aircraft containing several hundred bays. An alternative is to use classical rate feedback with matched, collocated transducer pairs to achieve active damping. Unfortunately, due to the conservative nature of this strategy, stability is guaranteed at the expense of achievable noise reduction. Therefore, this paper describes the development of a combined control strategy using robust active damping in addition to a high-authority controller based on linear quadratic Gaussian (LQG) theory. The combined control system is evaluated on a tensioned, two-bay model using piezoceramic actuators and ideal point velocity sensors. Transducer placement on the two-bay structure is discussed, and the advantages of a combined control strategy are presented.

Decentralized Control of Sound Radiation using a High-Authority/Low-Authority Control Strategy with Anisotropic Actuators

NASA Technical Reports Server (NTRS)

Schiller, Noah H.; Cabell, Randolph H.; Fuller, Chris R.

2008-01-01

This paper describes a combined control strategy designed to reduce sound radiation from stiffened aircraft-style panels. The control architecture uses robust active damping in addition to high-authority linear quadratic Gaussian (LQG) control. Active damping is achieved using direct velocity feedback with triangularly shaped anisotropic actuators and point velocity sensors. While active damping is simple and robust, stability is guaranteed at the expense of performance. Therefore the approach is often referred to as low-authority control. In contrast, LQG control strategies can achieve substantial reductions in sound radiation. Unfortunately, the unmodeled interaction between neighboring control units can destabilize decentralized control systems. Numerical simulations show that combining active damping and decentralized LQG control can be beneficial. In particular, augmenting the in-bandwidth damping supplements the performance of the LQG control strategy and reduces the destabilizing interaction between neighboring control units.

Numerical Simulation And Experimental Investigation Of The Lift-Off And Blowout Of Enclosed Laminar Flames

NASA Technical Reports Server (NTRS)

Venuturmilli, Rajasekhar; Zhang, Yong; Chen, Lea-Der

2003-01-01

Enclosed flames are found in many industrial applications such as power plants, gas-turbine combustors and jet engine afterburners. A better understanding of the burner stability limits can lead to development of combustion systems that extend the lean and rich limits of combustor operations. This paper reports a fundamental study of the stability limits of co-flow laminar jet diffusion flames. A numerical study was conducted that used an adaptive mesh refinement scheme in the calculation. Experiments were conducted in two test rigs with two different fuels and diluted with three inert species. The numerical stability limits were compared with microgravity experimental data. Additional normal-gravity experimental results were also presented.

Micromachined silicon acoustic delay line with 3D-printed micro linkers and tapered input for improved structural stability and acoustic directivity

NASA Astrophysics Data System (ADS)

Cho, Y.; Kumar, A.; Xu, S.; Zou, J.

2016-10-01

Recent studies have shown that micromachined silicon acoustic delay lines can provide a promising solution to achieve real-time photoacoustic tomography without the need for complex transducer arrays and data acquisition electronics. To achieve deeper imaging depth and wider field of view, a longer delay time and therefore delay length are required. However, as the length of the delay line increases, it becomes more vulnerable to structural instability due to reduced mechanical stiffness. In this paper, we report the design, fabrication, and testing of a new silicon acoustic delay line enhanced with 3D printed polymer micro linker structures. First, mechanical deformation of the silicon acoustic delay line (with and without linker structures) under gravity was simulated by using finite element method. Second, the acoustic crosstalk and acoustic attenuation caused by the polymer micro linker structures were evaluated with both numerical simulation and ultrasound transmission testing. The result shows that the use of the polymer micro linker structures significantly improves the structural stability of the silicon acoustic delay lines without creating additional acoustic attenuation and crosstalk. In addition, the improvement of the acoustic acceptance angle of the silicon acoustic delay lines was also investigated to better suppress the reception of unwanted ultrasound signals outside of the imaging plane. These two improvements are expected to provide an effective solution to eliminate current limitations on the achievable acoustic delay time and out-of-plane imaging resolution of micromachined silicon acoustic delay line arrays.

Guidelines for LTS magnet design based on transient stability

NASA Astrophysics Data System (ADS)

Seo, Kazutaka; Morita, Masao

2006-05-01

Stabilities of low critical temperature superconducting (LTS) magnets and their designs are studied and discussed. There are two contradictory necessities; those are low cost and high performance, in the other words, high magnetic field and large current density. Especially, the maximum magnetic fields of the latest high performance Nb 3Sn magnets are around 20 T. Mentioned necessities result in the small stability margins. Needless to say, the superconducting magnet must produce its nominal field reliably. Therefore, maintaining adequate stability margin, the magnet design to draw out the high potential of the superconductor is required. The transient stability of the superconducting magnet is determined by the relationship between mechanical disturbance energy and stability margin. The minimum quench energy (MQE) is one of the index of stability margin and it is defined as the minimum energy to trigger quenching of a superconductor. MQE should be beyond any possible disturbance energy during the operation. It is difficult to identify the mechanical disturbance energy quantitatively. On the contrary, MQE had been evaluated precisely by means of our developed resistive carbon paste heater (CPH). At the same time, we can predict MQE by numerical simulations. Because the magnet comes to quench if the mechanical disturbance exceeds the MQE, the disturbance energies are suspected to be equivalent to MQEs during the magnet-training. When we achieved somewhat larger MQE, we may exclude numbers of training quenches. In this paper, we discuss the guidelines of LTS magnet design from the standpoint of MQE. We represent some case studies for various superconducting magnets and/or some different winding methods.

An optimal implicit staggered-grid finite-difference scheme based on the modified Taylor-series expansion with minimax approximation method for elastic modeling

NASA Astrophysics Data System (ADS)

Yang, Lei; Yan, Hongyong; Liu, Hong

2017-03-01

Implicit staggered-grid finite-difference (ISFD) scheme is competitive for its great accuracy and stability, whereas its coefficients are conventionally determined by the Taylor-series expansion (TE) method, leading to a loss in numerical precision. In this paper, we modify the TE method using the minimax approximation (MA), and propose a new optimal ISFD scheme based on the modified TE (MTE) with MA method. The new ISFD scheme takes the advantage of the TE method that guarantees great accuracy at small wavenumbers, and keeps the property of the MA method that keeps the numerical errors within a limited bound at the same time. Thus, it leads to great accuracy for numerical solution of the wave equations. We derive the optimal ISFD coefficients by applying the new method to the construction of the objective function, and using a Remez algorithm to minimize its maximum. Numerical analysis is made in comparison with the conventional TE-based ISFD scheme, indicating that the MTE-based ISFD scheme with appropriate parameters can widen the wavenumber range with high accuracy, and achieve greater precision than the conventional ISFD scheme. The numerical modeling results also demonstrate that the MTE-based ISFD scheme performs well in elastic wave simulation, and is more efficient than the conventional ISFD scheme for elastic modeling.

Stability of Nonlinear Principal Components Analysis: An Empirical Study Using the Balanced Bootstrap

ERIC Educational Resources Information Center

Linting, Marielle; Meulman, Jacqueline J.; Groenen, Patrick J. F.; van der Kooij, Anita J.

2007-01-01

Principal components analysis (PCA) is used to explore the structure of data sets containing linearly related numeric variables. Alternatively, nonlinear PCA can handle possibly nonlinearly related numeric as well as nonnumeric variables. For linear PCA, the stability of its solution can be established under the assumption of multivariate…

Numerical Boundary Condition Procedures

NASA Technical Reports Server (NTRS)

1981-01-01

Topics include numerical procedures for treating inflow and outflow boundaries, steady and unsteady discontinuous surfaces, far field boundaries, and multiblock grids. In addition, the effects of numerical boundary approximations on stability, accuracy, and convergence rate of the numerical solution are discussed.

Do biodegradable magnesium alloy intramedullary interlocking nails prematurely lose fixation stability in the treatment of tibial fracture? A numerical simulation.

PubMed

Wang, Haosen; Hao, Zhixiu; Wen, Shizhu

2017-01-01

Intramedullary interlocking nailing is an effective technique used to treat long bone fractures. Recently, biodegradable metals have drawn increased attention as an intramedullary interlocking nailing material. In this study, numerical simulations were implemented to determine whether the degradation rate of magnesium alloy makes it a suitable material for manufacturing biodegradable intramedullary interlocking nails. Mechano-regulatory and bone-remodeling models were used to simulate the fracture healing process, and a surface corrosion model was used to simulate intramedullary rod degradation. The results showed that magnesium alloy intramedullary rods exhibited a satisfactory degradation rate; the fracture healed and callus enhancement was observed before complete dissolution of the intramedullary rod. Delayed magnesium degradation (using surface coating techniques) did not confer a significant advantage over the non-delayed degradation process; immediate degradation also achieved satisfactory healing outcomes. However, delayed degradation had no negative effect on callus enhancement, as it did not cause signs of stress shielding. To avoid risks of individual differences such as delayed union, delayed degradation is recommended. Although the magnesium intramedullary rod did not demonstrate rapid degradation, its ability to provide high fixation stiffness to achieve earlier load bearing was inferior to that of the conventional titanium alloy and stainless steel rods. Therefore, light physiological loads should be ensured during the early stages of healing to achieve bony healing; otherwise, with increased loading and degraded intramedullary rods, the fracture may ultimately fail to heal. Copyright © 2016 Elsevier Ltd. All rights reserved.

Experimental Results of OH Regime Investigation in Globus-M Spherical Torus

NASA Astrophysics Data System (ADS)

Golant, Victor; Gusev, Vasily; Levin, Roman; Petrov, Yuriy; Sakharov, Nikolay

2001-10-01

Plasma parameters were measured in novel spherical torus Globus-M in highly shaped plasmas with aspect ratio, A > 1.5, elongation, k < 1.9, triangularity < 0.5. Plasma column was created by direct induction method with the currents up to Ip 0.3 MA in the magnetic field, Bt - 0.08 - 0.5 T. In Globus-M spherical torus plasma column is closely fitted into the vacuum vessel and wall conditioning technology described in [1] was used to achieve good plasma performance. Plasma experiments were focused around achievement of ultimate OH regimes allowed by power supplies. The operational limits of the device were investigated. In the regime with extreme low q(cy1) < 1 and high normalized current > 4, the plasma current of almost 100kA was sustained transiently in low magnetic field 800 Gs. The first results on stability analysis with numerical code are presented. The runaway electrons behavior was studied in spherical tokamak conditions. Influence of plasma current and density ramp-up speeds, MHD events on plasma performance and stability was demonstrated. Magnetic reconstruction was performed with EFIT version adopted for PC simulations. Plans for auxiliary heating and current drive are discussed. 1. V.K. Gusev, …, V.E. Golant, et al., Nucl. Fusion 41, No 7, (2001), to be published

Optimal Control-Based Adaptive NN Design for a Class of Nonlinear Discrete-Time Block-Triangular Systems.

PubMed

Liu, Yan-Jun; Tong, Shaocheng

2016-11-01

In this paper, we propose an optimal control scheme-based adaptive neural network design for a class of unknown nonlinear discrete-time systems. The controlled systems are in a block-triangular multi-input-multi-output pure-feedback structure, i.e., there are both state and input couplings and nonaffine functions to be included in every equation of each subsystem. The design objective is to provide a control scheme, which not only guarantees the stability of the systems, but also achieves optimal control performance. The main contribution of this paper is that it is for the first time to achieve the optimal performance for such a class of systems. Owing to the interactions among subsystems, making an optimal control signal is a difficult task. The design ideas are that: 1) the systems are transformed into an output predictor form; 2) for the output predictor, the ideal control signal and the strategic utility function can be approximated by using an action network and a critic network, respectively; and 3) an optimal control signal is constructed with the weight update rules to be designed based on a gradient descent method. The stability of the systems can be proved based on the difference Lyapunov method. Finally, a numerical simulation is given to illustrate the performance of the proposed scheme.

Instability of a solidifying binary mixture

NASA Technical Reports Server (NTRS)

Antar, B. N.

1982-01-01

An analysis is performed on the stability of a solidifying binary mixture due to surface tension variation of the free liquid surface. The basic state solution is obtained numerically as a nonstationary function of time. Due to the time dependence of the basic state, the stability analysis is of the global type which utilizes a variational technique. Also due to the fact that the basic state is a complex function of both space and time, the stability analysis is performed through numerical means.

An extended macro model accounting for acceleration changes with memory and numerical tests

NASA Astrophysics Data System (ADS)

Cheng, Rongjun; Ge, Hongxia; Sun, Fengxin; Wang, Jufeng

2018-09-01

Considering effect of acceleration changes with memory, an improved continuum model of traffic flow is proposed in this paper. By applying the linear stability theory, we derived the new model's linear stability condition. Through nonlinear analysis, the KdV-Burgers equation is derived to describe the propagating behavior of traffic density wave near the neutral stability line. Numerical simulation is carried out to study the extended traffic flow model, which explores how acceleration changes with memory affected each car's velocity, density and fuel consumption and exhaust emissions. Numerical results demonstrate that acceleration changes with memory have significant negative effect on dynamic characteristic of traffic flow. Furthermore, research results verify that the effect of acceleration changes with memory will deteriorate the stability of traffic flow and increase cars' total fuel consumptions and emissions during the whole evolution of small perturbation.

Stability of the discretization of the electron avalanche phenomenon

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

Villa, Andrea, E-mail: andrea.villa@rse-web.it; Barbieri, Luca, E-mail: luca.barbieri@rse-web.it; Gondola, Marco, E-mail: marco.gondola@rse-web.it

2015-09-01

The numerical simulation of the discharge inception is an active field of applied physics with many industrial applications. In this work we focus on the drift-reaction equation that describes the electron avalanche. This phenomenon is one of the basic building blocks of the streamer model. The main difficulty of the electron avalanche equation lies in the fact that the reaction term is positive when a high electric field is applied. It leads to exponentially growing solutions and this has a major impact on the behavior of numerical schemes. We analyze the stability of a reference finite volume scheme applied tomore » this latter problem. The stability of the method may impose a strict mesh spacing, therefore a proper stabilized scheme, which is stable whatever spacing is used, has been developed. The convergence of the scheme is treated as well as some numerical experiments.« less

Multiple-correction hybrid k-exact schemes for high-order compressible RANS-LES simulations on fully unstructured grids

NASA Astrophysics Data System (ADS)

Pont, Grégoire; Brenner, Pierre; Cinnella, Paola; Maugars, Bruno; Robinet, Jean-Christophe

2017-12-01

A Godunov's type unstructured finite volume method suitable for highly compressible turbulent scale-resolving simulations around complex geometries is constructed by using a successive correction technique. First, a family of k-exact Godunov schemes is developed by recursively correcting the truncation error of the piecewise polynomial representation of the primitive variables. The keystone of the proposed approach is a quasi-Green gradient operator which ensures consistency on general meshes. In addition, a high-order single-point quadrature formula, based on high-order approximations of the successive derivatives of the solution, is developed for flux integration along cell faces. The proposed family of schemes is compact in the algorithmic sense, since it only involves communications between direct neighbors of the mesh cells. The numerical properties of the schemes up to fifth-order are investigated, with focus on their resolvability in terms of number of mesh points required to resolve a given wavelength accurately. Afterwards, in the aim of achieving the best possible trade-off between accuracy, computational cost and robustness in view of industrial flow computations, we focus more specifically on the third-order accurate scheme of the family, and modify locally its numerical flux in order to reduce the amount of numerical dissipation in vortex-dominated regions. This is achieved by switching from the upwind scheme, mostly applied in highly compressible regions, to a fourth-order centered one in vortex-dominated regions. An analytical switch function based on the local grid Reynolds number is adopted in order to warrant numerical stability of the recentering process. Numerical applications demonstrate the accuracy and robustness of the proposed methodology for compressible scale-resolving computations. In particular, supersonic RANS/LES computations of the flow over a cavity are presented to show the capability of the scheme to predict flows with shocks, vortical structures and complex geometries.

Numerical simulations of recent proton acceleration experiments with sub-100 TW laser systems

NASA Astrophysics Data System (ADS)

Sinigardi, Stefano

2016-09-01

Recent experiments carried out at the Italian National Research Center, National Optics Institute Department in Pisa, are showing interesting results regarding maximum proton energies achievable with sub-100 TW laser systems. While laser systems are being continuously upgraded in laboratories around the world, at the same time a new trend on stabilizing and making ion acceleration results reproducible is growing in importance. Almost all applications require a beam with fixed performance, so that the energy spectrum and the total charge exhibit moderate shot to shot variations. This result is surely far from being achieved, but many paths are being explored in order to reach it. Some of the reasons for this variability come from fluctuations in laser intensity and focusing, due to optics instability. Other variation sources come from small differences in the target structure. The target structure can vary substantially, when it is impacted by the main pulse, due to the prepulse duration and intensity, the shape of the main pulse and the total energy deposited. In order to qualitatively describe the prepulse effect, we will present a two dimensional parametric scan of its relevant parameters. A single case is also analyzed with a full three dimensional simulation, obtaining reasonable agreement between the numerical and the experimental energy spectrum.

The stability of freak waves with regard to external impact and perturbation of initial data

NASA Astrophysics Data System (ADS)

Smirnova, Anna; Shamin, Roman

2014-05-01

We investigate solutions of the equations, describing freak waves, in perspective of stability with regard to external impact and perturbation of initial data. The modeling of freak waves is based on numerical solution of equations describing a non-stationary potential flow of the ideal fluid with a free surface. We consider the two-dimensional infinitely deep flow. For waves modeling we use the equations in conformal variables. The variant of these equations is offered in [1]. Mathematical correctness of these equations was discussed in [2]. These works establish the uniqueness of solutions, offer the effective numerical solution calculation methods, prove the numerical convergence of these methods. The important aspect of numerical modeling of freak waves is the stability of solutions, describing these waves. In this work we study the questions of stability with regards to external impact and perturbation of initial data. We showed the stability of freak waves numerical model, corresponding to the external impact. We performed series of computational experiments with various freak wave initial data and random external impact. This impact means the power density on free surface. In each experiment examine two waves: the wave that was formed by external impact and without one. In all the experiments we see the stability of equation`s solutions. The random external impact practically does not change the time of freak wave formation and its form. Later our work progresses to the investigation of solution's stability under perturbations of initial data. We take the initial data that provide a freak wave and get the numerical solution. In common we take the numerical solution of equation with perturbation of initial data. The computing experiments showed that the freak waves equations solutions are stable under perturbations of initial data.So we can make a conclusion that freak waves are stable relatively external perturbation and perturbation of initial data both. 1. Zakharov V.E., Dyachenko A.I., Vasilyev O.A. New method for numerical simulation of a nonstationary potential flow of incompressible fluid with a free surface// Eur. J.~Mech. B Fluids. 2002. V. 21. P. 283-291. 2. R.V. Shamin. Dynamics of an Ideal Liquid with a Free Surface in Conformal Variables // Journal of Mathematical Sciences, Vol. 160, No. 5, 2009. P. 537-678. 3. R.V. Shamin, V.E. Zakharov, A.I. Dyachenko. How probability for freak wave formation can be found // THE EUROPEAN PHYSICAL JOURNAL - SPECIAL TOPICS Volume 185, Number 1, 113-124, DOI: 10.1140/epjst/e2010-01242-y

A control problem for Burgers' equation with bounded input/output

NASA Technical Reports Server (NTRS)

Burns, John A.; Kang, Sungkwon

1990-01-01

A stabilization problem for Burgers' equation is considered. Using linearization, various controllers are constructed which minimize certain weighted energy functionals. These controllers produce the desired degree of stability for the closed-loop nonlinear system. A numerical scheme for computing the feedback gain functional is developed and several numerical experiments are performed to show the theoretical results.

Stope Stability Assessment and Effect of Horizontal to Vertical Stress Ratio on the Yielding and Relaxation Zones Around Underground Open Stopes Using Empirical and Finite Element Methods

NASA Astrophysics Data System (ADS)

Sepehri, Mohammadali; Apel, Derek; Liu, Wei

2017-09-01

Predicting the stability of open stopes can be a challenging task for underground mine engineers. For decades, the stability graph method has been used as the first step of open stope design around the world. However, there are some shortcomings with this method. For instance, the stability graph method does not account for the relaxation zones around the stopes. Another limitation of the stability graph is that this method cannot to be used to evaluate the stability of the stopes with high walls made of backfill materials. However, there are several analytical and numerical methods that can be used to overcome these limitations. In this study, both empirical and numerical methods have been used to assess the stability of an open stope located between mine levels N9225 and N9250 at Diavik diamond underground mine. It was shown that the numerical methods can be used as complementary methods along with other analytical and empirical methods to assess the stability of open stopes. A three dimensional elastoplastic finite element model was constructed using Abaqus software. In this paper a sensitivity analysis was performed to investigate the impact of the stress ratio "k" on the extent of the yielding and relaxation zones around the hangingwall and footwall of the understudy stope.

Comparison of eigenvectors for coupled seismo-electromagnetic layered-Earth modelling

NASA Astrophysics Data System (ADS)

Grobbe, N.; Slob, E. C.; Thorbecke, J. W.

2016-07-01

We study the accuracy and numerical stability of three eigenvector sets for modelling the coupled poroelastic and electromagnetic layered-Earth response. We use a known eigenvector set, its flux-normalized version and a newly derived flux-normalized set. The new set is chosen such that the system is properly uncoupled when the coupling between the poroelastic and electromagnetic fields vanishes. We carry out two different numerical stability tests: the first test focuses on the internal system, eigenvector and eigenvalue consistency; the second test investigates the stability and preciseness of the flux-normalized systems by looking at identity relations. We find that the known set shows the largest deviation for both tests, whereas the new set performs best. In two additional numerical modelling experiments, these numerical inaccuracies are shown to generate numerical noise levels comparable to small signals, such as signals coming from the important interface conversion responses, especially when the coupling coefficient is small. When coupling vanishes completely, the known set does not produce proper results. The new set produces numerically stable and accurate results in all situations. We therefore strongly recommend to use this newly derived set for future layered-Earth seismo-electromagnetic modelling experiments.

Surface mediated assembly of small, metastable gold nanoclusters.

PubMed

Pettibone, John M; Osborn, William A; Rykaczewski, Konrad; Talin, A Alec; Bonevich, John E; Hudgens, Jeffrey W; Allendorf, Mark D

2013-07-21

The unique properties of metallic nanoclusters are attractive for numerous commercial and industrial applications but are generally less stable than nanocrystals. Thus, developing methodologies for stabilizing nanoclusters and retaining their enhanced functionality is of great interest. We report the assembly of PPh3-protected Au9 clusters from a heterogeneous mixture into films consisting of sub 3 nm nanocluster assemblies. The depositing nanoclusters are metastable in solution, but the resulting nanocluster assemblies are stabilized indefinitely in air or fresh solvent. The films exhibit distinct structure from Au nanoparticles observed by X-ray diffraction, and film dissolution data support the preservation of small nanoclusters. UV-Vis spectroscopy, electrospray ionization mass spectrometry, X-ray photoelectron spectroscopy and electron microscopy are used to elucidate information regarding the nanocluster formation and assembly mechanism. Preferential deposition of nanocluster assemblies can be achieved on multiple substrates, including polymer, Cr, Si, SiO2, SiNx, and metal-organic frameworks (MOFs). Unlike other vapor phase coating processes, nanocluster assembly on the MIL-68(In) MOF crystal is capable of preferentially coating the external surface and stabilizing the crystal structure in hydrothermal conditions, which should enhance their storage, separation and delivery capabilities.

Spontaneous nucleation and topological stabilization of skyrmions in magnetic nanodisks with the interfacial Dzyaloshinskii-Moriya interaction

NASA Astrophysics Data System (ADS)

Kolesnikov, A. G.; Samardak, A. S.; Stebliy, M. E.; Ognev, A. V.; Chebotkevich, L. A.; Sadovnikov, A. V.; Nikitov, S. A.; Kim, Yong Jin; Cha, In Ho; Kim, Young Keun

2017-05-01

One of the major societal challenges is reducing the power consumption of information technology (IT) devices and numerous data centers. Distinct from the current approaches based on switching of magnetic single-domain nanostructures or on movement of domain walls under high currents, an original magnetic skyrmion technology offers ultra-low power, fast, high-density, and scalable spintronic devices, including non-volatile random access memory. Using data-driven micromagnetic simulations, we demonstrate the possibility of spontaneous nucleation and stabilization of different skyrmionic states, such as skyrmions, merons, and meron-like configurations, in heavy metal/ferromagnetic nanodisks with the interfacial Dzyaloshinskii-Moriya interaction (iDMI) as a result of quasi-static magnetization reversal only. Since iDMI is not easily modulated in real systems, we show that skyrmion stabilization is easily achievable by manipulating magnetic anisotropy, saturation magnetization, and the diameters of nanodisks. The state diagrams, presented in terms of the topological charge, allow to explicitly distinguish the intermediate states between skyrmions and merons and can be used for developing a skyrmionic medium, which has been recently proposed to be a building block for future spin-orbitronic devices.

Dynamic response and stability of a gas-lubricated Rayleigh-step pad

NASA Technical Reports Server (NTRS)

Cheng, C.; Cheng, H. S.

1973-01-01

The quasi-static, pressure characteristics of a gas-lubricated thrust bearing with shrouded, Rayleigh-step pads are determined for a time-varying film thickness. The axial response of the thrust bearing to an axial forcing function or an axial rotor disturbance is investigated by treating the gas film as a spring having nonlinear restoring and damping forces. These forces are related to the film thickness by a power relation. The nonlinear equation of motion in the axial mode is solved by the Ritz-Galerkin method as well as the direct, numerical integration. Results of the nonlinear response by both methods are compared with the response based on the linearized equation. Further, the gas-film instability of an infinitely wide Rayleigh step thrust pad is determined by solving the transient Reynolds equation coupled with the equation of the motion of the pad. Results show that the Rayleigh-step geometry is very stable for bearing number A up to 50. The stability threshold is shown to exist only for ultrahigh values of Lambda equal to or greater than 100, where the stability can be achieved by making the mass heavier than the critical mass.

Comparison of numerical results and multicavity purge and rim seal data with extensions to dynamics

NASA Astrophysics Data System (ADS)

Athavale, Mahesh; Przekwas, Andrzej J.; Hendricks, Robert C.; Steinetz, Bruce M.

1995-05-01

The computation of flows within interconnected, multiple-disk cavities shows strong interaction between the cavities and the power stream. For this reason, simulations of single cavities in such cases are not realistic; the complete, linked configuration must be considered. Unsteady flow fields affect engine stability and can engender power-stream-driven secondary flows that produce local hot spotting or general cavity heating. Further, a concentric whirling rotor produces a circumferential pressure wave, but a statically eccentric whirling rotor produces a radial wave; both waves affect cavity ingestion and the stability of the entire engine. It is strongly suggested that seals be used to enhance turbojet engine stability. Simple devices, such as swirl brakes, honeycomb inserts, and new seal configurations, should be considered. The cost effectiveness of the NASA Lewis Research Center seals program can be expressed in terms of program goals (e.g., the Integrated High-pressure/Temperature Engine Technology (IHPTET) cannot be achieved without such a program), cost (savings to $250 million/1-percent decrease in specific fuel consumption), and indirect benefits (reduction of atmospheric NO(x) and CO2 and reduction of powerplant downtime).

Comparison of Numerical Results and Multicavity Purge and Rim Seal Data with Extensions to Dynamics

NASA Technical Reports Server (NTRS)

Athavale, Mahesh; Przekwas, Andrzej J.; Hendricks, Robert C.; Steinetz, Bruce M.

1995-01-01

The computation of flows within interconnected, multiple-disk cavities shows strong interaction between the cavities and the power stream. For this reason, simulations of single cavities in such cases are not realistic; the complete, linked configuration must be considered. Unsteady flow fields affect engine stability and can engender power-stream-driven secondary flows that produce local hot spotting or general cavity heating. Further, a concentric whirling rotor produces a circumferential pressure wave, but a statically eccentric whirling rotor produces a radial wave; both waves affect cavity ingestion and the stability of the entire engine. It is strongly suggested that seals be used to enhance turbojet engine stability. Simple devices, such as swirl brakes, honeycomb inserts, and new seal configurations, should be considered. The cost effectiveness of the NASA Lewis Research Center seals program can be expressed in terms of program goals (e.g., the Integrated High-pressure/Temperature Engine Technology (IHPTET) cannot be achieved without such a program), cost (savings to $250 million/1-percent decrease in specific fuel consumption), and indirect benefits (reduction of atmospheric NO(x) and CO2 and reduction of powerplant downtime).

Common-path digital holographic microscopy based on a beam displacer unit

NASA Astrophysics Data System (ADS)

Di, Jianglei; Zhang, Jiwei; Song, Yu; Wang, Kaiqiang; Wei, Kun; Zhao, Jianlin

2018-02-01

Digital holographic microscopy (DHM) has become a novel tool with advantages of full field, non-destructive, high-resolution and 3D imaging, which captures the quantitative amplitude and phase information of microscopic specimens. It's a well-established method for digital recording and numerical reconstructing the full complex field of wavefront of the samples with a diffraction-limited lateral resolution down to 0.3 μm depending on the numerical aperture of microscope objective. Meanwhile, its axial resolution through axial direction is less than 10 nm due to the interferometric nature in phase imaging. Compared with the typical optical configurations such as Mach-Zehnder interferometer and Michelson interferometer, the common-path DHM has the advantages of simple and compact configuration, high stability, and so on. Here, a simple, compact, and low-cost common-path DHM based on a beam displacer unit is proposed for quantitative phase imaging of biological cells. The beam displacer unit is completely compatible with commercial microscope and can be easily set up in the output port of the microscope as a compact independent device. This technique can be used to achieve the quantitative phase measurement of biological cells with an excellent temporal stability of 0.51 nm, which makes it having a good prospect in the fields of biological and medical science. Living mouse osteoblastic cells are quantitatively measured with the system to demonstrate its capability and applicability.

Numerical and arithmetical cognition: a longitudinal study of process and concept deficits in children with learning disability.

PubMed

Geary, D C; Hamson, C O; Hoard, M K

2000-11-01

Based on the stability and level of performance on standard achievement tests in first and second grade (mean age in first grade = 82 months), children with IQ scores in the low-average to high-average range were classified as learning disabled (LD) in mathematics (MD), reading (RD), or both (MD/RD). These children (n = 42), a group of children who showed variable achievement test performance across grades (n = 16), and a control group of academically normal peers (n = 35) were administered a series of experimental and psychometric tasks. The tasks assessed number comprehension and production skills, counting knowledge, arithmetic skills, working memory, the ease of activation of phonetic representations of words and numbers, and spatial abilities. The children with variable achievement test performance did not differ from the academically normal children in any cognitive domain, whereas the children in the LD groups showed specific patterns of cognitive deficit, above and beyond the influence of IQ. Discussion focuses on the similarities and differences across the groups of LD children. Copyright 2000 Academic Press.

Integrated modeling of high βN steady state scenario on DIII-D

DOE PAGES

Park, Jin Myung; Ferron, J. R.; Holcomb, Christopher T.; ...

2018-01-10

Theory-based integrated modeling validated against DIII-D experiments predicts that fully non-inductive DIII-D operation with β N > 4.5 is possible with certain upgrades. IPS-FASTRAN is a new iterative numerical procedure that integrates models of core transport, edge pedestal, equilibrium, stability, heating, and current drive self-consistently to find steady-state ( d/dt = 0) solutions and reproduces most features of DIII-D high β N discharges with a stationary current profile. Projecting forward to scenarios possible on DIII-D with future upgrades, the high q min > 2 scenario achieves stable operation at β N as high as 5 by using a very broadmore » current density profile to improve the ideal-wall stabilization of low- n instabilities along with confinement enhancement from low magnetic shear. This modeling guides the necessary upgrades of the heating and current drive system to realize reactor-relevant high β N steady-state scenarios on DIII-D by simultaneous optimization of the current and pressure profiles.« less

Approximate message passing for nonconvex sparse regularization with stability and asymptotic analysis

NASA Astrophysics Data System (ADS)

Sakata, Ayaka; Xu, Yingying

2018-03-01

We analyse a linear regression problem with nonconvex regularization called smoothly clipped absolute deviation (SCAD) under an overcomplete Gaussian basis for Gaussian random data. We propose an approximate message passing (AMP) algorithm considering nonconvex regularization, namely SCAD-AMP, and analytically show that the stability condition corresponds to the de Almeida-Thouless condition in spin glass literature. Through asymptotic analysis, we show the correspondence between the density evolution of SCAD-AMP and the replica symmetric (RS) solution. Numerical experiments confirm that for a sufficiently large system size, SCAD-AMP achieves the optimal performance predicted by the replica method. Through replica analysis, a phase transition between replica symmetric and replica symmetry breaking (RSB) region is found in the parameter space of SCAD. The appearance of the RS region for a nonconvex penalty is a significant advantage that indicates the region of smooth landscape of the optimization problem. Furthermore, we analytically show that the statistical representation performance of the SCAD penalty is better than that of \

The interaction between thermodynamic stability and buried free cysteines in regulating the functional half-life of fibroblast growth factor-1.

PubMed

Lee, Jihun; Blaber, Michael

2009-10-16

Protein biopharmaceuticals are an important and growing area of human therapeutics; however, the intrinsic property of proteins to adopt alternative conformations (such as during protein unfolding and aggregation) presents numerous challenges, limiting their effective application as biopharmaceuticals. Using fibroblast growth factor-1 as model system, we describe a cooperative interaction between the intrinsic property of thermostability and the reactivity of buried free-cysteine residues that can substantially modulate protein functional half-life. A mutational strategy that combines elimination of buried free cysteines and secondary mutations that enhance thermostability to achieve a substantial gain in functional half-life is described. Furthermore, the implementation of this design strategy utilizing stabilizing mutations within the core region resulted in a mutant protein that is essentially indistinguishable from wild type as regard protein surface and solvent structure, thus minimizing the immunogenic potential of the mutations. This design strategy should be generally applicable to soluble globular proteins containing buried free-cysteine residues.

Hierarchical Porous Interlocked Polymeric Microcapsules: Sulfonic Acid Functionalization as Acid Catalysts

NASA Astrophysics Data System (ADS)

Wang, Xiaomei; Gu, Jinyan; Tian, Lei; Zhang, Xu

2017-03-01

Owing to their unique structural and surface properties, mesoporous microspheres are widely applied in the catalytic field. Generally, increasing the surface area of the specific active phase of the catalyst is a good method, which can achieve a higher catalytic activity through the fabrication of the corresponding catalytic microspheres with the smaller size and hollow structure. However, one of the major challenges in the use of hollow microspheres (microcapsules) as catalysts is their chemical and structural stability. Herein, the grape-like hypercrosslinked polystyrene hierarchical porous interlocked microcapsule (HPIM-HCL-PS) is fabricated by SiO2 colloidal crystals templates, whose structure is the combination of open mouthed structure, mesoporous nanostructure and interlocked architecture. Numerous microcapsules assembling together and forming the roughly grape-like microcapsule aggregates can enhance the structural stability and recyclability of these microcapsules. After undergoing the sulfonation, the sulfonated HPIM-HCL-PS is served as recyclable acid catalyst for condensation reaction between benzaldehyde and ethylene glycol (TOF = 793 h-1), moreover, exhibits superior activity, selectivity and recyclability.

A new smooth robust control design for uncertain nonlinear systems with non-vanishing disturbances

NASA Astrophysics Data System (ADS)

Xian, Bin; Zhang, Yao

2016-06-01

In this paper, we consider the control problem for a general class of nonlinear system subjected to uncertain dynamics and non-varnishing disturbances. A smooth nonlinear control algorithm is presented to tackle these uncertainties and disturbances. The proposed control design employs the integral of a nonlinear sigmoid function to compensate the uncertain dynamics, and achieve a uniformly semi-global practical asymptotic stable tracking control of the system outputs. A novel Lyapunov-based stability analysis is employed to prove the convergence of the tracking errors and the stability of the closed-loop system. Numerical simulation results on a two-link robot manipulator are presented to illustrate the performance of the proposed control algorithm comparing with the layer-boundary sliding mode controller and the robust of integration of sign of error control design. Furthermore, real-time experiment results for the attitude control of a quadrotor helicopter are also included to confirm the effectiveness of the proposed algorithm.

Multi-disciplinary optimization of aeroservoelastic systems

NASA Technical Reports Server (NTRS)

Karpel, Mordechay

1990-01-01

Efficient analytical and computational tools for simultaneous optimal design of the structural and control components of aeroservoelastic systems are presented. The optimization objective is to achieve aircraft performance requirements and sufficient flutter and control stability margins with a minimal weight penalty and without violating the design constraints. Analytical sensitivity derivatives facilitate an efficient optimization process which allows a relatively large number of design variables. Standard finite element and unsteady aerodynamic routines are used to construct a modal data base. Minimum State aerodynamic approximations and dynamic residualization methods are used to construct a high accuracy, low order aeroservoelastic model. Sensitivity derivatives of flutter dynamic pressure, control stability margins and control effectiveness with respect to structural and control design variables are presented. The performance requirements are utilized by equality constraints which affect the sensitivity derivatives. A gradient-based optimization algorithm is used to minimize an overall cost function. A realistic numerical example of a composite wing with four controls is used to demonstrate the modeling technique, the optimization process, and their accuracy and efficiency.

Multidisciplinary optimization of aeroservoelastic systems using reduced-size models

NASA Technical Reports Server (NTRS)

Karpel, Mordechay

1992-01-01

Efficient analytical and computational tools for simultaneous optimal design of the structural and control components of aeroservoelastic systems are presented. The optimization objective is to achieve aircraft performance requirements and sufficient flutter and control stability margins with a minimal weight penalty and without violating the design constraints. Analytical sensitivity derivatives facilitate an efficient optimization process which allows a relatively large number of design variables. Standard finite element and unsteady aerodynamic routines are used to construct a modal data base. Minimum State aerodynamic approximations and dynamic residualization methods are used to construct a high accuracy, low order aeroservoelastic model. Sensitivity derivatives of flutter dynamic pressure, control stability margins and control effectiveness with respect to structural and control design variables are presented. The performance requirements are utilized by equality constraints which affect the sensitivity derivatives. A gradient-based optimization algorithm is used to minimize an overall cost function. A realistic numerical example of a composite wing with four controls is used to demonstrate the modeling technique, the optimization process, and their accuracy and efficiency.

Integrated modeling of high βN steady state scenario on DIII-D

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

Park, Jin Myung; Ferron, J. R.; Holcomb, Christopher T.

Theory-based integrated modeling validated against DIII-D experiments predicts that fully non-inductive DIII-D operation with β N > 4.5 is possible with certain upgrades. IPS-FASTRAN is a new iterative numerical procedure that integrates models of core transport, edge pedestal, equilibrium, stability, heating, and current drive self-consistently to find steady-state ( d/dt = 0) solutions and reproduces most features of DIII-D high β N discharges with a stationary current profile. Projecting forward to scenarios possible on DIII-D with future upgrades, the high q min > 2 scenario achieves stable operation at β N as high as 5 by using a very broadmore » current density profile to improve the ideal-wall stabilization of low- n instabilities along with confinement enhancement from low magnetic shear. This modeling guides the necessary upgrades of the heating and current drive system to realize reactor-relevant high β N steady-state scenarios on DIII-D by simultaneous optimization of the current and pressure profiles.« less

Integrated modeling of high βN steady state scenario on DIII-D

NASA Astrophysics Data System (ADS)

Park, J. M.; Ferron, J. R.; Holcomb, C. T.; Buttery, R. J.; Solomon, W. M.; Batchelor, D. B.; Elwasif, W.; Green, D. L.; Kim, K.; Meneghini, O.; Murakami, M.; Snyder, P. B.

2018-01-01

Theory-based integrated modeling validated against DIII-D experiments predicts that fully non-inductive DIII-D operation with βN > 4.5 is possible with certain upgrades. IPS-FASTRAN is a new iterative numerical procedure that integrates models of core transport, edge pedestal, equilibrium, stability, heating, and current drive self-consistently to find steady-state (d/dt = 0) solutions and reproduces most features of DIII-D high βN discharges with a stationary current profile. Projecting forward to scenarios possible on DIII-D with future upgrades, the high qmin > 2 scenario achieves stable operation at βN as high as 5 by using a very broad current density profile to improve the ideal-wall stabilization of low-n instabilities along with confinement enhancement from low magnetic shear. This modeling guides the necessary upgrades of the heating and current drive system to realize reactor-relevant high βN steady-state scenarios on DIII-D by simultaneous optimization of the current and pressure profiles.

Feasibility of potable water generators to meet vessel numeric ballast water discharge limits.

PubMed

Albert, Ryan J; Viveiros, Edward; Falatko, Debra S; Tamburri, Mario N

2017-07-15

Ballast water is taken on-board vessels into ballast water tanks to maintain vessel draft, buoyancy, and stability. Unmanaged ballast water contains aquatic organisms that, when transported and discharged to non-native waters, may establish as invasive species. Technologies capable of achieving regulatory limits designed to decrease the likelihood of invasion include onboard ballast water management systems. However, to date, the treatment development and manufacturing marketplace is limited to large vessels with substantial ballast requirements. For smaller vessels or vessels with reduced ballast requirements, we evaluated the feasibility of meeting the discharge limits by generating ballast water using onboard potable water generators. Case studies and parametric analyses demonstrated the architectural feasibility of installing potable water generators onboard actual vessels with minimal impacts for most vessel types evaluated. Furthermore, land-based testing of a potable water generator demonstrated capability to meet current numeric discharge limits for living organisms in all size classes. Published by Elsevier Ltd.

Supratransmission in a metastable modular metastructure for tunable non-reciprocal wave transmission

NASA Astrophysics Data System (ADS)

Wu, Zhen; Wang, K. W.

2018-03-01

In this research, we numerically and analytically investigate the nonlinear energy transmission phenomenon in a metastable modular metastructure. Numerical studies on a 1D metastable chain provide clear evidence that when driving frequency is within the stopband of the periodic structure, there exists a threshold for the driving amplitude, above which sudden increase in the energy transmission can be observed. This onset of transmission is due to nonlinear instability and is known as supratransmission. We discover that due to spatial asymmetry of strategically configured constituents, such transmission thresholds are considerably different when structure is excited from different ends and this discrepancy creates a region of non-reciprocal energy transmission. We demonstrate that when the loss of stability is due to saddlenode bifurcation, the transmission threshold can be predicted analytically using a localized nonlinear-linear system model, and analyzed via combining harmonic balancing and transfer matrix methods. These investigations elucidate the rich and complex dynamics achievable by nonlinearity and metastabilities, and provide synthesize tools for tunable bandgaps and non-reciprocal wave transmissions.

Vortex sheet modeling with higher order curved panels. Ph.D Thesis Final Technical Report

NASA Technical Reports Server (NTRS)

Nagati, M. G.

1985-01-01

A numerical technique is presented for modeling the vortex sheet with a deformable surface definition, along which a continuous vortex strength distribution in the spanwise direction is applied, so that by repeatedly modifying its shape, its true configuration is approached, in the proximity of its generating wing. Design problems requiring the inclusion of a realistic configuration of the vortex sheet are numerous. Examples discussed include: control effectiveness and stability derivatives, longitudinal stability, lateral stability, canards, propellers and helicopter rotors, and trailing vortex hazards.

Stability of sequences generated by nonlinear differential systems. [for analysis of glider jet aircraft motion

NASA Technical Reports Server (NTRS)

Brown, R. L.

1979-01-01

A local stability analysis is presented for both the analytic and numerical solutions of the initial value problem for a system of ordinary differential equations. It is shown that, using a proper choice of Liapunov function, a connected region of stable initial values of both the analytic solution and the one-leg k-step numerical solution can be approximated. Attention is given to the example of the two-dimensional problem involving the stability of the longitudinal equations of motion of a gliding jet aircraft.

Analysis of the Effects of Streamwise Lift Distribution on Sonic Boom Signature

NASA Technical Reports Server (NTRS)

Yoo, Paul

2013-01-01

Investigation of sonic boom has been one of the major areas of study in aeronautics due to the benefits a low-boom aircraft has in both civilian and military applications. This work conducts a numerical analysis of the effects of streamwise lift distribution on the shock coalescence characteristics. A simple wing-canard-stabilator body model is used in the numerical simulation. The streamwise lift distribution is varied by fixing the canard at a deflection angle while trimming the aircraft with the wing and the stabilator at the desired lift coefficient. The lift and the pitching moment coefficients are computed using the Missile DATCOM v. 707. The flow field around the wing-canard- stabilator body model is resolved using the OVERFLOW-2 flow solver. Overset/ chimera grid topology is used to simplify the grid generation of various configurations representing different streamwise lift distributions. The numerical simulations are performed without viscosity unless it is required for numerical stability. All configurations are simulated at Mach 1.4, angle-of-attack of 1.50, lift coefficient of 0.05, and pitching moment coefficient of approximately 0. Four streamwise lift distribution configurations were tested.

Relations of different types of numerical magnitude representations to each other and to mathematics achievement.

PubMed

Fazio, Lisa K; Bailey, Drew H; Thompson, Clarissa A; Siegler, Robert S

2014-07-01

We examined relations between symbolic and non-symbolic numerical magnitude representations, between whole number and fraction representations, and between these representations and overall mathematics achievement in fifth graders. Fraction and whole number symbolic and non-symbolic numerical magnitude understandings were measured using both magnitude comparison and number line estimation tasks. After controlling for non-mathematical cognitive proficiency, both symbolic and non-symbolic numerical magnitude understandings were uniquely related to mathematics achievement, but the relation was much stronger for symbolic numbers. A meta-analysis of 19 published studies indicated that relations between non-symbolic numerical magnitude knowledge and mathematics achievement are present but tend to be weak, especially beyond 6 years of age. Copyright © 2014 Elsevier Inc. All rights reserved.

Overview of ten-year operation of the superconducting linear accelerator at the Spallation Neutron Source

NASA Astrophysics Data System (ADS)

Kim, S.-H.; Afanador, R.; Barnhart, D. L.; Crofford, M.; Degraff, B. D.; Doleans, M.; Galambos, J.; Gold, S. W.; Howell, M. P.; Mammosser, J.; McMahan, C. J.; Neustadt, T. S.; Peters, C.; Saunders, J. W.; Strong, W. H.; Vandygriff, D. J.; Vandygriff, D. M.

2017-04-01

The Spallation Neutron Source (SNS) has acquired extensive operational experience of a pulsed proton superconducting linear accelerator (SCL) as a user facility. Numerous lessons have been learned in its first 10 years operation to achieve a stable and reliable operation of the SCL. In this paper, an overview of the SNS SCL design, qualification of superconducting radio frequency (SRF) cavities and ancillary subsystems, an overview of the SNS cryogenic system, the SCL operation including SCL output energy history and downtime statistics, performance stability of the SRF cavities, efforts for SRF cavity performance recovery and improvement at the SNS, and maintenance activities for cryomodules are introduced.

Hybrid Chaos Synchronization of Four-Scroll Systems via Active Control

NASA Astrophysics Data System (ADS)

Karthikeyan, Rajagopal; Sundarapandian, Vaidyanathan

2014-03-01

This paper investigates the hybrid chaos synchronization of identical Wang four-scroll systems (Wang, 2009), identical Liu-Chen four-scroll systems (Liu and Chen, 2004) and non-identical Wang and Liu-Chen four-scroll systems. Active control method is the method adopted to achieve the hybrid chaos synchronization of the four-scroll chaotic systems addressed in this paper and our synchronization results are established using Lyapunov stability theory. Since the Lyapunov exponents are not required for these calculations, the active control method is effective and convenient to hybrid synchronize identical and different Wang and Liu-Chen four-scroll chaotic systems. Numerical simulations are also shown to illustrate and validate the hybrid synchronization results derived in this paper.

A multi-layer discrete-ordinate method for vector radiative transfer in a vertically-inhomogeneous, emitting and scattering atmosphere. I - Theory. II - Application

NASA Technical Reports Server (NTRS)

Weng, Fuzhong

1992-01-01

A theory is developed for discretizing the vector integro-differential radiative transfer equation including both solar and thermal radiation. A complete solution and boundary equations are obtained using the discrete-ordinate method. An efficient numerical procedure is presented for calculating the phase matrix and achieving computational stability. With natural light used as a beam source, the Stokes parameters from the model proposed here are compared with the analytical solutions of Chandrasekhar (1960) for a Rayleigh scattering atmosphere. The model is then applied to microwave frequencies with a thermal source, and the brightness temperatures are compared with those from Stamnes'(1988) radiative transfer model.

Overview of ten-year operation of the superconducting linear accelerator at the Spallation Neutron Source

DOE PAGES

Kim, Sang-Ho; Afanador, Ralph; Barnhart, Debra L.; ...

2017-02-04

The Spallation Neutron Source (SNS) has acquired extensive operational experience of a pulsed proton superconducting linear accelerator (SCL) as a user facility. Numerous lessons have been learned in its first 10 years operation to achieve a stable and reliable operation of the SCL. In this paper, an overview of the SNS SCL design, qualification of superconducting radio frequency (SRF) cavities and ancillary subsystems, an overview of the SNS cryogenic system, the SCL operation including SCL output energy history and downtime statistics, performance stability of the SRF cavities, efforts for SRF cavity performance recovery and improvement at the SNS, and maintenancemore » activities for cryomodules are introduced.« less

Ultra-high brightness wavelength-stabilized kW-class fiber coupled diode laser

NASA Astrophysics Data System (ADS)

Huang, Robin K.; Chann, Bien; Glenn, John D.

2011-03-01

TeraDiode has produced a fiber-coupled direct diode laser with a power level of 1,040 W from a 200 μm core diameter, 0.18 numerical aperture (NA) output fiber at a single center wavelength. This was achieved with a novel beam combining and shaping technique using COTS diode lasers. The fiber-coupled output corresponds to a Beam Parameter Product (BPP) of 18 mm-mrad and is the lowest BPP kW-class direct diode laser yet reported. The laser has been used to demonstrate laser cutting and welding of steel sheet metal up to 6.65 mm thick. Further advances of these ultra-bright lasers are also projected.

Overview of ten-year operation of the superconducting linear accelerator at the Spallation Neutron Source

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

Kim, Sang-Ho; Afanador, Ralph; Barnhart, Debra L.

The Spallation Neutron Source (SNS) has acquired extensive operational experience of a pulsed proton superconducting linear accelerator (SCL) as a user facility. Numerous lessons have been learned in its first 10 years operation to achieve a stable and reliable operation of the SCL. In this paper, an overview of the SNS SCL design, qualification of superconducting radio frequency (SRF) cavities and ancillary subsystems, an overview of the SNS cryogenic system, the SCL operation including SCL output energy history and downtime statistics, performance stability of the SRF cavities, efforts for SRF cavity performance recovery and improvement at the SNS, and maintenancemore » activities for cryomodules are introduced.« less

[Marketing in health service--a necessity or a temporary whim?].

PubMed

Dworzański, Wojciech; Dworzańska, Anna; Niezabitowska, Ewa; Czerwonka, Maciej; Madej, Barbara; Burdan, Franciszek

2012-01-01

Permanent changes occurring in economy and politics, open borders, a greater number of people with higher education as well as growing social awareness concerning taking care of health are followed by changes in healthcare sector. In response to an increased social demand, new medical institutions are being founded. They offer their patients a vast array of services ranging from therapeutic to preventive and highly specialized procedures improving their standard of living. In order to attract regular patients, and therefore receive regular income, it is necessary to stabilize such institutions on services market. It can be achieved by numerous marketing tools which, when applied appropriately, enable to create a positive image of a medical institution.

Local bounds preserving stabilization for continuous Galerkin discretization of hyperbolic systems

NASA Astrophysics Data System (ADS)

Mabuza, Sibusiso; Shadid, John N.; Kuzmin, Dmitri

2018-05-01

The objective of this paper is to present a local bounds preserving stabilized finite element scheme for hyperbolic systems on unstructured meshes based on continuous Galerkin (CG) discretization in space. A CG semi-discrete scheme with low order artificial dissipation that satisfies the local extremum diminishing (LED) condition for systems is used to discretize a system of conservation equations in space. The low order artificial diffusion is based on approximate Riemann solvers for hyperbolic conservation laws. In this case we consider both Rusanov and Roe artificial diffusion operators. In the Rusanov case, two designs are considered, a nodal based diffusion operator and a local projection stabilization operator. The result is a discretization that is LED and has first order convergence behavior. To achieve high resolution, limited antidiffusion is added back to the semi-discrete form where the limiter is constructed from a linearity preserving local projection stabilization operator. The procedure follows the algebraic flux correction procedure usually used in flux corrected transport algorithms. To further deal with phase errors (or terracing) common in FCT type methods, high order background dissipation is added to the antidiffusive correction. The resulting stabilized semi-discrete scheme can be discretized in time using a wide variety of time integrators. Numerical examples involving nonlinear scalar Burgers equation, and several shock hydrodynamics simulations for the Euler system are considered to demonstrate the performance of the method. For time discretization, Crank-Nicolson scheme and backward Euler scheme are utilized.

Projecting High Beta Steady-State Scenarios from DIII-D Advanced Tokamk Discharges

NASA Astrophysics Data System (ADS)

Park, J. M.

2013-10-01

Fusion power plant studies based on steady-state tokamak operation suggest that normalized beta in the range of 4-6 is needed for economic viability. DIII-D is exploring a range of candidate high beta scenarios guided by FASTRAN modeling in a repeated cycle of experiment and modeling validation. FASTRAN is a new iterative numerical procedure coupled to the Integrated Plasma Simulator (IPS) that integrates models of core transport, heating and current drive, equilibrium and stability self-consistently to find steady state (d / dt = 0) solutions, and reproduces most features of DIII-D high beta discharges with a stationary current profile. Separately, modeling components such as core transport (TGLF) and off-axis neutral beam current drive (NUBEAM) show reasonable agreement with experiment. Projecting forward to scenarios possible on DIII-D with future upgrades, two self-consistent noninductive scenarios at βN > 4 are found: high qmin and high internal inductance li. Both have bootstrap current fraction fBS > 0 . 5 and rely on the planned addition of a second off-axis neutral beamline and increased electron cyclotron heating. The high qmin > 2 scenario achieves stable operation at βN as high as 5 by a very broad current density profile to improve the ideal-wall stabilization of low-n instabilities along with confinement enhancement from low magnetic shear. The li near 1 scenario does not depend on ideal-wall stabilization. Improved confinement from strong magnetic shear makes up for the lower pedestal needed to maintain li high. The tradeoff between increasing li and reduced edge pedestal determines the achievable βN (near 4) and fBS (near 0.5). This modeling identifies the necessary upgrades to achieve target scenarios and clarifies the pros and cons of particular scenarios to better inform the development of steady-state fusion. Supported by the US Department of Energy under DE-AC05-00OR22725 & DE-FC02-04ER54698.

Numerical modeling of carrier gas flow in atomic layer deposition vacuum reactor: A comparative study of lattice Boltzmann models

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

Pan, Dongqing; Chien Jen, Tien; Li, Tao

2014-01-15

This paper characterizes the carrier gas flow in the atomic layer deposition (ALD) vacuum reactor by introducing Lattice Boltzmann Method (LBM) to the ALD simulation through a comparative study of two LBM models. Numerical models of gas flow are constructed and implemented in two-dimensional geometry based on lattice Bhatnagar–Gross–Krook (LBGK)-D2Q9 model and two-relaxation-time (TRT) model. Both incompressible and compressible scenarios are simulated and the two models are compared in the aspects of flow features, stability, and efficiency. Our simulation outcome reveals that, for our specific ALD vacuum reactor, TRT model generates better steady laminar flow features all over the domainmore » with better stability and reliability than LBGK-D2Q9 model especially when considering the compressible effects of the gas flow. The LBM-TRT is verified indirectly by comparing the numerical result with conventional continuum-based computational fluid dynamics solvers, and it shows very good agreement with these conventional methods. The velocity field of carrier gas flow through ALD vacuum reactor was characterized by LBM-TRT model finally. The flow in ALD is in a laminar steady state with velocity concentrated at the corners and around the wafer. The effects of flow fields on precursor distributions, surface absorptions, and surface reactions are discussed in detail. Steady and evenly distributed velocity field contribute to higher precursor concentration near the wafer and relatively lower particle velocities help to achieve better surface adsorption and deposition. The ALD reactor geometry needs to be considered carefully if a steady and laminar flow field around the wafer and better surface deposition are desired.« less

Numerical analysis of the angular motion of a neutrally buoyant spheroid in shear flow at small Reynolds numbers.

PubMed

Rosén, T; Einarsson, J; Nordmark, A; Aidun, C K; Lundell, F; Mehlig, B

2015-12-01

We numerically analyze the rotation of a neutrally buoyant spheroid in a shear flow at small shear Reynolds number. Using direct numerical stability analysis of the coupled nonlinear particle-flow problem, we compute the linear stability of the log-rolling orbit at small shear Reynolds number Re(a). As Re(a)→0 and as the box size of the system tends to infinity, we find good agreement between the numerical results and earlier analytical predictions valid to linear order in Re(a) for the case of an unbounded shear. The numerical stability analysis indicates that there are substantial finite-size corrections to the analytical results obtained for the unbounded system. We also compare the analytical results to results of lattice Boltzmann simulations to analyze the stability of the tumbling orbit at shear Reynolds numbers of order unity. Theory for an unbounded system at infinitesimal shear Reynolds number predicts a bifurcation of the tumbling orbit at aspect ratio λ(c)≈0.137 below which tumbling is stable (as well as log rolling). The simulation results show a bifurcation line in the λ-Re(a) plane that reaches λ≈0.1275 at the smallest shear Reynolds number (Re(a)=1) at which we could simulate with the lattice Boltzmann code, in qualitative agreement with the analytical results.

Generalization of von Neumann analysis for a model of two discrete half-spaces: The acoustic case

USGS Publications Warehouse

Haney, M.M.

2007-01-01

Evaluating the performance of finite-difference algorithms typically uses a technique known as von Neumann analysis. For a given algorithm, application of the technique yields both a dispersion relation valid for the discrete time-space grid and a mathematical condition for stability. In practice, a major shortcoming of conventional von Neumann analysis is that it can be applied only to an idealized numerical model - that of an infinite, homogeneous whole space. Experience has shown that numerical instabilities often arise in finite-difference simulations of wave propagation at interfaces with strong material contrasts. These interface instabilities occur even though the conventional von Neumann stability criterion may be satisfied at each point of the numerical model. To address this issue, I generalize von Neumann analysis for a model of two half-spaces. I perform the analysis for the case of acoustic wave propagation using a standard staggered-grid finite-difference numerical scheme. By deriving expressions for the discrete reflection and transmission coefficients, I study under what conditions the discrete reflection and transmission coefficients become unbounded. I find that instabilities encountered in numerical modeling near interfaces with strong material contrasts are linked to these cases and develop a modified stability criterion that takes into account the resulting instabilities. I test and verify the stability criterion by executing a finite-difference algorithm under conditions predicted to be stable and unstable. ?? 2007 Society of Exploration Geophysicists.

Spherical gyroscopic moment stabilizer for attitude control of microsatellites

NASA Astrophysics Data System (ADS)

Keshtkar, Sajjad; Moreno, Jaime A.; Kojima, Hirohisa; Uchiyama, Kenji; Nohmi, Masahiro; Takaya, Keisuke

2018-02-01

This paper presents a new and improved concept of recently proposed two-degrees of freedom spherical stabilizer for triaxial orientation of microsatellites. The analytical analysis of the advantages of the proposed mechanism over the existing inertial attitude control devices are introduced. The extended equations of motion of the stabilizing satellite including the spherical gyroscope, for control law design and numerical simulations, are studied in detail. A new control algorithm based on continuous high-order sliding mode algorithms, for managing the torque produced by the stabilizer and therefore the attitude control of the satellite in the presence of perturbations/uncertainties, is presented. Some numerical simulations are carried out to prove the performance of the proposed mechanism and control laws.

Attitude stability of a spinning spacecraft during appendage deployment/retraction

NASA Technical Reports Server (NTRS)

Fitz-Coy, Norman; Fullerton, Wayne

1994-01-01

The work presented is motivated by the need for a national satellite rescue policy, not the ad hoc policy now in place. In studying different approaches for a national policy, the issue of capture and stabilization of a tumbling spacecraft must be addressed. For a rescue mission involving a tumbling spacecraft, it may be advantageous to have a rescue vehicle which is compact and 'rigid' during the rendezvous/capture phase. After capture, passive stabilization techniques could be utilized as an efficient means of detumbling the resulting system (i.e., both the rescue vehicle and captures spacecraft). Since the rescue vehicle is initially compact and 'rigid,' significant passive stabilization through energy dissipation can only be achieved through the deployment of flexible appendages. Once stabilization is accomplished, retraction of the appendages before maneuvering the system to its final destination may also prove advantageous. It is therefore of paramount interest that we study the effect of appendage deployment/retraction on the attitude stability of a spacecraft. Particular interest should be paid to appendage retraction, since if this process is destabilizing, passive stabilization as proposed may not be useful. Over the past three decades, it has been an 'on-again-off-again affair' with the problem of spacecraft appendage deployment. In most instances, these studies have been numerical simulations of specific spacecraft configurations for which there were specific concerns. The primary focus of these studies was the behavior of the appendage during deployment; the effects of appendage retraction was considered only in one of these studies. What is missing in the literature is a thorough study of the effects of appendage deployment/retraction on the attitude stability of a spacecraft. This paper presents a rigorous analysis of the stability of a spinning spacecraft during the deployment or the retraction of an appendage. The analysis is simplified such that meaningful insights into the problem can be inferred; it is not overly simplified such that critical dynamical behavior is neglected. The system is analyzed assuming that the spacecraft hub is rigid. The appendage deployment mechanism is modeled as a point mass on a massless rod whose length undergoes prescribed changes. Simplified flexibility effects of the appendage are included. The system is examined for stability by linearizing the equations in terms of small deviations from steady, noninterfering coning motion. Routh's procedure for analyzing small deviations from steady motion in dynamical systems is utilized in the analysis. The system of equations are nondimensionalized to facilitate parametric studies. The results are presented in terms of a reduced number of nondimensional parameters so that some general conclusions may be drawn. Verification of the linear analysis is presented through numerical simulations of the complete nonlinear, nonautonomous, coupled equations.

Computation of the stability derivatives via CFD and the sensitivity equations

NASA Astrophysics Data System (ADS)

Lei, Guo-Dong; Ren, Yu-Xin

2011-04-01

The method to calculate the aerodynamic stability derivates of aircrafts by using the sensitivity equations is extended to flows with shock waves in this paper. Using the newly developed second-order cell-centered finite volume scheme on the unstructured-grid, the unsteady Euler equations and sensitivity equations are solved simultaneously in a non-inertial frame of reference, so that the aerodynamic stability derivatives can be calculated for aircrafts with complex geometries. Based on the numerical results, behavior of the aerodynamic sensitivity parameters near the shock wave is discussed. Furthermore, the stability derivatives are analyzed for supersonic and hypersonic flows. The numerical results of the stability derivatives are found in good agreement with theoretical results for supersonic flows, and variations of the aerodynamic force and moment predicted by the stability derivatives are very close to those obtained by CFD simulation for both supersonic and hypersonic flows.

Analysis of Time Filters in Multistep Methods

NASA Astrophysics Data System (ADS)

Hurl, Nicholas

Geophysical ow simulations have evolved sophisticated implicit-explicit time stepping methods (based on fast-slow wave splittings) followed by time filters to control any unstable models that result. Time filters are modular and parallel. Their effect on stability of the overall process has been tested in numerous simulations, but never analyzed. Stability is proven herein for the Crank-Nicolson Leapfrog (CNLF) method with the Robert-Asselin (RA) time filter and for the Crank-Nicolson Leapfrog method with the Robert-Asselin-Williams (RAW) time filter for systems by energy methods. We derive an equivalent multistep method for CNLF+RA and CNLF+RAW and stability regions are obtained. The time step restriction for energy stability of CNLF+RA is smaller than CNLF and CNLF+RAW time step restriction is even smaller. Numerical tests find that RA and RAW add numerical dissipation. This thesis also shows that all modes of the Crank-Nicolson Leap Frog (CNLF) method are asymptotically stable under the standard timestep condition.

Experimental and numerical analysis of the influence of tyres' properties on the straight running stability of a sport-touring motorcycle

NASA Astrophysics Data System (ADS)

Cossalter, Vittore; Doria, Alberto; Formentini, Matteo; Peretto, Martino

2012-03-01

The behaviour of a motorcycle on the road is largely governed by tyre properties. This paper presents experimental and numerical analyses dealing with the influence of tyre properties on the stability of weave and wobble in straight running. The final goal is to find optimal sets of tyre properties that improve the stability of a motorcycle. The investigation is based on road tests carried out on a sport-touring motorcycle equipped with sensors. Three sets of tyres are tested at different speeds in the presence of weave and wobble. The analysis of telemetry data highlights significant differences in the trends of frequency and damping of weave and wobble against speed. The experimental analysis is integrated by a parametric numerical analysis. Tyre properties are varied according to the design of experiments method, in order to highlight the single effects on stability of lateral and cornering coefficient of front and rear tyres.

Theoretical stability in coefficient inverse problems for general hyperbolic equations with numerical reconstruction

NASA Astrophysics Data System (ADS)

Yu, Jie; Liu, Yikan; Yamamoto, Masahiro

2018-04-01

In this article, we investigate the determination of the spatial component in the time-dependent second order coefficient of a hyperbolic equation from both theoretical and numerical aspects. By the Carleman estimates for general hyperbolic operators and an auxiliary Carleman estimate, we establish local Hölder stability with either partial boundary or interior measurements under certain geometrical conditions. For numerical reconstruction, we minimize a Tikhonov functional which penalizes the gradient of the unknown function. Based on the resulting variational equation, we design an iteration method which is updated by solving a Poisson equation at each step. One-dimensional prototype examples illustrate the numerical performance of the proposed iteration.

A stability analysis on forced convection boundary layer stagnation-point slip flow in Darcy-Forchheimer porous medium towards a shrinking sheet

NASA Astrophysics Data System (ADS)

Bakar, Shahirah Abu; Arifin, Norihan Md; Ali, Fadzilah Md; Bachok, Norfifah; Nazar, Roslinda

2017-08-01

The stagnation-point flow over a shrinking sheet in Darcy-Forchheimer porous medium is numerically studied. The governing partial differential equations are transformed into ordinary differential equations using a similarity transformation, and then solved numerically by using shooting technique method with Maple implementation. Dual solutions are observed in a certain range of the shrinking parameter. Regarding on numerical solutions, we prepared stability analysis to identify which solution is stable between non-unique solutions by bvp4c solver in Matlab. Further we obtain numerical results or each solution, which enable us to discuss the features of the respective solutions.

Stabilization of the Rayleigh-Taylor instability in quantum magnetized plasmas

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

Wang, L. F.; Ye, W. H.; He, X. T.

2012-07-15

In this research, stabilization of the Rayleigh-Taylor instability (RTI) due to density gradients, magnetic fields, and quantum effects, in an ideal incompressible plasma, is studied analytically and numerically. A second-order ordinary differential equation (ODE) for the RTI including quantum corrections, with a continuous density profile, in a uniform external magnetic field, is obtained. Analytic expressions of the linear growth rate of the RTI, considering modifications of density gradients, magnetic fields, and quantum effects, are presented. Numerical approaches are performed to solve the second-order ODE. The analytical model proposed here agrees with the numerical calculation. It is found that the densitymore » gradients, the magnetic fields, and the quantum effects, respectively, have a stabilizing effect on the RTI (reduce the linear growth of the RTI). The RTI can be completely quenched by the magnetic field stabilization and/or the quantum effect stabilization in proper circumstances leading to a cutoff wavelength. The quantum effect stabilization plays a central role in systems with large Atwood number and small normalized density gradient scale length. The presence of external transverse magnetic fields beside the quantum effects will bring about more stability on the RTI. The stabilization of the linear growth of the RTI, for parameters closely related to inertial confinement fusion and white dwarfs, is discussed. Results could potentially be valuable for the RTI treatment to analyze the mixing in supernovas and other RTI-driven objects.« less

Boundary control of anti-symmetric vibration of satellite with flexible appendages in planar motion with exponential stability

NASA Astrophysics Data System (ADS)

Rad, Hossein Kaviani; Salarieh, Hassan; Alasty, Aria; Vatankhah, Ramin

2018-06-01

In this research, we have investigated the planar maneuver of a flexible satellite with appendages anti-symmetric vibration. The hybrid governing equations are comprised of coupled partial and ordinary differential equations which are derived by employing Hamilton's principle. In this paper, control goals are the tracking desired pitch angle along with the flexible appendages vibration suppression simultaneously by using only one control torque which is applied to the central hub. The boundary control is proposed to fulfill these control aims; furthermore, this boundary control ensures that spillover instability phenomenon is eliminated, and in-domain sensors and actuators implement are excluded. Indeed, the proposed boundary control is able to stabilize an infinite number of vibration modes, which is one of the important benefits of the proposed control when it is considered that different factors including external disturbances and even the satellite maneuver can excite the various vibration modes of the flexible appendages and consequently the excitement of the high order vibration modes will be possible. Lyapunov's direct method is used to prove the exponential stability; moreover, this Proof is achieved in absence of any damping effect in modeling the vibrations of flexible appendages. In addition, the procedure for finding the boundary control coefficients which ensures the exponential stability is provided. Eventually, numerical simulations are presented to illustrate the effectiveness of the proposed boundary control.

Flow past a rotating cylinder

NASA Astrophysics Data System (ADS)

Mittal, Sanjay; Kumar, Bhaskar

2003-02-01

Flow past a spinning circular cylinder placed in a uniform stream is investigated via two-dimensional computations. A stabilized finite element method is utilized to solve the incompressible Navier Stokes equations in the primitive variables formulation. The Reynolds number based on the cylinder diameter and free-stream speed of the flow is 200. The non-dimensional rotation rate, [alpha] (ratio of the surface speed and freestream speed), is varied between 0 and 5. The time integration of the flow equations is carried out for very large dimensionless time. Vortex shedding is observed for [alpha] < 1.91. For higher rotation rates the flow achieves a steady state except for 4.34 < [alpha] < 4:70 where the flow is unstable again. In the second region of instability, only one-sided vortex shedding takes place. To ascertain the instability of flow as a function of [alpha] a stabilized finite element formulation is proposed to carry out a global, non-parallel stability analysis of the two-dimensional steady-state flow for small disturbances. The formulation and its implementation are validated by predicting the Hopf bifurcation for flow past a non-rotating cylinder. The results from the stability analysis for the rotating cylinder are in very good agreement with those from direct numerical simulations. For large rotation rates, very large lift coefficients can be obtained via the Magnus effect. However, the power requirement for rotating the cylinder increases rapidly with rotation rate.

Dynamical approach study of spurious steady-state numerical solutions of nonlinear differential equations. Part 1: The ODE connection and its implications for algorithm development in computational fluid dynamics

NASA Technical Reports Server (NTRS)

Yee, H. C.; Sweby, P. K.; Griffiths, D. F.

1990-01-01

Spurious stable as well as unstable steady state numerical solutions, spurious asymptotic numerical solutions of higher period, and even stable chaotic behavior can occur when finite difference methods are used to solve nonlinear differential equations (DE) numerically. The occurrence of spurious asymptotes is independent of whether the DE possesses a unique steady state or has additional periodic solutions and/or exhibits chaotic phenomena. The form of the nonlinear DEs and the type of numerical schemes are the determining factor. In addition, the occurrence of spurious steady states is not restricted to the time steps that are beyond the linearized stability limit of the scheme. In many instances, it can occur below the linearized stability limit. Therefore, it is essential for practitioners in computational sciences to be knowledgeable about the dynamical behavior of finite difference methods for nonlinear scalar DEs before the actual application of these methods to practical computations. It is also important to change the traditional way of thinking and practices when dealing with genuinely nonlinear problems. In the past, spurious asymptotes were observed in numerical computations but tended to be ignored because they all were assumed to lie beyond the linearized stability limits of the time step parameter delta t. As can be seen from the study, bifurcations to and from spurious asymptotic solutions and transitions to computational instability not only are highly scheme dependent and problem dependent, but also initial data and boundary condition dependent, and not limited to time steps that are beyond the linearized stability limit.

The design and pre-clinical evaluation of knee replacements for osteoarthritis.

PubMed

Walker, Peter S

2015-03-18

One of the concepts that Rik Huiskes promoted was that implants such as knee and hip replacements could be analyzed and optimized using numerical models such as finite element analysis, or by experimental testing, an area he called pre-clinical testing. The design itself could be formulated or improved by defining a specific goal or asking a key question. These propositions are examined in the light of almost five decades of experience with knee implants. Achieving the required laxity and stability, was achieved by attempting to reproduce anatomical values by suitable radii of curvature and selective ligament retention. Obtaining durable fixation was based on testing many configurations to obtain the most uniform stress distribution at the implant-bone interface. Achieving the best overall kinematics has yet to be fully solved due to the variations in activities and patients. These and many other factors have usually been addressed individually rather than as a composite, although as time has gone on, successful features have gradually been assimilated into most designs. But even a systematic approach has been flawed because some unrecognized response was not accounted for in the pre-clinical model, a limitation of models in general. In terms of the design process, so far no method has emerged for systematically reaching an optimal solution from all aspects, although this is possible in principle. Overall however, predictive numerical or physical models should be an essential element in the design of new or improved knee replacements, a part of the design process itself. Copyright © 2015. Published by Elsevier Ltd.

Continuous uniformly finite time exact disturbance observer based control for fixed-time stabilization of nonlinear systems with mismatched disturbances

PubMed Central

Liu, Chongxin; Liu, Hang

2017-01-01

This paper presents a continuous composite control scheme to achieve fixed-time stabilization for nonlinear systems with mismatched disturbances. The composite controller is constructed in two steps: First, uniformly finite time exact disturbance observers are proposed to estimate and compensate the disturbances. Then, based on adding a power integrator technique and fixed-time stability theory, continuous fixed-time stable state feedback controller and Lyapunov functions are constructed to achieve global fixed-time system stabilization. The proposed control method extends the existing fixed-time stable control results to high order nonlinear systems with mismatched disturbances and achieves global fixed-time system stabilization. Besides, the proposed control scheme improves the disturbance rejection performance and achieves performance recovery of nominal system. Simulation results are provided to show the effectiveness, the superiority and the applicability of the proposed control scheme. PMID:28406966

Evaluating the Effect of Rainfall Infiltration on the Slope Stability of T16 tower of Taipei Mao-kong Gondola by Numerical Methods

NASA Astrophysics Data System (ADS)

RUNG, J.

2013-12-01

In this study, a series of rainfall-stability analyses were performed to simulate the failure mechanism and the function of remediation works of the down slope of T-16 tower pier, Mao-Kong gondola (or T-16 Slope) at the hillside of Taipei City using two-dimensional finite element method. The failure mechanism of T-16 Slope was simulated using the rainfall hyetograph of Jang-Mi typhoon in 2008 based on the field investigation data, monitoring data, soil/rock mechanical testing data and detail design plots of remediation works. Eventually, the numerical procedures and various input parameters in the analysis were verified by comparing the numerical results with the field observations. In addition, 48 hrs design rainfalls corresponding to 5, 10, 25 and 50 years return periods were prepared using the 20 years rainfall data of Mu-Zha rainfall observation station, Central Weather Bureau for the rainfall-stability analyses of T-16 Slope to inspect the effect of the compound stabilization works on the overall stability of the slope. At T-16 Slope, without considering the longitudinal and transverse drainages on the ground surface, there totally 4 types of stabilization works were installed to stabilize the slope. From the slope top to the slope toe, the stabilization works of T-16 Slope consists of RC-retaining wall with micro-pile foundation at the up-segment, earth anchor at the up-middle-segment, soil nailing at the middle-segment and retaining pile at the down-segment of the slope. The effect of each individual stabilization work on the slope stability under rainfall condition was examined and evaluated by raising field groundwater level.

Two-Relaxation-Time Lattice Boltzmann Method for Advective-Diffusive-Reactive Transport

NASA Astrophysics Data System (ADS)

Yan, Z.; Hilpert, M.

2016-12-01

The lattice Boltzmann method (LBM) has been applied to study a wide range of reactive transport in porous and fractured media. The single-relaxation-time (SRT) LBM, employing single relaxation time, is the most popular LBM due to its simplicity of understanding and implementation. Nevertheless, the SRT LBM may suffer from numerical instability for small value of the relaxation time. By contrast, the multiple-relaxation-time (MRT) LBM, employing multiple relaxation times, can improve the numerical stability through tuning the multiple relaxation times, but the complexity of implementing this method restricts its applications. The two-relaxation-time (TRT) LBM, which employs two relaxation times, combines the advantages of SRT and MRT LBMs. The TRT LBM can produce simulations with better accuracy and stability than the SRT one, and is easier to implement than the MRT one. This work evaluated the numerical accuracy and stability of the TRT method by comparing the simulation results with analytical solutions of Gaussian hill transport and Taylor dispersion under different advective velocities. The accuracy generally increased with the tunable relaxation time τ, and the stability first increased and then decreased as τ increased, showing an optimal TRT method emerging the best numerical stability. The free selection of τ enabled the TRT LBM to simulate the Gaussian hill transport and Taylor dispersion under relatively high advective velocity, under which the SRT LBM suffered from numerical instability. Finally, the TRT method was applied to study the contaminant degradation by chemotactic microorganisms in porous media, which acted as a reprehensive of reactive transport in this study, and well predicted the evolution of microorganisms and degradation of contaminants for different transport scenarios. To sum up, the TRT LBM produced simulation results with good accuracy and stability for various advective-diffusive-reactive transport through tuning the relaxation time τ, illustrating its potential to study various biogeochemical behaviors in the subsurface environment.

An extended continuum model accounting for the driver's timid and aggressive attributions

NASA Astrophysics Data System (ADS)

Cheng, Rongjun; Ge, Hongxia; Wang, Jufeng

2017-04-01

Considering the driver's timid and aggressive behaviors simultaneously, a new continuum model is put forwarded in this paper. By applying the linear stability theory, we presented the analysis of new model's linear stability. Through nonlinear analysis, the KdV-Burgers equation is derived to describe density wave near the neutral stability line. Numerical results verify that aggressive driving is better than timid act because the aggressive driver will adjust his speed timely according to the leading car's speed. The key improvement of this new model is that the timid driving deteriorates traffic stability while the aggressive driving will enhance traffic stability. The relationship of energy consumption between the aggressive and timid driving is also studied. Numerical results show that aggressive driver behavior can not only suppress the traffic congestion but also reduce the energy consumption.

Interpreting Popov criteria in Lure´ systems with complex scaling stability analysis

NASA Astrophysics Data System (ADS)

Zhou, J.

2018-06-01

The paper presents a novel frequency-domain interpretation of Popov criteria for absolute stability in Lure´ systems by means of what we call complex scaling stability analysis. The complex scaling technique is developed for exponential/asymptotic stability in LTI feedback systems, which dispenses open-loop poles distribution, contour/locus orientation and prior frequency sweeping. Exploiting the technique for alternatively revealing positive realness of transfer functions, re-interpreting Popov criteria is explicated. More specifically, the suggested frequency-domain stability conditions are conformable both in scalar and multivariable cases, and can be implemented either graphically with locus plotting or numerically without; in particular, the latter is suitable as a design tool with auxiliary parameter freedom. The interpretation also reveals further frequency-domain facts about Lure´ systems. Numerical examples are included to illustrate the main results.

From h to p efficiently: optimal implementation strategies for explicit time-dependent problems using the spectral/hp element method

PubMed Central

Bolis, A; Cantwell, C D; Kirby, R M; Sherwin, S J

2014-01-01

We investigate the relative performance of a second-order Adams–Bashforth scheme and second-order and fourth-order Runge–Kutta schemes when time stepping a 2D linear advection problem discretised using a spectral/hp element technique for a range of different mesh sizes and polynomial orders. Numerical experiments explore the effects of short (two wavelengths) and long (32 wavelengths) time integration for sets of uniform and non-uniform meshes. The choice of time-integration scheme and discretisation together fixes a CFL limit that imposes a restriction on the maximum time step, which can be taken to ensure numerical stability. The number of steps, together with the order of the scheme, affects not only the runtime but also the accuracy of the solution. Through numerical experiments, we systematically highlight the relative effects of spatial resolution and choice of time integration on performance and provide general guidelines on how best to achieve the minimal execution time in order to obtain a prescribed solution accuracy. The significant role played by higher polynomial orders in reducing CPU time while preserving accuracy becomes more evident, especially for uniform meshes, compared with what has been typically considered when studying this type of problem.© 2014. The Authors. International Journal for Numerical Methods in Fluids published by John Wiley & Sons, Ltd. PMID:25892840

Numerical simulations and linear stability analysis of a boundary layer developed on wavy surfaces

NASA Astrophysics Data System (ADS)

Siconolfi, Lorenzo; Camarri, Simone; Fransson, Jens H. M.

2015-11-01

The development of passive methods leading to a laminar to turbulent transition delay in a boundary layer (BL) is a topic of great interest both for applications and academic research. In literature it has been shown that a proper and stable spanwise velocity modulation can reduce the growth rate of Tollmien-Schlichting (TS) waves and delay transition. In this study, we investigate numerically the possibility of obtaining a stabilizing effect of the TS waves through the use of a spanwise sinusoidal modulation of a flat plate. This type of control has been already successfully investigated experimentally. An extensive set of direct numerical simulations is carried out to study the evolution of a BL flow developed on wavy surfaces with different geometric characteristics, and the results will be presented here. Moreover, since this configuration is characterized by a slowly-varying flow field in streamwise direction, a local stability analysis is applied to define the neutral stability curves for the BL flow controlled by this type of wall modifications. These results give the possibility of investigating this control strategy and understanding the effect of the free parameters on the stabilization mechanism.

Guidelines for Computing Longitudinal Dynamic Stability Characteristics of a Subsonic Transport

NASA Technical Reports Server (NTRS)

Thompson, Joseph R.; Frank, Neal T.; Murphy, Patrick C.

2010-01-01

A systematic study is presented to guide the selection of a numerical solution strategy for URANS computation of a subsonic transport configuration undergoing simulated forced oscillation about its pitch axis. Forced oscillation is central to the prevalent wind tunnel methodology for quantifying aircraft dynamic stability derivatives from force and moment coefficients, which is the ultimate goal for the computational simulations. Extensive computations are performed that lead in key insights of the critical numerical parameters affecting solution convergence. A preliminary linear harmonic analysis is included to demonstrate the potential of extracting dynamic stability derivatives from computational solutions.

The stability issues in problems of mathematical modeling

NASA Astrophysics Data System (ADS)

Mokin, A. Yu.; Savenkova, N. P.; Udovichenko, N. S.

2018-03-01

In the paper it is briefly considered various aspects of stability concepts, which are used in physics, mathematics and numerical methods of solution. The interrelation between these concepts is described, the questions of preliminary stability research before the numerical solution of the problem and the correctness of the mathematical statement of the physical problem are discussed. Examples of concrete mathematical statements of individual physical problems are given: a nonlocal problem for the heat equation, the Korteweg-de Fries equation with boundary conditions at infinity, the sine-Gordon equation, the problem of propagation of femtosecond light pulses in an area with a cubic nonlinearity.

Linear stability analysis of detonations via numerical computation and dynamic mode decomposition

NASA Astrophysics Data System (ADS)

Kabanov, Dmitry I.; Kasimov, Aslan R.

2018-03-01

We introduce a new method to investigate linear stability of gaseous detonations that is based on an accurate shock-fitting numerical integration of the linearized reactive Euler equations with a subsequent analysis of the computed solution via the dynamic mode decomposition. The method is applied to the detonation models based on both the standard one-step Arrhenius kinetics and two-step exothermic-endothermic reaction kinetics. Stability spectra for all cases are computed and analyzed. The new approach is shown to be a viable alternative to the traditional normal-mode analysis used in detonation theory.

Blade-Cast Nonfullerene Organic Solar Cells in Air with Excellent Morphology, Efficiency, and Stability.

PubMed

Zhang, Lin; Lin, Baojun; Hu, Bo; Xu, Xianbin; Ma, Wei

2018-04-17

Blade-coating serving as a prototype tool for slot-die coating can be very compatible with large-area roll-to-roll coating. Using blade-coating in an ambient environment, an average power conversion efficiency (PCE) of 10.03% is achieved in nonfullerene organic solar cells, which is higher than that of the optimal spin-coated device with a PCE of 9.41%. It is demonstrated that blade-coating can induce a higher degree of molecular packing for both conjugated polymer donors and small-molecular acceptors as it helps to produce a seeding film containing numerous crystal grains, subsequently providing nucleation sites for the residual solution when the motion of the blade exposes a liquid front. Due to this effect, blade-coating can partially replace the role of the additive 1,8-diiodooctane (DIO) and thus achieves the optimized morphology with fewer additives. Moreover, it is found that the blade-coated film with 0.25% DIO possesses not only a smaller domain size but also higher domain purity, suggesting more D/A (donor/acceptor) interfaces and a purer phase domain as compared to the spin-coated film with 1% DIO. Encouragingly, the blade-coated device with less DIO (0.25%) exhibits much better stability than the spin-coated device with 1% DIO, showing excellent prospects. © 2018 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim.

Enhanced methodology of focus control and monitoring on scanner tool

NASA Astrophysics Data System (ADS)

Chen, Yen-Jen; Kim, Young Ki; Hao, Xueli; Gomez, Juan-Manuel; Tian, Ye; Kamalizadeh, Ferhad; Hanson, Justin K.

2017-03-01

As the demand of the technology node shrinks from 14nm to 7nm, the reliability of tool monitoring techniques in advanced semiconductor fabs to achieve high yield and quality becomes more critical. Tool health monitoring methods involve periodic sampling of moderately processed test wafers to detect for particles, defects, and tool stability in order to ensure proper tool health. For lithography TWINSCAN scanner tools, the requirements for overlay stability and focus control are very strict. Current scanner tool health monitoring methods include running BaseLiner to ensure proper tool stability on a periodic basis. The focus measurement on YIELDSTAR by real-time or library-based reconstruction of critical dimensions (CD) and side wall angle (SWA) has been demonstrated as an accurate metrology input to the control loop. The high accuracy and repeatability of the YIELDSTAR focus measurement provides a common reference of scanner setup and user process. In order to further improve the metrology and matching performance, Diffraction Based Focus (DBF) metrology enabling accurate, fast, and non-destructive focus acquisition, has been successfully utilized for focus monitoring/control of TWINSCAN NXT immersion scanners. The optimal DBF target was determined to have minimized dose crosstalk, dynamic precision, set-get residual, and lens aberration sensitivity. By exploiting this new measurement target design, 80% improvement in tool-to-tool matching, >16% improvement in run-to-run mean focus stability, and >32% improvement in focus uniformity have been demonstrated compared to the previous BaseLiner methodology. Matching <2.4 nm across multiple NXT immersion scanners has been achieved with the new methodology of set baseline reference. This baseline technique, with either conventional BaseLiner low numerical aperture (NA=1.20) mode or advanced illumination high NA mode (NA=1.35), has also been evaluated to have consistent performance. This enhanced methodology of focus control and monitoring on multiple illumination conditions, opens an avenue to significantly reduce Focus-Exposure Matrix (FEM) wafer exposure for new product/layer best focus (BF) setup.

Stabilization of ring dark solitons in Bose-Einstein condensates

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

Wang, Wenlong; Kevrekidis, P. G.; Carretero-González, R.

Earlier work has shown that ring dark solitons in two-dimensional Bose-Einstein condensates are generically unstable. In this work, we propose a way of stabilizing the ring dark soliton via a radial Gaussian external potential. We investigate the existence and stability of the ring dark soliton upon variations of the chemical potential and also of the strength of the radial potential. Numerical results show that the ring dark soliton can be stabilized in a suitable interval of external potential strengths and chemical potentials. Furthermore, we also explore different proposed particle pictures considering the ring as a moving particle and find, wheremore » appropriate, results in very good qualitative and also reasonable quantitative agreement with the numerical findings.« less

Stabilization of ring dark solitons in Bose-Einstein condensates

DOE PAGES

Wang, Wenlong; Kevrekidis, P. G.; Carretero-González, R.; ...

2015-09-14

Earlier work has shown that ring dark solitons in two-dimensional Bose-Einstein condensates are generically unstable. In this work, we propose a way of stabilizing the ring dark soliton via a radial Gaussian external potential. We investigate the existence and stability of the ring dark soliton upon variations of the chemical potential and also of the strength of the radial potential. Numerical results show that the ring dark soliton can be stabilized in a suitable interval of external potential strengths and chemical potentials. Furthermore, we also explore different proposed particle pictures considering the ring as a moving particle and find, wheremore » appropriate, results in very good qualitative and also reasonable quantitative agreement with the numerical findings.« less

Linear stability and nonlinear analyses of traffic waves for the general nonlinear car-following model with multi-time delays

NASA Astrophysics Data System (ADS)

Sun, Dihua; Chen, Dong; Zhao, Min; Liu, Weining; Zheng, Linjiang

2018-07-01

In this paper, the general nonlinear car-following model with multi-time delays is investigated in order to describe the reactions of vehicle to driving behavior. Platoon stability and string stability criteria are obtained for the general nonlinear car-following model. Burgers equation and Korteweg de Vries (KdV) equation and their solitary wave solutions are derived adopting the reductive perturbation method. We investigate the properties of typical optimal velocity model using both analytic and numerical methods, which estimates the impact of delays about the evolution of traffic congestion. The numerical results show that time delays in sensing relative movement is more sensitive to the stability of traffic flow than time delays in sensing host motion.

Equilibrium points, stability and numerical solutions of fractional-order predator-prey and rabies models

NASA Astrophysics Data System (ADS)

Ahmed, E.; El-Sayed, A. M. A.; El-Saka, H. A. A.

2007-01-01

In this paper we are concerned with the fractional-order predator-prey model and the fractional-order rabies model. Existence and uniqueness of solutions are proved. The stability of equilibrium points are studied. Numerical solutions of these models are given. An example is given where the equilibrium point is a centre for the integer order system but locally asymptotically stable for its fractional-order counterpart.

Oscillations and stability of numerical solutions of the heat conduction equation

NASA Technical Reports Server (NTRS)

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

1976-01-01

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

Origin of vertical orientation in two-dimensional metal halide perovskites and its effect on photovoltaic performance.

PubMed

Chen, Alexander Z; Shiu, Michelle; Ma, Jennifer H; Alpert, Matthew R; Zhang, Depei; Foley, Benjamin J; Smilgies, Detlef-M; Lee, Seung-Hun; Choi, Joshua J

2018-04-06

Thin films based on two-dimensional metal halide perovskites have achieved exceptional performance and stability in numerous optoelectronic device applications. Simple solution processing of the 2D perovskite provides opportunities for manufacturing devices at drastically lower cost compared to current commercial technologies. A key to high device performance is to align the 2D perovskite layers, during the solution processing, vertical to the electrodes to achieve efficient charge transport. However, it is yet to be understood how the counter-intuitive vertical orientations of 2D perovskite layers on substrates can be obtained. Here we report a formation mechanism of such vertically orientated 2D perovskite in which the nucleation and growth arise from the liquid-air interface. As a consequence, choice of substrates can be liberal from polymers to metal oxides depending on targeted application. We also demonstrate control over the degree of preferential orientation of the 2D perovskite layers and its drastic impact on device performance.

The miniature parachute of the dandelion fruit

NASA Astrophysics Data System (ADS)

Cummins, Cathal; Viola, Ignazio Maria; Seale, Madeleine; Mastropaolo, Enrico; Nakayama, Naomi

2017-11-01

At the low Reynolds number at which small plant fruit (the seed-bearing structure in flowering plants) fly, there are a variety of modes of flight available: from parachuting to gliding and autorotation. Here we will explore the aerodynamics of small plumed fruit (dandelions) that utilise the parachuting mode of flight. If a parachute-type fruit is picked up by the breeze, it can be carried over formidable distances. Incredibly, these parachutes are mostly empty space, yet they are effectively impervious to the airflow as they descend. In addition, the fruit can become more or less streamlined depending on the environmental conditions. In this talk, we will present results from our numerical and physical modelling that clarify how these tiny parachutes achieve such impermeability despite their high porosity. We reveal that the dandelion's parachute tunes its permeability to achieve the aerodynamic stability as it flies, which helps confer the fruit's incredible flight capacity. This work was supported by the Leverhulme Trust [RPG-2015-255].

USP7/HAUSP: A SUMO deubiquitinase at the heart of DNA replication.

PubMed

Smits, Veronique A J; Freire, Raimundo

2016-09-01

DNA replication is both highly conserved and controlled. Problematic DNA replication can lead to genomic instability and therefore carcinogenesis. Numerous mechanisms work together to achieve this tight control and increasing evidence suggests that post-translational modifications (phosphorylation, ubiquitination, SUMOylation) of DNA replication proteins play a pivotal role in this process. Here we discuss such modifications in the light of a recent article that describes a novel role for the deubiquitinase (DUB) USP7/HAUSP in the control of DNA replication. USP7 achieves this function by an unusual and novel mechanism, namely deubiquitination of SUMOylated proteins at the replication fork, making USP7 also a SUMO DUB (SDUB). This work extends previous observations of increased levels of SUMO and low levels of ubiquitin at the on-going replication fork. Here, we discuss this novel study, its contribution to the DNA replication and genomic stability field and what questions arise from this work. © 2016 WILEY Periodicals, Inc.

Implant Insertion Torque: Its Role in Achieving Primary Stability of Restorable Dental Implants.

PubMed

Greenstein, Gary; Cavallaro, John

2017-02-01

A literature review was conducted to determine the role of insertion torque in attaining primary stability of dental implants. The review is comprised of articles that discussed the amount of torque needed to achieve primary implant stability in healed ridges and fresh extraction sockets prior to immediate implant loading. Studies were appraised that addressed the effects of minimum and maximum forces that can be used to successfully place implants. The minimum torque that can be employed to attain primary stability is undefined. Forces ≥30 Ncm are routinely used to place implants into healed ridges and fresh extraction sockets prior to immediate loading of implants. Increased insertion torque (≥50 Ncm) reduces micromotion and does not appear to damage bone. In general, the healing process after implant insertion provides a degree of biologic stability that is similar whether implants are placed with high or low initial insertion torque. Primary stability is desirable when placing implants, but the absence of micromotion is what facilitates predictable implant osseointegration. Increased insertion torque helps achieve primary stability by reducing implant micromotion. Furthermore, tactile information provided by the first surgical twist drill can aid in selecting the initial insertion torque to achieve predictable stability of inserted dental implants.

Dynamical approach study of spurious steady-state numerical solutions of nonlinear differential equations. I - The dynamics of time discretization and its implications for algorithm development in computational fluid dynamics

NASA Technical Reports Server (NTRS)

Yee, H. C.; Sweby, P. K.; Griffiths, D. F.

1991-01-01

Spurious stable as well as unstable steady state numerical solutions, spurious asymptotic numerical solutions of higher period, and even stable chaotic behavior can occur when finite difference methods are used to solve nonlinear differential equations (DE) numerically. The occurrence of spurious asymptotes is independent of whether the DE possesses a unique steady state or has additional periodic solutions and/or exhibits chaotic phenomena. The form of the nonlinear DEs and the type of numerical schemes are the determining factor. In addition, the occurrence of spurious steady states is not restricted to the time steps that are beyond the linearized stability limit of the scheme. In many instances, it can occur below the linearized stability limit. Therefore, it is essential for practitioners in computational sciences to be knowledgeable about the dynamical behavior of finite difference methods for nonlinear scalar DEs before the actual application of these methods to practical computations. It is also important to change the traditional way of thinking and practices when dealing with genuinely nonlinear problems. In the past, spurious asymptotes were observed in numerical computations but tended to be ignored because they all were assumed to lie beyond the linearized stability limits of the time step parameter delta t. As can be seen from the study, bifurcations to and from spurious asymptotic solutions and transitions to computational instability not only are highly scheme dependent and problem dependent, but also initial data and boundary condition dependent, and not limited to time steps that are beyond the linearized stability limit.

Frequency-tunable Pre-stabilized Lasers for LISA via Sideband-locking

NASA Technical Reports Server (NTRS)

Livas, Jeffrey C.; Thorpe, James I.; Numata, Kenji; Mitryk, Shawn; Mueller, Guido; Wand, Vinzenz

2008-01-01

Laser frequency noise mitigation is one of the most challenging aspects of the LISA interferometric measurement system. The unstabilized frequency fluctuations must be suppressed by roughly twelve orders of magnitude in order to achieve stability sufficient for gravitational wave detection. This enormous suppression will be achieved through a combination of stabilization and common-mode rejection. The stabilization component will itself be achieved in two stages: pre-stabilization to a local optical cavity followed by arm-locking to some combination of the inter-spacecraft distances. In order for these two stabilization stages to work simultaneously, the lock-point of the pre-stabilization loop must be frequency tunable. The current baseline stabilization technique, locking to an optical cavity, does not provide tunability between cavity resonances, which are typically spaced by 100s of MHz. Here we present a modification to the traditional Pound-Drever-Hall cavity locking technique that allows the laser to be locked to a cavity resonance with an adjustable frequency offset. This technique requires no modifications to the optical cavity itself, thus preserving the stability of the frequency reference. We present measurements of the system performance and demonstrate that we can meet implement the first two stages of stabilization.

On the stability of dust orbits in mean-motion resonances perturbed by from an interstellar wind

NASA Astrophysics Data System (ADS)

Pástor, Pavol

2014-09-01

Circumstellar dust particles can be captured in a mean-motion resonance (MMR) with a planet and simultaneously be affected by non-gravitational effects. It is possible to describe the secular variations of a particle orbit in the MMR analytically using averaged resonant equations. We derive the averaged resonant equations from the equations of motion in near-canonical form. The secular variations of the particle orbit depending on the orientation of the orbit in space are taken into account. The averaged resonant equations can be derived/confirmed also from Lagrange's planetary equations. We apply the derived theory to the case when the non-gravitational effects are the Poynting-Robertson effect, the radial stellar wind, and an interstellar wind. The analytical and numerical results obtained are in excellent agreement. We found that the types of orbits correspond to libration centers of the conservative problem. The averaged resonant equations can lead to a system of equations which holds for stationary points in a subset of resonant variables. Using this system we show analytically that for the considered non-gravitational effects, all stationary points should correspond to orbits which are stationary in interplanetary space after an averaging over a synodic period. In an exact resonance, the stationary orbits are stable. The stability is achieved by a periodic repetition of the evolution during the synodic period. Numerical solutions of this system show that there are no stationary orbits for either the exact or non-exact resonances.

SHARP: A Spatially Higher-order, Relativistic Particle-in-cell Code

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

Shalaby, Mohamad; Broderick, Avery E.; Chang, Philip

Numerical heating in particle-in-cell (PIC) codes currently precludes the accurate simulation of cold, relativistic plasma over long periods, severely limiting their applications in astrophysical environments. We present a spatially higher-order accurate relativistic PIC algorithm in one spatial dimension, which conserves charge and momentum exactly. We utilize the smoothness implied by the usage of higher-order interpolation functions to achieve a spatially higher-order accurate algorithm (up to the fifth order). We validate our algorithm against several test problems—thermal stability of stationary plasma, stability of linear plasma waves, and two-stream instability in the relativistic and non-relativistic regimes. Comparing our simulations to exact solutionsmore » of the dispersion relations, we demonstrate that SHARP can quantitatively reproduce important kinetic features of the linear regime. Our simulations have a superior ability to control energy non-conservation and avoid numerical heating in comparison to common second-order schemes. We provide a natural definition for convergence of a general PIC algorithm: the complement of physical modes captured by the simulation, i.e., those that lie above the Poisson noise, must grow commensurately with the resolution. This implies that it is necessary to simultaneously increase the number of particles per cell and decrease the cell size. We demonstrate that traditional ways for testing for convergence fail, leading to plateauing of the energy error. This new PIC code enables us to faithfully study the long-term evolution of plasma problems that require absolute control of the energy and momentum conservation.« less

On the impact of atmospheric thermal stability on the characteristics of nocturnal downslope flows

NASA Astrophysics Data System (ADS)

Ye, Z. J.; Garratt, J. R.; Segal, M.; Pielke, R. A.

1990-04-01

The impacts of background (or ambient) and local atmospheric thermal stabilities, and slope steepness, on nighttime thermally induced downslope flow in meso-β domains (i.e., 20 200 km horizontal extent) have been investigated using analytical and numerical model approaches. Good agreement between the analytical and numerical evaluations was found. It was concluded that: (i) as anticipated, the intensity of the downslope flow increases with increased slope steepness, although the depth of the downslope flow was found to be insensitive to slope steepness in the studied situations; (ii) the intensity of the downslope flow is generally independent of background atmospheric thermal stability; (iii) for given integrated nighttime cooling across the nocturnal boundary layer (NBL), Q s the local atmospheric thermal stability exerts a strong influence on downslope flow behavior: the downslope flow intensity increases when local atmospheric thermal stability increases; and (iv) the downslope flow intensity is proportional to Q s 1/2.

Stability analysis of nonlinear Roesser-type two-dimensional systems via a homogenous polynomial technique

NASA Astrophysics Data System (ADS)

Zhang, Tie-Yan; Zhao, Yan; Xie, Xiang-Peng

2012-12-01

This paper is concerned with the problem of stability analysis of nonlinear Roesser-type two-dimensional (2D) systems. Firstly, the fuzzy modeling method for the usual one-dimensional (1D) systems is extended to the 2D case so that the underlying nonlinear 2D system can be represented by the 2D Takagi—Sugeno (TS) fuzzy model, which is convenient for implementing the stability analysis. Secondly, a new kind of fuzzy Lyapunov function, which is a homogeneous polynomially parameter dependent on fuzzy membership functions, is developed to conceive less conservative stability conditions for the TS Roesser-type 2D system. In the process of stability analysis, the obtained stability conditions approach exactness in the sense of convergence by applying some novel relaxed techniques. Moreover, the obtained result is formulated in the form of linear matrix inequalities, which can be easily solved via standard numerical software. Finally, a numerical example is also given to demonstrate the effectiveness of the proposed approach.

A multi-scale Q1/P0 approach to langrangian shock hydrodynamics.

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

Shashkov, Mikhail; Love, Edward; Scovazzi, Guglielmo

A new multi-scale, stabilized method for Q1/P0 finite element computations of Lagrangian shock hydrodynamics is presented. Instabilities (of hourglass type) are controlled by a stabilizing operator derived using the variational multi-scale analysis paradigm. The resulting stabilizing term takes the form of a pressure correction. With respect to currently implemented hourglass control approaches, the novelty of the method resides in its residual-based character. The stabilizing residual has a definite physical meaning, since it embeds a discrete form of the Clausius-Duhem inequality. Effectively, the proposed stabilization samples and acts to counter the production of entropy due to numerical instabilities. The proposed techniquemore » is applicable to materials with no shear strength, for which there exists a caloric equation of state. The stabilization operator is incorporated into a mid-point, predictor/multi-corrector time integration algorithm, which conserves mass, momentum and total energy. Encouraging numerical results in the context of compressible gas dynamics confirm the potential of the method.« less

A discrete time-varying internal model-based approach for high precision tracking of a multi-axis servo gantry.

PubMed

Zhang, Zhen; Yan, Peng; Jiang, Huan; Ye, Peiqing

2014-09-01

In this paper, we consider the discrete time-varying internal model-based control design for high precision tracking of complicated reference trajectories generated by time-varying systems. Based on a novel parallel time-varying internal model structure, asymptotic tracking conditions for the design of internal model units are developed, and a low order robust time-varying stabilizer is further synthesized. In a discrete time setting, the high precision tracking control architecture is deployed on a Voice Coil Motor (VCM) actuated servo gantry system, where numerical simulations and real time experimental results are provided, achieving the tracking errors around 3.5‰ for frequency-varying signals. Copyright © 2014 ISA. Published by Elsevier Ltd. All rights reserved.

Design and control strategies for CELSS - Integrating mechanistic paradigms and biological complexities

NASA Technical Reports Server (NTRS)

Moore, B., III; Kaufmann, R.; Reinhold, C.

1981-01-01

Systems analysis and control theory consideration are given to simulations of both individual components and total systems, in order to develop a reliable control strategy for a Controlled Ecological Life Support System (CELSS) which includes complex biological components. Because of the numerous nonlinearities and tight coupling within the biological component, classical control theory may be inadequate and the statistical analysis of factorial experiments more useful. The range in control characteristics of particular species may simplify the overall task by providing an appropriate balance of stability and controllability to match species function in the overall design. The ultimate goal of this research is the coordination of biological and mechanical subsystems in order to achieve a self-supporting environment.

Distributed Secure Coordinated Control for Multiagent Systems Under Strategic Attacks.

PubMed

Feng, Zhi; Wen, Guanghui; Hu, Guoqiang

2017-05-01

This paper studies a distributed secure consensus tracking control problem for multiagent systems subject to strategic cyber attacks modeled by a random Markov process. A hybrid stochastic secure control framework is established for designing a distributed secure control law such that mean-square exponential consensus tracking is achieved. A connectivity restoration mechanism is considered and the properties on attack frequency and attack length rate are investigated, respectively. Based on the solutions of an algebraic Riccati equation and an algebraic Riccati inequality, a procedure to select the control gains is provided and stability analysis is studied by using Lyapunov's method.. The effect of strategic attacks on discrete-time systems is also investigated. Finally, numerical examples are provided to illustrate the effectiveness of theoretical analysis.

CFD simulation of liquid-liquid dispersions in a stirred tank bioreactor

NASA Astrophysics Data System (ADS)

Gelves, R.

2013-10-01

In this paper simulations were developed in order to allow the examinations of drop sizes in liquid-liquid dispersions (oil-water) in a stirred tank bioreactor using CFD simulations (Computational Fluid Dynamics). The effects of turbulence, rotating flow, drop breakage were simulated by using the k-e, MRF (Multiple Reference Frame) and PBM (Population Balance Model), respectively. The numerical results from different operational conditions are compared with experimental data obtained from an endoscope technique and good agreement is achieved. Motivated by these simulated and experimental results CFD simulations are qualified as a very promising tool for predicting hydrodynamics and drop sizes especially useful for liquid-liquid applications which are characterized by the challenging problem of emulsion stability due to undesired drop sizes.

Adaptive PID formation control of nonholonomic robots without leader's velocity information.

PubMed

Shen, Dongbin; Sun, Weijie; Sun, Zhendong

2014-03-01

This paper proposes an adaptive proportional integral derivative (PID) algorithm to solve a formation control problem in the leader-follower framework where the leader robot's velocities are unknown for the follower robots. The main idea is first to design some proper ideal control law for the formation system to obtain a required performance, and then to propose the adaptive PID methodology to approach the ideal controller. As a result, the formation is achieved with much more enhanced robust formation performance. The stability of the closed-loop system is theoretically proved by Lyapunov method. Both numerical simulations and physical vehicle experiments are presented to verify the effectiveness of the proposed adaptive PID algorithm. Copyright © 2013 ISA. Published by Elsevier Ltd. All rights reserved.

Exponential H ∞ Synchronization of Chaotic Cryptosystems Using an Improved Genetic Algorithm

PubMed Central

Hsiao, Feng-Hsiag

2015-01-01

This paper presents a systematic design methodology for neural-network- (NN-) based secure communications in multiple time-delay chaotic (MTDC) systems with optimal H ∞ performance and cryptography. On the basis of the Improved Genetic Algorithm (IGA), which is demonstrated to have better performance than that of a traditional GA, a model-based fuzzy controller is then synthesized to stabilize the MTDC systems. A fuzzy controller is synthesized to not only realize the exponential synchronization, but also achieve optimal H ∞ performance by minimizing the disturbance attenuation level. Furthermore, the error of the recovered message is stated by using the n-shift cipher and key. Finally, a numerical example with simulations is given to demonstrate the effectiveness of our approach. PMID:26366432

Feedback control of combustion instabilities from within limit cycle oscillations using H∞ loop-shaping and the ν-gap metric

PubMed Central

Morgans, Aimee S.

2016-01-01

Combustion instabilities arise owing to a two-way coupling between acoustic waves and unsteady heat release. Oscillation amplitudes successively grow, until nonlinear effects cause saturation into limit cycle oscillations. Feedback control, in which an actuator modifies some combustor input in response to a sensor measurement, can suppress combustion instabilities. Linear feedback controllers are typically designed, using linear combustor models. However, when activated from within limit cycle, the linear model is invalid, and such controllers are not guaranteed to stabilize. This work develops a feedback control strategy guaranteed to stabilize from within limit cycle oscillations. A low-order model of a simple combustor, exhibiting the essential features of more complex systems, is presented. Linear plane acoustic wave modelling is combined with a weakly nonlinear describing function for the flame. The latter is determined numerically using a level set approach. Its implication is that the open-loop transfer function (OLTF) needed for controller design varies with oscillation level. The difference between the mean and the rest of the OLTFs is characterized using the ν-gap metric, providing the minimum required ‘robustness margin’ for an H∞ loop-shaping controller. Such controllers are designed and achieve stability both for linear fluctuations and from within limit cycle oscillations. PMID:27493558

Optimal Recursive Digital Filters for Active Bending Stabilization

NASA Technical Reports Server (NTRS)

Orr, Jeb S.

2013-01-01

In the design of flight control systems for large flexible boosters, it is common practice to utilize active feedback control of the first lateral structural bending mode so as to suppress transients and reduce gust loading. Typically, active stabilization or phase stabilization is achieved by carefully shaping the loop transfer function in the frequency domain via the use of compensating filters combined with the frequency response characteristics of the nozzle/actuator system. In this paper we present a new approach for parameterizing and determining optimal low-order recursive linear digital filters so as to satisfy phase shaping constraints for bending and sloshing dynamics while simultaneously maximizing attenuation in other frequency bands of interest, e.g. near higher frequency parasitic structural modes. By parameterizing the filter directly in the z-plane with certain restrictions, the search space of candidate filter designs that satisfy the constraints is restricted to stable, minimum phase recursive low-pass filters with well-conditioned coefficients. Combined with optimal output feedback blending from multiple rate gyros, the present approach enables rapid and robust parametrization of autopilot bending filters to attain flight control performance objectives. Numerical results are presented that illustrate the application of the present technique to the development of rate gyro filters for an exploration-class multi-engined space launch vehicle.

Experimental and numerical investigation on the ignition and combustion stability in solid fuel ramjet with swirling flow

NASA Astrophysics Data System (ADS)

Musa, Omer; Xiong, Chen; Changsheng, Zhou

2017-08-01

The present article investigates experimentally and numerically the ignition and flame stability of high-density polyethylene solid fuel with incoming swirling air through a solid fuel ramjet (SFRJ). A new design of swirler is proposed and used in this work. Experiments on connected pipes test facility were performed for SFRJ with and without swirl. An in-house code has been developed to simulate unsteady, turbulent, reacting, swirling flow in the SFRJ. Four different swirl intensities are utilized to study experimentally and numerically the effect of swirl number on the transient regression, ignition of the solid fuel in a hot-oxidizing flow and combustion phenomenon in the SFRJ. The results showed that using swirl flow decreases the ignition time delay, recirculation zone length, and the distance between the flame and the wall, meanwhile, increases the residence time, heat transfer, regression rate and mixing degree, thus, improving the combustion efficiency and stability.

Dynamic State Estimation for Multi-Machine Power System by Unscented Kalman Filter With Enhanced Numerical Stability

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

Qi, Junjian; Sun, Kai; Wang, Jianhui

In this paper, in order to enhance the numerical stability of the unscented Kalman filter (UKF) used for power system dynamic state estimation, a new UKF with guaranteed positive semidifinite estimation error covariance (UKFGPS) is proposed and compared with five existing approaches, including UKFschol, UKF-kappa, UKFmodified, UKF-Delta Q, and the squareroot UKF (SRUKF). These methods and the extended Kalman filter (EKF) are tested by performing dynamic state estimation on WSCC 3-machine 9-bus system and NPCC 48-machine 140-bus system. For WSCC system, all methods obtain good estimates. However, for NPCC system, both EKF and the classic UKF fail. It is foundmore » that UKFschol, UKF-kappa, and UKF-Delta Q do not work well in some estimations while UKFGPS works well in most cases. UKFmodified and SRUKF can always work well, indicating their better scalability mainly due to the enhanced numerical stability.« less

Mixed Arlequin method for multiscale poromechanics problems: Mixed Arlequin method for multiscale poromechanics problems

DOE PAGES

Sun, WaiChing; Cai, Zhijun; Choo, Jinhyun

2016-11-18

An Arlequin poromechanics model is introduced to simulate the hydro-mechanical coupling effects of fluid-infiltrated porous media across different spatial scales within a concurrent computational framework. A two-field poromechanics problem is first recast as the twofold saddle point of an incremental energy functional. We then introduce Lagrange multipliers and compatibility energy functionals to enforce the weak compatibility of hydro-mechanical responses in the overlapped domain. Here, to examine the numerical stability of this hydro-mechanical Arlequin model, we derive a necessary condition for stability, the twofold inf–sup condition for multi-field problems, and establish a modified inf–sup test formulated in the product space ofmore » the solution field. We verify the implementation of the Arlequin poromechanics model through benchmark problems covering the entire range of drainage conditions. Finally, through these numerical examples, we demonstrate the performance, robustness, and numerical stability of the Arlequin poromechanics model.« less

Stabilized high-order Galerkin methods based on a parameter-free dynamic SGS model for LES

NASA Astrophysics Data System (ADS)

Marras, Simone; Nazarov, Murtazo; Giraldo, Francis X.

2015-11-01

The high order spectral element approximation of the Euler equations is stabilized via a dynamic sub-grid scale model (Dyn-SGS). This model was originally designed for linear finite elements to solve compressible flows at large Mach numbers. We extend its application to high-order spectral elements to solve the Euler equations of low Mach number stratified flows. The major justification of this work is twofold: stabilization and large eddy simulation are achieved via one scheme only. Because the diffusion coefficients of the regularization stresses obtained via Dyn-SGS are residual-based, the effect of the artificial diffusion is minimal in the regions where the solution is smooth. The direct consequence is that the nominal convergence rate of the high-order solution of smooth problems is not degraded. To our knowledge, this is the first application in atmospheric modeling of a spectral element model stabilized by an eddy viscosity scheme that, by construction, may fulfill stabilization requirements, can model turbulence via LES, and is completely free of a user-tunable parameter. From its derivation, it will be immediately clear that Dyn-SGS is independent of the numerical method; it could be implemented in a discontinuous Galerkin, finite volume, or other environments alike. Preliminary discontinuous Galerkin results are reported as well. The straightforward extension to non-linear scalar problems is also described. A suite of 1D, 2D, and 3D test cases is used to assess the method, with some comparison against the results obtained with the most known Lilly-Smagorinsky SGS model.

An ACC Design Method for Achieving Both String Stability and Ride Comfort

NASA Astrophysics Data System (ADS)

Yamamura, Yoshinori; Seto, Yoji; Nishira, Hikaru; Kawabe, Taketoshi

An investigation was made of a method for designing adaptive cruise control (ACC) so as to achieve a headway distance response that feels natural to the driver while at the same time obtaining high levels of both string stability and ride comfort. With this design method, the H∞ norm is adopted as the index of string stability. Additionally, two norms are introduced for evaluating ride comfort and natural vehicle behavior. The relationship between these three norms and headway distance response characteristics was analyzed, and an evaluation method was established for achieving high levels of the various performance characteristics required of ACC. An ACC system designed with this method was evaluated in driving tests conducted on a proving ground course, and the results confirmed that it achieved the targeted levels of string stability, ride comfort and natural vehicle behavior.

Stability of large DC power systems using switching converters, with application to the international space station

NASA Technical Reports Server (NTRS)

Manners, B.; Gholdston, E. W.; Karimi, K.; Lee, F. C.; Rajagopalan, J.; Panov, Y.

1996-01-01

As space direct current (dc) power systems continue to grow in size, switching power converters are playing an ever larger role in power conditioning and control. When designing a large dc system using power converters of this type, special attention must be placed on the electrical stability of the system and of the individual loads on the system. In the design of the electric power system (EPS) of the International Space Station (ISS), the National Aeronautics and Space Administration (NASA) and its contractor team led by Boeing Defense & Space Group has placed a great deal of emphasis on designing for system and load stability. To achieve this goal, the team has expended considerable effort deriving a dear concept on defining system stability in both a general sense and specifically with respect to the space station. The ISS power system presents numerous challenges with respect to system stability, such as high power, complex sources and undefined loads. To complicate these issues, source and load components have been designed in parallel by three major subcontractors (Boeing, Rocketdyne, and McDonnell Douglas) with interfaces to both sources and loads being designed in different countries (Russia, Japan, Canada, Europe, etc.). These issues, coupled with the program goal of limiting costs, have proven a significant challenge to the program. As a result, the program has derived an impedance specification approach for system stability. This approach is based on the significant relationship between source and load impedances and the effect of this relationship on system stability. This approach is limited in its applicability by the theoretical and practical limits on component designs as presented by each system segment. As a result, the overall approach to system stability implemented by the ISS program consists of specific hardware requirements coupled with extensive system analysis and hardware testing. Following this approach, the ISS program plans to begin construction of the world's largest orbiting power system in 1997.

Relations of Different Types of Numerical Magnitude Representations to Each Other and to Mathematics Achievement

ERIC Educational Resources Information Center

Fazio, Lisa K.; Bailey, Drew H.; Thompson, Clarissa A.; Siegler, Robert S.

2014-01-01

We examined relations between symbolic and non-symbolic numerical magnitude representations, between whole number and fraction representations, and between these representations and overall mathematics achievement in fifth graders. Fraction and whole number symbolic and non-symbolic numerical magnitude understandings were measured using both…

Multiply scaled constrained nonlinear equation solvers. [for nonlinear heat conduction problems

NASA Technical Reports Server (NTRS)

Padovan, Joe; Krishna, Lala

1986-01-01

To improve the numerical stability of nonlinear equation solvers, a partitioned multiply scaled constraint scheme is developed. This scheme enables hierarchical levels of control for nonlinear equation solvers. To complement the procedure, partitioned convergence checks are established along with self-adaptive partitioning schemes. Overall, such procedures greatly enhance the numerical stability of the original solvers. To demonstrate and motivate the development of the scheme, the problem of nonlinear heat conduction is considered. In this context the main emphasis is given to successive substitution-type schemes. To verify the improved numerical characteristics associated with partitioned multiply scaled solvers, results are presented for several benchmark examples.

Stability of numerical integration techniques for transient rotor dynamics

NASA Technical Reports Server (NTRS)

Kascak, A. F.

1977-01-01

A finite element model of a rotor bearing system was analyzed to determine the stability limits of the forward, backward, and centered Euler; Runge-Kutta; Milne; and Adams numerical integration techniques. The analysis concludes that the highest frequency mode determines the maximum time step for a stable solution. Thus, the number of mass elements should be minimized. Increasing the damping can sometimes cause numerical instability. For a uniform shaft, with 10 mass elements, operating at approximately the first critical speed, the maximum time step for the Runge-Kutta, Milne, and Adams methods is that which corresponds to approximately 1 degree of shaft movement. This is independent of rotor dimensions.

The stagnation-point flow towards a shrinking sheet with homogeneous - heterogeneous reactions effects: A stability analysis

NASA Astrophysics Data System (ADS)

Ismail, Nurul Syuhada; Arifin, Norihan Md.; Bachok, Norfifah; Mahiddin, Norhasimah

2017-01-01

A numerical study is performed to evaluate the problem of stagnation - point flow towards a shrinking sheet with homogeneous - heterogeneous reaction effects. By using non-similar transformation, the governing equations be able to reduced to an ordinary differential equation. Then, results of the equations can be obtained numerically by shooting method with maple implementation. Based on the numerical results obtained, the velocity ratio parameter λ< 0, the dual solutions do exist. Then, the stability analysis is carried out to determine which solution is more stable between both of the solutions by bvp4c solver in Matlab.

An energy-saving nonlinear position control strategy for electro-hydraulic servo systems.

PubMed

Baghestan, Keivan; Rezaei, Seyed Mehdi; Talebi, Heidar Ali; Zareinejad, Mohammad

2015-11-01

The electro-hydraulic servo system (EHSS) demonstrates numerous advantages in size and performance compared to other actuation methods. Oftentimes, its utilization in industrial and machinery settings is limited by its inferior efficiency. In this paper, a nonlinear backstepping control algorithm with an energy-saving approach is proposed for position control in the EHSS. To achieve improved efficiency, two control valves including a proportional directional valve (PDV) and a proportional relief valve (PRV) are used to achieve the control objectives. To design the control algorithm, the state space model equations of the system are transformed to their normal form and the control law through the PDV is designed using a backstepping approach for position tracking. Then, another nonlinear set of laws is derived to achieve energy-saving through the PRV input. This control design method, based on the normal form representation, imposes internal dynamics on the closed-loop system. The stability of the internal dynamics is analyzed in special cases of operation. Experimental results verify that both tracking and energy-saving objectives are satisfied for the closed-loop system. Copyright © 2015 ISA. Published by Elsevier Ltd. All rights reserved.

Calculation of the lateral-dynamic stability of aircraft

NASA Technical Reports Server (NTRS)

Raikh, A

1952-01-01

Graphs and formulas are given with the aid of which all the aerodynamic coefficients required for computing the lateral dynamic stability can be determined. A number of numerical examples are given for obtaining the stability derivatives and solving the characteristic-stability equation. Approximate formulas are derived with the aid of which rapid preliminary computations may be made and the stability coefficients corrected for certain modifications of the airplane. A derivation of the lateral-dynamic-stability equations is included.

Numerical and Permeability Constraints on Simulation of Sill-Driven Hydrothermal Convection

NASA Astrophysics Data System (ADS)

Carr, P. M.; Cathles, L. M.; Barrie, C. T.; Manhardt, P.

2004-05-01

Volcanic-associated massive sulfide deposits are formed where seawater, heated to ~350oC by subsurface magma intrusions, is quenched by cold water at or near the seafloor. Many VMS districts, like the one at Matagami, Quebec, contain their zinc, lead, and copper in about a dozen discrete ore bodies, with one or two deposits containing more than half of the district's resources. We construct numerical models to investigate the causes of variations in deposit size. These models show that a process which stabilizes the location of hydrothermal venting plumes is required to numerically generate discrete VMS ore bodies by sill-driven hydrothermal convection. This is achieved in our models by increasing rock permeability in a fashion that makes vent plumes more permeable than their surroundings. Maintaining the Courant number ≤1 (so that a thermal anomaly traverses only one grid cell in one timestep of the simulation) is shown to be crucial to numerical convergence. If this rule is violated, visually compelling but incorrect hydrothermal vents result. Small hydrothermal convection cells over the interior of an areally-extensive sill with a tabular edge are smaller than those formed at the sill edge. However, for a sill with the geometry of that at Matagami, numerical simulations indicate that large ore deposits should form near the thickest part of the sill where metals extracted from the underside of the still-hot portions of the sill can optimally contribute. Thus it is essential to construct a model of the entire domain rather than slicing a portion local to the deposition. The numerical models replicate the ten-fold range in deposit size variation, and predict the largest deposits at Matagami will be discovered at 5 to 8 km depth between currently known deposits in the South Flank and Phelps Dodge areas.

Neural network-based sliding mode control for atmospheric-actuated spacecraft formation using switching strategy

NASA Astrophysics Data System (ADS)

Sun, Ran; Wang, Jihe; Zhang, Dexin; Shao, Xiaowei

2018-02-01

This paper presents an adaptive neural networks-based control method for spacecraft formation with coupled translational and rotational dynamics using only aerodynamic forces. It is assumed that each spacecraft is equipped with several large flat plates. A coupled orbit-attitude dynamic model is considered based on the specific configuration of atmospheric-based actuators. For this model, a neural network-based adaptive sliding mode controller is implemented, accounting for system uncertainties and external perturbations. To avoid invalidation of the neural networks destroying stability of the system, a switching control strategy is proposed which combines an adaptive neural networks controller dominating in its active region and an adaptive sliding mode controller outside the neural active region. An optimal process is developed to determine the control commands for the plates system. The stability of the closed-loop system is proved by a Lyapunov-based method. Comparative results through numerical simulations illustrate the effectiveness of executing attitude control while maintaining the relative motion, and higher control accuracy can be achieved by using the proposed neural-based switching control scheme than using only adaptive sliding mode controller.

Compatible diagonal-norm staggered and upwind SBP operators

NASA Astrophysics Data System (ADS)

Mattsson, Ken; O'Reilly, Ossian

2018-01-01

The main motivation with the present study is to achieve a provably stable high-order accurate finite difference discretisation of linear first-order hyperbolic problems on a staggered grid. The use of a staggered grid makes it non-trivial to discretise advective terms. To overcome this difficulty we discretise the advective terms using upwind Summation-By-Parts (SBP) operators, while the remaining terms are discretised using staggered SBP operators. The upwind and staggered SBP operators (for each order of accuracy) are compatible, here meaning that they are based on the same diagonal norms, allowing for energy estimates to be formulated. The boundary conditions are imposed using a penalty (SAT) technique, to guarantee linear stability. The resulting SBP-SAT approximations lead to fully explicit ODE systems. The accuracy and stability properties are demonstrated for linear hyperbolic problems in 1D, and for the 2D linearised Euler equations with constant background flow. The newly derived upwind and staggered SBP operators lead to significantly more accurate numerical approximations, compared with the exclusive usage of (previously derived) central-difference first derivative SBP operators.

Modeling of Feedback Stabilization of External MHD Modes in Toroidal Geometry

NASA Astrophysics Data System (ADS)

Chu, M. S.; Chance, M. S.; Okabayashi, M.

2000-10-01

The intelligent shell feedback scheme(C.M. Bishop, Plasma Phys. Contr. Nucl. Fusion 31), 1179 (1989). seeks to utilize external coils to suppress the unstable MHD modes slowed down by the resistive shell. We present a new formulation and numerical results of the interaction between the plasma and its outside vacuum region, with complete plasma response and the inclusion of a resistive vessel in general toroidal geometry. This is achieved by using the Green's function technique, which is a generalization of that previously used for the VACUUM(M.S. Chance, Phys. Plasmas 4), 2161 (1997). code and coupled with the ideal MHD code GATO. The effectiveness of different realizations of the intelligent shell concept is gauged by their ability to minimize the available free energy to drive the MHD mode. Computations indicate poloidal coverage of 30% of the total resistive wall surface area and 6 or 7 segments of ``intelligent coil'' arrays superimposed on the resistive wall will allow recovery of up to 90% the effectiveness of the ideal shell in stabilizing the ideal external kink.

Nonlinear Parametric Excitation Effect Induces Stability Transitions in Swimming Direction of Flexible Superparamagnetic Microswimmers.

PubMed

Harduf, Yuval; Jin, Dongdong; Or, Yizhar; Zhang, Li

2018-04-05

Microscopic artificial swimmers have recently become highly attractive due to their promising potential for biomedical microrobotic applications. Previous pioneering work has demonstrated the motion of a robotic microswimmer with a flexible chain of superparamagnetic beads, which is actuated by applying an oscillating external magnetic field. Interestingly, they have shown that the microswimmer's orientation undergoes a 90°-transition when the magnetic field's oscillation amplitude is increased above a critical value. This unexpected transition can cause severe problems in steering and manipulation of flexible magnetic microrobotic swimmers. Thus, theoretical understanding and analysis of the physical origins of this effect are of crucial importance. In this work, we investigate this transition both theoretically and experimentally by using numerical simulations and presenting a novel flexible microswimmer with an anisotropic superparamagnetic head. We prove that this effect depends on both frequency and amplitude of the oscillating magnetic field, and demonstrate existence of an optimal amplitude achieving maximal swimming speed. Asymptotic analysis of a minimal two-link model reveals that the changes in the swimmer's direction represent stability transitions, which are induced by a nonlinear parametric excitation.

A Riccati solution for the ideal MHD plasma response with applications to real-time stability control

NASA Astrophysics Data System (ADS)

Glasser, Alexander; Kolemen, Egemen; Glasser, A. H.

2018-03-01

Active feedback control of ideal MHD stability in a tokamak requires rapid plasma stability analysis. Toward this end, we reformulate the δW stability method with a Hamilton-Jacobi theory, elucidating analytical and numerical features of the generic tokamak ideal MHD stability problem. The plasma response matrix is demonstrated to be the solution of an ideal MHD matrix Riccati differential equation. Since Riccati equations are prevalent in the control theory literature, such a shift in perspective brings to bear a range of numerical methods that are well-suited to the robust, fast solution of control problems. We discuss the usefulness of Riccati techniques in solving the stiff ordinary differential equations often encountered in ideal MHD stability analyses—for example, in tokamak edge and stellarator physics. We demonstrate the applicability of such methods to an existing 2D ideal MHD stability code—DCON [A. H. Glasser, Phys. Plasmas 23, 072505 (2016)]—enabling its parallel operation in near real-time, with wall-clock time ≪1 s . Such speed may help enable active feedback ideal MHD stability control, especially in tokamak plasmas whose ideal MHD equilibria evolve with inductive timescale τ≳ 1s—as in ITER.

A note on a corrector formula for the numerical solution of ordinary differential equations

NASA Technical Reports Server (NTRS)

Chien, Y.-C.; Agrawal, K. M.

1979-01-01

A new corrector formula for predictor-corrector methods for numerical solutions of ordinary differential equations is presented. Two considerations for choosing corrector formulas are given: (1) the coefficient in the error term and (2) its stability properties. The graph of the roots of an equation plotted against its stability region, of different values, is presented along with the tables that correspond to various corrector equations, including Hamming's and Milne and Reynolds'.

A note on the effects of viscosity on the stability of a trailing-line vortex

NASA Technical Reports Server (NTRS)

Duck, Peter W.; Khorrami, Mehdi R.

1992-01-01

The linear stability of the Batchelor (1964) vortex is examined with emphasis on new viscous modes recently found numerically by Khorrami (1991). Unlike the previously reported inviscid modes of instability, these modes are destabilized by viscosity and exhibit small growth rates at large Reynolds numbers. The analysis presented here uses a combination of asymptotic and numerical techniques. The results confirm the existence of the additional modes of instability due to viscosity.

The use of the modified Cholesky decomposition in divergence and classification calculations

NASA Technical Reports Server (NTRS)

Van Rooy, D. L.; Lynn, M. S.; Snyder, C. H.

1973-01-01

This report analyzes the use of the modified Cholesky decomposition technique as applied to the feature selection and classification algorithms used in the analysis of remote sensing data (e.g., as in LARSYS). This technique is approximately 30% faster in classification and a factor of 2-3 faster in divergence, as compared with LARSYS. Also numerical stability and accuracy are slightly improved. Other methods necessary to deal with numerical stability problems are briefly discussed.

p21 stability: linking chaperones to a cell cycle checkpoint.

PubMed

Liu, Geng; Lozano, Guillermina

2005-02-01

Progression through the cell cycle is regulated by numerous proteins, one of which is the cyclin-dependent kinase inhibitor, p21. A new study identifies a novel protein complex that stabilizes p21. The stability of this complex is critical in effecting the p53-mediated cell cycle checkpoint.

Three-dimensional Cascaded Lattice Boltzmann Model for Thermal Convective Flows

NASA Astrophysics Data System (ADS)

Hajabdollahi, Farzaneh; Premnath, Kannan

2017-11-01

Fluid motion driven by thermal effects, such as due to buoyancy in differentially heated enclosures arise in several natural and industrial settings, whose understanding can be achieved via numerical simulations. Lattice Boltzmann (LB) methods are efficient kinetic computational approaches for coupled flow physics problems. In this study, we develop three-dimensional (3D) LB models based on central moments and multiple relaxation times for D3Q7 and D3Q15 lattices to solve the energy transport equations in a double distribution function approach. Their collision operators lead to a cascaded structure involving higher order terms resulting in improved stability. This is coupled to a central moment based LB flow solver with source terms. The new 3D cascaded LB models for the convective flows are first validated for natural convection of air driven thermally on two vertically opposite faces in a cubic cavity at different Rayleigh numbers against prior numerical and experimental data, which show good quantitative agreement. Then, the detailed structure of the 3D flow and thermal fields and the heat transfer rates at different Rayleigh numbers are analyzed and interpreted.

Dikin-type algorithms for dextrous grasping force optimization

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

Buss, M.; Faybusovich, L.; Moore, J.B.

1998-08-01

One of the central issues in dextrous robotic hand grasping is to balance external forces acting on the object and at the same time achieve grasp stability and minimum grasping effort. A companion paper shows that the nonlinear friction-force limit constraints on grasping forces are equivalent to the positive definiteness of a certain matrix subject to linear constraints. Further, compensation of the external object force is also a linear constraint on this matrix. Consequently, the task of grasping force optimization can be formulated as a problem with semidefinite constraints. In this paper, two versions of strictly convex cost functions, onemore » of them self-concordant, are considered. These are twice-continuously differentiable functions that tend to infinity at the boundary of possible definiteness. For the general class of such cost functions, Dikin-type algorithms are presented. It is shown that the proposed algorithms guarantee convergence to the unique solution of the semidefinite programming problem associated with dextrous grasping force optimization. Numerical examples demonstrate the simplicity of implementation, the good numerical properties, and the optimality of the approach.« less

Recent developments of nano-structured materials as the catalysts for oxygen reduction reaction

NASA Astrophysics Data System (ADS)

Kang, SungYeon; Kim, HuiJung; Chung, Yong-Ho

2018-04-01

Developments of high efficient materials for electrocatalyst are significant topics of numerous researches since a few decades. Recent global interests related with energy conversion and storage lead to the expansion of efforts to find cost-effective catalysts that can substitute conventional catalytic materials. Especially, in the field of fuel cell, novel materials for oxygen reduction reaction (ORR) have been noticed to overcome disadvantages of conventional platinum-based catalysts. Various approaching methods have been attempted to achieve low cost and high electrochemical activity comparable with Pt-based catalysts, including reducing Pt consumption by the formation of hybrid materials, Pt-based alloys, and not-Pt metal or carbon based materials. To enhance catalytic performance and stability, numerous methods such as structural modifications and complex formations with other functional materials are proposed, and they are basically based on well-defined and well-ordered catalytic active sites by exquisite control at nanoscale. In this review, we highlight the development of nano-structured catalytic materials for ORR based on recent findings, and discuss about an outlook for the direction of future researches.

Stabilized finite element methods to simulate the conductances of ion channels

NASA Astrophysics Data System (ADS)

Tu, Bin; Xie, Yan; Zhang, Linbo; Lu, Benzhuo

2015-03-01

We have previously developed a finite element simulator, ichannel, to simulate ion transport through three-dimensional ion channel systems via solving the Poisson-Nernst-Planck equations (PNP) and Size-modified Poisson-Nernst-Planck equations (SMPNP), and succeeded in simulating some ion channel systems. However, the iterative solution between the coupled Poisson equation and the Nernst-Planck equations has difficulty converging for some large systems. One reason we found is that the NP equations are advection-dominated diffusion equations, which causes troubles in the usual FE solution. The stabilized schemes have been applied to compute fluids flow in various research fields. However, they have not been studied in the simulation of ion transport through three-dimensional models based on experimentally determined ion channel structures. In this paper, two stabilized techniques, the SUPG and the Pseudo Residual-Free Bubble function (PRFB) are introduced to enhance the numerical robustness and convergence performance of the finite element algorithm in ichannel. The conductances of the voltage dependent anion channel (VDAC) and the anthrax toxin protective antigen pore (PA) are simulated to validate the stabilization techniques. Those two stabilized schemes give reasonable results for the two proteins, with decent agreement with both experimental data and Brownian dynamics (BD) simulations. For a variety of numerical tests, it is found that the simulator effectively avoids previous numerical instability after introducing the stabilization methods. Comparison based on our test data set between the two stabilized schemes indicates both SUPG and PRFB have similar performance (the latter is slightly more accurate and stable), while SUPG is relatively more convenient to implement.

Enhancement of numeric cognition in children with low achievement in mathematic after a non-instrumental musical training.

PubMed

Ribeiro, Fabiana Silva; Santos, Flávia H

2017-03-01

Studies suggest that musical training enhances spatial-temporal reasoning and leads to greater learning of mathematical concepts. The aim of this prospective study was to verify the efficacy of a Non-Instrumental Musical Training (NIMT) on the Numerical Cognition systems in children with low achievement in math. For this purpose, we examined, with a cluster analysis, whether children with low scores on Numerical Cognition would be grouped in the same cluster at pre and post-NIMT. Participants were primary school children divided into two groups according to their scores on an Arithmetic test. Results with a specialized battery of Numerical Cognition revealed improvements for Cluster 2 (children with low achievement in math) especially for number production capacity compared to normative data. Besides, the number of children with low scores in Numerical Cognition decreased at post-NIMT. These findings suggest that NIMT enhances Numerical Cognition and seems to be a useful tool for rehabilitation of children with low achievement in math. Copyright © 2016 Elsevier Ltd. All rights reserved.

Stability Analysis of Finite Difference Schemes for Hyperbolic Systems, and Problems in Applied and Computational Linear Algebra.

DTIC Science & Technology

FINITE DIFFERENCE THEORY, * LINEAR ALGEBRA , APPLIED MATHEMATICS, APPROXIMATION(MATHEMATICS), BOUNDARY VALUE PROBLEMS, COMPUTATIONS, HYPERBOLAS, MATHEMATICAL MODELS, NUMERICAL ANALYSIS, PARTIAL DIFFERENTIAL EQUATIONS, STABILITY.

High- β equilibrium and ballooning stability of the low aspect ratio CNT stellarator

DOE PAGES

Hammond, K. C.; Lazerson, S. A.; Volpe, F. A.

2017-04-07

In the paper, the existence and ballooning-stability of low aspect ratio stellarator equilibria is predicted for the Columbia Neutral Torus (CNT) with the aid of 3D numerical tools. In addition to having a low aspect ratio, CNT is characterized by a low magnetic field and small plasma volume. Also, highly overdense plasmas were recently heated in CNT by means of microwaves. These characteristics suggest that CNT might attain relatively high values of plasma beta and thus be of use in the experimental study of stellarator stability to high-beta instabilities such as ballooning modes. As a first step in that direction,more » here the ballooning stability limit is found numerically. Depending on the particular magnetic configuration we expect volume-averaged β limits in the range 0.9%–3.0%, and possibly higher, and observe indications of a second region of ballooning stability. As the aspect ratio is reduced, stability is found to increase in some configurations and decrease in others. Energy-balance estimates using stellarator scaling laws indicate that the lower β limit may be attainable with overdense heating at powers of 40 to 100 kW. The present study serves the additional purpose of testing VMEC and other stellarator codes at high values of β and at low aspect ratios. For this reason, the study was carried out both for free boundary, for maximum fidelity to experiment, as well as with a fixed boundary, as a numerical test.« less

Control of nonlinear systems using periodic parametric perturbations with application to a reversed field pinch

NASA Astrophysics Data System (ADS)

Mirus, Kevin Andrew

In this thesis, the possibility of controlling low- and high-dimensional chaotic systems by periodically driving an accessible system parameter is examined. This method has been carried out on several numerical systems and the MST Reversed Field Pinch. The numerical systems investigated include the logistic equation, the Lorenz equations, the Rossler equations, a coupled lattice of logistic equations, a coupled lattice of Lorenz equations, the Yoshida equations, which model tearing mode fluctuations in a plasma, and a neural net model for magnetic fluctuations on MST. This method was tested on the MST by sinusoidally driving a magnetic flux through the toroidal gap of the device. Numerically, periodic drives were found to be most effective at producing limit cycle behavior or significantly reducing the dimension of the system when the perturbation frequency was near natural frequencies of unstable periodic orbits embedded in the attractor of the unperturbed system. Several different unstable periodic orbits have been stabilized in this way for the low-dimensional numerical systems, sometimes with perturbation amplitudes that were less than 5% of the nominal value of the parameter being perturbed. In high- dimensional systems, limit cycle behavior and significant decreases in the system dimension were also achieved using perturbations with frequencies near the natural unstable periodic orbit frequencies. Results for the MST were not this encouraging, most likely because of an insufficient drive amplitude, the extremely high dimension of the plasma behavior, large amounts of noise, and a lack of stationarity in the transient plasma pulses.

Applications of Laplace transform methods to airfoil motion and stability calculations

NASA Technical Reports Server (NTRS)

Edwards, J. W.

1979-01-01

This paper reviews the development of generalized unsteady aerodynamic theory and presents a derivation of the generalized Possio integral equation. Numerical calculations resolve questions concerning subsonic indicial lift functions and demonstrate the generation of Kutta waves at high values of reduced frequency, subsonic Mach number, or both. The use of rational function approximations of unsteady aerodynamic loads in aeroelastic stability calculations is reviewed, and a reformulation of the matrix Pade approximation technique is given. Numerical examples of flutter boundary calculations for a wing which is to be flight tested are given. Finally, a simplified aerodynamic model of transonic flow is used to study the stability of an airfoil exposed to supersonic and subsonic flow regions.

Model of a Frame of Dynamic Routing and Its Equilibrium

NASA Astrophysics Data System (ADS)

Zhang, Shu; Yuan, Yuan; Xu, Jian

Dynamic routing algorithm based on the shortest path principle is criticized due to the oscillation induced by such routing scheme. In the present work, we propose the model of TCP/RED algorithm by a new frame of dynamic routing, based on the measurement of occupation ratio of router buffer for different links, which only requires the information of the queue size at the buffer of the router, to stabilize the system. We classify several types of equilibrium and employ the numerical method to study the stability of the steady state. Our numerical results show that the careful selection of the parameters characterizing the dynamic routing algorithm can stabilize the system in some cases.

Numerical implementation of complex orthogonalization, parallel transport on Stiefel bundles, and analyticity

NASA Astrophysics Data System (ADS)

Avitabile, Daniele; Bridges, Thomas J.

2010-06-01

Numerical integration of complex linear systems of ODEs depending analytically on an eigenvalue parameter are considered. Complex orthogonalization, which is required to stabilize the numerical integration, results in non-analytic systems. It is shown that properties of eigenvalues are still efficiently recoverable by extracting information from a non-analytic characteristic function. The orthonormal systems are constructed using the geometry of Stiefel bundles. Different forms of continuous orthogonalization in the literature are shown to correspond to different choices of connection one-form on the Stiefel bundle. For the numerical integration, Gauss-Legendre Runge-Kutta algorithms are the principal choice for preserving orthogonality, and performance results are shown for a range of GLRK methods. The theory and methods are tested by application to example boundary value problems including the Orr-Sommerfeld equation in hydrodynamic stability.

On the effects of viscosity on the stability of a trailing-line vortex

NASA Technical Reports Server (NTRS)

Duck, Peter W.; Khorrami, Mehdi R.

1991-01-01

The linear stability of the Batchelor (1964) vortex is investigated. Particular emphasis is placed on modes found recently in a numerical study by Khorrami (1991). These modes have a number of features very distinct from those found previously for this vortex, including exhibiting small growth rates at large Reynolds numbers and susceptibility to destabilization by viscosity. These modes are described using asymptotic techniques, producing results which compare favorably with fully numerical results at large Reynolds numbers.

On dynamics in a Keynesian model of monetary stabilization policy with debt effect

NASA Astrophysics Data System (ADS)

Asada, Toichiro; Demetrian, Michal; Zimka, Rudolf

2018-05-01

In this paper, a four-dimensional model of flexible prices with the central bank's stabilization policy, describing the development of the firms' private debt, the output, the expected rate of inflation and the rate of interest is analyzed. Questions concerning the existence of limit cycles around its normal equilibrium point are investigated. The bifurcation equation is found. The formulae for the calculation of its coefficients are gained. A numerical example is presented by means of numerical simulations.

A numerical study of variable density flow and mixing in porous media

NASA Astrophysics Data System (ADS)

Fan, Yin; Kahawita, René

1994-10-01

A numerical study of a negatively buoyant plume intruding into a neutrally stratified porous medium has been undertaken using finite different methods. Of particular interest has been to ascertain whether the experimentally observed gravitational instabilities that form along the lower edge of the plume are reproduced in the numerical model. The model has been found to faithfully reproduce the mean flow as well as the gravitational instabilities in the intruding plume. A linear stability analysis has confirmed the fact that the negatively buoyant plume is in fact gravitationally unstable and that the stability depends on two parameters: a concentration Rayleigh number and a characteristic length scale which is dependent on the transverse dispersivity.

Problem-Based Instructional Strategy and Numerical Ability as Determinants of Senior Secondary Achievement in Mathematics

ERIC Educational Resources Information Center

Badru, Ademola K.

2016-01-01

The study investigated Problem-based Instructional Strategy and Numerical ability as determinants of Senior Secondary Achievement in Mathematics. This study used 4 x 2 x 2 non-randomised control group Pretest-Posttest Quasi-experimental Factorial design. It consisted of two independent variables (treatment and Numerical ability) and one moderating…

Micromachined silicon acoustic delay line with improved structural stability and acoustic directivity for real-time photoacoustic tomography

NASA Astrophysics Data System (ADS)

Cho, Young; Kumar, Akhil; Xu, Song; Zou, Jun

2017-03-01

Recent studies have shown that micromachined silicon acoustic delay lines can provide a promising solution to achieve real-time photoacoustic tomography without the need for complex transducer arrays and data acquisition electronics. However, as its length increases to provide longer delay time, the delay line becomes more vulnerable to structural instability due to reduced mechanical stiffness. In addition, the small cross-section area of the delay line results in a large acoustic acceptance angle and therefore poor directivity. To address these two issues, this paper reports the design, fabrication, and testing of a new silicon acoustic delay line enhanced with 3D printed polymer micro linker structures. First, mechanical deformation of the silicon acoustic delay line (with and without linker structures) under gravity was simulated by using finite element method. Second, the acoustic crosstalk and acoustic attenuation caused by the polymer micro linker structures were evaluated with both numerical simulation and ultrasound transmission testing. The result shows that the use of the polymer micro linker structures significantly improves the structural stability of the silicon acoustic delay lines without creating additional acoustic attenuation and crosstalk. In addition, a new tapered design for the input terminal of the delay line was also investigate to improve its acoustic directivity by reducing the acoustic acceptance angle. These two improvements are expected to provide an effective solution to eliminate current limitations on the achievable acoustic delay time and out-of-plane imaging resolution of micromachined silicon acoustic delay line arrays.

The Measurement Structure, Stability and Mediating Effects of Achievement Goals in Math with Middle-School Student Data

ERIC Educational Resources Information Center

Chiang, Yu-Tzu; Lin, Sunny S. J.

2014-01-01

This study examined the measurement structure, cross-year stability of achievement goals, and mediating effects of achievement goals between self-efficacy and math grades in a national sample of Taiwan middle school students. The measurement model with factorial structure showed good fit to the data. In the panel data (N?=?343), four achievement…

Numerical simulation of a two-sex human papillomavirus (HPV) vaccination model

NASA Astrophysics Data System (ADS)

Suryani, I.; Adi-Kusumo, F.

2014-02-01

Human Papillomavirus (HPV) is a major cause of cervical cancer, precancerous lesions, cancer and other disease. HPV is the most common sexually transmitted infection. Although HPV virus primarily affects woman but it can also affects man because it cause of cancer of the anus, vulva, vagina, penis and some other cancers. HPV vaccines now used to prevent cervical cancer and genital warts because the vaccine protect against four types of HPV that most commonly cause disease are types 6, 11, 16, and 18. This paper is sequel work of Elbasha (2008). Difference with Elbasha (2008) are give alternative proof global stability, numerical simulation and interpretation. Global stability of the equilibrium on the model of a two-sex HPV vaccination were explored by using Lyapunov. Although we use the same lyapunov function, we use the largest invariant set to proof the global stability. The result show that the global stability of the equilibrium depends on the effective reproduction number (R). If R < 1 then the infection-free equilibrium is asymptotically stable globally. If R > 1 then endemic equilibrium have globally asymptotically stable properties. Then equilibrium proceed with the interpretation of numerical simulation.

Basic numerical competences in large-scale assessment data: Structure and long-term relevance.

PubMed

Hirsch, Stefa; Lambert, Katharina; Coppens, Karien; Moeller, Korbinian

2018-03-01

Basic numerical competences are seen as building blocks for later numerical and mathematical achievement. The current study aimed at investigating the structure of early numeracy reflected by different basic numerical competences in kindergarten and its predictive value for mathematical achievement 6 years later using data from large-scale assessment. This allowed analyses based on considerably large sample sizes (N > 1700). A confirmatory factor analysis indicated that a model differentiating five basic numerical competences at the end of kindergarten fitted the data better than a one-factor model of early numeracy representing a comprehensive number sense. In addition, these basic numerical competences were observed to reliably predict performance in a curricular mathematics test in Grade 6 even after controlling for influences of general cognitive ability. Thus, our results indicated a differentiated view on early numeracy considering basic numerical competences in kindergarten reflected in large-scale assessment data. Consideration of different basic numerical competences allows for evaluating their specific predictive value for later mathematical achievement but also mathematical learning difficulties. Copyright © 2017 Elsevier Inc. All rights reserved.

Formation flying design and applications in weak stability boundary regions.

PubMed

Folta, David

2004-05-01

Weak stability regions serve as superior locations for interferomertric scientific investigations. These regions are often selected to minimize environmental disturbances and maximize observation efficiency. Designs of formations in these regions are becoming ever more challenging as more complex missions are envisioned. The development of algorithms to enable the capability for formation design must be further enabled to incorporate better understanding of weak stability boundary solution space. This development will improve the efficiency and expand the capabilities of current approaches. The Goddard Space Flight Center (GSFC) is currently supporting multiple formation missions in weak stability boundary regions. This end-to-end support consists of mission operations, trajectory design, and control. It also includes both algorithm and software development. The Constellation-X, Maxim, and Stellar Imager missions are examples of the use of improved numeric methods to attain constrained formation geometries and control their dynamical evolution. This paper presents a survey of formation missions in the weak stability boundary regions and a brief description of formation design using numerical and dynamical techniques.

Long-time stability effects of quadrature and artificial viscosity on nodal discontinuous Galerkin methods for gas dynamics

NASA Astrophysics Data System (ADS)

Durant, Bradford; Hackl, Jason; Balachandar, Sivaramakrishnan

2017-11-01

Nodal discontinuous Galerkin schemes present an attractive approach to robust high-order solution of the equations of fluid mechanics, but remain accompanied by subtle challenges in their consistent stabilization. The effect of quadrature choices (full mass matrix vs spectral elements), over-integration to manage aliasing errors, and explicit artificial viscosity on the numerical solution of a steady homentropic vortex are assessed over a wide range of resolutions and polynomial orders using quadrilateral elements. In both stagnant and advected vortices in periodic and non-periodic domains the need arises for explicit stabilization beyond the numerical surface fluxes of discontinuous Galerkin spectral elements. Artificial viscosity via the entropy viscosity method is assessed as a stabilizing mechanism. It is shown that the regularity of the artificial viscosity field is essential to its use for long-time stabilization of small-scale features in nodal discontinuous Galerkin solutions of the Euler equations of gas dynamics. Supported by the Department of Energy Predictive Science Academic Alliance Program Contract DE-NA0002378.

On the manifestation of coexisting nontrivial equilibria leading to potential well escapes in an inhomogeneous floating body

NASA Astrophysics Data System (ADS)

Sequeira, Dane; Wang, Xue-She; Mann, B. P.

2018-02-01

This paper examines the bifurcation and stability behavior of inhomogeneous floating bodies, specifically a rectangular prism with asymmetric mass distribution. A nonlinear model is developed to determine the stability of the upright and tilted equilibrium positions as a function of the vertical position of the center of mass within the prism. These equilibria positions are defined by an angle of rotation and a vertical position where rotational motion is restricted to a two dimensional plane. Numerical investigations are conducted using path-following continuation methods to determine equilibria solutions and evaluate stability. Bifurcation diagrams and basins of attraction that illustrate the stability of the equilibrium positions as a function of the vertical position of the center of mass within the prism are generated. These results reveal complex stability behavior with many coexisting solutions. Static experiments are conducted to validate equilibria orientations against numerical predictions with results showing good agreement. Dynamic experiments that examine potential well hopping behavior in a waveflume for various wave conditions are also conducted.

Stability of soliton families in nonlinear Schrödinger equations with non-parity-time-symmetric complex potentials

NASA Astrophysics Data System (ADS)

Yang, Jianke; Nixon, Sean

2016-11-01

Stability of soliton families in one-dimensional nonlinear Schrödinger equations with non-parity-time (PT)-symmetric complex potentials is investigated numerically. It is shown that these solitons can be linearly stable in a wide range of parameter values both below and above phase transition. In addition, a pseudo-Hamiltonian-Hopf bifurcation is revealed, where pairs of purely-imaginary eigenvalues in the linear-stability spectra of solitons collide and bifurcate off the imaginary axis, creating oscillatory instability, which resembles Hamiltonian-Hopf bifurcations of solitons in Hamiltonian systems even though the present system is dissipative and non-Hamiltonian. The most important numerical finding is that, eigenvalues of linear-stability operators of these solitons appear in quartets (λ , - λ ,λ* , -λ*), similar to conservative systems and PT-symmetric systems. This quartet eigenvalue symmetry is very surprising for non- PT-symmetric systems, and it has far-reaching consequences on the stability behaviors of solitons.

Langley Stability and Transition Analysis Code (LASTRAC) Version 1.2 User Manual

NASA Technical Reports Server (NTRS)

Chang, Chau-Lyan

2004-01-01

LASTRAC is a general-purposed, physics-based transition prediction code released by NASA for Laminar Flow Control studies and transition research. The design and development of the LASTRAC code is aimed at providing an engineering tool that is easy to use and yet capable of dealing with a broad range of transition related issues. It was written from scratch based on the state-of-the-art numerical methods for stability analysis and modern software technologies. At low fidelity, it allows users to perform linear stability analysis and N-factor transition correlation for a broad range of flow regimes and configurations by using either the linear stability theory or linear parabolized stability equations method. At high fidelity, users may use nonlinear PSE to track finite-amplitude disturbances until the skin friction rise. This document describes the governing equations, numerical methods, code development, detailed description of input/output parameters, and case studies for the current release of LASTRAC.

Orthogonal optimization of a water hydraulic pilot-operated pressure-reducing valve

NASA Astrophysics Data System (ADS)

Mao, Xuyao; Wu, Chao; Li, Bin; Wu, Di

2017-12-01

In order to optimize the comprehensive characteristics of a water hydraulic pilot-operated pressure-reducing valve, numerical orthogonal experimental design was adopted. Six parameters of the valve, containing diameters of damping plugs, volume of spring chamber, half cone angle of main spool, half cone angle of pilot spool, mass of main spool and diameter of main spool, were selected as the orthogonal factors, and each factor has five different levels. An index of flowrate stability, pressure stability and pressure overstrike stability (iFPOS) was used to judge the merit of each orthogonal attempt. Embedded orthogonal process turned up and a final optimal combination of these parameters was obtained after totally 50 numerical orthogonal experiments. iFPOS could be low to a fairly low value which meant that the valve could have much better stabilities. During the optimization, it was also found the diameters of damping plugs and main spool played important roles in stability characteristics of the valve.

Numerical investigation of stability of breather-type solutions of the nonlinear Schrödinger equation

NASA Astrophysics Data System (ADS)

Calini, A.; Schober, C. M.

2013-09-01

In this article we present the results of a broad numerical investigation on the stability of breather-type solutions of the nonlinear Schrödinger (NLS) equation, specifically the one- and two-mode breathers for an unstable plane wave, which are frequently used to model rogue waves. The numerical experiments involve large ensembles of perturbed initial data for six typical random perturbations. Ensemble estimates of the "closeness", A(t), of the perturbed solution to an element of the respective unperturbed family indicate that the only neutrally stable breathers are the ones of maximal dimension, that is: given an unstable background with N unstable modes, the only neutrally stable breathers are the N-dimensional ones (obtained as a superimposition of N simple breathers via iterated Backlund transformations). Conversely, breathers which are not fully saturated are sensitive to noisy environments and are unstable. Interestingly, A(t) is smallest for the coalesced two-mode breather indicating the coalesced case may be the most robust two-mode breather in a laboratory setting. The numerical simulations confirm and provide a realistic realization of the stability behavior established analytically by the authors.