#### Sample records for achieve numerical stability

1. Examining the Stability of Achievement Goal Orientation

ERIC Educational Resources Information Center

Muis, Krista R.; Edwards, Ordene

2009-01-01

We examined stability and change in students' achievement goal orientations over varying tasks. Two naturalistic longitudinal studies were conducted in undergraduate courses. Students completed self-reports designed to measure their achievement goals. Achievement goals were measured four times: prior to two assignments and two exams. Four…

2. Numerical Tribute to Achievement of Euler

Figueroa-Navarro, Carlos; Molinar-Tabares, Martín Eduardo; Castro-Arce, Lamberto; Campos-García, Julio Cesar

2014-03-01

This work aims to make a tribute to one of the world's brightest personalities as it was the mathematical physicist Leonhard Euler (1707-1783). Some results where the influence of Euler persists with the novelty of applying numerical analysis using Matlab are here exposed. A first analysis was done with the series that defines Euler numbers and polynomials of Frobenius-Euler; another result is the characterization of the functions that carry to Euler-Macheroni constant. In hydrodynamics is also feasible to evaluate graphically the relationship between dimensions in diameter and the exit angle of the height of Euler for turbomachines. In differential equations of Cauchy-Euler solutions for the cases of distinct real roots and complex roots are generated. Furthermore we report the generation of the Fourier series and the Fourier transform calculated by using Direct Commands of Matlab. In variational calculus it is possible to obtain plots from a problem of the Euler Lagrange equations. Finally, the Euler function is analyzed. Our purpose is to present a tribute to this giant of science also it could be an excuse to study his legacy by utilizing modern computational techniques.

3. Numerical experiments on unstructured PIC stability.

SciTech Connect

Day, David Minot

2011-04-01

Particle-In-Cell (PIC) is a method for plasmas simulation. Particles are pushed with Verlet time integration. Fields are modeled using finite differences on a tensor product mesh (cells). The Unstructured PIC methods studied here use instead finite element discretizations on unstructured (simplicial) meshes. PIC is constrained by stability limits (upper bounds) on mesh and time step sizes. Numerical evidence (2D) and analysis will be presented showing that similar bounds constrain unstructured PIC.

4. Numerical robust stability estimation in milling process

Zhang, Xiaoming; Zhu, Limin; Ding, Han; Xiong, Youlun

2012-09-01

The conventional prediction of milling stability has been extensively studied based on the assumptions that the milling process dynamics is time invariant. However, nominal cutting parameters cannot guarantee the stability of milling process at the shop floor level since there exists many uncertain factors in a practical manufacturing environment. This paper proposes a novel numerical method to estimate the upper and lower bounds of Lobe diagram, which is used to predict the milling stability in a robust way by taking into account the uncertain parameters of milling system. Time finite element method, a milling stability theory is adopted as the conventional deterministic model. The uncertain dynamics parameters are dealt with by the non-probabilistic model in which the parameters with uncertainties are assumed to be bounded and there is no need for probabilistic distribution densities functions. By doing so, interval instead of deterministic stability Lobe is obtained, which guarantees the stability of milling process in an uncertain milling environment. In the simulations, the upper and lower bounds of Lobe diagram obtained by the changes of modal parameters of spindle-tool system and cutting coefficients are given, respectively. The simulation results show that the proposed method is effective and can obtain satisfying bounds of Lobe diagrams. The proposed method is helpful for researchers at shop floor to making decision on machining parameters selection.

5. Numerical stability in problems of linear algebra.

NASA Technical Reports Server (NTRS)

Babuska, I.

1972-01-01

Mathematical problems are introduced as mappings from the space of input data to that of the desired output information. Then a numerical process is defined as a prescribed recurrence of elementary operations creating the mapping of the underlying mathematical problem. The ratio of the error committed by executing the operations of the numerical process (the roundoff errors) to the error introduced by perturbations of the input data (initial error) gives rise to the concept of lambda-stability. As examples, several processes are analyzed from this point of view, including, especially, old and new processes for solving systems of linear algebraic equations with tridiagonal matrices. In particular, it is shown how such a priori information can be utilized as, for instance, a knowledge of the row sums of the matrix. Information of this type is frequently available where the system arises in connection with the numerical solution of differential equations.

6. RELAP-7 Numerical Stabilization: Entropy Viscosity Method

SciTech Connect

R. A. Berry; M. O. Delchini; J. Ragusa

2014-06-01

The RELAP-7 code is the next generation nuclear reactor system safety analysis code being developed at the Idaho National Laboratory (INL). The code is based on the INL's modern scientific software development framework, MOOSE (Multi-Physics Object Oriented Simulation Environment). The overall design goal of RELAP-7 is to take advantage of the previous thirty years of advancements in computer architecture, software design, numerical integration methods, and physical models. The end result will be a reactor systems analysis capability that retains and improves upon RELAP5's capability and extends the analysis capability for all reactor system simulation scenarios. RELAP-7 utilizes a single phase and a novel seven-equation two-phase flow models as described in the RELAP-7 Theory Manual (INL/EXT-14-31366). The basic equation systems are hyperbolic, which generally require some type of stabilization (or artificial viscosity) to capture nonlinear discontinuities and to suppress advection-caused oscillations. This report documents one of the available options for this stabilization in RELAP-7 -- a new and novel approach known as the entropy viscosity method. Because the code is an ongoing development effort in which the physical sub models, numerics, and coding are evolving, so too must the specific details of the entropy viscosity stabilization method. Here the fundamentals of the method in their current state are presented.

7. Mapping numerical magnitudes onto symbols: the numerical distance effect and individual differences in children's mathematics achievement.

PubMed

Holloway, Ian D; Ansari, Daniel

2009-05-01

Although it is often assumed that abilities that reflect basic numerical understanding, such as numerical comparison, are related to children's mathematical abilities, this relationship has not been tested rigorously. In addition, the extent to which symbolic and nonsymbolic number processing play differential roles in this relationship is not yet understood. To address these questions, we collected mathematics achievement measures from 6- to 8-year-olds as well as reaction times from a numerical comparison task. Using the reaction times, we calculated the size of the numerical distance effect exhibited by each child. In a correlational analysis, we found that the individual differences in the distance effect were related to mathematics achievement but not to reading achievement. This relationship was found to be specific to symbolic numerical comparison. Implications for the role of basic numerical competency and the role of accessing numerical magnitude information from Arabic numerals for the development of mathematical skills and their impairment are discussed. PMID:18513738

8. Achieving stabilization in interferometric logic operations.

PubMed

Zavalin, Andrey I; Shamir, Joseph; Vikram, Chandra S; Caulfield, H John

2006-01-10

Interferometric systems with amplitude beam splitters can implement reversible operations that, on detection, become Boolean operators. Being passive, they consume no energy, do not limit the operating bandwidth, and have negligible latency. Unfortunately, conventional interferometric systems are notoriously sensitive to uncontrolled disturbances. Here the use of polarization in a common-path interferometric logic gate with and without polarization beam splitters is explored as an attractive alternative to overcome those difficulties. Two of three device configurations considered offer significant stability and lower drive modulator voltage as advantages over the previous systems. The first experimental tests of such a system are reported. Common-path interferometry lends itself to even more stability and robustness by compatibility with no-air-gap, solid optics. PMID:16422166

9. Achieving stabilization in interferometric logic operations

Zavalin, Andrey I.; Shamir, Joseph; Vikram, Chandra S.; Caulfield, H. John

2006-01-01

Interferometric systems with amplitude beam splitters can implement reversible operations that, on detection, become Boolean operators. Being passive, they consume no energy, do not limit the operating bandwidth, and have negligible latency. Unfortunately, conventional interferometric systems are notoriously sensitive to uncontrolled disturbances. Here the use of polarization in a common-path interferometric logic gate with and without polarization beam splitters is explored as an attractive alternative to overcome those difficulties. Two of three device configurations considered offer significant stability and lower drive modulator voltage as advantages over the previous systems. The first experimental tests of such a system are reported. Common-path interferometry lends itself to even more stability and robustness by compatibility with no-air-gap, solid optics.

10. Achieving better cooling of turbine blades using numerical simulation methods

Inozemtsev, A. A.; Tikhonov, A. S.; Sendyurev, C. I.; Samokhvalov, N. Yu.

2013-02-01

A new design of the first-stage nozzle vane for the turbine of a prospective gas-turbine engine is considered. The blade's thermal state is numerically simulated in conjugate statement using the ANSYS CFX 13.0 software package. Critical locations in the blade design are determined from the distribution of heat fluxes, and measures aimed at achieving more efficient cooling are analyzed. Essentially lower (by 50-100°C) maximal temperature of metal has been achieved owing to the results of the performed work.

11. Numerical computation of homogeneous slope stability.

PubMed

Xiao, Shuangshuang; Li, Kemin; Ding, Xiaohua; Liu, Tong

2015-01-01

To simplify the computational process of homogeneous slope stability, improve computational accuracy, and find multiple potential slip surfaces of a complex geometric slope, this study utilized the limit equilibrium method to derive expression equations of overall and partial factors of safety. This study transformed the solution of the minimum factor of safety (FOS) to solving of a constrained nonlinear programming problem and applied an exhaustive method (EM) and particle swarm optimization algorithm (PSO) to this problem. In simple slope examples, the computational results using an EM and PSO were close to those obtained using other methods. Compared to the EM, the PSO had a small computation error and a significantly shorter computation time. As a result, the PSO could precisely calculate the slope FOS with high efficiency. The example of the multistage slope analysis indicated that this slope had two potential slip surfaces. The factors of safety were 1.1182 and 1.1560, respectively. The differences between these and the minimum FOS (1.0759) were small, but the positions of the slip surfaces were completely different than the critical slip surface (CSS). PMID:25784927

12. Numerical stability in multifluid gas dynamics with implicit drag forces

Ramshaw, J. D.; Chang, C. H.

2015-10-01

The numerical stability of a conventional explicit numerical scheme for solving the inviscid multifluid dynamical equations describing a multicomponent gas mixture is investigated both analytically and computationally. Although these equations do not explicitly contain diffusion terms, it is well known that they reduce to a single-fluid diffusional description when the drag coefficients in the species momentum equations are large. The question then arises as to whether their numerical solution is subject to a diffusional stability restriction on the time step in addition to the usual Courant sound-speed stability condition. An analytical stability analysis is performed for the special case of a quiescent binary gas mixture with equal sound speeds and temperatures. It is found that the Courant condition is always sufficient to ensure stability, so that no additional diffusional stability restriction arises for any value of the drag coefficient, however large. This result is confirmed by one-dimensional computational results for binary and ternary mixtures with unequal sound speeds, which remain stable even when the time step exceeds the usual diffusional limit by factors of order 100.

13. Verbal, Numerical and Perceptual Skills Related to Chemistry Achievement.

ERIC Educational Resources Information Center

Bodner, George M.; And Others

The relationship between students' relative ability in visual-spatial tasks as well as their verbal and numerical skills to their performance in an introductory college chemistry course was investigated. For 700 subjects, verbal and mathematics Scholastic Aptitude Test scores (SAT-V) and (SAT-M) and the following four perceptual tests were…

14. Busted Butte: Achieving the Objectives and Numerical Modeling Results

SciTech Connect

W.E. Soll; M. Kearney; P. Stauffer; P. Tseng; H.J. Turin; Z. Lu

2002-10-07

The Unsaturated Zone Transport Test (UZTT) at Busted Butte is a mesoscale field/laboratory/modeling investigation designed to address uncertainties associated with flow and transport in the UZ site-process models for Yucca Mountain. The UZTT test facility is located approximately 8 km southeast of the potential Yucca Mountain repository area. The UZTT was designed in two phases, to address five specific objectives in the UZ: the effect of heterogeneities, flow and transport (F&T) behavior at permeability contrast boundaries, migration of colloids , transport models of sorbing tracers, and scaling issues in moving from laboratory scale to field scale. Phase 1A was designed to assess the influence of permeability contrast boundaries in the hydrologic Calico Hills. Visualization of fluorescein movement , mineback rock analyses, and comparison with numerical models demonstrated that F&T are capillary dominated with permeability contrast boundaries distorting the capillary flow. Phase 1B was designed to assess the influence of fractures on F&T and colloid movement. The injector in Phase 1B was located at a fracture, while the collector, 30 cm below, was placed at what was assumed to be the same fracture. Numerical simulations of nonreactive (Br) and reactive (Li) tracers show the experimental data are best explained by a combination of molecular diffusion and advective flux. For Phase 2, a numerical model with homogeneous unit descriptions was able to qualitatively capture the general characteristics of the system. Numerical simulations and field observations revealed a capillary dominated flow field. Although the tracers showed heterogeneity in the test block, simulation using heterogeneous fields did not significantly improve the data fit over homogeneous field simulations. In terms of scaling, simulations of field tracer data indicate a hydraulic conductivity two orders of magnitude higher than measured in the laboratory. Simulations of Li, a weakly sorbing tracer

15. Association between Basic Numerical Abilities and Mathematics Achievement

ERIC Educational Resources Information Center

Sasanguie, Delphine; De Smedt, Bert; Defever, Emmy; Reynvoet, Bert

2012-01-01

Various measures have been used to investigate number processing in children, including a number comparison or a number line estimation task. The present study aimed to examine whether and to which extent these different measures of number representation are related to performance on a curriculum-based standardized mathematics achievement test in…

16. Determining the Numerical Stability of Quantum Chemistry Algorithms.

PubMed

Knizia, Gerald; Li, Wenbin; Simon, Sven; Werner, Hans-Joachim

2011-08-01

We present a simple, broadly applicable method for determining the numerical properties of quantum chemistry algorithms. The method deliberately introduces random numerical noise into computations, which is of the same order of magnitude as the floating point precision. Accordingly, repeated runs of an algorithm give slightly different results, which can be analyzed statistically to obtain precise estimates of its numerical stability. This noise is produced by automatic code injection into regular compiler output, so that no substantial programming effort is required, only a recompilation of the affected program sections. The method is applied to investigate: (i) the numerical stability of the three-center Obara-Saika integral evaluation scheme for high angular momenta, (ii) if coupled cluster perturbative triples can be evaluated with single precision arithmetic, (iii) how to implement the density fitting approximation in Møller-Plesset perturbation theory (MP2) most accurately, and (iv) which parts of density fitted MP2 can be safely evaluated with single precision arithmetic. In the integral case, we find a numerical instability in an equation that is used in almost all integral programs. Due to the results of (ii) and (iv), we conjecture that single precision arithmetic can be applied whenever a calculation is done in an orthogonal basis set and excessively long linear sums are avoided. PMID:26606614

17. 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.

18. Numerical experiments on the stability of preplanetary disks

NASA Technical Reports Server (NTRS)

Cassen, P. M.; Smith, B. F.; Reynolds, R. T.; Miller, R. H.

1981-01-01

Gravitational stability of gaseous protostellar disks is relevant to theories of planetary formation. Stable gas disks favor formation of planetesimals by the accumulation of solid material; unstable disks allow the possibility of direct condensation of gaseous protoplanets. This paper presents the results of numerical experiments designed to test the stability of thin disks against large-scale, self-gravitational disruption. It is found that a disk as massive as 1 solar mass, surrounding a 1 solar mass protostar, can be stable against long-wavelength gravitational disruption if its temperature is about 300 K or greater. Stability of a cooler disk requires that it be less massive, but even at 100 K a stable disk can have an appreciable fraction (about 1/3) of a solar mass.

19. Mapping Numerical Magnitudes onto Symbols: The Numerical Distance Effect and Individual Differences in Children's Mathematics Achievement

ERIC Educational Resources Information Center

Holloway, Ian D.; Ansari, Daniel

2009-01-01

Although it is often assumed that abilities that reflect basic numerical understanding, such as numerical comparison, are related to children's mathematical abilities, this relationship has not been tested rigorously. In addition, the extent to which symbolic and nonsymbolic number processing play differential roles in this relationship is not yet…

20. Stability of Numerical Interface Conditions for Fluid/Structure Interaction

SciTech Connect

Banks, J W; Sjogreen, B

2009-08-13

In multi physics computations, where a compressible fluid is coupled with a linearly elastic solid, it is standard to enforce continuity of the normal velocities and of the normal stresses at the interface between the fluid and the solid. In a numerical scheme, there are many ways that the velocity- and stress-continuity can be enforced in the discrete approximation. This paper performs a normal mode analysis to investigate the stability of different numerical interface conditions for a model problem approximated by upwind type of finite difference schemes. The analysis shows that depending on the ratio of densities between the solid and the fluid, some numerical interface conditions are stable up to the maximal CFL-limit, while other numerical interface conditions suffer from a severe reduction of the stable CFL-limit. The paper also presents a new interface condition, obtained as a simplified charcteristic boundary condition, that is proved to not suffer from any reduction of the stable CFL-limit. Numerical experiments in one space dimension show that the new interface condition is stable also for computations with the non-linear Euler equations of compressible fluid flow coupled with a linearly elastic solid.

1. Stability of Inviscid Flow over Airfoils Admitting Multiple Numerical Solutions

Liu, Ya; Xiong, Juntao; Liu, Feng; Luo, Shijun

2012-11-01

Multiple numerical solutions at the same flight condition are found of inviscid transonic flow over certain airfoils (Jameson et al., AIAA 2011-3509) within some Mach number range. Both symmetric and asymmetric solutions exist for a symmetric airfoil at zero angle of attack. Global linear stability analysis of the multiple solutions is conducted. Linear perturbation equations of the Euler equations around a steady-state solution are formed and discretized numerically. An eigenvalue problem is then constructed using the modal analysis approach. Only a small portion of the eigen spectrum is needed and thus can be found efficiently by using Arnoldi's algorithm. The least stable or unstable mode corresponds to the eigenvalue with the largest real part. Analysis of the NACA 0012 airfoil indicates stability of symmetric solutions of the Euler equations at conditions where buffet is found from unsteady Navier-Stokes equations. Euler solutions of the same airfoil but modified to include the displacement thickness of the boundary layer computed from the Navier-Stokes equations, however, exhibit instability based on the present linear stability analysis. Graduate Student.

2. Numerical studies of transverse curvature effects on transonic flow stability

NASA Technical Reports Server (NTRS)

Macaraeg, M. G.; Daudpota, Q. I.

1992-01-01

A numerical study of transverse curvature effects on compressible flow temporal stability for transonic to low supersonic Mach numbers is presented for axisymmetric modes. The mean flows studied include a similar boundary-layer profile and a nonsimilar axisymmetric boundary-layer solution. The effect of neglecting curvature in the mean flow produces only small quantitative changes in the disturbance growth rate. For transonic Mach numbers (1-1.4) and aerodynamically relevant Reynolds numbers (5000-10,000 based on displacement thickness), the maximum growth rate is found to increase with curvature - the maximum occurring at a nondimensional radius (based on displacement thickness) between 30 and 100.

3. Stability constraints for oceanic numerical models: implications for the formulation of time and space discretizations

Lemarié, F.; Debreu, L.; Madec, G.; Demange, J.; Molines, J. M.; Honnorat, M.

2015-08-01

Except for vertical diffusion (and possibly the external mode and bottom drag), oceanic models usually rely on explicit time-stepping algorithms subject to Courant-Friedrichs-Lewy (CFL) stability criteria. Implicit methods could be unconditionally stable, but an algebraic system must be solved at each time step and other considerations such as accuracy and efficiency are less straightforward to achieve. Depending on the target application, the process limiting the maximum allowed time-step is generally different. In this paper, we introduce offline diagnostics to predict stability limits associated with internal gravity waves, advection, diffusion, and rotation. This suite of diagnostics is applied to a set of global, regional and coastal numerical simulations with several horizontal/vertical resolutions and different numerical models. We show that, for resolutions finer that 1/2°, models with an Eulerian vertical coordinate are generally constrained by vertical advection in a few hot spots and that numerics must be extremely robust to changes in Courant number. Based on those results, we review the stability and accuracy of existing numerical kernels in vogue in primitive equations oceanic models with a focus on advective processes and the dynamics of internal waves. We emphasize the additional value of studying the numerical kernel of oceanic models in the light of coupled space-time approaches instead of studying the time schemes independently from spatial discretizations. From this study, we suggest some guidelines for the development of temporal schemes in future generation multi-purpose oceanic models.

4. Achieving Stability Requirements for Nanoprobe and Long Beam Lines at NSLS II. A Comprehensive Study

SciTech Connect

Simos,N.; Fallier, M.; Hill, J.; Berman, L.; Evans-Lutterodt, K.; Broadbent, A.

2008-06-23

Driven by beam stability requirements at the NSLS II synchrotron, such that the desired small beam sizes and high brightness are both realized and stable, a comprehensive study has been launched seeking to provide assurances that stability at the nanometer level at critical x-ray beam-lines, is achievable, given the natural and cultural vibration environment at the selected site. The study consists of (a) an extensive investigation of the site to evaluate the existing ground vibration, in terms of amplitude, frequency content and coherence, and (b) of a numerical study of wave propagation and interaction with the infrastructure of the sensitive lines. The paper presents results from both aspects of the study.

5. Stability of vegetated slopes in unsaturated conditions: a numerical study

Battista Chirico, Giovanni; Borga, Marco; Tarolli, Paolo; Rigon, Riccardo; Preti, Federico

2014-05-01

Extreme rainfall events can trigger shallow landslides with failure planes located in soils far from saturated conditions. The stability of shallow soils on very steep slopes under unsaturated conditions can be highly influenced by the vegetation, according to both geo-mechanical and soil-hydrological factors, particularly in regions characterized by a strong climatic seasonality. The root structure of the vegetation reinforces the shallow soils, by providing additional apparent cohesion to the soil. The root water uptake enhances the stability by increasing the frequency of high suction pressure heads in the soil layers explored by the roots. In water controlled eco-systems, such as Mediterranean areas, these two factors are mutually related. Plants develop their root structure in order to optimize the uptake of the water available in the soil, since water availability is limited during the growing season. In this study we present the results of some numerical experiments with the aim to assess the relative importance of these two factors. We simulated the soil water dynamics within homogeneous loamy-sand soils, assuming climatic conditions and root structures typically observed in a deciduous forest of central and southern Italy. An infinite slope stability model is employed for assessing the temporal evolution of the contribute of the soil suction regime to the slope stability, as compared with the contribute of the soil root reinforcement. The results suggest that, during the wet season, the effect of the soil suction state on slope stability is much smaller than that attributable to the mechanical reinforcement provided by the root structure, at least within soil depths explored by the plant roots. Instead, during the growing and dry summer seasons, the soil suction state is far more relevant than the mechanical reinforcement. Thus, accounting for the antecedent soil suction state can be relevant for an appropriate prediction of shallow landslide hazards in

6. a Numerical Method for Stability Analysis of Pinned Flexible Mechanisms

Beale, D. G.; Lee, S. W.

1996-05-01

A technique is presented to investigate the stability of mechanisms with pin-jointed flexible members. The method relies on a special floating frame from which elastic link co-ordinates are defined. Energies are easily developed for use in a Lagrange equation formulation, leading to a set of non-linear and mixed ordinary differential-algebraic equations of motion with constraints. Stability and bifurcation analysis is handled using a numerical procedure (generalized co-ordinate partitioning) that avoids the tedious and difficult task of analytically reducing the system of equations to a number equalling the system degrees of freedom. The proposed method was then applied to (1) a slider-crank mechanism with a flexible connecting rod and crank of constant rotational speed, and (2) a four-bar linkage with a flexible coupler with a constant speed crank. In both cases, a single pinned-pinned beam bending mode is employed to develop resonance curves and stability boundaries in the crank length-crank speed parameter plane. Flip and fold bifurcations are common occurrences in both mechanisms. The accuracy of the proposed method was also verified by comparison with previous experimental results [1].

7. The numerical stroop effect in primary school children: a comparison of low, normal, and high achievers.

PubMed

Heine, Angela; Tamm, Sascha; De Smedt, Bert; Schneider, Michael; Thaler, Verena; Torbeyns, Joke; Stern, Elsbeth; Verschaffel, Lieven; Jacobs, Arthur

2010-01-01

Sixty-six primary school children were selected, of which 21 scored low on a standardized math achievement test, 23 were normal, and 22 high achievers. In a numerical Stroop experiment, children were asked to make numerical and physical size comparisons on digit pairs. The effects of congruity and numerical distance were determined. All children exhibited congruity and distance effects in the numerical comparison. In the physical comparison, children of all performance groups showed Stroop effects when the numerical distance between the digits was large but failed to show them when the distance was small. Numerical distance effects depended on the congruity condition, with a typical effect of distance in the congruent, and a reversed distance effect in the incongruent condition. Our results are hard to reconcile with theories that suggest that deficits in the automaticity of numerical processing can be related to differential math achievement levels. Immaturity in the precision of mappings between numbers and their numerical magnitudes might be better suited to explain the Stroop effects in children. However, as the results for the high achievers demonstrate, in addition to numerical processing capacity per se, domain-general functions might play a crucial role in Stroop performance, too. PMID:20437281

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

ERIC Educational Resources Information Center

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…

9. 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…

10. Stability analysis and numerical simulation of simplified solid rocket motors

Boyer, G.; Casalis, G.; Estivalèzes, J.-L.

2013-08-01

This paper investigates the Parietal Vortex Shedding (PVS) instability that significantly influences the Pressure Oscillations of the long and segmented solid rocket motors. The eigenmodes resulting from the stability analysis of a simplified configuration, namely, a cylindrical duct with sidewall injection, are presented. They are computed taking into account the presence of a wall injection defect, which is shown to induce hydrodynamic instabilities at discrete frequencies. These instabilities exhibit eigenfunctions in good agreement with the measured PVS vortical structures. They are successfully compared in terms of temporal evolution and frequencies to the unsteady hydrodynamic fluctuations computed by numerical simulations. In addition, this study has shown that the hydrodynamic instabilities associated with the PVS are the driving force of the flow dynamics, since they are responsible for the emergence of pressure waves propagating at the same frequency.

11. Liquid falling films: linear stability and direct numerical simulation

Schmidt, Patrick; O'Naraigh, Lennon; Valluri, Prashant; Lucquiaud, Mathieu

2013-11-01

Interfacial instability of falling liquid films in counter-current contact with a turbulent gas phase is investigated by means of an Orr-Sommerfeld analysis. This study is complemented by a full energy budget analysis, identifying the key mechanisms of the instability. This gives first insight into the dynamic behaviour of the two-phase system, which is relevant for a wide range of technical applications, such as absorption and distillation. The linear stability analysis is also used to identify the operating limits of a counter-current operation i.e. the so-called loading and flooding limits. In addition, the results of this analysis are benchmark for direct numerical simulations using the newly launched Two-Phase Level Set (http://sourceforge.net/projects/tpls/) solver. High resolution DNS is used to obtain detailed knowledge of important mechanisms at play, especially with regard to interfacial instability and transient system behaviour, which can help to design more efficient mass transfer equipment such as structured packings. Sulzer Chemtech Ltd, EPSRC, Energy Technology Partnership.

12. The Influence of Pre-University Students' Mathematics Test Anxiety and Numerical Anxiety on Mathematics Achievement

ERIC Educational Resources Information Center

Seng, Ernest Lim Kok

2015-01-01

This study examines the relationship between mathematics test anxiety and numerical anxiety on students' mathematics achievement. 140 pre-university students who studied at one of the institutes of higher learning were being investigated. Gender issue pertaining to mathematics anxieties was being addressed besides investigating the magnitude of…

13. Numerical implications of stabilization by the use of integrals

NASA Technical Reports Server (NTRS)

Beaudet, P. R.

1975-01-01

Liapunov or energy restraint methods for dynamic stabilization in two body motion perturbation problems are considered. Results of computerized orbital stabilization estimates show that the application of energy restraint prevents the occurrence of consistent timing errors in the stepwise integration of equations of motion for a nearly circular orbit.

14. Recommendations for Achieving Accurate Numerical Simulation of Tip Clearance Flows in Transonic Compressor Rotors

NASA Technical Reports Server (NTRS)

VanZante, Dale E.; Strazisar, Anthony J.; Wood, Jerry R,; Hathaway, Michael D.; Okiishi, Theodore H.

2000-01-01

The tip clearance flows of transonic compressor rotors are important because they have a significant impact on rotor and stage performance. While numerical simulations of these flows are quite sophisticated. they are seldom verified through rigorous comparisons of numerical and measured data because these kinds of measurements are rare in the detail necessary to be useful in high-speed machines. In this paper we compare measured tip clearance flow details (e.g. trajectory and radial extent) with corresponding data obtained from a numerical simulation. Recommendations for achieving accurate numerical simulation of tip clearance flows are presented based on this comparison. Laser Doppler Velocimeter (LDV) measurements acquired in a transonic compressor rotor, NASA Rotor 35, are used. The tip clearance flow field of this transonic rotor was simulated using a Navier-Stokes turbomachinery solver that incorporates an advanced k-epsilon turbulence model derived for flows that are not in local equilibrium. Comparison between measured and simulated results indicates that simulation accuracy is primarily dependent upon the ability of the numerical code to resolve important details of a wall-bounded shear layer formed by the relative motion between the over-tip leakage flow and the shroud wall. A simple method is presented for determining the strength of this shear layer.

15. Advanced material modelling in numerical simulation of primary acetabular press-fit cup stability.

PubMed

Souffrant, R; Zietz, C; Fritsche, A; Kluess, D; Mittelmeier, W; Bader, R

2012-01-01

Primary stability of artificial acetabular cups, used for total hip arthroplasty, is required for the subsequent osteointegration and good long-term clinical results of the implant. Although closed-cell polymer foams represent an adequate bone substitute in experimental studies investigating primary stability, correct numerical modelling of this material depends on the parameter selection. Material parameters necessary for crushable foam plasticity behaviour were originated from numerical simulations matched with experimental tests of the polymethacrylimide raw material. Experimental primary stability tests of acetabular press-fit cups consisting of static shell assembly with consecutively pull-out and lever-out testing were subsequently simulated using finite element analysis. Identified and optimised parameters allowed the accurate numerical reproduction of the raw material tests. Correlation between experimental tests and the numerical simulation of primary implant stability depended on the value of interference fit. However, the validated material model provides the opportunity for subsequent parametric numerical studies. PMID:22817471

16. A numerical study of the thermal stability of solar loops

NASA Technical Reports Server (NTRS)

Klimchuk, J. A.; Antiochos, S. K.; Mariska, J. T.

1987-01-01

An important property of all loops is their thermal stability. If low lying hot loops were thermally unstable, for example, a great majority of the low loops on the Sun might be expected to be cool. How small perturbations evolve in low lying, linearly unstable hot loops was determined and how high lying, linearly stable hot loops respond to large amplitude disturbances such as might be expected on the Sun were examined. Only general descriptions and results are given.

17. Numerical stability of pseudo-spectral PIC code generalizations

Godfrey, Brendan B.; Vay, Jean-Luc

2014-10-01

Laser Plasma Accelerator (LPA) particle-in-cell (PIC) simulations are computationally demanding, because they require beam transport over times and distances long compared with the natural scales of the acceleration mechanism and because they are prone to numerical instabilities. To provide greater flexibility in LPA PIC simulations, we have generalized the Pseudo-Spectral Time Domain (PSTD) algorithm to accommodate arbitrary order spatial derivative approximations and substantially longer time steps. Here, we show that, by extending approaches developed by us for other PIC algorithms, numerical Cherenkov instabilities can be suppressed for the generalized PSTD algorithm. We also illustrate the relationships between the generalized PSTD and other PIC algorithms, such as Finite Difference Time Domain (FDTD) and Pseudo-Spectral Analytical Time Domain (PSATD) algorithms. Background information can be found at http://hifweb.lbl.gov/public/BLAST/Godfrey/. Work supported in part by DOE under Contract DE-AC02-05CH11231.

18. Stabilization of numerical interchange in spectral-element magnetohydrodynamics

Sovinec, C. R.

2016-08-01

Auxiliary numerical projections of the divergence of flow velocity and vorticity parallel to magnetic field are developed and tested for the purpose of suppressing unphysical interchange instability in magnetohydrodynamic simulations. The numerical instability arises with equal-order C0 finite- and spectral-element expansions of the flow velocity, magnetic field, and pressure and is sensitive to behavior at the limit of resolution. The auxiliary projections are motivated by physical field-line bending, and coercive responses to the projections are added to the flow-velocity equation. Their incomplete expansions are limited to the highest-order orthogonal polynomial in at least one coordinate of the spectral elements. Cylindrical eigenmode computations show that the projections induce convergence from the stable side with first-order ideal-MHD equations during h-refinement and p-refinement. Hyperbolic and parabolic projections and responses are compared, together with different methods for avoiding magnetic divergence error. The projections are also shown to be effective in linear and nonlinear time-dependent computations with the NIMROD code Sovinec et al. [17], provided that the projections introduce numerical dissipation.

19. Numerical simulation of fluid-structure interactions with stabilized finite element method

Sváček, Petr

2016-03-01

This paper is interested to the interactions of the incompressible flow with a flexibly supported airfoil. The bending and the torsion modes are considered. The problem is mathematically described. The numerical method is based on the finite element method. A combination of the streamline-upwind/Petrov-Galerkin and pressure stabilizing/Petrov-Galerkin method is used for the stabilization of the finite element method. The numerical results for a three-dimensional problem of flow over an airfoil are shown.

20. A numerical model for stability considerations in HTS magnets

Lehtonen, Jorma; Mikkonen, Risto; Paasi, Jaakko

2000-03-01

We propose that in an HTS application, stability is lost more likely because of a global increase in temperature caused by heat generation distributed over the whole coil than because of a local normal zone which starts to propagate. For consideration of stability in HTS magnets, we present a computational model based on the heat conduction equation coupled with Maxwell's equations, whereby analysis can be performed by using commercial software packages for computational electromagnetics and thermodynamics. For temperature distribution inside the magnet, we derive the magnetic field dependent effective values of thermal conductivity, specific heat, and heat generated by electromagnetic phenomena for the composite structure of the magnet, while cooling conditions and external heat sources are described as boundary conditions. Our model enables the magnet designer to estimate a safe level of the operation current before a thermal runaway. Finally, as examples, we present some calculations of the HTS magnet with ac to review the effects of slanted electric field-current density E (J ) characteristics and high critical temperature of HTS materials.

1. Mechanical stability of propped hydraulic fractures: A numerical study

SciTech Connect

Asgian, M.I.; Cundall, P.A.; Brady, B.H.

1995-03-01

Proppant is sometimes produced along with hydrocarbons in hydraulically fractured petroleum wells. Sometimes 10% to 20% of the proppant is backproduced, which can lead to damaged equipment and downtime. Furthermore, proppant flowback can lead to a substantial loss of fracture conductivity. A numerical study was conducted to help understand what conditions are likely to lead to proppant flowback. In the simulations, the mechanical interaction of a larger number (several thousand) individual proppant grains was modeled with a distinct-element-type code. The numerical simulations show that hydraulic fractures propped with cohesionless, unbonded proppant fail under closure stress at a critical ratio of mean grain diameter to fracture width. This is consistent with published laboratory studies. The simulations identify the mechanism (arch failure) that triggers the mechanical instability and also show that the primary way that drawdowns (less than {approx} 75 psi/ft) affect proppant flowback is to transport loose proppant grains in front of the stable arch to the wellbore. Drawdowns > 75 psi/ft are sufficient to destabilize the arch and to cause progressive failure of the propped fractures.

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

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.

3. Accurate numerical solution of compressible, linear stability equations

NASA Technical Reports Server (NTRS)

Malik, M. R.; Chuang, S.; Hussaini, M. Y.

1982-01-01

The present investigation is concerned with a fourth order accurate finite difference method and its application to the study of the temporal and spatial stability of the three-dimensional compressible boundary layer flow on a swept wing. This method belongs to the class of compact two-point difference schemes discussed by White (1974) and Keller (1974). The method was apparently first used for solving the two-dimensional boundary layer equations. Attention is given to the governing equations, the solution technique, and the search for eigenvalues. A general purpose subroutine is employed for solving a block tridiagonal system of equations. The computer time can be reduced significantly by exploiting the special structure of two matrices.

4. Stability and Change in Achievement Goal Orientations: A Person-Centered Approach

ERIC Educational Resources Information Center

Tuominen-Soini, Heta; Salmela-Aro, Katariina; Niemivirta, Markku

2011-01-01

Two longitudinal studies used a person-centered approach to examine the stability and change in students' achievement goal orientations within a school year (i.e., during 9th grade; measurement period 4 months, N = 530) and between school years (i.e., across 11th and 12th grade; measurement period 12 months, N = 519). Distinct groups of students…

5. Numerical Stability of Partitioned Approach in Fluid-Structure Interaction for a Deformable Thin-Walled Vessel

PubMed Central

Wong, Kelvin K. L.; Thavornpattanapong, Pongpat; Cheung, Sherman C. P.; Tu, Jiyuan

2013-01-01

Added-mass instability is known to be an important issue in the partitioned approach for fluid-structure interaction (FSI) solvers. Despite the implementation of the implicit approach, convergence of solution can be difficult to achieve. Relaxation may be applied to improve this implicitness of the partitioned algorithm, but this commonly leads to a significant increase in computational time. This is because the critical relaxation factor that allows stability of the coupling tends to be impractically small. In this study, a mathematical analysis for optimizing numerical performance based on different time integration schemes that pertain to both the fluid and solid accelerations is presented. The aim is to determine the most efficient configuration for the FSI architecture. Both theoretical and numerical results suggest that the choice of time integration schemes has a significant influence on the stability of FSI coupling. This concludes that, in addition to material and its geometric properties, the choice of time integration schemes is important in determining the stability of the numerical computation. A proper selection of the associated parameters can improve performance considerably by influencing the condition of coupling stability. PMID:24222785

6. Stability of numerous novel potassium chlorides at high pressure.

PubMed

Zhang, Weiwei; Oganov, Artem R; Zhu, Qiang; Lobanov, Sergey S; Stavrou, Elissaios; Goncharov, Alexander F

2016-01-01

K-Cl is a simple system displaying all four main types of bonding, as it contains (i) metallic potassium, (ii) elemental chlorine made of covalently bonded Cl2 molecules held together by van der Waals forces, and (iii) an archetypal ionic compound KCl. The charge balance rule, assigning classical charges of "+1" to K and "-1" to Cl, predicts that no compounds other than KCl are possible. However, our quantum-mechanical variable-composition evolutionary simulations predict an extremely complex phase diagram, with new thermodynamically stable compounds K3Cl, K2Cl, K3Cl2, K4Cl3, K5Cl4, K3Cl5, KCl3 and KCl7. Of particular interest are 2D-metallic homologs Kn+1Cln, the presence of positively charged Cl atoms in KCl7, and the predicted stability of KCl3 already at nearly ambient pressures at zero Kelvin. We have synthesized cubic -KCl3 at 40-70 GPa and trigonal -KCl3 at 20-40 GPa in a laser-heated diamond anvil cell (DAC) at temperature exceeding 2000 K from KCl and Cl2. These phases were identified using in situ synchrotron X-ray diffraction and Raman spectroscopy. Upon unloading to 10 GPa, -KCl3 transforms to a yet unknown structure before final decomposition to KCl and Cl2 at near-ambient conditions. PMID:27211847

7. Stability of numerous novel potassium chlorides at high pressure

Zhang, Weiwei; Oganov, Artem R.; Zhu, Qiang; Lobanov, Sergey S.; Stavrou, Elissaios; Goncharov, Alexander F.

2016-05-01

K-Cl is a simple system displaying all four main types of bonding, as it contains (i) metallic potassium, (ii) elemental chlorine made of covalently bonded Cl2 molecules held together by van der Waals forces, and (iii) an archetypal ionic compound KCl. The charge balance rule, assigning classical charges of “+1” to K and “‑1” to Cl, predicts that no compounds other than KCl are possible. However, our quantum-mechanical variable-composition evolutionary simulations predict an extremely complex phase diagram, with new thermodynamically stable compounds K3Cl, K2Cl, K3Cl2, K4Cl3, K5Cl4, K3Cl5, KCl3 and KCl7. Of particular interest are 2D-metallic homologs Kn+1Cln, the presence of positively charged Cl atoms in KCl7, and the predicted stability of KCl3 already at nearly ambient pressures at zero Kelvin. We have synthesized cubic -KCl3 at 40–70 GPa and trigonal -KCl3 at 20–40 GPa in a laser-heated diamond anvil cell (DAC) at temperature exceeding 2000 K from KCl and Cl2. These phases were identified using in situ synchrotron X-ray diffraction and Raman spectroscopy. Upon unloading to 10 GPa, -KCl3 transforms to a yet unknown structure before final decomposition to KCl and Cl2 at near-ambient conditions.

8. Stability of numerous novel potassium chlorides at high pressure

PubMed Central

Zhang, Weiwei; Oganov, Artem R.; Zhu, Qiang; Lobanov, Sergey S.; Stavrou, Elissaios; Goncharov, Alexander F.

2016-01-01

K-Cl is a simple system displaying all four main types of bonding, as it contains (i) metallic potassium, (ii) elemental chlorine made of covalently bonded Cl2 molecules held together by van der Waals forces, and (iii) an archetypal ionic compound KCl. The charge balance rule, assigning classical charges of “+1” to K and “−1” to Cl, predicts that no compounds other than KCl are possible. However, our quantum-mechanical variable-composition evolutionary simulations predict an extremely complex phase diagram, with new thermodynamically stable compounds K3Cl, K2Cl, K3Cl2, K4Cl3, K5Cl4, K3Cl5, KCl3 and KCl7. Of particular interest are 2D-metallic homologs Kn+1Cln, the presence of positively charged Cl atoms in KCl7, and the predicted stability of KCl3 already at nearly ambient pressures at zero Kelvin. We have synthesized cubic -KCl3 at 40–70 GPa and trigonal -KCl3 at 20–40 GPa in a laser-heated diamond anvil cell (DAC) at temperature exceeding 2000 K from KCl and Cl2. These phases were identified using in situ synchrotron X-ray diffraction and Raman spectroscopy. Upon unloading to 10 GPa, -KCl3 transforms to a yet unknown structure before final decomposition to KCl and Cl2 at near-ambient conditions. PMID:27211847

9. Effect of virtual mass on the characteristics and the numerical stability in two-phase flows

SciTech Connect

No, H.C.; Kazimi, M.S.

1981-04-01

It is known that the typical six equation two-fluid model of the two-phase flow possesses complex characteristics, exhibits unbounded instabilities in the short-wavelength limit and constitutes an ill-posed initial value problem. Among the suggestions to overcome these difficulties, one model for the virtual mass force terms were studied here, because the virtual mass represents real physical effects to accomplish the dissipation for numerical stability. It was found that the virtual mass has a profound effect upon the mathematical characteristic and numerical stability. Here a quantitative bound on the coefficient of the virtual mass terms was suggested for mathematical hyperbolicity and numerical stability. It was concluded that the finite difference scheme with the virtual mass model is restricted only by the convective stability conditions with the above suggested value.

10. 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. PMID:25080896

11. Stability of numerous novel potassium chlorides at high pressure

DOE PAGESBeta

Zhang, Weiwei; Oganov, Artem R.; Lobanov, Sergey S.; Stavrou, Elissaios; Goncharov, Alexander F.; Zhu, Qiang

2016-05-23

K-Cl is a simple system displaying all four main types of bonding, as it contains (i) metallic potassium, (ii) elemental chlorine made of covalently bonded Cl2 molecules held together by van der Waals forces, and (iii) an archetypal ionic compound KCl. The charge balance rule, assigning classical charges of “+1” to K and “–1” to Cl, predicts that no compounds other than KCl are possible. However, our quantum-mechanical variable-composition evolutionary simulations predict an extremely complex phase diagram, with new thermodynamically stable compounds K3Cl, K2Cl, K3Cl2, K4Cl3, K5Cl4, K3Cl5, KCl3 and KCl7. Of particular interest are 2D-metallic homologs Kn+1Cln, the presencemore » of positively charged Cl atoms in KCl7, and the predicted stability of KCl3 already at nearly ambient pressures at zero Kelvin. We have synthesized cubic Pm3¯n -KCl3 at 40–70 GPa and trigonal P3¯cl -KCl3 at 20–40 GPa in a laser-heated diamond anvil cell (DAC) at temperature exceeding 2000 K from KCl and Cl2. These phases were identified using in situ synchrotron X-ray diffraction and Raman spectroscopy. Lastly, upon unloading to 10 GPa, P3¯cl -KCl3 transforms to a yet unknown structure before final decomposition to KCl and Cl2 at near-ambient conditions.« less

12. Numerical simulations and stability of magnetic structures in the heliosheath

Cox, S.; Avinash, K.; Shaikh, D.; Zank, G. P.

2008-12-01

We extend the three fluid model of Avinash and Zank [2007] for magnetic structures in the heliosheath to a four fluid model consisting of electrons, pick-up ions (PUIs), solar wind ions (SWI), and neutral hydrogen. The PUIs are generated by neutrals via charge exchange with SWI. Since the kinetic pressure of PUIs is nearly three to four times the pressure of SWI, these are more suited to mediate small scale structures in the heliosheath such as magnetic holes/humps etc. The constant energy exchange between these two fluids drives them non-adiabatic. The PUIs are isothermal while SWI are non adiabatic with an index ~1.25. The four fluid model captures these effects via a modified equation of state for PUI and SWI. The phase space of time independent solutions in terms of the Mach numbers of PUI and SWI is constructed to delineate the parameter space which allows structure formation in the heliosheath. We examine the stability of the time independent solutions by evolving them via a full system of Hall -MHD equations. The simulation results show that these solutions are not quite stable. As the structure propagates it develops growing oscillations in the wings. Concomitantly, there are changes in the amplitude and width of the structure. This instability could be due to local changes in the velocity of the structure and reflects an exchange between the kinetic and magnetic parts of the total energy. Our results showing the presence of growing oscillations in the wings of solitary wave solutions are consistent with the recent analysis of magnetic holes in the heliosheth by Burlaga et al [2007].

13. Ponderomotive stabilization of flute modes in mirrors Feedback control and numerical results

NASA Technical Reports Server (NTRS)

Similon, P. L.

1987-01-01

Ponderomotive stabilization of rigid plasma flute modes is numerically investigated by use of a variational principle, for a simple geometry, without eikonal approximation. While the near field of the studied antenna can be stabilizing, the far field has a small contribution only, because of large cancellation by quasi mode-coupling terms. The field energy for stabilization is evaluated and is a nonnegligible fraction of the plasma thermal energy. A new antenna design is proposed, and feedback stabilization is investigated. Their use drastically reduces power requirements.

14. Numerical Stability and Accuracy of Temporally Coupled Multi-Physics Modules in Wind-Turbine CAE Tools

SciTech Connect

Gasmi, A.; Sprague, M. A.; Jonkman, J. M.; Jones, W. B.

2013-02-01

In this paper we examine the stability and accuracy of numerical algorithms for coupling time-dependent multi-physics modules relevant to computer-aided engineering (CAE) of wind turbines. This work is motivated by an in-progress major revision of FAST, the National Renewable Energy Laboratory's (NREL's) premier aero-elastic CAE simulation tool. We employ two simple examples as test systems, while algorithm descriptions are kept general. Coupled-system governing equations are framed in monolithic and partitioned representations as differential-algebraic equations. Explicit and implicit loose partition coupling is examined. In explicit coupling, partitions are advanced in time from known information. In implicit coupling, there is dependence on other-partition data at the next time step; coupling is accomplished through a predictor-corrector (PC) approach. Numerical time integration of coupled ordinary-differential equations (ODEs) is accomplished with one of three, fourth-order fixed-time-increment methods: Runge-Kutta (RK), Adams-Bashforth (AB), and Adams-Bashforth-Moulton (ABM). Through numerical experiments it is shown that explicit coupling can be dramatically less stable and less accurate than simulations performed with the monolithic system. However, PC implicit coupling restored stability and fourth-order accuracy for ABM; only second-order accuracy was achieved with RK integration. For systems without constraints, explicit time integration with AB and explicit loose coupling exhibited desired accuracy and stability.

15. Numerical stability of the electromagnetic quasinormal and quasibound modes of Kerr black holes

Staicova, Denitsa; Fiziev, Plamen

2015-07-01

The proper understanding of the electromagnetic counterpart of gravity-waves emitters is of serious interest to the multimessenger astronomy. In this article, we study the numerical stability of the quasinormal modes (QNM) and quasibound modes (QBM) obtained as solutions of the Teukolsky Angular Equation and the Teukolsky Radial Equation with appropriate boundary conditions. We use the epsilon-method for the system featuring the confluent Heun functions to study the stability of the spectra with respect to changes in the radial variable. We find that the QNM and QBM are stable in certain regions of the complex plane, just as expected, while the third ``spurious'' spectrum was found to be numerically unstable and thus unphysical. This analysis shows the importance of understanding the numerical results in the framework of the theory of the confluent Heun functions, in order to be able to distinguish the physical spectra from the numerical artifacts.

16. The Role of a Steepness Parameter in the Exponential Stability of a Model Problem. Numerical Aspects

Todorovic, N.

2011-06-01

The Nekhoroshev theorem considers quasi integrable Hamiltonians providing stability of actions in exponentially long times. One of the hypothesis required by the theorem is a mathematical condition called steepness. Nekhoroshev conjectured that different steepness properties should imply numerically observable differences in the stability times. After a recent study on this problem (Guzzo et al. 2011, Todorovic et al. 2011) we show some additional numerical results on the change of resonances and the diffusion laws produced by the increasing effect of steepness. The experiments are performed on a 4-dimensional steep symplectic map designed in a way that a parameter smoothly regulates the steepness properties in the model.

17. Comparisons of numerical magnitudes in children with different levels of mathematical achievement. An ERP study.

PubMed

Gómez-Velázquez, Fabiola Reveca; Berumen, Gustavo; González-Garrido, Andrés Antonio

2015-11-19

The ability to map between non-symbolic and symbolic magnitude representations is crucial in the development of mathematics and this map is disturbed in children with math difficulties. In addition, positive parietal ERPs have been found to be sensitive to the number distance effect and skills solving arithmetic problems. Therefore we aimed to contrast the behavioral and ERP responses in children with different levels of mathematical achievement: low (LA), average (AA) and high (HA), while comparing symbolic and non-symbolic magnitudes. The results showed that LA children repeatedly failed when comparing magnitudes, particularly the symbolic ones. In addition, a positive correlation between correct responses while analyzing symbolic quantities and WRAT-4 scores emerged. The amplitude of N200 was significantly larger during non-symbolic comparisons. In addition, P2P amplitude was consistently smaller in LA children while comparing both symbolic and non-symbolic quantities, and correlated positively with the WRAT-4 scores. The latency of P3 seemed to be sensitive to the type of numerical comparison. The results suggest that math difficulties might be related to a more general magnitude representation problem, and that ERP are useful to study its timecourse in children with different mathematical skills. PMID:26385418

18. Numerical study on the stabilization of neoclassical tearing modes by electron cyclotron current drive

SciTech Connect

Wang, Xiaoguang; Zhang, Xiaodong; Wu, Bin; Zhu, Sizheng; Hu, Yemin

2015-02-15

It is well known that electron cyclotron current drive (ECCD) around the o-point of magnetic island along the plasma current direction can stabilize neoclassical tearing modes (NTMs) in tokamak devices. The effects of the radial misalignment between the island and the driven current, the phase misalignment, and the on-duty ratio for modulated current drive on NTM stabilization are studied numerically in this paper. A small radial misalignment is found to significantly decrease the stabilizing effect. When a sufficiently large phase misalignment occurs for the modulated ECCD, the stabilization effect is also reduced a lot. The optimal on-duty ratio of modulated ECCD to stabilize NTMs is found to be in the range of 60%–70%. A larger on-duty ratio than 50% could also mitigate the effect of phase misalignment if it is not too large. There is no benefit from modulation if the phase misalignment is larger than a threshold.

19. Numerical study on the stabilization of neoclassical tearing modes by electron cyclotron current drive

Wang, Xiaoguang; Zhang, Xiaodong; Wu, Bin; Zhu, Sizheng; Hu, Yemin

2015-02-01

It is well known that electron cyclotron current drive (ECCD) around the o-point of magnetic island along the plasma current direction can stabilize neoclassical tearing modes (NTMs) in tokamak devices. The effects of the radial misalignment between the island and the driven current, the phase misalignment, and the on-duty ratio for modulated current drive on NTM stabilization are studied numerically in this paper. A small radial misalignment is found to significantly decrease the stabilizing effect. When a sufficiently large phase misalignment occurs for the modulated ECCD, the stabilization effect is also reduced a lot. The optimal on-duty ratio of modulated ECCD to stabilize NTMs is found to be in the range of 60%-70%. A larger on-duty ratio than 50% could also mitigate the effect of phase misalignment if it is not too large. There is no benefit from modulation if the phase misalignment is larger than a threshold.

20. On the stability of numerical integration routines for ordinary differential equations.

NASA Technical Reports Server (NTRS)

Glover, K.; Willems, J. C.

1973-01-01

Numerical integration methods for the solution of initial value problems for ordinary vector differential equations may be modelled as discrete time feedback systems. The stability criteria discovered in modern control theory are applied to these systems and criteria involving the routine, the step size and the differential equation are derived. Linear multistep, Runge-Kutta, and predictor-corrector methods are all investigated.

1. The stability of numerical methods for second order ordinary differential equations

NASA Technical Reports Server (NTRS)

Gear, C. W.

1978-01-01

An important characterization of a numerical method for first order ODE's is the region of absolute stability. If all eigenvalues of the linear problem dy/dt = Ay are inside this region, the numerical method is stable. If the second order system d/dt(dy/dt) = 2Ady/dt - By is solved as a first order system, the same result applies to the eigenvalues of the generalized eigenvalue problem (lambda-squared)I 2(lambda)A + B. No such region exists for general methods for second order equations, but in some cases a region of absolute stability can be defined for methods for the single second order equation d/dt(dy/dt) = 2ady/dt - by. The absence of a region of absolute stability can occur when different members of a system of first order equations are solved by different methods.

2. Stability analysis of numerical boundary conditions and implicit difference approximations for hyperbolic equations

NASA Technical Reports Server (NTRS)

Beam, R. M.; Warming, R. F.; Yee, H. C.

1981-01-01

Implicit, noniterative, finite difference schemes were recently developed by several authors for multidimensional systems of nonlinear hyperbolic partial differential equations. When applied to linear model equations with periodic boundary conditions those schemes are unconditionally stable (A-stable). As applied in practice the algorithms often face a severe time step restriction. A major source of the difficulty is the treatment of the numerical boundary conditions. One conjecture was that unconditional stability requires implicit numerical boundary conditions. An apparent counter example was the space time extrapolation considered by Gustafsson, Kreiss, and Sunstrom. Spatial (implicit) and space time (explicit) extrapolation using normal mode analysis for a finite and infinite number of spatial mesh intervals are examined. The results indicate that for unconditional stability with a finite number of spatial mesh intervals, the numerical boundary conditions must be implicit.

3. The role of opacity and transparency in achieving strategic stability in South Asia.

SciTech Connect

Rajain, Arpit; Ashraf, Tariq Mahmud

2005-08-01

According to international relations theory, deterrence can be used as a tool to achieve stability between potentially hostile nations. India and Pakistan's long history of periodic crises raises the question of how they can achieve deterrence stability. 'Transparency' describes the flow of information between parties and plays a key role in establishing a deterrence relationship. This paper studies the balance needed between opacity and transparency in nuclear topics for the maintenance of deterrence stability between India and Pakistan. States with nuclear weapons are postulated to implement transparency in four categories: potential, capability, intent, and resolve. The study applies these categories to the nuclear components of the ongoing India-Pakistan Composite Dialogue Working Group for Peace and Security including CBMs. To focus our efforts, we defined four scenarios to characterize representative strategic/military/political conditions. The scenarios are combinations of these two sets of opposite poles: competition - cooperation; extremism - moderation (to be understood primarily in a religious/nationalistic sense). We describe each scenario in terms of select focal areas (nuclear doctrine, nuclear command and control, nuclear stockpile, nuclear delivery/defensive systems, and conventional force posture). The scenarios help frame the realm of possibilities, and have been described in terms of expected conditions for the focal areas. We then use the conditions in each scenario to prescribe a range of information-sharing actions that the two countries could take to increase stability. We also highlight the information that should not be shared. These actions can be political (e.g., declarations), procedural (e.g., advance notice of certain military activities), or technologically based (e.g., seismic monitoring of the nuclear test moratorium).

4. Numerical stability analysis of the pseudo-spectral analytical time-domain PIC algorithm

SciTech Connect

Godfrey, Brendan B.; Vay, Jean-Luc; Haber, Irving

2014-02-01

The pseudo-spectral analytical time-domain (PSATD) particle-in-cell (PIC) algorithm solves the vacuum Maxwell's equations exactly, has no Courant time-step limit (as conventionally defined), and offers substantial flexibility in plasma and particle beam simulations. It is, however, not free of the usual numerical instabilities, including the numerical Cherenkov instability, when applied to relativistic beam simulations. This paper derives and solves the numerical dispersion relation for the PSATD algorithm and compares the results with corresponding behavior of the more conventional pseudo-spectral time-domain (PSTD) and finite difference time-domain (FDTD) algorithms. In general, PSATD offers superior stability properties over a reasonable range of time steps. More importantly, one version of the PSATD algorithm, when combined with digital filtering, is almost completely free of the numerical Cherenkov instability for time steps (scaled to the speed of light) comparable to or smaller than the axial cell size.

5. Numerical Study of the Transverse Stability of Compressive and Rarefactive Alfven Solitons

Hamilton, R.; Haneberg, C.

2015-12-01

A numerical study of the stability of DNLS bright and dark solitons subject to oblique perturbations is reported. The DNLS equation is a weakly nonlinear, weakly dispersive and one dimensional limiting form of MHD with the inclusion of Hall dispersion which has been shown to remain valid for plane wave propagation parallel, as well as quasiparallel, to the ambient magnetic field. Related analytic work has dealt with the transverse stability of circularly polarized Alfven waves [E. Mjolhus, T. Hada, J. Plasma Phys., 43, 257 - 268 (1990)] describing stability in relation to the propagation angle of the perturbation and the wave's amplitude and wavenumber. The amplitude and wavenumber relation for transverse stability has a striking similarity to the criterion for modulational instability. A prior analytic work [M. S. Ruderman, Fluid Dyn. 22, 299, (1987)] found the dark soliton to be unstable. Our numerical results are established in the context of these analytic results. Additionally, the transverse stability properties of dark solitons will be addressed as they relate to their role in representing magnetic decreases observed in interplanetary space.

6. Tailings dams stability analysis using numerical modelling of geotechnical and geophysical data

Mihai, S.; Zlagnean, M.; Oancea, I.; Petrescu, A.

2009-04-01

dam's state of safety. This study considered the SSR (Shear Strength Reduction) technique for slope stability numerical modelling. In the SSR finite element technique, elasto-plastic strength is assumed for dam's materials and shear strengths are progressively reduced until collapse occurs. Numerical modelling was performed on the most critical profile choosed through analysis of geophysical and geotechnical informational volume achieved by insitu or in laboratory tests. Finite element analysis were considered in two situations: first, before geophysical investigations and second considering the whole informational of data achieved. Both situations were analysed in static and pseudo-static conditions. The factor of safety before geophysical investigations is high enough to describe a stable state of stability even for the seismic load. The total displacement distributions were modified by the presence of internal erosional element giving a high state of instability, especially for the pseudo-static case. These analysis using the finite element method prove the importance of structural disturbance elements that may occure inside the dam body produced by internal erosional processes.

7. Numerical Simulations of Power Law Heating Functions for Quiescent Loops: Stability and Observables

Martens, P. C.; Winter, H. D.; Munetsi-Mugomba, K.

2007-12-01

We present the numerical simulations of quiescent coronal loops with heating functions that are power law functions of pressure and temperature. These simulations are made using a time-dependent, 1D hydrodynamics code with heating functions that are treated as dynamic variables which are constantly re- evaluated during the loops' lifetimes. These numerical simulations provide a stability test for the analytical solutions formulated by Martens (2007, submitted) for the same heating functions. TRACE and XRT datasets are simulated to determine if present observables can provide adequate information to discriminate between power law heating functions.

8. Numerical solution of wave equations for the stability of the inner cometo-sheath

NASA Technical Reports Server (NTRS)

Srivastava, Krishna M.; Tsurutani, Bruce T.; Goldstein, Bruce E.

1993-01-01

Numerical solution of the MHD wave equations for stability of the cometary sheath determined by the balance between the inward Lorentz body force and the outward ion-neutral drag force is obtained by using a two-point boundary value method. The eigenvalues and the eigenfunctions are obtained numerically by treating the cometary inner sheath as a layer of finite thickness, bounded by the contact surface, i.e., the diamagnetic cavity boundary. The magnetic field structure discovered in the ionosphere of Comets Halley and Giacobini-Zinner is found to be unstable. The effects of finite plasma pressure, dissociative recombination, and mass loading due to photoionization are found to be stabilizing but are unable to quench the instability completely. It is also found that the higher the neutral production rate the lesser is the growth rate for the instability.

9. A numerical model based on closed form solution for elastic stability of thin plates

Ciaramella, S.; Migliore, M.; Minutolo, V.; Ruocco, E.

2010-06-01

An analytical approach for studying the elastic stability of thin rectangular plates under arbitrary boundary conditions is presented. Because the solution is given in closed-form, the approach can be regarded as "exact" under the Kirchhoff-Love assumption. The proposed procedure allows us to obtain the buckling load and modal displacements that do not depend on the number of elements adopted in the numerical discretization using, say, the finite element method.

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

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.

11. Numerical study on the stability of weakly collisional plasma in E×B fields

SciTech Connect

Horký, M.

2015-02-15

Plasma stability in weakly collisional plasmas in the presence of E×B fields is studied with numerical simulations. Different types of ion-neutral collisions are considered in a fully magnetized regime. We study the influence of ion-neutral collisions and the role of collision types on the stability of plasma. It is found that the stability of plasma depends on the type of ion-neutral collisions, with the plasma being unstable for charge exchange collisions, and stable for the elastic scattering. The analysis focuses on the temporal evolution of the velocity phase space, RMS values of the potential fluctuations, and coherent structures in potential densities. For the unstable case, we observe growth and propagation of electrostatic waves. Simulations are performed with a three-dimensional electrostatic particle in cell code.

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

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.

13. Achieving Vibration Stability of the NSLS-II Hard X-ray Nanoprobe Beamline

SciTech Connect

Simos, N.; Chu, Y. N.; Broadbent, A.; Nazaretski, E.; Margulies, L.; Dyling, O.; Shen, Q.; Fallier, M.

2010-08-30

The Hard X-ray Nanoprobe (HXN) Beamline of National Synchrotron Light Source II (NSLS-lI) requires high levels of stability in order to achieve the desired instrument resolution. To ensure that the design of the endstation helps meet the stringent criteria and that natural and cultural vibration is mitigated both passively and actively, a comprehensive study complimentary to the design process has been undertaken. Vibration sources that have the potential to disrupt sensitive experiments such as wind, traffic and NSLS II operating systems have been studied using state of the art simulations and an array of field data. Further, final stage vibration isolation principles have been explored in order to be utilized in supporting endstation instruments. This paper presents results of the various study aspects and their influence on the HXN design optimization.

14. 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…

15. Combined analytical and numerical approaches in Dynamic Stability analyses of engineering systems

Náprstek, Jiří

2015-03-01

Dynamic Stability is a widely studied area that has attracted many researchers from various disciplines. Although Dynamic Stability is usually associated with mechanics, theoretical physics or other natural and technical disciplines, it is also relevant to social, economic, and philosophical areas of our lives. Therefore, it is useful to occasionally highlight the general aspects of this amazing area, to present some relevant examples and to evaluate its position among the various branches of Rational Mechanics. From this perspective, the aim of this study is to present a brief review concerning the Dynamic Stability problem, its basic definitions and principles, important phenomena, research motivations and applications in engineering. The relationships with relevant systems that are prone to stability loss (encountered in other areas such as physics, other natural sciences and engineering) are also noted. The theoretical background, which is applicable to many disciplines, is presented. In this paper, the most frequently used Dynamic Stability analysis methods are presented in relation to individual dynamic systems that are widely discussed in various engineering branches. In particular, the Lyapunov function and exponent procedures, Routh-Hurwitz, Liénard, and other theorems are outlined together with demonstrations. The possibilities for analytical and numerical procedures are mentioned together with possible feedback from experimental research and testing. The strengths and shortcomings of these approaches are evaluated together with examples of their effective complementing of each other. The systems that are widely encountered in engineering are presented in the form of mathematical models. The analyses of their Dynamic Stability and post-critical behaviour are also presented. The stability limits, bifurcation points, quasi-periodic response processes and chaotic regimes are discussed. The limit cycle existence and stability are examined together with their

16. Continuation of equilibria and stability of slender elastic rods using an asymptotic numerical method

Lazarus, A.; Miller, J. T.; Reis, P. M.

2013-08-01

We present a theoretical and numerical framework to compute bifurcations of equilibria and stability of slender elastic rods. The 3D kinematics of the rod is treated in a geometrically exact way by parameterizing the position of the centerline and making use of quaternions to represent the orientation of the material frame. The equilibrium equations and the stability of their solutions are derived from the mechanical energy which takes into account the contributions due to internal moments (bending and twist), external forces and torques. Our use of quaternions allows for the equilibrium equations to be written in a quadratic form and solved efficiently with an asymptotic numerical continuation method. This finite element perturbation method gives interactive access to semi-analytical equilibrium branches, in contrast with the individual solution points obtained from classical minimization or predictor-corrector techniques. By way of example, we apply our numerics to address the specific problem of a naturally curved and heavy rod under extreme twisting and perform a detailed comparison against our own precision model experiments of this system. Excellent quantitative agreement is found between experiments and simulations for the underlying 3D buckling instabilities and the characterization of the resulting complex configurations. We believe that our framework is a powerful alternative to other methods for the computation of nonlinear equilibrium 3D shapes of rods in practical scenarios.

17. Landfast ice affects the stability of the Arctic halocline: Evidence from a numerical model

Itkin, Polona; Losch, Martin; Gerdes, Rüdiger

2015-04-01

Landfast ice covers large surface areas of the winter Siberian Seas. The immobile landfast ice cover inhibits divergent and convergent motion, hence dynamical sea ice growth and redistribution, decouples winter river plumes in coastal seas from the atmosphere, and positions polynyas at the landfast ice edge offshore. In spite of the potentially large effects, state-of-the-art numerical models usually do not represent landfast ice in its correct extent. A simple parametrization of landfast ice based on bathymetry and internal sea ice strength is introduced and its effects on the simulated Arctic Ocean are demonstrated. The simulations suggest that the Siberian landfast ice impacts the Arctic halocline stability through enhanced brine production in polynyas located closer to the shelf break and by redirecting river water to the Canadian Basin. These processes strengthen the halocline in the Canadian Basin, but erode its stability in the Makarov and Eurasian Basins.

18. Three-Dimensional Numerical Analysis for Posture Stability of Laser Propulsion Vehicle

Takahashi, Masayuki; Ohnishi, Naofumi

2011-11-01

We have developed a three-dimensional hydrodynamics code coupling equation of motion of a rigid body for analyzing posture stability of laser propulsion vehicle through numerical simulations of flowfield interacting with unsteady motion of the vehicle. Asymmetric energy distribution is initially added around the focal spot (ring) in order to examine posture stability against an asymmetric blast wave resulting from a laser offset for a lightcraft-type vehicle. The vehicle moves to cancel out the offset from initial offset. However, the Euler angle grows and never returns to zero in a time scale of laser pulse. Also, we found that the vehicle moves to cancel tipping angle when the laser is irradiated to the vehicle with initial tipping angle over the wide angle range, through the vehicle cannot get sufficient restoring force in particular angle, and the tipping angle does not decrease from the initial value for that case.

19. Stability analysis and numerical simulation of 1 prey - 2 predator system

2015-03-01

In this paper, we study an ecological system that consists of 1 prey and 2 predators populations. The prey population grows logistically while Holling type II functional response is applied for both predators . The first predator preys on the prey and the second predator preys on the first one. The study starts with the stability analysis of critical points of the systems. Then, by using normal form and centre manifold method the information about other nontrivial solutions due to bifurcation including possible limit cycles appearance is obtained. The results are confirmed by numerical simulation using MatCont and biological interpretation of the results is also presented.

20. Investigation of Nozzle Stability for the First Ovalization Mode by Numerical Solution of the Fluid Structure Interaction Problem

Schwane, R.; Zia, Y.

2005-02-01

The present paper validates results from numerical simulations for side load generation in rocket nozzles against related data from analytical models that are presently used for rocket engine nozzle design activities. Key words: Nozzle stability; nozzle ovalization; flow structure interaction.

1. Time-resolved photoemission apparatus achieving sub-20-meV energy resolution and high stability

SciTech Connect

Ishida, Y.; Togashi, T.; Yamamoto, K.; Tanaka, M.; Kiss, T.; Otsu, T.; Kobayashi, Y.; Shin, S.

2014-12-15

The paper describes a time- and angle-resolved photoemission apparatus consisting of a hemispherical analyzer and a pulsed laser source. We demonstrate 1.48-eV pump and 5.92-eV probe measurements at the ⩾10.5-meV and ⩾240-fs resolutions by use of fairly monochromatic 170-fs pulses delivered from a regeneratively amplified Ti:sapphire laser system operating typically at 250 kHz. The apparatus is capable to resolve the optically filled superconducting peak in the unoccupied states of a cuprate superconductor, Bi{sub 2}Sr{sub 2}CaCu{sub 2}O{sub 8+δ}. A dataset recorded on Bi(111) surface is also presented. Technical descriptions include the followings: A simple procedure to fine-tune the spatio-temporal overlap of the pump-and-probe beams and their diameters; achieving a long-term stability of the system that enables a normalization-free dataset acquisition; changing the repetition rate by utilizing acoustic optical modulator and frequency-division circuit.

2. A new free-surface stabilization algorithm for geodynamical modelling: Theory and numerical tests

Andrés-Martínez, Miguel; Morgan, Jason P.; Pérez-Gussinyé, Marta; Rüpke, Lars

2015-09-01

The surface of the solid Earth is effectively stress free in its subaerial portions, and hydrostatic beneath the oceans. Unfortunately, this type of boundary condition is difficult to treat computationally, and for computational convenience, numerical models have often used simpler approximations that do not involve a normal stress-loaded, shear-stress free top surface that is free to move. Viscous flow models with a computational free surface typically confront stability problems when the time step is bigger than the viscous relaxation time. The small time step required for stability (< 2 Kyr) makes this type of model computationally intensive, so there remains a need to develop strategies that mitigate the stability problem by making larger (at least ∼10 Kyr) time steps stable and accurate. Here we present a new free-surface stabilization algorithm for finite element codes which solves the stability problem by adding to the Stokes formulation an intrinsic penalization term equivalent to a portion of the future load at the surface nodes. Our algorithm is straightforward to implement and can be used with both Eulerian or Lagrangian grids. It includes α and β parameters to respectively control both the vertical and the horizontal slope-dependent penalization terms, and uses Uzawa-like iterations to solve the resulting system at a cost comparable to a non-stress free surface formulation. Four tests were carried out in order to study the accuracy and the stability of the algorithm: (1) a decaying first-order sinusoidal topography test, (2) a decaying high-order sinusoidal topography test, (3) a Rayleigh-Taylor instability test, and (4) a steep-slope test. For these tests, we investigate which α and β parameters give the best results in terms of both accuracy and stability. We also compare the accuracy and the stability of our algorithm with a similar implicit approach recently developed by Kaus et al. (2010). We find that our algorithm is slightly more accurate

3. Achieving high performance in numerical computations on RISC workstations and parallel systems

SciTech Connect

Goedecker, S.; Hoisie, A.

1997-08-20

The nominal peak speeds of both serial and parallel computers is raising rapidly. At the same time however it is becoming increasingly difficult to get out a significant fraction of this high peak speed from modern computer architectures. In this tutorial the authors give the scientists and engineers involved in numerically demanding calculations and simulations the necessary basic knowledge to write reasonably efficient programs. The basic principles are rather simple and the possible rewards large. Writing a program by taking into account optimization techniques related to the computer architecture can significantly speedup your program, often by factors of 10--100. As such, optimizing a program can for instance be a much better solution than buying a faster computer. If a few basic optimization principles are applied during program development, the additional time needed for obtaining an efficient program is practically negligible. In-depth optimization is usually only needed for a few subroutines or kernels and the effort involved is therefore also acceptable.

4. Diesel flame lift-off stabilization in the presence of laser-ignition: a numerical study

Gong, Cheng; Jangi, Mehdi; Bai, Xue-Song

2015-11-01

Diesel flame lift-off and stabilization in the presence of laser-ignition were numerically investigated with the method of Eulerian stochastic fields. The aim was to scrutinise the interaction between the lifted diesel flame and an ignition kernel upstream of the lifted flame. The numerical simulation was carried out in a constant-volume combustion vessel with n-heptane as fuel. The process was studied previously in an experiment employing Diesel #2 as the fuel in the same combustion vessel. In the experiment a lifted flame was first established at a position downstream of the nozzle. An ignition kernel was then initiated using a high-energy pulse laser at a position upstream of the natural lift-off position of the diesel flame. The laser-ignition kernel was modelled using a high-temperature (∼2000 K) hot spot. In both experiment and simulations the upstream front of the ignition kernel was shown to remain around the initial laser ignition site for a substantially long period of time, while the downstream front of the ignition kernel propagates rapidly towards the natural lift-off position downstream of the laser ignition site. The lift-off position oscillated before the final stabilization at the natural lift-off position. The structures and the propagation speed of the reaction fronts in the laser-ignition kernel and the main flame were analysed. Two different stabilization mechanisms, the auto-ignition mechanism and the flame propagation mechanism, were identified for the naturally lifted flame and the laser-induced reaction front, respectively. A mechanism was proposed to explain the oscillation of the lift-off position.

5. A numerical study of wind turbine wakes under various atmospheric stability conditions

Xie, Shengbai

The goal of this research is to investigate the properties of wind turbine wakes and their interactions with the atmospheric boundary layer (ABL) via large-eddy simulations (LES) with special emphasis on the effects of atmospheric stability. The ABL is considered stable when the ground surface is cooler than the air, unstable when the opposite happens, and neutral when the temperature effect is negligible. In the literature, neutral conditions have been studied extensively, whereas the effects of stability have not. A new LES code, named Wind Turbine and Turbulence Simulator (WiTTS), was developed based on finite-difference (FD) schemes. First, the code's sensitivity to numerous aspects of the FD LES, such as the subgrid-scale (SGS) model, resolution, numerical treatment of the convective term, and filter types, was analyzed by simulating a neutral ABL. It was found that the Lagrangian-averaged scale-dependent (LASD) SGS model performs better than other scale-invariant Smagorinsky-type models. Second, the WiTTS was used to study the wakes from a miniature wind turbine inside a wind tunnel, following the setup of past experimental and numerical studies. It was found that those wakes are spatially anisotropic, with lateral growth faster than the vertical. Based on this, a new wake model is proposed and the Gaussian-type self-similarity is obtained for this simplified scenario. Third, to study a more realistic ABL, the stability conditions have been considered by the Boussinesq approximation and by varying thermal conditions on the ground surface, together with a constant Coriolis force. The LES results indicate that the properties of utility-scale wind turbine wakes are strongly correlated to the stability conditions. The wake recovery is enhanced by the increased turbulence due to buoyant convection in the unstable ABL, while in the stable ABL the spreading of the wake is significantly larger in the lateral direction than in the vertical direction. The stability

6. Stabilized numerical solution for transient dynamic contact of inelastic solids on rough surfaces

Suwannachit, A.; Nackenhorst, U.; Chiarello, R.

2012-06-01

A simulation technique to deal with transient dynamic contact of tire rubber compounds on rough road surfaces is presented. The segment-to-surface approach is used for modeling the contact between tire tread rubber and road track. While the rubber components are deformable and described by a sophisticated viscoelastic damage constitutive model, the road surface is assumed to be rigid and characterized by an analytical function. A spectral approach based on an inverse computation of the 2D-Fast Fourier transform has been suggested for the reconstruction of rough surface profiles. The Newmark time-stepping method is used for the integration of transient dynamic equations. With the so-called contact-stabilized Newmark method the spurious oscillation at contact boundary has been removed. The detailed investigation on the dynamic contact of inelastic rubber block with rough road surfaces has been made. The robustness of the contact-stabilized Newmark method within a finite deformation framework is underlined by numerical studies, in which it is compared with several dissipation-based stabilization techniques selected from literature.

7. von Neumann Stability Analysis of Numerical Solution Schemes for 1D and 2D Euler Equations

Konangi, Santosh; Palakurthi, Nikhil Kumar; Ghia, Urmila

2014-11-01

A von Neumann stability analysis is conducted for numerical schemes for the full system of coupled, density-based 1D and 2D Euler equations, closed by an isentropic equation of state. The governing equations are discretized on a staggered grid, which permits equivalence to finite-volume discretization. Presently, first-order accurate spatial and temporal finite-difference techniques are analyzed. The momentum convection term is treated as explicit, semi-implicit or implicit. Density upwind bias is included in the spatial operator of the continuity equation. By combining the discretization techniques, ten solution schemes are formulated. For each scheme, unstable and stable regimes are identified through the stability analysis, and the maximum allowable CFL number is predicted. The predictions are verified for selected schemes, using the Riemann problem at incompressible and compressible Mach numbers. Very good agreement is obtained between the analytically predicted and ``experimentally'' observed CFL values for all cases, thereby validating the analysis. The demonstrated analysis provides an accurate indication of stability conditions for the Euler equations, in contrast to the simplistic conditions arising from model equations, such as the wave equation.

8. Peaks, plateaus, numerical instabilities, and achievable accuracy in Galerkin and norm minimizing procedures for solving Ax=b

SciTech Connect

Cullum, J.

1994-12-31

Plots of the residual norms generated by Galerkin procedures for solving Ax = b often exhibit strings of irregular peaks. At seemingly erratic stages in the iterations, peaks appear in the residual norm plot, intervals of iterations over which the norms initially increase and then decrease. Plots of the residual norms generated by related norm minimizing procedures often exhibit long plateaus, sequences of iterations over which reductions in the size of the residual norm are unacceptably small. In an earlier paper the author discussed and derived relationships between such peaks and plateaus within corresponding Galerkin/Norm Minimizing pairs of such methods. In this paper, through a set of numerical experiments, the author examines connections between peaks, plateaus, numerical instabilities, and the achievable accuracy for such pairs of iterative methods. Three pairs of methods, GMRES/Arnoldi, QMR/BCG, and two bidiagonalization methods are studied.

9. Numerical solvers to the stabilizing solution of perturbed algebraic Riccati equations in LQ zero-sum games

Ivanov, I. G.; Netov, N. C.; Bogdanova, B. C.

2015-10-01

This paper addresses the problem of solving a generalized algebraic Riccati equation with an indefinite sign of its quadratic term. We extend the approach introduced by Lanzon, Feng, Anderson and Rotkowitz (2008) for solving similar Riccati equations. We numerically investigate two types of iterative methods for computing the stabilizing solution. The first type of iterative methods constructs two matrix sequences, where the sum of them converges to the stabilizing solution. The second type of methods defines one matrix sequence which converges to the stabilizing solution. Computer realizations of the presented methods are numerically tested and compared on the test of family examples. Based on the experiments some conclusions are derived.

10. Numerical Calculations Demonstrating Complete Stabilization of the Ideal Magnetohydrodynamic Resistive Wall Mode by Longitudinal Flow

SciTech Connect

Smith, S.; Jardin, S. C.; Freidberg, J. P.; Guazzotto, L.

2009-05-20

The cylindrical ideal magnetohydrodynamic (MHD) stability problem, including ow and a resistive wall, is cast in the standard mathematical form, ωA∙x = B∙x, without discretizing the vacuum regions surrounding the plasma. This is accomplished by means of a finite element expansion for the plasma perturbations, by coupling the plasma surface perturbations to the resistive wall using a Green's function approach, and by expanding the unknown vector, x, to include the perturbed current in the resistive wall as an additional degree of freedom. The ideal MHD resistive wall mode (RWM) can be stabilized when the plasma has a uniform equilibrium ow such that the RWM frequency resonates with the plasma's Doppler-shifted sound continuum modes. The resonance induces a singularity in the parallel component of the plasma perturbations, which must be adequately resolved. Complete stabilization within the ideal MHD model (i.e. without parallel damping being added) is achieved as the grid spacing in the region of the resonance is extrapolated to 0 step size

11. Comparative Study of Vibration Stability at Operating Light Source Facilities and Lessons Learned in Achieving NSLS II Stability Goals

SciTech Connect

Simos,N.; Fallier, M.; Amick, H.

2008-06-23

In an effort to ensure that the stability goals of the NSLS II will be met once the accelerator structure is set on the selected BNL site a comprehensive evaluation of the ground vibration observed at existing light source facilities has been undertaken. The study has relied on measurement data collected and reported by the operating facilities as well as on new data collected in the course of this study. The primary goal of this comprehensive effort is to compare the green-field conditions that exist in the various sites both in terms of amplitude as well as frequency content and quantify the effect of the interaction of these accelerator facilities with the green-field vibration. The latter represents the ultimate goal of this effort where the anticipated motion of the NSLS II ring is estimated prior to its construction and compared with the required stability criteria.

12. Numerical stability of relativistic beam multidimensional PIC simulations employing the Esirkepov algorithm

SciTech Connect

Godfrey, Brendan B.; Vay, Jean-Luc

2013-09-01

Rapidly growing numerical instabilities routinely occur in multidimensional particle-in-cell computer simulations of plasma-based particle accelerators, astrophysical phenomena, and relativistic charged particle beams. Reducing instability growth to acceptable levels has necessitated higher resolution grids, high-order field solvers, current filtering, etc. except for certain ratios of the time step to the axial cell size, for which numerical growth rates and saturation levels are reduced substantially. This paper derives and solves the cold beam dispersion relation for numerical instabilities in multidimensional, relativistic, electromagnetic particle-in-cell programs employing either the standard or the Cole–Karkkainnen finite difference field solver on a staggered mesh and the common Esirkepov current-gathering algorithm. Good overall agreement is achieved with previously reported results of the WARP code. In particular, the existence of select time steps for which instabilities are minimized is explained. Additionally, an alternative field interpolation algorithm is proposed for which instabilities are almost completely eliminated for a particular time step in ultra-relativistic simulations.

13. Is Cu involved in prion oligopeptide stability? Experiments and numerical simulations

Minicozzi, V.; Morante, S.

The high-sociological impact of neurodegenerative diseases (like Alzheimer disease, Transmissible Spongiform Encephalopathies, Parkinson disease, etc.) has renewed the interest of researchers in the study of misfolding processes and in particular of the rôle played by metals in plaque formation as their unbalanced concentration can be regarded as a possible concurrent cause of protein aggregation. Metals are essential players in many of the fundamental activities of cells. Storing, metabolism, and trafficking of metals through the cellular membrane and within the cytoplasm are mediated by many proteins via well-tuned mechanisms because of the toxicity of free ions. In this review article, we summarize the results of the most recent experimental and numerical investigations aimed at understanding the possible rôle of Cu in stabilizing the Prion protein structure and in the formation of protein polymers.

14. A numerical study of the nonlinear thermal stability of solar loops

NASA Technical Reports Server (NTRS)

Klimchuk, J. A.; Antiochos, S. K.; Mariska, J. T.

1987-01-01

A time-dependent numerical model is used to investigate the nonlinear thermal stability of static loops of various heights. Simulations show that the instability of a hot state with loop heights of less than about 1000 km is physically significant, with an initially hot atmosphere in low-lying compact loops evolving to an extended atmosphere with temperatures far below 100,000 K. Results also show that high-lying loops are stable to all reasonable perturbations, including those of large initial amplitude and long wavelength. The simulation results suggest that low-lying compact loops should not be common to the sun, and that cool loops with temperatures near 100,000 K must be formed in the cool state initially and cannot evolve from preexisiting loops.

15. Experimental and numerical studies of a microfluidic device with compliant chambers for flow stabilization

Iyer, V.; Raj, A.; Annabattula, R. K.; Sen, A. K.

2015-07-01

This paper reports experimental and numerical studies of a passive microfluidic device that stabilizes a pulsating incoming flow and delivers a steady flow at the outlet. The device employs a series of chambers along the flow direction with a thin polymeric membrane (of thickness 75-250 µm) serving as the compliant boundary. The deformation of the membrane allows accumulation of fluid during an overflow and discharge of fluid during an underflow for flow stabilization. Coupled fluid-structure simulations are performed using Mooney-Rivlin formulations to account for a thin hyperelastic membrane material undergoing large deformations to accurately predict the device performance. The device was fabricated with PDMS as the substrate material and thin PDMS membrane as the compliant boundary. The performance of the device is defined in terms of a parameter called ‘Attenuation Factor (AF)’. The effect of various design parameters including membrane thickness, elastic modulus, chamber size and number of chambers in series as well as operating conditions including the outlet pressure, mean input flow rate, fluctuation amplitude and frequency on the device performance were studied using experiments and simulations. The simulation results successfully confront the experimental data (within 10%) which validates the numerical simulations. The device was used at the exit of a PZT actuated valveless micropump to take pulsating flow at the upstream and deliver steady flow downstream. The amplitude of the pulsating flow delivered by the micropump was significantly reduced (AF = 0.05 for a device with three 4 mm chambers) but at the expense of a reduction in the pressure capability (<20%). The proposed device could potentially be used for reducing flow pulsations in practical microfluidic circuits.

16. Linear stability analysis and direct numerical simulation of two layer channel flow

Sahu, Kirti; Govindarajan, Rama; Tripathi, Manoj

2015-11-01

We study the stability of two-fluid flow through a plane channel at Reynolds numbers of a hundred to a thousand. The two fluids have the same density but different viscosities. The fluids, when miscible, are separated from each other by a mixed layer of small but finite thickness, across which viscosity changes from that of one fluid to that of the other. When immiscible, the interface is sharp. Our study spans a range of Schmidt numbers, viscosity ratios and location and thickness of the mixed layer. Our two-dimensional linear stability results predict well the behaviour displayed by our three-dimensional direct numerical simulations at early times. In both linear and non-linear regimes, the miscible flow is more unstable than the corresponding immiscible one, and the miscible flow breaks spanwise symmetry more readily to go into three-dimensionality. We show that the miscible flow over our range of parameters is always significantly more unstable than the corresponding immiscible case.

17. A numerical investigation of coorbital stability and libration in three dimensions

Morais, M. H. M.; Namouni, F.

2016-05-01

Motivated by the dynamics of resonance capture, we study numerically the coorbital resonance for inclination 0≤ I≤ 180° in the circular restricted three-body problem. We examine the similarities and differences between planar and three dimensional coorbital resonance capture and seek their origin in the stability of coorbital motion at arbitrary inclination. After we present stability maps of the planar prograde and retrograde coorbital resonances, we characterize the new coorbital modes in three dimensions. We see that retrograde mode I (R1) and mode II (R2) persist as we change the relative inclination, while retrograde mode III (R3) seems to exist only in the planar problem. A new coorbital mode (R4) appears in 3D which is a retrograde analogue to an horseshoe-orbit. The Kozai-Lidov resonance is active for retrograde orbits as well as prograde orbits and plays a key role in coorbital resonance capture. Stable coorbital modes exist at all inclinations, including retrograde and polar obits. This result confirms the robustness the coorbital resonance at large inclination and encourages the search for retrograde coorbital companions of the solar system's planets.

18. Weakly nonlinear stability analysis of magnetohydrodynamic channel flow using an efficient numerical approach

Hagan, Jonathan; Priede, Jānis

2013-12-01

We analyze weakly nonlinear stability of a flow of viscous conducting liquid driven by pressure gradient in the channel between two parallel walls subject to a transverse magnetic field. Using a non-standard numerical approach, we compute the linear growth rate correction and the first Landau coefficient, which in a sufficiently strong magnetic field vary with the Hartmann number as μ 1˜ (0.814-i19.8)× 10^{-3}textit {Ha} and μ 2˜ (2.73-i1.50)× 10^{-5}textit {Ha}^{-4}. These coefficients describe a subcritical transverse velocity perturbation with the equilibrium amplitude |A|2=Re [μ 1]/Re [μ 2](textit {Re}c-textit {Re})˜ 29.8textit {Ha}5(textit {Re}c-textit {Re}), which exists at Reynolds numbers below the linear stability threshold textit {Re}c˜ 4.83× 104textit {Ha}. We find that the flow remains subcritically unstable regardless of the magnetic field strength. Our method for computing Landau coefficients differs from the standard one by the application of the solvability condition to the discretized rather than continuous problem. This allows us to bypass both the solution of the adjoint problem and the subsequent evaluation of the integrals defining the inner products, which results in a significant simplification of the method.

19. Impact of hydrothermal alteration on lava dome stability: a numerical modelling approach

Detienne, Marie; Delmelle, Pierre

2016-04-01

Lava domes are a common feature of many volcanoes worldwide. They represent a serious volcanic hazard as they are prone to repeated collapses, generating devastating debris avalanches and pyroclastic flows. While it has long been known that hydrothermal alteration degrades rock properties and weakens rock mass cohesion and strength, there is still little quantitative information allowing the description of this effect and its consequences for assessing the stability of a volcanic rock mass such as a lava dome. In this study, we use the finite difference numerical model FLAC 3D to investigate the impact of hydrothermal alteration on the stability of a volcanic dome lying on a flat surface. Different hydrothermal alteration distributions were tested to encompass the variability observed in natural lava domes. Rock shear strength parameters (minimum, maximum and mean cohesion "c" and friction angle "φ" values) representative of various degrees of hydrothermal rock alteration were used in the simulations. The model predicts that reduction of the basement rock's shear strength decreases the factor of safety significantly. A similar result is found by increasing the vertical and horizontal extension of hydrothermal alteration in the basement rocks. In addition, pervasive hydrothermal alteration within the lava dome is predicted to exert a strong negative influence on the factor of safety. Through reduction of rock porosity and permeability, hydrothermal alteration may also affect pore fluid pressure within a lava dome. The results of new FLAC 3D runs which simulate the effect of hydrothermal alteration-induced pore pressure changes on lava dome stability will be presented.

20. Experimental and Numerical Investigation of a Swirl Stabilized Premixed Combustor under Cold Flow Conditions

SciTech Connect

P.A. Strakey; M.J. Yip

2007-07-01

Planar velocity measurements under cold-flow conditions in a swirl-stabilized dump combustor typical of land-based gas turbine combustors were carried out using two-dimensional particle image velocimetry (PIV). Axial, radial, and tangential velocity components were measured sequentially using two experimental configurations. Mean and root-mean-squared velocity components are presented along with instantaneous realizations of the flowfield. A numerical study of the flowfield using large-eddy simulation (LES) and Reynolds-averaged Navier-Stokes (RANS) techniques was conducted in an effort to help understand the complex hydrodynamics observed in the experiments. The agreement between the experimental data and LES simulation was good with both showing evidence of a precessing vortex core. The results of the RANS simulation were not as encouraging. The results provide a fundamental understanding of the complex flowfield associated with the relatively simple geometry and also serve as a baseline validation dataset for further numerical simulations of the current geometry. Validation of LES models in a highly swirled, nonreacting flowfield such as the work presented here is an essential step towards more accurate prediction in a reacting environment.

1. A numerical experiment on the equilibrium and stability of a rotating galactic bar

NASA Technical Reports Server (NTRS)

Miller, R. H.; Vandervoort, P. O.; Welty, D. E.; Smith, B. F.

1982-01-01

A self-consistent, three-dimensional numerical experiment is performed on an N-body system whose initial state is a realization of a certain theoretical model of a rotating triaxial galaxy. The model is a stellar-dynamical counterpart of a uniformly rotating polytrope of index equal to 0.5. The aim of the experiment is to study the equilibrium of the system and, in particular, to test its stability. The experimental system behaves in the mean like a realization of the theoretical model for at least seven crossing times. The principal departure of the system from equilibrium is an oscillation which is identified as a radial pulsation. There is no indication in its behavior that the system is unstable with respect to anu mode with an e-folding time shorter than or of the order of two crossing times. Certain changes that occur in the state of the system are interpreted, with the aid of the theoretical model, as secular changes which result from a slight failure of our numerical methods to conserve the mass, energy, and angular momentum of the system; these effects are small enough that they do not vitiate the experiment on a dynamical time scale.

2. The importance of trochanteric lag screws to achieve primary stability in cementless fixation of the RM hip prosthesis.

PubMed

Heitemeyer, U; Hierholzer, G; Haines, J

1987-01-01

To allow the bony incorporation of a cementless prosthesis it is important to achieve stability at the time of operation. To neutralize tension and torsional stresses the RM-shaft prosthesis is fixed with two lag screws in the trochanteric part of the femur. By measuring the applied torque intraoperatively we could demonstrate that the threads of the screws found a better grip when inserted from the bone to the prosthesis. Thus, the stronger fixation of the screws enhanced the primary stability of the cementless prosthesis. PMID:3566504

3. Explicit numerical formulas of improved stability and accuracy for the solution of parabolic equations

NASA Technical Reports Server (NTRS)

1979-01-01

A class of explicit numerical formulas which involve next nearest neighbor as well as nearest neighbor points are explored in this paper. These formulas are formal approximations to the linear parabolic partial-differential equation of first order in time and second order in distance. It was found that some of these formulas can employ time steps as much as four times that for the conventional explicit technique without becoming unstable. Others showed improved accuracy for a given time step and spatial grid spacing. One formula achieved a steady-state solution of specified accuracy for an example problem in less than 4 percent of the total computational time required by the conventional explicit technique.

4. Numerical simulation of compressible Couette flow stability: Non-linear and three-dimensional effects

Saati, Abdulmannan Abdulhamid

1991-02-01

The direct numerical simulation of the stability and transition of compressible Couette flow is studied. The effects of a constant body force along the vertical direction are also studied. Cartesian geometry is adopted to approximate Couette flow produced in the gap between two coaxial cylinders rotating at high-speed, with the body force representing the effects of the centrifugal force. A new, compressible flow solver for two- and three-dimensional, time dependent Navier-Stokes equations, using both the MacCormack and the high-order Two-Four methods was developed. In order to facilitate the simulations with greater detail and accuracy, a high-speed supercomputer with large core memory is required. Thus, the computer code was written in FORTRAN for its execution on the CRAY2, at NASA Langley. In a concurrent effort, in order to study the feasibility and efficiency of massively parallel super-computers and to speed up the computations, the work was further extended by rewriting the computer code in both C* and PARIS languages, for execution on the massively parallel Connection Machine CM 2 at the University of Colorado. Extensive testing of this new computer code was performed using wave propagation problems involving small- and large-amplitude two- and three-dimensional disturbances. Numerical simulations on the stability of compressible Couette flow between two infinite, parallel plates, with the inclusion of (1) a sudden body force, and (2) a body force in equilibrium, were performed. First, two-dimensional disturbances were considered and then the work was extended by considering three-dimensional disturbances on the rectangular Couette flow problem. Effects of body force magnitude, Mach number, and Reynolds number were also investigated. The simulations provide excellent agreement with the linear theory, thus documenting the phase and amplitude accuracy of the computed results; the overall amplitude error remains less than one percent. The results show that

5. School Stability: Improving Academic Achievement for NJ Foster Children. Policy Brief

ERIC Educational Resources Information Center

Bernard-Rance, Kourtney; Parello, Nancy

2014-01-01

Children in New Jersey's foster care system are more likely to remain in their home school when they enter foster care, thanks to a law passed in 2010, giving these fragile children improved educational stability. The law allows children to remain in their "school of origin" when they are placed in foster care, even if the foster…

6. Investigating motion and stability of particles in flows using numerical models

Khurana, Nidhi

The phenomenon of transport of particles in a fluid is ubiquitous in nature and a detailed understanding of its mechanism continues to remain a fundamental question for physicists. In this thesis, we use numerical methods to study the dynamics and stability of particles advected in flows. First, we investigate the dynamics of a single, motile particle advected in a two-dimensional chaotic flow. The particle can be either spherical or ellipsoidal. Particle activity is modeled as a constant intrinsic swimming velocity and stochastic fluctuations in both the translational and rotational motions are also taken into account. Our results indicate that interaction of swimming with flow structures causes a reduction in long-term transport at low speeds. Swimmers can get trapped at the transport barriers of the flow. We show that elongated swimmers respond more strongly to the dynamical structures of the flow field. At low speeds, their macroscopic transport is reduced even further than in the case of spherical swimmers. However, at high speeds these elongated swimmers tend to get attracted to the stable manifolds of hyperbolic fixed points, leading to increased transport. We then investigate the collective dynamics of a system of particles. The particles may interact both with each other and with the background flow. We focus on two different cases. In the fist case, we examine the stability of aggregation models in a turbulent-like flow. We use a simple aggregation model in which a point-like particle moves with a constant intrinsic speed while its velocity vector is reoriented according to the average direction of motion of its neighbors. We generate a strongly fluctuating, spatially correlated background flow using Kinematic Simulation, and show that flocks are highly sensitive to this background flow and break into smaller clusters. Our results indicate that such environmental perturbations must be taken into account for models which aim to capture the collective

7. Numerical Modeling of a Vortex Stabilized Arcjet. Ph.D. Thesis, 1991 Final Report

NASA Technical Reports Server (NTRS)

Pawlas, Gary E.

1992-01-01

Arcjet thrusters are being actively considered for use in Earth orbit maneuvering applications. Experimental studies are currently the chief means of determining an optimal thruster configuration. Earlier numerical studies have failed to include all of the effects found in typical arcjets including complex geometries, viscosity, and swirling flow. Arcjet geometries are large area ratio converging nozzles with centerbodies in the subsonic portion of the nozzle. The nozzle walls serve as the anode while the centerbody functions as the cathode. Viscous effects are important because the Reynolds number, based on the throat radius, is typically less than 1,000. Experimental studies have shown that a swirl or circumferential velocity component stabilizes a constricted arc. This dissertation describes the equations governing flow through a constricted arcjet thruster. An assumption that the flowfield is in local thermodynamic equilibrium leads to a single fluid plasma temperature model. An order of magnitude analysis reveals the governing fluid mechanics equations are uncoupled from the electromagnetic field equations. A numerical method is developed to solve the governing fluid mechanics equations, the Thin Layer Navier-Stokes equations. A coordinate transformation is employed in deriving the governing equations to simplify the application of boundary conditions in complex geometries. An axisymmetric formulation is employed to include the swirl velocity component as well as the axial and radial velocity components. The numerical method is an implicit finite-volume technique and allows for large time steps to reach a converged steady-state solution. The inviscid fluxes are flux-split, and Gauss-Seidel line relaxation is used to accelerate convergence. Converging-diverging nozzles with exit-to-throat area ratios up to 100:1 and annular nozzles were examined. Quantities examined included Mach number and static wall pressure distributions, and oblique shock structures. As

8. Numerical Modeling of the Stability of Face-Centered Cubic Metals with High Vacancy Concentration

SciTech Connect

Brian P. Somerday; M. I. Baskes

1998-12-01

The objective of this research is to assess the possibility of forming an atomically porous structure in a low-density metal, e.g., Al with vacancies up to 0.20/lattice site; and to examine the effects of hydrogen and vacancy concentration on the stability of an atomically porous structure that has been experimentally produced in nickel. The approach involves numerical modeling using the Embedded-Atom Method (EAM). High vacancy concentrations cause the Al lattice to disorder at 300K. In contrast, Ni retains the face-centered-cubic structure at 300K for vacancy concentrations up to 0.15 Vac/lattice site. Unexpectedly, the lattice with 0.15 Vac/lattice site is more stable than the lattice with 0.10 or 0.20 Vac/lattice site. The Ni systems with 0.10 and 0.15 Vac/lattice site exhibit domains consisting of uniform lattice rotations. The Ni lattice with 0.15 Vac/lattice site is more stable with an initial distribution of random vacancies compared to ordered vacancies. The equilibrium lattice structures of Ni a d Al containing vacancies and H are less ordered to structures with vacancies only at 300K.

9. Prevalence, Stability, and Functionality of Achievement Goal Profiles in Mathematics from Third to Seventh Grade

ERIC Educational Resources Information Center

Schwinger, Malte; Wild, Elke

2012-01-01

This study examined the longitudinal trajectories of achievement goal profiles in mathematics from third to seventh grade in a sample of 302 German students. Latent profile analyses were conducted separately for each school year and revealed three subgroups of students with distinct goal profiles labeled high multiple goals, moderate multiple…

10. Stability of Self-Esteem Ratings and Their Relation to Academic Achievement: A Longitudinal Study

ERIC Educational Resources Information Center

Rubin, Rosalyn A.

1978-01-01

In a sample of 380 children, self-esteem as measured by the Coopersmith Self-Esteem Inventory became more stable and more highly correlated with school achievement as the children grew older. Children's self-perceptions appear less firmly established, and therefore may be more responsive to intervention, at earlier ages. (Author)

11. Study on sub-cycling algorithm for flexible multi-body system: stability analysis and numerical examples

Miao, J. C.; Zhu, P.; Shi, G. L.; Chen, G. L.

2008-01-01

Numerical stability is an important issue for any integral procedure. Since sub-cycling algorithm has been presented by Belytschko et al. (Comput Methods Appl Mech Eng 17/18: 259-275, 1979), various kinds of these integral procedures were developed in later 20 years and their stability were widely studied. However, on how to apply the sub-cycling to flexible multi-body dynamics (FMD) is still a lack of investigation up to now. A particular sub-cycling algorithm for the FMD based on the central difference method was introduced in detail in part I (Miao et al. in Comp Mech doi: 10.1007/s00466-007-0183-9) of this paper. Adopting an integral approximation operator method, stability of the presented algorithm is transformed to a generalized eigenvalue problem in the paper and is discussed by solving the problem later. Numerical examples are performed to verify the availability and efficiency of the algorithm further.

12. Nonlinear effects in Paul traps operated in the second stability region: analytical analysis and numerical verification.

PubMed

Xiong, Caiqiao; Zhou, Xiaoyu; Zhang, Ning; Zhan, Lingpeng; Chen, Suming; Nie, Zongxiu

2014-11-01

Paul trap working in the second stability region has long been recognized as a possible approach for achieving high-resolution mass spectrometry (MS), which however is still far away from the experimental implementations because of the narrow working area and inefficient ion trapping. Full understanding of the ion motional behavior is helpful for solving the problem. In this article, the ion motion in a superimposed octopole field, which was characterized by the nonlinear Mathieu equation, was solved analytically using Poincare-Lighthill-Kuo (PLK) method. This method equivalently described the nonlinear disturbance by an effective quadrupole field with perturbed Mathieu parameters, a(u) and q(u), which would bring huge convenience in the studies of nonlinear ion dynamics and was, therefore, used for rapid evaluation of the nonlinear effects of ion motion. Fourth-order Runge-Kutta method (4th R-K) indicated the error of PLK for characterizing the frequency shift of ion motion was within 15%. PMID:25183226

13. An experimental and numerical study on the stability and propagation of laminar premixed flames

Vagelopoulos, Christina Maria

The laminar flame speed is a very important property of laminar premixed flames, especially for the validation of chemical kinetics and modeling of turbulent combustion. The counterflow technique is one of the best approaches for the experimental determination of this property because it allows for the establishment of planar, nearly adiabatic, steady, quasi-one dimensional flames that are subjected to well-defined aerodynamic strain rate. However non-linear effects as the strain rate goes to zero lead to overprediction of the laminar flame speed. In the present study these non-linear effects were investigated experimentally and numerically and significant overprediction was verified, particularly for weakly-burning hydrogen/air flames. Subsequently effort was made to establish and study flame properties at a very-low strain rate regime and qualitative and quantitative conclusions were drawn for the stability of the flame surface subjected to very low aerodynamic strain rate, coupled with the effect of gravity and preferential diffusion. A new experimental technique was developed, based on the observation that if a laminar premixed flame undergoes a transition from planar to Bunsen the strain rate changes from positive to negative values and a near-zero strain-rate regime is established. Flame speed measurements were conducted by using LDV for this regime; the flame speed measured is the true laminar flame speed and this is the first time that this property is directly and experimentally measured. Particle Streak velocimetry was developed to evaluate the strain-rates for near-zero strain-rate regime. The laminar flame speed was measured for atmospheric methane/air, ethane/air and propane/air mixtures for the whole range of equivalence ratios; the new data are lower when compared to previous ones and the overprediction is at the order of 15%.

14. Anion control as a strategy to achieve high-mobility and high-stability oxide thin-film transistors

PubMed Central

Kim, Hyun-Suk; Jeon, Sang Ho; Park, Joon Seok; Kim, Tae Sang; Son, Kyoung Seok; Seon, Jong-Baek; Seo, Seok-Jun; Kim, Sun-Jae; Lee, Eunha; Chung, Jae Gwan; Lee, Hyungik; Han, Seungwu; Ryu, Myungkwan; Lee, Sang Yoon; Kim, Kinam

2013-01-01

Ultra-definition, large-area displays with three-dimensional visual effects represent megatrend in the current/future display industry. On the hardware level, such a “dream” display requires faster pixel switching and higher driving current, which in turn necessitate thin-film transistors (TFTs) with high mobility. Amorphous oxide semiconductors (AOS) such as In-Ga-Zn-O are poised to enable such TFTs, but the trade-off between device performance and stability under illumination critically limits their usability, which is related to the hampered electron-hole recombination caused by the oxygen vacancies. Here we have improved the illumination stability by substituting oxygen with nitrogen in ZnO, which may deactivate oxygen vacancies by raising valence bands above the defect levels. Indeed, the stability under illumination and electrical bias is superior to that of previous AOS-based TFTs. By achieving both mobility and stability, it is highly expected that the present ZnON TFTs will be extensively deployed in next-generation flat-panel displays. PMID:23492854

15. Anion control as a strategy to achieve high-mobility and high-stability oxide thin-film transistors.

PubMed

Kim, Hyun-Suk; Jeon, Sang Ho; Park, Joon Seok; Kim, Tae Sang; Son, Kyoung Seok; Seon, Jong-Baek; Seo, Seok-Jun; Kim, Sun-Jae; Lee, Eunha; Chung, Jae Gwan; Lee, Hyungik; Han, Seungwu; Ryu, Myungkwan; Lee, Sang Yoon; Kim, Kinam

2013-01-01

Ultra-definition, large-area displays with three-dimensional visual effects represent megatrend in the current/future display industry. On the hardware level, such a "dream" display requires faster pixel switching and higher driving current, which in turn necessitate thin-film transistors (TFTs) with high mobility. Amorphous oxide semiconductors (AOS) such as In-Ga-Zn-O are poised to enable such TFTs, but the trade-off between device performance and stability under illumination critically limits their usability, which is related to the hampered electron-hole recombination caused by the oxygen vacancies. Here we have improved the illumination stability by substituting oxygen with nitrogen in ZnO, which may deactivate oxygen vacancies by raising valence bands above the defect levels. Indeed, the stability under illumination and electrical bias is superior to that of previous AOS-based TFTs. By achieving both mobility and stability, it is highly expected that the present ZnON TFTs will be extensively deployed in next-generation flat-panel displays. PMID:23492854

16. Enhanced erythrocyte suspension layer stability achieved by surface tension lowering additives

NASA Technical Reports Server (NTRS)

Omenyi, S. N.; Snyder, R. S.; Absolom, D. R.; Van Oss, C. J.; Neumann, A. W.

1982-01-01

In connection with a fractionation procedure involving the separation of particles, a dilute suspension of these particles in a liquid is carefully layered on a dense liquid. Under ideal conditions, the suspension forms a zone of finite thickness with a 'sharp' interface between the suspension layer and the supporting liquid. Under an applied field, e.g., gravitational or electrical, the particles in the suspension layer migrate to form different layers according to their size and/or density or according to their electrophoretic mobilities. However, in many cases the ideal conditions necessary for the fractionation process are not obtained. Many studies have been conducted to explore the reasons for suspension layer 'instability'. The present investigation represents an extension of a study conducted by Omenyi et al. (1981). An electrostatic repulsion-van der Waals mechanism was used to study the stability of fixed erythrocyte suspensions layered on a D2O cushion.

17. Body temperature stability achieved by the large body mass of sea turtles.

PubMed

Sato, Katsufumi

2014-10-15

To investigate the thermal characteristics of large reptiles living in water, temperature data were continuously recorded from 16 free-ranging loggerhead turtles, Caretta caretta, during internesting periods using data loggers. Core body temperatures were 0.7-1.7°C higher than ambient water temperatures and were kept relatively constant. Unsteady numerical simulations using a spherical thermodynamic model provided mechanistic explanations for these phenomena, and the body temperature responses to fluctuating water temperature can be simply explained by a large body mass with a constant thermal diffusivity and a heat production rate rather than physiological thermoregulation. By contrast, body temperatures increased 2.6-5.1°C in 107-152 min during their emergences to nest on land. The estimated heat production rates on land were 7.4-10.5 times the calculated values in the sea. The theoretical prediction that temperature difference between body and water temperatures would increase according to the body size was confirmed by empirical data recorded from several species of sea turtles. Comparing previously reported data, the internesting intervals of leatherback, green and loggerhead turtles were shorter when the body temperatures were higher. Sea turtles seem to benefit from a passive thermoregulatory strategy, which depends primarily on the physical attributes of their large body masses. PMID:25147244

18. Achieving synergy between chemical oxidation and stabilization in a contaminated soil.

PubMed

Srivastava, Vipul J; Hudson, Jeffrey Michael; Cassidy, Daniel P

2016-07-01

Eight in situ solidification/stabilization (ISS) amendments were tested to promote in situ chemical oxidation (ISCO) with activated persulfate (PS) in a contaminated soil. A 3% (by weight) dose of all ISS amendments selected for this study completely activated a 1.5% dose of PS within 3 h by raising temperatures above 30 °C (heat activation) and/or increasing pH above 10.5 (alkaline activation). Heat is released by the reaction of CaO with water, and pH increases because this reaction produces Ca(OH)2. Heat activation is preferred because it generates 2 mol of oxidizing radicals per mole of PS, whereas alkaline activation releases only 1. The relative contribution of heat vs. alkaline activation increased with CaO content of the ISS amendment, which was reflected by enhanced contaminant oxidation with increasing CaO content, and was confirmed by comparing to controls promoting purely heat or alkaline (NaOH) activation. The test soil was contaminated with benzene, toluene, ethylbenzene, and xylenes (BTEX) and polycyclic aromatic hydrocarbons (PAH), particularly naphthalene (NAP). ISS-activated PS oxidized between 47% and 84% of the BTEX & NAP, and between 13% and 33% of the higher molecular weight PAH. ISS-activated PS reduced the leachability of BTEX & NAP by 76%-91% and of the 17 PAH by 83%-96%. Combined ISCO/ISS reduced contaminant leachability far than ISCO or ISS treatments alone, demonstrating the synergy that is possible with combined remedies. PMID:27088536

19. Numerical Models of Slow Slip and Dynamic Rupture Including Dilatant Stabilization and Thermal Pressurization

2009-12-01

Seismic and geodetic observations indicate that slow-slip events (SSE) occur down-dip of locked megathrusts, in areas of high pore-pressure, p. We suggest that at low effective stress (σ -p) dilatancy stabilizes rate-weakening faults, whereas at higher (σ -p) thermal pressurization overwhelms dilatancy leading to dynamic slip. 2D simulations include rate-state (slip-law) friction with Linker-Dieterich normal stress effect, Segall-Rice dilatancy linked to state evolution, and heat and pore-fluid flow normal to the fault. The fault is loaded by down-dip slip at v∞ . We discretize the fault normal direction with log spacing, and employ explicit-implicit time integration to improve speed and accuracy. The governing equations involve numerous physical parameters, but relatively few non-dimensional groups. Ep={ɛ h}/ [2 {β (σ - p∞)}√ {{v{∞ }}/ {chyd}dc] and ET=({f_0Λ }/{2ρ cp})√ {{chydd_c v∞}}/{cth}2 represent dilatancy and shear heating efficiency, respectively. For a nominal set of parameters (given below), spatially uniform properties, and σ -p=1MPa, (Ep = 1.5 × 10-3, ET = 3 × 10-5), we find a series of propagating SSE, that are stabilized by dilatancy-induced drops in p at the rupture tips. For a broad range of parameters we observe slow-slip events driven by down-dip slip, with negative stress drop, as well as faster (but quasi-static) events that relax the accumulated stress. At 10 MPa effective stress, the models exhibit both SSE and dynamic ruptures. Following dynamic stress drops, a sequence of slow slip events is driven from the down-dip end of the fault, with generally increasing maximum slip-speeds. We also consider spatially variable (σ -p), ranging from 2 MPa down-dip to 10 MPa up-dip (with arctangent distribution such that 80% of the variation occurs across 20% of the fault), and uniform material properties. The models exhibit both SSE and dynamic events. Following a dynamic rupture there are initially no slow events, and the

20. Analysis of numerical stability and amplification matrices: Fourth-order Runge-Kutta methods

NASA Technical Reports Server (NTRS)

Kennedy, E. W.

1979-01-01

Amplification matrices, numerical kernels, stable, and exponentially stable numerical solutions are examined. The various techniques involved in these concepts are applied to certain systems that have Jordan forms, which are nondiagonal, with particular interest in the case of imaginary or zero eigenvalues.

1. 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

2. Sensitivity of free bar morphology in rivers to secondary flow modeling: Linear stability analysis and numerical simulation

Iwasaki, Toshiki; Shimizu, Yasuyuki; Kimura, Ichiro

2016-06-01

A number of numerical models have been proposed to understand and simulate fluvial river morphodynamics; however, it is somewhat unclear whether all the models are able to consistently simulate flow-bed instability phenomena. This study investigates the sensitivity of free bar morphology in rivers to secondary flow models used in depth-averaged models using linear stability analyses and numerical simulations. Both the linear analyses and numerical simulations suggest that under certain hydraulic conditions, an equilibrium-type secondary flow model, which has been widely used in river morphodynamic models, fails to generate a finite wavelength and bar mode, allowing the inception of bars of infinitely short scale and infinitely high mode. Using a nonequilibrium-type secondary flow model avoids the unphysical formation of these incipient free bars, and gives better solutions regarding finite amplitude bars. Since free bars are essential, intrinsic river morphological features, the findings of this study can be applied to a wide range of river morphodynamic calculations.

3. Stability analysis of Western flank of Cumbre Vieja volcano (La Palma) using numerical modelling

Bru, Guadalupe; Gonzalez, Pablo J.; Fernandez-Merodo, Jose A.; Fernandez, Jose

2016-04-01

La Palma volcanic island is one of the youngest of the Canary archipelago, being a composite volcano formed by three overlapping volcanic centers. There are clear onshore and offshore evidences of past giant landslides that have occurred during its evolution. Currently, the active Cumbre Vieja volcano is in an early development state (Carracedo et al., 2001). The study of flank instability processes aim to assess, among other hazards, catastrophic collapse and potential tsunami generation. Early studies of the potential instability of Cumbre Vieja volcano western flank have focused on the use of sparse geodetic networks (Moss et al. 1999), surface geological mapping techniques (Day et al. 1999) and offshore bathymetry (Urgeles et al. 1999). Recently, a dense GNSS network and satellite radar interferometry results indicate ground motion consistent with deep-seated creeping processes (Prieto et al. 2009, Gonzalez et al. 2010). In this work, we present a geomechanical advanced numerical model that captures the ongoing deformation processes at Cumbre Vieja. We choose the Finite Elements Method (FEM) which is based in continuum mechanics and is the most used for geotechnical applications. FEM has the ability of using arbitrary geometry, heterogeneities, irregular boundaries and different constitutive models representative of the geotechnical units involved. Our main contribution is the introduction of an inverse approach to constrain the geomechanical parameters using satellite radar interferometry displacements. This is the first application of such approach on a large volcano flank study. We suggest that the use of surface displacements and inverse methods to rigorously constrain the geomechanical model parameter space is a powerful tool to understand volcano flank instability. A particular important result of the studied case is the estimation of displaced rock volume, which is a parameter of critical importance for simulations of Cumbre Vieja tsunamigenic hazard

4. The Effects of Teaching Numerical Control Concepts Via Simulator Versus Non-Simulator Activities on the Achievement, Programming Proficiency and Attitude of High School Students.

ERIC Educational Resources Information Center

Pine, Douglas Taylor

This study utilized 120 metalworking students and six teachers from Columbus, Ohio area high schools to ascertain the effects of teaching numerical control to industrial arts students by means of simulator-aided activities versus nonsimulator aided activities. Scores obtained from an achievement test, attitude inventory, and word address…

5. Stability analysis of numerical solutions of wave propagation when pseudoviscosity is replaced by real viscosity

Fanget, A.

1996-05-01

The artificial viscocity method was originally introduced by Von Neuman and Richtmeyer [R.D. Richtmeyer and K. W. Morton, in Difference methods for Initial Value Problem (Interscience, New York, 1967)] to spread the shock over several cells of the mesh. This method is useful if the physical phenomena present in the shock front are not of interest. In some cases, such as for polymers, the amount of viscosity is large enough to spread the front shock naturally. In our problem, pseudoviscosity is ignored. A stability analysis shows that time step restrictions for stability can be less stringent than with pseudoviscosity.

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

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.

7. Numerical stability analysis and flow simulation of lid-driven cavity subjected to high magnetic field

Marioni, L.; Bay, F.; Hachem, E.

2016-05-01

In this work, we study the flow of a conducting fluid inside a two-dimensional square domain. The problem is solved by using a variational multiscale finite element approach. The study focuses on a high magnetic interaction parameter range and high Reynolds number. Under the imposition of a high magnetic field, the flow gets regularized, but fast transient phenomena take place, which could lead to numerical errors. An expression to compute the maximum time step that guarantees convergence in explicit schemes is proposed and validated through numerical tests.

8. Achieving ligament stability and correct rotational alignment of the femur in knee arthroplasty: a study using the Medial Pivot knee.

PubMed

Shakespeare, David; Kinzel, Vera; Ledger, Michael

2005-12-01

In a series of 90 Medial Pivot arthroplasties rotational alignment of the femur was achieved by provisionally reconstructing the lateral side of the joint and tensioning the medial side with feeler gauges. Axial CT scans were employed to measure the rotational alignment relative to surgical epicondylar axis. In valgus knees the cutting block was externally rotated to adjust for posterolateral bone loss. The mean rotational alignment of the femur was 0.6 degrees of external rotation (S.D. 1.3, range 3 degrees of ER to 4 degrees of IR). The mean laxity of the medial ligament was 1 mm in flexion (SD 1, range 0-5 mm) and 0.5 mm in flexion (S.D. 0.5, range 0-2 mm) In those knees in which the medial ligament had been released the CT alignment was perfect, but when internally rotated against the hip 3-4 mm of gapping was noted. In valgus knees the mean rotation of the femoral component was 0.8 degrees of internal rotation (S.D. 1.5, range 1 degrees of IR to 4 degrees of ER). In spite of externally rotating the cutting block there was still a tendency to internally rotate the femur in some knees. This simple technique achieves the two goals of ligament stability and correct rotational alignment in a high proportion of cases. It may be applicable to any instrument system which employs posterior referencing. PMID:15967668

9. Numerical computation of the linear stability of the diffusion model for crystal growth simulation

SciTech Connect

Yang, C.; Sorensen, D.C.; Meiron, D.I.; Wedeman, B.

1996-12-31

We consider a computational scheme for determining the linear stability of a diffusion model arising from the simulation of crystal growth. The process of a needle crystal solidifying into some undercooled liquid can be described by the dual diffusion equations with appropriate initial and boundary conditions. Here U{sub t} and U{sub a} denote the temperature of the liquid and solid respectively, and {alpha} represents the thermal diffusivity. At the solid-liquid interface, the motion of the interface denoted by r and the temperature field are related by the conservation relation where n is the unit outward pointing normal to the interface. A basic stationary solution to this free boundary problem can be obtained by writing the equations of motion in a moving frame and transforming the problem to parabolic coordinates. This is known as the Ivantsov parabola solution. Linear stability theory applied to this stationary solution gives rise to an eigenvalue problem of the form.

10. Application of numerical simulation using a progressive failure approach to underground-coal-mine stability analysis

SciTech Connect

Ash, N.F.

1987-01-01

Stability in underground coal mines is of major concern to the coal industry due to its effect on both safety and productivity. Consequently, this can have a great influence on the design of efficient mine systems. In this work a progressive failure approach was used to simulate underground coal mine stability at two different mines. The two mines considered have different characteristics. Two- and three- dimensional finite element models were created to model different areas of a longwall mine. Different chain pillar configurations were considered and the resulting stress distributions were comparable to field measurements. A complete mine section was successfully modeled taking into consideration face advancement. The roof above entry intersections was also modeled using laminated composite simulation and the finite element method. The results showed trends similar to field observations. In addition, the progressive development of subsidence for the two different mines was simulated. The same variation in subsidence behavior recorded at the mine was realized in the finite element simulation. The progressive failure approach used in this work can successfully simulate underground coal mine stability. It can also be a helpful tool in the design of more efficient mine systems which can increase productivity and maintain a high level of safety.

11. A numerical study of global axisymmetric states and their stability to quasi-geostrophic disturbances

NASA Technical Reports Server (NTRS)

Henderson, H. W.; Dutton, J. A.

1985-01-01

A steady, axisymmetric model of the general circulation is developed as a basis for climate stability studies. The model includes the effects of heating, rotation, and internal friction, but neglects topography. It is assumed that the axisymmetric flow may be modeled by making the Boussinesq and deep convection approximations. The hydrostatic assumption is not made, thus permitting the advective terms to be included in the vertical equation of motion. The initial set of five primitive equations is reduced to three equations in terms of the zonal velocity, meridional streamfunction, and the potential temperature perturbation. The application of the Boussinesq, deep convection, and quasi-geostrophic assumptions limits the ranges of the heating and rotation rates. For values not too far from typical atmospheric values, the model produces a stability boundary separating Hadley from Rossby flow. The boundary is characterized by a particular value of vertical wind shear, which suggests that baroclinic instability is the primary mechanism for the loss of stability. The initial growth rates are largest for longitudinal waves 4-7, also in agreement with studies of baroclinic instability.

12. Experimental and Numerical Investigation of Vortical Structures in Lean Premixed Swirl-Stabilized Combustion

Taamallah, Soufien; Chakroun, Nadim; Shanbhogue, Santosh; Kewlani, Gaurav; Ghoniem, Ahmed

2015-11-01

A combined experimental and LES investigation is performed to identify the origin of major flow dynamics and vortical structures in a model gas turbine's swirl-stabilized turbulent combustor. Swirling flows in combustion lead to the formation of complex flow dynamics and vortical structures that can interact with flames and influence its stabilization. Our experimental results for non-reacting flow show the existence of large scale precession motion. The precessing vortex core (PVC) dynamics disappears with combustion but only above a threshold of equivalence ratio. In addition, large scale vortices along the inner shear layer (ISL) are observed. These structures interact with the ISL stabilized flame and contribute to its wrinkling. Next, the LES setup is validated against the flow field's low-order statistics and point temperature measurement in relevant areas of the chamber. Finally, we show that LES is capable of predicting the precession motion as well as the ISL vortices in the reacting case: we find that ISL vortices originate from a vortex core that is formed right downstream of the swirler's centerbody. The vortex core has a conical spiral shape resembling a corkscrew that interacts - as it winds out - with the flame when it reaches the ISL.

13. Numerical Analysis of Large Telescopes in Terms of Induced Loads and Resulting Geometrical Stability

Upnere, S.; Jekabsons, N.; Joffe, R.

2013-03-01

Comprehensive numerical studies, involving structural and Computational Fluid Dynamics (CFD) analysis, have been carried out at the Engineering Research Institute "Ventspils International Radio Astron- omy Center" (VIRAC) of the Ventspils University College to investigate the gravitational and wind load effects on large, ground-based radio tele- scopes RT-32 performance. Gravitational distortions appear to be the main limiting factor for the reflector performance in everyday operation. Random loads caused by wind gusts (unavoidable at zenith) contribute to the fatigue accumulation.

14. Individual differences in children's mathematics achievement: The roles of symbolic numerical magnitude processing and domain-general cognitive functions.

PubMed

Vanbinst, K; De Smedt, B

2016-01-01

This contribution reviewed the available evidence on the domain-specific and domain-general neurocognitive determinants of children's arithmetic development, other than nonsymbolic numerical magnitude processing, which might have been overemphasized as a core factor of individual differences in mathematics and dyscalculia. We focused on symbolic numerical magnitude processing, working memory, and phonological processing, as these determinants have been most researched and their roles in arithmetic can be predicted against the background of brain imaging data. Our review indicates that symbolic numerical magnitude processing is a major determinant of individual differences in arithmetic. Working memory, particularly the central executive, also plays a role in learning arithmetic, but its influence appears to be dependent on the learning stage and experience of children. The available evidence on phonological processing suggests that it plays a more subtle role in children's acquisition of arithmetic facts. Future longitudinal studies should investigate these factors in concert to understand their relative contribution as well as their mediating and moderating roles in children's arithmetic development. PMID:27339010

15. Slope stability analysis for Valles Marineris, Mars: a numerical analysis of controlling conditions and failure types

Crosta, G.; Castellanza, R.; De Blasio, F.; Utili, S.

2012-04-01

Valles Marineris (VM hereafter) in the equatorial area of Mars exhibits several gravitative failures often involving the whole 6-8 km thickness of the valley walls. The failures have resulted in a series of long-runout landslides up to several hundred cubic kilometres in volume (Quantin et al., 2004), and the formation of sub-circular alcoves perched on the top. Several questions arise as to forces at play in the stability of the walls of VM, the geometrical shape of the alcoves and the shape and long-runout of the landslides (see for example Lucas et al., 2011). In this work, we concentrate on the stability analysis of the walls of VM with two precise questions in mind starting from past studies (Bigot-Cormier and Montgomery, 2006; Neuffer and Schultz, 2006, Schultz, 2002). The first concerns the properties of the materials that give origin to instability. We performed several finite element and discrete element calculations tailored to slope stability analysis based on the genuine shape of the walls of VM taken from the MOLA topographic data. We considered stratified and differently altered/degraded materials to define the range of physical mechanical properties required for failure to occur and to explain the discrete distribution of failures along the VM valley flanks. A second question addressed in this work is the geometrical shape of the sub-circular alcoves. Normally, these shapes are commonplace for slopes made of uniform and isotropic properties, and are also observed in subaqueous environment. We performed calculations taking into consideration the progressive failure in the slope showing the final results in terms of surface failure geometry. Bigot-Cormier, F., Montgomery, D.R. (2007) Valles Marineris landslides: Evidence for a strength limit to Martian relief? Earth and Planetary Science Letters, 260, 1-2, 15, 179-186 Lucas, A., Mangeney, A., Mège, D., and Bouchut, F., 2011. Influence of the scar geometry on landslide dynamics and deposits

16. Numerical assessment of efficiency and control stability of an HTS synchronous motor

Xian, Wei; Yuan, Weijia; Coombs, T. A.

2010-06-01

A high temperature superconducting (HTS) permanent magnet synchronous motor (PMSM) is designed and developed in Cambridge University. It is expected to become cost competitive with the conventional PMSM owing to its high efficiency, high power density, high torque density, etc. The structure and parameters of HTS PMSM are detailed. Both AC losses by transport current and applied filed in stator armature winding of HTS PMSM are also analyzed. Computed and simulated results of the characteristics of the HTS PMSM and conventional PMSM are compared. The improvement on stability of direct torque control (DTC) on the HTS PMSM is estimated, and proved by simulation on Matlab/Simulink.

17. A novel numerical approach and stability analysis of thermo-acoustic phenomenon in the Rijke tube problem

2013-11-01

The modeling of thermo-acoustic coupling in reactive flows represents a challenging task. In this study, we focus on the Rijke tube problem, which includes relevant features such as a compact acoustic source, an empirical modeling of the heat source, and non-linearities. This thermo-acoustic system features a complex dynamical behavior, which renders the characterization of the different encountered flow regimes difficult. In order to synthesize accurate time series, we tackle this problem from a numerical point-of-view, and start by proposing a dedicated solver designed for dealing with the underlying stiffness, in particular, the retarded time and the discontinuity at the location of the heat source. Convergence and parametric studies are carried out to assess the accuracy of the discretization, hence laying a foundation for a stability analysis of the semi-discrete system. This stability analysis is performed by means of the projection method proposed by Jarlebring, which alleviates the linearization of the retarded term, and is used to validate the numerical results. Finally, the focus is set on the application of the dynamic mode decomposition technique to study bifurcations.

18. A numerical study of the thermal stability of low-lying coronal loops

NASA Technical Reports Server (NTRS)

Klimchuk, J. A.; Antiochos, S. K.; Mariska, J. T.

1986-01-01

The nonlinear evolution of loops that are subjected to a variety of small but finite perturbations was studied. Only the low-lying loops are considered. The analysis was performed numerically using a one-dimensional hydrodynamical model developed at the Naval Research Laboratory. The computer codes solve the time-dependent equations for mass, momentum, and energy transport. The primary interest is the active region filaments, hence a geometry appropriate to those structures was considered. The static solutions were subjected to a moderate sized perturbation and allowed to evolve. The results suggest that both hot and cool loops of the geometry considered are thermally stable against amplitude perturbations of all kinds.

19. Development and stability of bed forms: a numerical analysis of dune pattern coarsening and giant dunes

Gao, Xin; Narteau, Clement; Rozier, Olivier

2015-04-01

We investigate the development and stability of transverse dunes for ranges of flow depths and velocities using a cellular automaton dune model. Subsequent to the initial bed instability, dune pattern coarsening is driven by bed form interactions. Collisions lead to two types of coalescence associated with upstream or downstream dominant dunes. In addition, a single collision-ejection mechanism enhances the exchange of mass between two consecutive bed forms (through-passing dunes). The power-law increases in wavelength and amplitude exhibit the same exponents, which are independent of flow properties. Contrary to the wavelength, dune height is not only limited by flow depth but also by the strength of the flow. Superimposed bedforms may propagate and continuously destabilize the largest dunes. Then, we identify three classes of steady-state transverse dune fields according to the periodicity in crest-to-crest spacing and the mechanism of size limitation. In all cases, the steady state is reached when the bed shear stress in the dune trough regions is close to its critical value for motion inception. Such a critical shear stress value is reached and maintained through the dynamic equilibrium between flow strength and dune aspect ratio. Comparisons with natural dune fields show that many of them may have reached such a steady state. Finally, we infer that the sedimentary patterns in the model may be used to bring new constraints on the stability of modern and ancient dune fields.

20. Convective flows generated by evaporation: experiments, linear stability analysis and numerical simulations

Dunstan, Jocelyn; Lee, Kyoung Jin; Park, Simon; Goldstein, Raymond E.

A novel form of convection was observed in a suspension of non-motile Photobacterium phosphoreum bacteria. The pattern resembles classical bioconvection, however this strain has limited if any motility, which excludes this possible explanation. After performing a series of control experiments we found that the convection was actually driven by the evaporation of the salty bacterial medium, and the same kind of plumes were observed using polystyrene beads suspended in water with salt added. A mathematical model was formulated for the process and studied using a linear stability analysis and finite element method simulations, reproducing most of the observed experimental features. From the linear stability analysis, a threshold in salt concentration to observe convective motion was obtained, as well as the wavelength of the pattern at the onset of the instability. This was complemented by finite element simulations, which produced plume dynamics remarkably similar to the experimental observations. Evaporation-driven convection on the millimeter scale has not been studied extensively, and its effect may have been underestimated in other experiments.

1. Numerical prediction of turbulent flame stability in premixed/prevaporized (HSCT) combustors

NASA Technical Reports Server (NTRS)

Winowich, Nicholas S.

1990-01-01

A numerical analysis of combustion instabilities that induce flashback in a lean, premixed, prevaporized dump combustor is performed. KIVA-II, a finite volume CFD code for the modeling of transient, multidimensional, chemically reactive flows, serves as the principal analytical tool. The experiment of Proctor and T'ien is used as a reference for developing the computational model. An experimentally derived combustion instability mechanism is presented on the basis of the observations of Proctor and T'ien and other investigators of instabilities in low speed (M less than 0.1) dump combustors. The analysis comprises two independent procedures that begin from a calculated stable flame: The first is a linear increase of the equivalence ratio and the second is the linear decrease of the inflow velocity. The objective is to observe changes in the aerothermochemical features of the flow field prior to flashback. It was found that only the linear increase of the equivalence ratio elicits a calculated flashback result. Though this result did not exhibit large scale coherent vortices in the turbulent shear layer coincident with a flame flickering mode as was observed experimentally, there were interesting acoustic effects which were resolved quite well in the calculation. A discussion of the k-e turbulence model used by KIVA-II is prompted by the absence of combustion instabilities in the model as the inflow velocity is linearly decreased. Finally, recommendations are made for further numerical analysis that may improve correlation with experimentally observed combustion instabilities.

2. Numerical simulation of a crack in the cement stabilized stone using cohesive zone models

Liu, Jing-hui; Wang, Duan-yi

2008-11-01

Reflective cracking arising from cracks in base materials has been a major distress of semi-rigid asphalt concrete road. Previous studies in base mixture cracking have typically considered the materials homogeneity. Adopting Digital Image Processing techniques and Finite Element Method, in the micro scale, the cement and aggregates are treated as distinct materials with different materials parameters. The potential crack zones are simulated by Cohesive Zone Model. The initiation and propagation of the crack in the cylindrical specimen under conventional Indirect Tensile Test (IDT) are modeled. The numerical results from the micromechanical analysis match well with the results from the macro experiment. Even though this study only presented a attempt to a numerical simulation of a simple IDT test, the theory and methods adopted by this study can be applied to the fatigue damage scenario under complicated loading conditions including material heterogeneity. It effectively allows researchers to link the micro-scale damage observed on the local scale with the real pavements failing on the global scale.

3. Numerical proof of stability of roll waves in the small-amplitude limit for inclined thin film flow

Barker, Blake

2014-10-01

We present a rigorous numerical proof based on interval arithmetic computations categorizing the linearized and nonlinear stability of periodic viscous roll waves of the KdV-KS equation modeling weakly unstable flow of a thin fluid film on an incline in the small-amplitude KdV limit. The argument proceeds by verification of a stability condition derived by Bar-Nepomnyashchy and Johnson-Noble-Rodrigues-Zumbrun involving inner products of various elliptic functions arising through the KdV equation. One key point in the analysis is a bootstrap argument balancing the extremely poor sup norm bounds for these functions against the extremely good convergence properties for analytic interpolation in order to obtain a feasible computation time. Another is the way of handling analytic interpolation in several variables by a two-step process carving up the parameter space into manageable pieces for rigorous evaluation. These and other general aspects of the analysis should serve as blueprints for more general analyses of spectral stability.

4. Numerical analysis of two-ring flame stabilizer flows in aircraft afterburners

SciTech Connect

Ravichandran, M.; Ganesan, V.; Sampath, S.

1995-07-01

A computational procedure is described for the calculation of isothermal flow fields of two-ring flame stabilizers in practical afterburners and also for the reacting flow parameters in a research afterburner. The predictions have been obtained using a finite volume solution procedure for the steady three-dimensional elliptic equations of fluid flow. The physical models include the {kappa}-{epsilon} turbulence model, eddy breakup model, two-step reaction model, droplet vaporization and combustion model, and six-flux radiation model. The presence of radial gutters and staggering the distance between two-ring gutters alter the afterburner flow patterns considerably, including the elimination of diffuser stall and flow separation losses near the diffuser. The predicted and measured emissions in the research afterburner geometry agree qualitatively. The effect of the fuel to air ratio and location of the fuel injector on the temperature and species concentration distribution are presented in detail.

5. Numerical Simulations of the Onset and Stability of Dynamical Mass Transfer in Binaries

D'Souza, Mario C. R.; Motl, Patrick M.; Tohline, Joel E.; Frank, Juhan

2006-05-01

Hydrodynamical simulations of semidetached, polytropic binary stars are presented in an effort to study the onset and stability of dynamical mass transfer events. Initial, synchronously rotating equilibrium models are constructed using a self-consistent field technique and then evolved with an Eulerian hydrodynamics code in a fully self-consistent manner. We describe code improvements introduced over the past few years that permit us to follow dynamical mass transfer events through more than 30 orbits. Mass transfer evolutions are presented for two different initial configurations: a dynamically unstable binary with initial mass ratio (donor/accretor) q0=1.3 that leads to a complete merger in ~10 orbits, and a double-degenerate binary with initial mass ratio q0=0.5 that, after some initial unstable growth of mass transfer, tends to separate as the mass transfer rate levels off.

6. Experimental and numerical investigation on the motorcycle front frame flexibility and its effect on stability

Cossalter, V.; Doria, A.; Massaro, M.; Taraborrelli, L.

2015-08-01

It is well known that front fork flexibility may have a significant effect on motorcycle stability. This work addresses the problem of developing lumped element models of the front fork from experimental results. The front forks of an enduro motorcycle and of a super sport motorcycle are characterized performing static, dynamic and modal tests by means of specific testing equipment. The concept of wheel twisting axis is proposed to characterize static and dynamic deformability of the front fork. Modal analysis results show the presence of two important modes of vibration of the front assembly in the low frequency range: the lateral mode and the longitudinal mode. Different lumped models are discussed and a new model that takes into account information obtained from static and dynamic tests is proposed. Simulations are carried out by means of a multibody code and show the effect of the front assembly deformability on the weave and wobble vibration modes.

7. A numerical and analytical study of nonlinear bifurcations associated with the morphological stability of two-dimensional single crystals

NASA Technical Reports Server (NTRS)

Brush, L. N.; Mcfadden, G. B.; Sekerka, R. F.

1990-01-01

The nonlinear stability of a two-dimensional single crystal of pure material in an undercooled melt is studied both analytically and numerically. The quasi-steady state approximation is used for the thermal fields, and the effects of different solid and liquid thermal conductivities, isotropic interfacial growth kinetics, and isotropic surface tension are included. The bifurcation analysis is performed by calculating the instantaneous value of the fundamental component of the local normal growth speed for an interface perturbed by a single Fourier shape component. Numerically, the fundamental component of the interfacial growth speed is found by Fourier analysis of the solution to an integrodifferential equation obeyed at the interface. Analytically, an expansion technique is used to derive a solvability condition defining each of these bifurcation points. The analytical and numerical results are in very close agreement. Almost all of the bifurcations are subcritical, and the results are presented by giving values of the Landau coefficient as a function of the different dimensionless parameters used in the model.

8. A numerical study of the stability of radiative shocks. [in accretion flows onto white dwarf stars

NASA Technical Reports Server (NTRS)

Imamura, J. N.; Wolff, M. T.; Durisen, R. H.

1984-01-01

Attention is given to the oscillatory instability of optically thin radiative shocks in time-dependent numerical calculations of accretion flows onto degenerate dwarfs. The present nonlinear calculations yield good quantitative agreement with the linear results obtained for oscillation frequencies, damping rates, and critical alpha-values. The fundamental mode and the first overtone in the shock radius and luminosity variations can be clearly identified, and evidence is sometimes seen for the second overtone. Time-dependent calculations are also performed which include additional physics relevant to degenerate dwarf accretion, such as electron thermal conduction, unequal electron and ion temperatures, Compton cooling, and relativistic corrections to the bremsstrahlung cooling law. All oscillatory modes are found to be damped, and hence stable, in the case of a 1-solar mass white dwarf accreting in spherical symmetry.

9. Stabilizing coherence with nested environments: a numerical study using kicked Ising models

González-Gutiérrez, C.; Villaseñor, E.; Pineda, C.; Seligman, T. H.

2016-08-01

We study a tripartite system of coupled spins, where a first set of one or two spins is our central system which is coupled to another set considered, the near environment, in turn coupled to the third set, the far environment. The dynamics considered are those of a generalized kicked spin chain in the regime of quantum chaotic dynamics. This allows us to test recent results that suggest that the presence of a far environment, coupled to the near environment, slows decoherence of the central system. After an extensive numerical study, we confirm previous results for extreme values and special cases. In particular, under a wide variety of circumstances an increasing coupling between near and far environment, slows decoherence, as measured by purity, and protects internal entanglement.

10. A numerical study of a marine subtropical stratus cloud layer and its stability

NASA Technical Reports Server (NTRS)

Moeng, C.-H.; Arakawa, A.

1980-01-01

A model for numerical simulation of stratus cloud layers is constructed by combining a second-order closure, turbulent transfer model with a thermal radiative transfer model. The turbulent transfer model allows water vapor saturation. The combined turbulence-radiation model is applied to both a horizontally uniform one-dimensional case and a horizontally nonuniform two-dimensional case. In the latter, the dynamics of mesoscale circulations are also incorporated. Results of the two-dimensional simulation show that the layer cloud instability occurs where the sea surface temperature is high and the large-scale subsidence is weak. The simulated instability is analyzed in view of an instability criterion, the eddy kinetic energy budget, and evaporative cooling near the cloud top.

11. The stability of Taylor bubbles in large-diameter tubes: direct numerical simulations

Dhanjal, Amanjalot; Saravan-Butler, Maya; Smith, Sydney; Yang, Junfeng; Matar, Omar

2015-11-01

Slug flow corresponds to intermittent Taylor bubbles and liquid slugs, and is widely observed in the oil-and-gas industry. The fluctuating flow rate caused by Taylor bubbles is problematical; thus, the destabilisation of this regime would be beneficial. To gain better understanding of this regime in vertical tubes, three-dimensional CFD simulations of Taylor air bubble rise in initially stagnant water and progressively larger diameter tubes, are carried out. Tubes with diameters in the range of 0.032m-0.290m and a height of 2m are considered. The topology of the Taylor bubbles and their rise velocity are predicted and validated against experimental results. Our results suggest that the wake of leading bubbles plays a key role in the deformation and break-up of trailing bubbles. Motivated by these results, the effect of bubble separation distance, and aspect ratio, on bubble stability and the slug flow regime is investigated. EPSRC Programme Grant, MEMPHIS, EP/K0039761/1.

12. Theoretical Hill-type muscle and stability: numerical model and application.

PubMed

Schmitt, S; Günther, M; Rupp, T; Bayer, A; Häufle, D

2013-01-01

The construction of artificial muscles is one of the most challenging developments in today's biomedical science. The application of artificial muscles is focused both on the construction of orthotics and prosthetics for rehabilitation and prevention purposes and on building humanoid walking machines for robotics research. Research in biomechanics tries to explain the functioning and design of real biological muscles and therefore lays the fundament for the development of functional artificial muscles. Recently, the hyperbolic Hill-type force-velocity relation was derived from simple mechanical components. In this contribution, this theoretical yet biomechanical model is transferred to a numerical model and applied for presenting a proof-of-concept of a functional artificial muscle. Additionally, this validated theoretical model is used to determine force-velocity relations of different animal species that are based on the literature data from biological experiments. Moreover, it is shown that an antagonistic muscle actuator can help in stabilising a single inverted pendulum model in favour of a control approach using a linear torque generator. PMID:24319495

13. Theoretical Hill-Type Muscle and Stability: Numerical Model and Application

PubMed Central

Schmitt, S.; Günther, M.; Rupp, T.; Bayer, A.; Häufle, D.

2013-01-01

The construction of artificial muscles is one of the most challenging developments in today's biomedical science. The application of artificial muscles is focused both on the construction of orthotics and prosthetics for rehabilitation and prevention purposes and on building humanoid walking machines for robotics research. Research in biomechanics tries to explain the functioning and design of real biological muscles and therefore lays the fundament for the development of functional artificial muscles. Recently, the hyperbolic Hill-type force-velocity relation was derived from simple mechanical components. In this contribution, this theoretical yet biomechanical model is transferred to a numerical model and applied for presenting a proof-of-concept of a functional artificial muscle. Additionally, this validated theoretical model is used to determine force-velocity relations of different animal species that are based on the literature data from biological experiments. Moreover, it is shown that an antagonistic muscle actuator can help in stabilising a single inverted pendulum model in favour of a control approach using a linear torque generator. PMID:24319495

14. Stability of submarine slopes in the northern South China Sea: a numerical approach

Zhang, Liang; Luan, Xiwu

2013-01-01

Submarine landslides occur frequently on most continental margins. They are effective mechanisms of sediment transfer but also a geological hazard to seafloor installations. In this paper, submarine slope stability is evaluated using a 2D limit equilibrium method. Considerations of slope, sediment, and triggering force on the factor of safety (FOS) were calculated in drained and undrained ( Φ=0) cases. Results show that submarine slopes are stable when the slope is <16° under static conditions and without a weak interlayer. With a weak interlayer, slopes are stable at <18° in the drained case and at <9° in the undrained case. Earthquake loading can drastically reduce the shear strength of sediment with increased pore water pressure. The slope became unstable at >13° with earthquake peak ground acceleration (PGA) of 0.5 g; whereas with a weak layer, a PGA of 0.2 g could trigger instability at slopes >10°, and >3° for PGA of 0.5 g. The northern slope of the South China Sea is geomorphologically stable under static conditions. However, because of the possibility of high PGA at the eastern margin of the South China Sea, submarine slides are likely on the Taiwan Bank slope and eastern part of the Dongsha slope. Therefore, submarine slides recognized in seismic profiles on the Taiwan Bank slope would be triggered by an earthquake, the most important factor for triggering submarine slides on the northern slope of the South China Sea. Considering the distribution of PGA, we consider the northern slope of the South China Sea to be stable, excluding the Taiwan Bank slope, which is tectonically active.

15. Numerical studies on convective stability and flow pattern in three-dimensional spherical mantle of terrestrial planets

Yanagisawa, Takatoshi; Kameyama, Masanori; Ogawa, Masaki

2016-06-01

We explore thermal convection of a fluid with a temperature-dependent viscosity in a basally heated three-dimensional spherical shell using linear stability analyses and numerical experiments, while considering the application of our results to terrestrial planets. The inner to outer radius ratio of the shell f assumed in the linear stability analyses is in the range of 0.11-0.88. The critical Rayleigh number Rc for the onset of thermal convection decreases by two orders of magnitude as f increases from 0.11 to 0.88, when the viscosity depends sensitively on the temperature, as is the case for real mantle materials. Numerical simulations carried out in the range of f = 0.11-0.55 show that a thermal boundary layer (TBL) develops both along the surface and bottom boundaries to induce cold and hot plumes, respectively, when f is 0.33 or larger. However, for smaller f values, a TBL develops only on the bottom boundary. Convection occurs in the stagnant-lid regime where the root mean square velocity on the surface boundary is less than 1% of its maximum at depth, when the ratio of the viscosity at the surface boundary to that at the bottom boundary exceeds a threshold that depends on f. The threshold decreases from 106.5 at f = 0.11 to 104 at f = 0.55. If the viscosity at the base of the convecting mantle is 1020-1021 Pa s, the Rayleigh number exceeds Rc for Mars, Venus and the Earth, but does not for the Moon and Mercury; convection is unlikely to occur in the latter planets unless the mantle viscosity is much lower than 1020 Pa s and/or the mantle contains a strong internal heat source.

16. Stability Assessment of the Gökgöl Karstic Cave (Zonguldak, Turkey) by Analytical and Numerical Methods

Geniş, Melih; Çolak, Berdan

2015-11-01

Turkey has many cave formations and some caves are important tourist attractions. Theses caves have natural and historical characteristics that need to be investigated and protected. An important rock engineering topic is the identification of possible rock falls and collapses in caves that are open to the public and planning preventive measures. In this study, the locations and types of instability problems were identified along a section of 875 m of the Gökgöl cave in Zonguldak, Turkey. The main instabilities were identified as block falls controlled by discontinuities, rock block sliding, and flexural toppling. To obtain the intact rock and rock mass properties of the stratified limestone which is the host rock of the cave, field and laboratory studies were conducted. Stability assessments for four different sections inside the cave were carried out using numerical and analytical methods, and the results were compared. The results of the analyses revealed that some countermeasures and precautions must be implemented to increase the stability in some areas of the Gökgöl cave, especially those which have potential for regional failure.

17. Numerical slope stability simulations of chasma walls in Valles Marineris/Mars using a distinct element method (dem).

Imre, B.

2003-04-01

NUMERICAL SLOPE STABILITY SIMULATIONS OF CHASMA WALLS IN VALLES MARINERIS/MARS USING A DISTINCT ELEMENT METHOD (DEM). B. Imre (1) (1) German Aerospace Center, Berlin Adlershof, bernd.imre@gmx.net The 8- to 10-km depths of Valles Marineris (VM) offer excellent views into the upper Martian crust. Layering, fracturing, lithology, stratigraphy and the content of volatiles have influenced the evolution of the Valles Marineris wallslopes. But these parameters also reflect the development of VM and its wall slopes. The scope of this work is to gain understanding in these parameters by back-simulating the development of wall slopes. For that purpose, the two dimensional Particle Flow Code PFC2D has been chosen (ITASCA, version 2.00-103). PFC2D is a distinct element code for numerical modelling of movements and interactions of assemblies of arbitrarily sized circular particles. Particles may be bonded together to represent a solid material. Movements of particles are unlimited. That is of importance because results of open systems with numerous unknown variables are non-unique and therefore highly path dependent. This DEM allows the simulation of whole development paths of VM walls what makes confirmation of the model more complete (e.g. Oreskes et al., Science 263, 1994). To reduce the number of unknown variables a proper (that means as simple as possible) field-site had to be selected. The northern wall of eastern Candor Chasma has been chosen. This wall is up to 8-km high and represents a significant outcrop of the upper Martian crust. It is quite uncomplex, well-aligned and of simple morphology. Currently the work on the model is at the stage of performing the parameter study. Results will be presented via poster by the EGS-Meeting.

18. A Combined Remote Sensing-Numerical Modelling Approach to the Stability Analysis of Delabole Slate Quarry, Cornwall, UK

Havaej, Mohsen; Coggan, John; Stead, Doug; Elmo, Davide

2016-04-01

Rock slope geometry and discontinuity properties are among the most important factors in realistic rock slope analysis yet they are often oversimplified in numerical simulations. This is primarily due to the difficulties in obtaining accurate structural and geometrical data as well as the stochastic representation of discontinuities. Recent improvements in both digital data acquisition and incorporation of discrete fracture network data into numerical modelling software have provided better tools to capture rock mass characteristics, slope geometries and digital terrain models allowing more effective modelling of rock slopes. Advantages of using improved data acquisition technology include safer and faster data collection, greater areal coverage, and accurate data geo-referencing far exceed limitations due to orientation bias and occlusion. A key benefit of a detailed point cloud dataset is the ability to measure and evaluate discontinuity characteristics such as orientation, spacing/intensity and persistence. This data can be used to develop a discrete fracture network which can be imported into the numerical simulations to study the influence of the stochastic nature of the discontinuities on the failure mechanism. We demonstrate the application of digital terrestrial photogrammetry in discontinuity characterization and distinct element simulations within a slate quarry. An accurately geo-referenced photogrammetry model is used to derive the slope geometry and to characterize geological structures. We first show how a discontinuity dataset, obtained from a photogrammetry model can be used to characterize discontinuities and to develop discrete fracture networks. A deterministic three-dimensional distinct element model is then used to investigate the effect of some key input parameters (friction angle, spacing and persistence) on the stability of the quarry slope model. Finally, adopting a stochastic approach, discrete fracture networks are used as input for 3D

19. Steric repulsion as a way to achieve the required stability for the preparation of ionic liquid-based ferrofluids.

PubMed

Rodríguez-Arco, Laura; López-López, Modesto T; González-Caballero, Fernando; Durán, Juan D G

2011-05-01

With this work we would like to emphasize the necessity of steric repulsion to stabilize novel ionic liquid-based ferrofluids. For this purpose, we prepared a suspension of magnetite nanoparticles coated with a double layer of oleic acid, dispersed in 1-ethyl-3-methylimidazolium ethylsulphate ([EMIM][EtSO(4)]). For comparison, a suspension of bare magnetite nanoparticles in [EMIM][EtSO(4)] was also prepared. The stability of these suspensions was checked by magnetic sedimentation and centrifugation processes. Furthermore, their yield stress was measured as a function of the applied magnetic field, which gave additional information on their stability. The results of these experiments showed that the suspension of bare nanoparticles was rather unstable, whereas the suspension of double layer coated nanoparticles gave rise to a true (stable) ferrofluid. PMID:21345446

20. A Longitudinal Study of Reading Comprehension Achievement from Grades 3 to 10: Investigating Models of Stability, Cumulative Growth, and Compensation

ERIC Educational Resources Information Center

Kwiatkowska-White, Bozena; Kirby, John R.; Lee, Elizabeth A.

2016-01-01

This longitudinal study of 78 Canadian English-speaking students examined the applicability of the stability, cumulative, and compensatory models in reading comprehension development. Archival government-mandated assessments of reading comprehension at Grades 3, 6, and 10, and the Canadian Test of Basic Skills measure of reading comprehension…

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

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.

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

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.

3. Spatial and temporal variations of glacial erosion in the European Alps: numerical models and implications for slope stability (Invited)

Sternai, P.; Herman, F.; Willett, S.; Champagnac, J.; Fox, M.; Valla, P.; Salcher, B.

2013-12-01

Glacial erosion in alpine landscapes can be highly variable in space and time and lead to significant morphologic modification and mass redistribution at virtually all scales. Because they affect the near-surface stress and strain distribution by producing cyclic variations of the surface load, removing and abrading rocks, storing/releasing sediments and affecting the surface and subsurface hydrology, glaciations have multiple effects on slope stability. Understanding how glacial erosion evolves in space and time is thus important for investigating potential feedbacks between glacial erosion and deep-seated gravitational slope deformation (DSGSD). The present-day topography of the European Alps shows evidence of intense glacial erosion. However, significant questions regarding Alpine landscape evolution during glaciations still persist. For example, large-scale topographic analyses suggest that glacial erosion is maximized at and above the glaciers' long-term Equilibrium Line Altitude. In contrast, measurements of long-term denudation rates from low-temperature thermochronology and reconstructions of the pre-glacial Alpine topography suggest high erosion towards low altitudes and formation of overdeepnenings, in turn indicating an increase of local relief in response to glacial processes. Based on sediment record, low-temperature thermochronology and burial cosmogenic nuclide dating, it has also been proposed that the mid-Pleistocene climatic transition from symmetric, 40kyr to asymmetric, 100kyr glacial/interglacial oscillations sets the onset of intense glacial erosion within the Alps. However, this climate threshold in glacial erosion has not been showed in other orogens, and positive feedbacks between climate periodicity and glacial erosion efficiency still remain to be proven. Numerical modeling provides estimates of the patterns and magnitudes of glacial erosion through time. Modeling results on an advanced reconstruction of the pre-glacial topography and the

4. Prediction of soil stability and erosion in semiarid regions using numerical hydrological model (MCAT) and airborne hyperspectral imagery

Brook, Anna; Wittenberg, Lea

2015-04-01

promising models is the MCAT, which is a MATLAB library of visual and numerical analysis tools for the evaluation of hydrological and environmental models. The model applied in this paper presents an innovative infrastructural system for predicting soil stability and erosion impacts. This integrated model is applicable to mixed areas with spatially varying soil properties, landscape, and land-cover characteristics. Data from a semiarid site in southern Israel was used to evaluate the model and analyze fundamental erosion mechanisms. The findings estimate the sensitivity of the suggested model to the physical parameters and encourage the use of hyperspectral remote sensing imagery (HSI). The proposed model is integrated according to the following stages: 1. The soil texture, aggregation, soil moisture estimated via airborne HSI data, including soil surface clay and calcium carbonate erosions; 2. The mechanical stability of soil assessed via pedo-transfer function corresponding to load dependent changes in soil physical properties due to pre-compression stress (set of equations study shear strength parameters take into account soil texture, aggregation, soil moisture and ecological soil variables); 3. The precipitation-related runoff model program (RMP) satisfactorily reproduces the observed seasonal mean and variation of surface runoff for the current climate simulation; 4. The Monte Carlo Analysis Toolbox (MCAT), a library of visual and numerical analysis tools for the evaluation of hydrological and environmental models, is proposed as a tool for integrate all the approaches to an applicable model. The presented model overcomes the limitations of existing modeling methods by integrating physical data produced via HSI and yet stays generic in terms of space and time independency.

5. Item Calibration Samples and the Stability of Achievement Estimates and System Rankings: Another Look at the PISA Model

ERIC Educational Resources Information Center

Rutkowski, Leslie; Rutkowski, David; Zhou, Yan

2016-01-01

Using an empirically-based simulation study, we show that typically used methods of choosing an item calibration sample have significant impacts on achievement bias and system rankings. We examine whether recent PISA accommodations, especially for lower performing participants, can mitigate some of this bias. Our findings indicate that standard…

6. The Stability of Teacher Effects on Student Math and Reading Achievement over Time: A Study of One Midsize District

ERIC Educational Resources Information Center

Bessolo, Jennifer

2013-01-01

Increased attention on teacher accountability via student achievement gains has brought proposed policy changes which recommend teachers receive performance pay as recognition for effective teaching. In response to growing consensus that teachers' contributions to student learning should be a part of the teacher evaluation process, education…

7. Influence of the fluoride ion in the stability of bone hydroxyapatite achieved in patients after the fluo-calcic therapy.

PubMed

Mongiorgi, R; Bertocchi, G; Gnudi, S; Maggi, G; Moroni, A; Coppola, G; Riva di Sanseverino, L

1991-04-01

It has been proposed that fluoride ion might actively enter into the basic structural framework of the bone, within its main constituent, the hydroxyapatite (HA), by substituting systematically the OH- ion in the crystal chemical formula Ca5[(PO4)3OH]. Accurate compositional studies have become necessary in order to explain an eventual stabilizing effect. Loss in weight (TG) and chemical reaction (DTG) while varying the temperature have been carried out for the first time on specimens under normal conditions and after fluo-calcic treatment, in parallel with accurate chemical compositional determinations by atomic absorption analysis. Our investigation shows that the new structure present after treatment is mechanically stable and proves more resistant to osteolytic processes. PMID:1910745

8. Using interlayer step-wise triplet transfer to achieve an efficient white organic light-emitting diode with high color-stability

SciTech Connect

Wang, Qi; Ma, Dongge Ding, Junqiao; Wang, Lixiang; Leo, Karl; Qiao, Qiquan; Jia, Huiping; Gnade, Bruce E.

2014-05-12

An efficient phosphorescent white organic light emitting-diode with a red-green-blue tri-emitting-layer structure is reported. The host of the red dopant possesses a lower triplet-energy than the green dye. An interlayer step-wise triplet transfer via blue dye → green dye → red host → red dye is achieved. This mechanism allows an efficient triplet harvesting by the three dopants, thus maintaining a balanced white light and reducing energy loss. Moreover, the color stability of the device is improved significantly. The white device not only achieves a peak external quantum efficiency of 21.1 ± 0.8% and power efficiency of 37.5 ± 1.4 lm/W but shows no color shift over a wide range of voltages.

9. Numerical stability of an explicit finite difference scheme for the solution of transient conduction in composite media

NASA Technical Reports Server (NTRS)

Campbell, W.

1981-01-01

A theoretical evaluation of the stability of an explicit finite difference solution of the transient temperature field in a composite medium is presented. The grid points of the field are assumed uniformly spaced, and media interfaces are either vertical or horizontal and pass through grid points. In addition, perfect contact between different media (infinite interfacial conductance) is assumed. A finite difference form of the conduction equation is not valid at media interfaces; therefore, heat balance forms are derived. These equations were subjected to stability analysis, and a computer graphics code was developed that permitted determination of a maximum time step for a given grid spacing.

10. Estimation of Geomechanical Deformation and Stability Change Due to Geologic Storage of Carbon Dioxide Using a Thermo-Hydro-Mechanical Numerical Model

Park, S.; Kihm, J.; Kim, J.; SNU CO2 GEO-SEQ TEAM

2011-12-01

The purpose of this study is to quantitatively estimate groundwater and carbon dioxide flow and geomechanical deformation and stability changes due to geologic storage of carbon dioxide. A hypothetical layered aquifer system with alternation of sandstone and shale layers is set up considering geologic characteristics of the Sindong Group in the Gyeongsang Basin, Korea. A series of numerical simulations is then performed for 1,000 years of geologic storage of carbon dioxide including 10 years of carbon dioxide injection period with various injection rates using a multiphase thermo-hydro-mechanical numerical model. The numerical simulation results show that geomechanical deformation and stability change occur as the fluid pressure builds up mainly near the screen interval of the injection well, and the injection rate of carbon dioxide significantly influences geomechanical deformation and stability change. For the base case with the injection rate of 0.2 Mton/year, the safety factor with respect to tensile failure is reduced to about 2.0 near the screen interval of the injection well with a sudden increase of the fluid pressure in the early stage of the carbon dioxide injection period. Such geomechanical deformation and instability near the screen interval of the injection well are then reduced as the fluid pressure build-up propagates and diminishes outward from the injection well. The vertical displacement on the ground surface right above the injection well is maximal, and it is reduced as the horizontal distance from the injection well increases. In addition, the maximum vertical displacement on the ground surface right above the injection well increases almost linearly as the injection rate of carbon dioxide increases. The results of this study suggest that integrated thermo-hydro-mechanical numerical simulation is required when the stability of injection wells and surface injection facilities has to be ensured for geologic storage of carbon dioxide. This work

11. Numerical investigation of CO{sub 2} emission and thermal stability of a convective and radiative stockpile of reactive material in a cylindrical pipe of variable thermal conductivity

SciTech Connect

Lebelo, Ramoshweu Solomon

2014-10-24

In this paper the CO{sub 2} emission and thermal stability in a long cylindrical pipe of combustible reactive material with variable thermal conductivity are investigated. It is assumed that the cylindrical pipe loses heat by both convection and radiation at the surface. The nonlinear differential equations governing the problem are tackled numerically using Runge-Kutta-Fehlberg method coupled with shooting technique method. The effects of various thermophysical parameters on the temperature and carbon dioxide fields, together with critical conditions for thermal ignition are illustrated and discussed quantitatively.

12. A tabulated chemistry approach for numerical modeling of diesel spray evaporation in a 'stabilized cool flame' environment

SciTech Connect

Kolaitis, D.I.; Founti, M.A.

2006-04-15

Droplet evaporation in a 'stabilized cool flame' environment leads to a homogeneous, heated air-fuel vapor mixture that can be subsequently either burnt or utilized in fuel-reforming applications for fuel cell systems. The paper investigates the locally occurring physico-chemical phenomena in an atmospheric pressure, diesel spray, stabilized cool flame reactor, utilizing a tabulated chemistry approach in conjunction with a two-phase, Eulerian-Lagrangian computational fluid dynamics code. Actual diesel oil physical properties are used to model spray evaporation in the two-phase simulations, whereas the corresponding chemistry is represented by n-heptane. A lookup table is constructed by performing a plethora of perfectly stirred reactor simulations, utilizing a semidetailed n-heptane oxidation chemical kinetics mechanism. The overall exothermicity of the preignition n-heptane oxidation chemistry and the fuel consumption rates are examined as a function of selected independent parameters, namely temperature, fuel concentration, and residence time; their influence on cool flame reactivity is thoroughly studied. It is shown that the tabulated chemistry approach allows accurate investigation of the chemical phenomena with low computational cost. The two-phase flow inside the stabilized cool flame reactor is simulated, utilizing the developed lookup table. Predictions are presented for a variety of test cases and are compared to available experimental data, with satisfactory agreement. Model validation tests indicate that prediction quality improves with increasing values of air temperature at the reactor's inlet. (author)

13. Mathematical modeling and numerical simulation of an actively stabilized beam-column with circular cross-section

Schaeffner, Maximilian; Enss, Georg C.; Platz, Roland

2014-04-01

Buckling of axially loaded beam-columns represents a critical design constraint for light-weight structures. Besides passive solutions to increase the critical buckling load, active buckling control provides a possibility to stabilize slender elements in structures. So far, buckling control by active forces or bending moments has been mostly investigated for beam-columns with rectangular cross-section and with a preferred direction of buckling. The proposed approach investigates active buckling control of a beam-column with circular solid cross-section which is fixed at its base and pinned at its upper end. Three controlled active lateral forces are applied near the fixed base with angles of 120° to each other to stabilize the beam-column and allow higher critical axial loads. The beam-column is subject to supercritical static axial loads and lateral disturbance forces with varying directions and offsets. Two independent modal state space systems are derived for the bending planes in the lateral y- and z-directions of the circular cross-section. These are used to design two linear-quadratic regulators (LQR) that determine the necessary control forces which are transformed into the directions of the active lateral forces. The system behavior is simulated with a finite element model using one-dimensional beam elements with six degrees of freedom at each node. With the implemented control, it is possible to actively stabilize a beam-column with circular cross-section in arbitrary buckling direction for axial loads significantly above the critical axial buckling load.

14. A predictive numerical model for potential mapping of the gas hydrate stability zone in the Gulf of Cadiz

Leon, R.; Somoza, L.

2009-04-01

This comunication presents a computational model for mapping the regional 3D distribution in which seafloor gas hydrates would be stable, that is carried out in a Geographical Information System (GIS) environment. The construction of the model is comprised of three primary steps, namely (1) the construction of surfaces for the various variables based on available 3D data (seafloor temperature, geothermal gradient and depth-pressure); (2) the calculation of the gas function equilibrium functions for the various hydrocarbon compositions reported from hydrate and sediment samples; and (3) the calculation of the thickness of the hydrate stability zone. The solution is based on a transcendental function, which is solved iteratively in a GIS environment. The model has been applied in the northernmost continental slope of the Gulf of Cadiz, an area where an abundant supply for hydrate formation, such as extensive hydrocarbon seeps, diapirs and fault structures, is combined with deep undercurrents and a complex seafloor morphology. In the Gulf of Cadiz, model depicts the distribution of the base of the gas hydrate stability zone for both biogenic and thermogenic gas compositions, and explains the geometry and distribution of geological structures derived from gas venting in the Tasyo Field (Gulf of Cadiz) and the generation of BSR levels on the upper continental slope.

15. Evaluation of Rainfall Impacts on Groundwater Flow and Land Deformation in an Unsaturated Heterogeneous Slope and Slope Stability Using a Fully Coupled Hydrogeomechanical Numerical Model

Kihm, J.; Kim, J.

2006-12-01

A series of numerical simulations using a fully coupled hydrogeomechanical numerical model, which is named COWADE123D, is performed to analyze groundwater flow and land deformation in an unsaturated heterogeneous slope and its stability under various rainfall rates. The slope is located along a dam lake in Republic of Korea. The slope consists of the Cretaceous granodiorite and can be subdivided into the four layers such as weathered soil, weathered rock, intermediate rock, and hard rock from its ground surface due to weathering process. The numerical simulation results show that both rainfall rate and heterogeneity play important roles in controlling groundwater flow and land deformation in the unsaturated slope. The slope becomes more saturated, and thus its overall hydrogeomechanical stability deteriorates, especially in the weathered rock and weathered soil layers, as the rainfall increases up to the maximum daily rainfall rate in the return period of one year. However, the slope becomes fully saturated, and thus its hydrogeomechanical responses are almost identical under more than such a critical rainfall rate. From the viewpoint of hydrogeology, the pressure head, and hence the hydraulic head increase as the rainfall rate increases. As a result, the groundwater table rises, the unsaturated zone reduces, the seepage face expands from the slope toe toward the slope crest, and the groundwater flow velocity increases along the seepage face. Particularly, the groundwater flow velocity increases significantly in the weathered soil and weathered rock layers as the rainfall rate increases. This is because their hydraulic conductivity is relatively higher than that of the intermediate rock and hard rock layers. From the viewpoint of geomechanics, the horizontal displacement increases, while the vertical displacement decreases toward the slope toe as the rainfall rate increases. This may result from the buoyancy effect associated with the groundwater table rise as the

16. Numerical analysis of stress distribution in Cu-stabilized GdBCO CC tapes during anvil tests for the evaluation of transverse delamination strength

Dizon, John Ryan C.; Gorospe, Alking B.; Shin, Hyung-Seop

2014-05-01

Rare-earth-Ba-Cu-O (REBCO) based coated conductors (CCs) are now being used for electric device applications. For coil-based applications such as motors, generators and magnets, the CC tape needs to have robust mechanical strength along both the longitudinal and transverse directions. The CC tape in these coils is subjected to transverse tensile stresses during cool-down and operation, which results in delamination within and between constituent layers. In this study, in order to explain the behaviour observed in the evaluation of c-axis delamination strength in Cu-stabilized GdBCO CC tapes by anvil tests, numerical analysis of the mechanical stress distribution within the CC tape has been performed. The upper anvil size was varied in the analysis to understand the effect of anvil size on stress distribution within the multilayered CC tape, which is closely related to the delamination strength, delamination mode and delamination sites that were experimentally observed. The numerical simulation results showed that, when an anvil size covering the whole tape width was used, the REBCO coating film was subjected to the largest stress, which could result in low mechanical delamination and electromechanical delamination strengths. Meanwhile, when smaller-sized anvils were used, the copper stabilizer layer would experience the largest stress among all the constituent layers of the CC tape, which could result in higher mechanical and electromechanical delamination strengths, as well as high scattering of both of these delamination strengths. As a whole, the numerical simulation results could explain the damage evolution observed in CC tapes tested under transverse tensile stress, as well as the transverse tensile stress response of the critical current, Ic.

17. Three-dimensional direct numerical simulation of a turbulent lifted hydrogen jet flame in a heated coflow: flame stabilization and structure

SciTech Connect

Chen, Jackie; Sankaran, Ramanan; Yoo, Chun S

2009-01-01

Direct numerical simulation (DNS) of the near field of a three-dimensional spatially developing turbulent lifted hydrogen jet flame in heated coflow is performed with a detailed mechanism to determine the stabilization mechanism and the flame structure. The DNS was performed at a jet Reynolds number of 11,000 with over 940 million grid points. The results show that auto-ignition in a fuel-lean mixture at the flame base is the main source of stabilization of the lifted jet flame. A chemical flux analysis shows the occurrence of near-isothermal chemical chain branching preceding thermal runaway upstream of the stabilization point, indicative of hydrogen auto-ignition in the second limit. The Damkoehler number and key intermediate-species behaviour near the leading edge of the lifted flame also verify that auto-ignition occurs at the flame base. At the lifted-flame base, it is found that heat release occurs predominantly through ignition in which the gradients of reactants are opposed. Downstream of the flame base, both rich-premixed and non-premixed flames develop and coexist with auto-ignition. In addition to auto-ignition, Lagrangian tracking of the flame base reveals the passage of large-scale flow structures and their correlation with the fluctuations of the flame base. In particular, the relative position of the flame base and the coherent flow structure induces a cyclic motion of the flame base in the transverse and axial directions about a mean lift-off height. This is confirmed by Lagrangian tracking of key scalars, heat release rate and velocity at the stabilization point.

18. Numerical studies of nitric oxide formation in nanosecond-pulsed discharge-stabilized flames of premixed methane/air.

PubMed

Bak, Moon Soo; Cappelli, Mark A

2015-08-13

A simulation is developed to investigate the kinetics of nitric oxide (NO) formation in premixed methane/air combustion stabilized by nanosecond-pulsed discharges. The simulation consists of two connected parts. The first part calculates the kinetics within the discharge while considering both plasma/combustion reactions and species diffusion, advection and thermal conduction to the surrounding flow. The second part calculates the kinetics of the overall flow after mixing the discharge flow with the surrounding flow to account for the effect that the discharge has on the overall kinetics. The simulation reveals that the discharge produces a significant amount of atomic oxygen (O) as a result of the high discharge temperature and dissociative quenching of excited state nitrogen by molecular oxygen. This atomic oxygen subsequently produces hydroxyl (OH) radicals. The fractions of these O and OH then undergo Zel'dovich reactions and are found to contribute to as much as 73% of the total NO that is produced. The post-discharge simulation shows that the NO survives within the flow once produced. PMID:26170428

19. Spreading of a surfactant-laden droplet down an inclined and pre-wetted plane - Numerics, Asymptotics and Linear Stability Analysis

Goddard, Joseph

2012-11-01

Non-uniformities in surfactant concentration result in a surface shear stress, known as the Marangoni stress. This stress, if sufficiently large, can influence the flow at the interface. Naturally occurring surfactants in the mammalian lung reduce the surface tension within the liquid lining the airways and help to prevent collapse of smaller airways. Premature infants produce insufficient surfactant because the lungs are under developed. Resulting Respiratory Distress Syndrome is treated by Surfactant Replacement Therapy. Motivated by this medical application we theoretically investigate a model problem involving the spreading of a drop laden with an insoluble surfactant down an inclined and pre-wetted plane. Our focus is on understanding the mechanisms behind a ``fingering'' instability observed experimentally by high-resolution numerics revealing a multi-region asymptotic structure of the spreading droplet. Approximate solutions for each region are then derived using asymptotic analysis. In particular, a quasi-steady similarity solution is obtained for the leading edge of the droplet and a linear stability analysis shows that the base state is linearly unstable to long-wavelength perturbations for all inclination angles. The Marangoni effect is shown to be behind this instability at small inclination angles. A stability criterion is derived at small wave number and it's implication in the onset of the instability will be discussed.

20. Validation of three-dimensional incompressible spatial direct numerical simulation code: A comparison with linear stability and parabolic stability equation theories for boundary-layer transition on a flat plate

NASA Technical Reports Server (NTRS)

Joslin, Ronald D.; Streett, Craig L.; Chang, Chau-Lyan

1992-01-01

Spatially evolving instabilities in a boundary layer on a flat plate are computed by direct numerical simulation (DNS) of the incompressible Navier-Stokes equations. In a truncated physical domain, a nonstaggered mesh is used for the grid. A Chebyshev-collocation method is used normal to the wall; finite difference and compact difference methods are used in the streamwise direction; and a Fourier series is used in the spanwise direction. For time stepping, implicit Crank-Nicolson and explicit Runge-Kutta schemes are used to the time-splitting method. The influence-matrix technique is used to solve the pressure equation. At the outflow boundary, the buffer-domain technique is used to prevent convective wave reflection or upstream propagation of information from the boundary. Results of the DNS are compared with those from both linear stability theory (LST) and parabolized stability equation (PSE) theory. Computed disturbance amplitudes and phases are in very good agreement with those of LST (for small inflow disturbance amplitudes). A measure of the sensitivity of the inflow condition is demonstrated with both LST and PSE theory used to approximate inflows. Although the DNS numerics are very different than those of PSE theory, the results are in good agreement. A small discrepancy in the results that does occur is likely a result of the variation in PSE boundary condition treatment in the far field. Finally, a small-amplitude wave triad is forced at the inflow, and simulation results are compared with those of LST. Again, very good agreement is found between DNS and LST results for the 3-D simulations, the implication being that the disturbance amplitudes are sufficiently small that nonlinear interactions are negligible.

1. Symbolic generation of elastic rotor blade equations using a FORTRAN processor and numerical study on dynamic inflow effects on the stability of helicopter rotors

NASA Technical Reports Server (NTRS)

Reddy, T. S. R.

1986-01-01

2. 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.

3. Slope stability and stabilization methods

SciTech Connect

Abramson, L.W.; Lee, T.S.; Boyce, G.M.; Sharma, S.S.

1995-12-01

Slope stability can be a major problem during the construction of surface facilities. Cutting into existing ground disturbs the mechanics of the surrounding area, which can result in landslides and rock falls. This practical reference gives you the comprehensive information you need for slope stability analysis, suitable methods of analysis with and without the use of computers, and examples of common stability problems and stabilization methods for cuts and fills. It includes detailed discussions of methods used in slope stability analysis, including the Ordinary Method of Slices, Simplified Janbu Method, Simplified Bishop Method, Spencer`s Method, other limit equilibrium methods, numerical methods, total stress analysis, effective stress analysis, and the use of computer programs to solve problems. Chapters include: General Slope Stability Concepts; Engineering Geology Principles; Groundwater Conditions; Geologic Site Exploration; Laboratory Testing Interpretation; Slope Stability Concepts; Slope Stabilization Methods; and Design, Construction and Maintenance.

4. Wechsler Individual Achievement Test.

ERIC Educational Resources Information Center

Taylor, Ronald L.

1999-01-01

This article describes the Wechsler Individual Achievement Test, a comprehensive measure of achievement for individuals in grades K-12. Eight subtests assess mathematics reasoning, spelling, reading comprehension, numerical operations, listening comprehension, oral expression, and written expression. Its administration, standardization,…

5. Numerical simulation of the insertion process of an uncemented hip prosthesis in order to evaluate the influence of residual stress and contact distribution on the stem initial stability.

PubMed

Monea, Aida Georgeta; Pastrav, Leonard Cezar; Mulier, Michiel; Van der Perre, Georges; Jaecques, Siegfried V

2014-01-01

The long-term success of a cementless total hip arthroplasty depends on the implant geometry and interface bonding characteristics (fit, coating and ingrowth) and on stem stiffness. This study evaluates the influence of stem geometry and fitting conditions on the evolution and distribution of the bone-stem contact, stress and strain during and after the hip stem insertion, by means of dynamic finite element techniques. Next, the influence of the mechanical state (bone-stem contact, stress and strain) resulted from the insertion process on the stem initial resistance to subsidence is investigated. In addition, a study on the influence of bone-stem interface conditions (friction) on the insertion process and on the initial stem stability under physiological loading is performed. The results indicate that for a stem with tapered shape the contact in the proximal part of the stem was improved, but contact in the calcar region was achieved only when extra press-fit conditions were considered. Changes in stem geometry towards a more tapered shape and extra press fit and variation in the bone-stem interface conditions (contact amount and high friction) led to a raise in the total insertion force. A direct positive relationship was found between the stem resistance to subsidence and stem geometry (tapering and press fit), bone-stem interface conditions (bone-stem contact and friction interface) and the mechanical status at the end of the insertion (residual stress and strain). Therefore, further studies on evaluating the initial performance of different stem types should consider the parameters describing the bone-stem interface conditions and the mechanical state resulted from the insertion process. PMID:22548243

6. Achieving extremely concentrated aqueous dispersions of graphene flakes and catalytically efficient graphene-metal nanoparticle hybrids with flavin mononucleotide as a high-performance stabilizer.

PubMed

Ayán-Varela, M; Paredes, J I; Guardia, L; Villar-Rodil, S; Munuera, J M; Díaz-González, M; Fernández-Sánchez, C; Martínez-Alonso, A; Tascón, J M D

2015-05-20

The stable dispersion of graphene flakes in an aqueous medium is highly desirable for the development of materials based on this two-dimensional carbon structure, but current production protocols that make use of a number of surfactants typically suffer from limitations regarding graphene concentration or the amount of surfactant required to colloidally stabilize the sheets. Here, we demonstrate that an innocuous and readily available derivative of vitamin B2, namely the sodium salt of flavin mononucleotide (FMNS), is a highly efficient dispersant in the preparation of aqueous dispersions of defect-free, few-layer graphene flakes. Most notably, graphene concentrations in water as high as ∼50 mg mL(-1) using low amounts of FMNS (FMNS/graphene mass ratios of about 0.04) could be attained, which facilitated the formation of free-standing graphene films displaying high electrical conductivity (∼52000 S m(-1)) without the need of carrying out thermal annealing or other types of post-treatment. The excellent performance of FMNS as a graphene dispersant could be attributed to the combined effect of strong adsorption on the sheets through the isoalloxazine moiety of the molecule and efficient colloidal stabilization provided by its negatively charged phosphate group. The FMNS-stabilized graphene sheets could be decorated with nanoparticles of several noble metals (Ag, Pd, and Pt), and the resulting hybrids exhibited a high catalytic activity in the reduction of nitroarenes and electroreduction of oxygen. Overall, the present results should expedite the processing and implementation of graphene in, e.g., conductive inks, composites, and hybrid materials with practical utility in a wide range of applications. PMID:25915172

7. Deposition of composite materials using a wire-bar coater for achieving processability and air-stability in Organic Field-Effect Transistors (OFETs)

Del Pozo, Freddy G.; Galindo, Sergi; Pfattner, Raphael; Rovira, Concepció; Mas-Torrent, Marta

2015-08-01

Organic thin films based on composite materials of semiconducting dibenzo-tetrathiafulvalene (DB-TTF) and insulating styrenic matrices (Polystyrene (PS10k) and Poly-alpha methylstyrene (PAMS10k) ) have been fabricated by the wire-bar coating technique in ambient conditions (air, light, humidity) and contrasted with the ones prepared by thermally evaporating the organic semiconductor. The transistors fabricated with DB-TTF:PS10k composites show a clear fieldeffect behavior with p-type characteristics, exhibiting charge carriers mobilities in the range of 0.01 cm2/Vs, fully comparable with the films obtained by thermal evaporation. However, while the thermally evaporated films show poor stability in air, the wire-bar coated composites films and devices are highly reproducible and exhibit lower threshold voltage values. Thus, we demonstrate the suitability of the wire-bar technique for manufacturing large area devices.

8. Improved electrochemical stability at the surface of La(0.8)Sr(0.2)CoO3 achieved by surface chemical modification.

PubMed

Tsvetkov, Nikolai; Lu, Qiyang; Yildiz, Bilge

2015-01-01

The degradation of the surface chemistry on perovskite (ABO3) oxides is a critical issue for their performance in energy conversion systems such as solid oxide fuel/electrolysis cells and in splitting of H2O and CO2 to produce fuels. This degradation is typically in the form of segregation and phase separation of dopant cations from the A-site, driven by elastic and electrostatic energy minimization and kinetic demixing. In this study, deposition of Ti at the surface was found to hinder the dopant segregation and the corresponding electrochemical degradation on a promising SOFC cathode material, La(0.8)Sr(0.2)CoO3 (LSC). The surface of the LSC films was modified by Ti (denoted as LSC-T) deposited from a TiCl4 solution. The LSC and LSC-T thin films were investigated by electrochemical impedance spectroscopy, nano-probe Auger electron spectroscopy, and X-ray photoelectron spectroscopy (XPS), upon annealing at 420-530 °C in air up to about 90 hours. The oxygen exchange coefficient, k(q), on LSC-T cathodes was found to be up to 8 times higher than that on LSC cathodes at 530 °C and retained its stability. Sr-rich insulating particles formed at the surface of the annealed LSC and LSC-T films, but with significantly less coverage of such particles on the LSC-T. From this result, it appears that modification of the LSC surface with Ti reduces the segregation of the blocking Sr-rich particles at the surface, and a larger area on LSC surface (with a higher Sr doping level in the lattice) is available for the oxygen reduction reaction. The stabilization of the LSC surface through Ti-deposition can open a new route for designing surface modifications on perovskite oxide electrodes for high temperature electro- and thermo-chemical applications. PMID:26227310

9. Stability of Bareiss algorithm

Bojanczyk, Adam W.; Brent, Richard P.; de Hoog, F. R.

1991-12-01

In this paper, we present a numerical stability analysis of Bareiss algorithm for solving a symmetric positive definite Toeplitz system of linear equations. We also compare Bareiss algorithm with Levinson algorithm and conclude that the former has superior numerical properties.

10. Numerical Stability and Control Analysis Towards Falling-Leaf Prediction Capabilities of Splitflow for Two Generic High-Performance Aircraft Models

NASA Technical Reports Server (NTRS)

Charlton, Eric F.

1998-01-01

Aerodynamic analysis are performed using the Lockheed-Martin Tactical Aircraft Systems (LMTAS) Splitflow computational fluid dynamics code to investigate the computational prediction capabilities for vortex-dominated flow fields of two different tailless aircraft models at large angles of attack and sideslip. These computations are performed with the goal of providing useful stability and control data to designers of high performance aircraft. Appropriate metrics for accuracy, time, and ease of use are determined in consultations with both the LMTAS Advanced Design and Stability and Control groups. Results are obtained and compared to wind-tunnel data for all six components of forces and moments. Moment data is combined to form a "falling leaf" stability analysis. Finally, a handful of viscous simulations were also performed to further investigate nonlinearities and possible viscous effects in the differences between the accumulated inviscid computational and experimental data.

11. Automatic stabilization

NASA Technical Reports Server (NTRS)

Haus, FR

1936-01-01

This report concerns the study of automatic stabilizers and extends it to include the control of the three-control system of the airplane instead of just altitude control. Some of the topics discussed include lateral disturbed motion, static stability, the mathematical theory of lateral motion, and large angles of incidence. Various mechanisms and stabilizers are also discussed. The feeding of Diesel engines by injection pumps actuated by engine compression, achieves the required high speeds of injection readily and permits rigorous control of the combustible charge introduced into each cylinder and of the peak pressure in the resultant cycle.

12. Numerical Study of Velocity Shear Stabilization of 3D and Theoretical Considerations for Centrifugally Confined Plasmas and Other Interchange-Limited Fusion Concepts

SciTech Connect

2015-09-21

We studied the feasibility of resonantly driving GAMs in tokamaks. A numerical simulation was carried out and showed the essential features and limitations. It was shown further that GAMs can damp by phase-mixing, from temperature gradients, or nonlinear detuning, thus broadening the resonance. Experimental implications of this were quantified. Theoretical support was provided for the Maryland Centrifugal Experiment, funded in a separate grant by DOE. Plasma diamagnetism from supersonic rotation was established. A theoretical model was built to match the data. Additional support to the experiment in terms of numerical simulation of the interchange turbulence was provided. Spectra from residual turbulence on account of velocity shear suppression were obtained and compared favorably to experiment. A new drift wave, driven solely by the thermal force, was identified.

13. Terrestrial Photogrammetry and Numerical Modelling for the Stability Analysis of Rock Slopes in High Mountain Areas: Aiguilles Marbrées case

Curtaz, M.; Ferrero, A. M.; Roncella, R.; Segalini, A.; Umili, G.

2014-03-01

Several high-altitude slope instability phenomena, involving rock blocks of different volumes, have been observed in recent years. The increase in these phenomena could be correlated to climatic variations and to a general increase in temperature that has induced both ice melting with consequent water seepage and glacial lowering, with a consequent loss of support of the rock face. The degradation of the high-altitude thermal layer, which is known as "permafrost", can determine the formation of highly fractured rock slopes where instabilities can concentrate. The present research has developed a methodology to improve the understanding and assessment of rock slope stability conditions in high mountain environments where access is difficult. The observed instabilities are controlled by the presence of discontinuities that can determine block detachments. Consequently, a detailed survey of the rock faces is necessary, both in terms of topography and geological structure, and in order to locate the discontinuities on the slope to obtain a better geometric reconstruction and subsequent stability analysis of the blocky rock mass. Photogrammetric surveys performed at different times allow the geostructure of the rock mass to be determined and the rock block volumes and detachment mechanisms to be estimated, in order to assess the stability conditions and potential triggering mechanisms. Photogrammetric surveys facilitate both the characterisation of the rock mass and the monitoring of slope instabilities over time. The methodology has been applied in a case study pertaining to the North Face of Aiguilles Marbrées in the Mont Blanc massif, which suffers from frequent instability phenomena. A slope failure that occurred in 2007 has been back-analysed using both the limit equilibrium method (LEM) and 3D distinct element modelling (DEM). The method has been supported and validated with traditional in situ surveys and measurements of the discontinuity orientation and other

14. Experimental and Numerical Optimization of a High-Lift System to Improve Low-Speed Performance, Stability, and Control of an Arrow-Wing Supersonic Transport

NASA Technical Reports Server (NTRS)

Hahne, David E.; Glaab, Louis J.

1999-01-01

15. Numerical Dissipation and Wrong Propagation Speed of Discontinuities for Stiff Source Terms

NASA Technical Reports Server (NTRS)

Yee, H. C.; Kotov, D. V.; Sjogreen, B.

2011-01-01

In compressible turbulent combustion/nonequilibrium flows, the constructions of numerical schemes for (a) stable and accurate simulation of turbulence with strong shocks, and (b) obtaining correct propagation speed of discontinuities for stiff reacting terms on coarse grids share one important ingredient - minimization of numerical dissipation while maintaining numerical stability. Here coarse grids means standard mesh density requirement for accurate simulation of typical non-reacting flows. This dual requirement to achieve both numerical stability and accuracy with zero or minimal use of numerical dissipation is most often conflicting for existing schemes that were designed for non-reacting flows. The goal of this paper is to relate numerical dissipations that are inherited in a selected set of high order shock-capturing schemes with the onset of wrong propagation speed of discontinuities for two representative stiff detonation wave problems.

16. Numerical Dissipation and Wrong Propagation Speed of Discontinuities for Stiff Source Terms

NASA Technical Reports Server (NTRS)

Yee, H. C.; Kotov, D. V.; Sjoegreen, B.

2012-01-01

In compressible turbulent combustion/nonequilibrium flows, the constructions of numerical schemes for (a) stable and accurate simulation of turbulence with strong shocks, and (b) obtaining correct propagation speed of discontinuities for stiff reacting terms on coarse grids share one important ingredient - minimization of numerical dissipation while maintaining numerical stability. Here coarse grids means standard mesh density requirement for accurate simulation of typical non-reacting flows. This dual requirement to achieve both numerical stability and accuracy with zero or minimal use of numerical dissipation is most often conflicting for existing schemes that were designed for non-reacting flows. The goal of this paper is to relate numerical dissipations that are inherited in a selected set of high order shock-capturing schemes with the onset of wrong propagation speed of discontinuities as a function of stiffness of the source term and the grid spacing.

17. Stochastic-deterministic modeling of bed load transport in shallow water flow over erodible slope: Linear stability analysis and numerical simulation

Bohorquez, Patricio; Ancey, Christophe

2015-09-01

In this article we propose a stochastic bed load transport formulation within the framework of the frictional shallow-water equations in which the sediment transport rate results from the difference between the entrainment and deposition of particles. First we show that the Saint-Venant-Exner equations are linearly unstable in most cases for a uniform base flow down an inclined erodible bed for Shields numbers in excess of the threshold for incipient sediment motion allowing us to compute noise-induced pattern formation for Froude numbers below 2. The wavelength of the bed forms are selected naturally due to the absolute character of the bed instability and the existence of a maximum growth rate at a finite wavelength when the particle diffusivity coefficient and the water eddy viscosity are present as for Turing-like instability. A numerical method is subsequently developed to analyze the performance of the model and theoretical results through three examples: the simulation of the fluctuations of the particle concentration using a stochastic Langevin equation, the deterministic simulation of anti-dunes formation over an erodible slope in full sediment-mobility conditions, and the computation of noise-induced pattern formation in hybrid stochastic-deterministic flows down a periodic flume. The full non-linear numerical simulations are in excellent agreement with the theoretical solutions. We conclude highlighting that the proposed depth-averaged formulation explains the developments of upstream migrating anti-dunes in straight flumes since the seminar experiments by Gilbert (1914).

18. An updated numerical simulation of the ground-water flow system for the Castle Lake debris dam, Mount St. Helens, Washington, and implications for dam stability against heave

USGS Publications Warehouse

Roeloffs, Evelyn A.

1994-01-01

A numerical simulation of the ground-water flow system in the Castle Lake debris dam, calibrated to data from the 1991 and 1992 water years, was used to estimate factors of safety against heave and internal erosion. The Castle Lake debris dam, 5 miles northwest of the summit of Mount St. Helens, impounds 19,000 acre-ft of water that could pose a flood hazard in the event of a lake breakout. A new topographic map of the Castle Lake area prior to the 1980 eruption of Mount St. Helens was prepared and used to calculate the thickness of the debris avalanche deposits that compose the dam. Water levels in 22 piezometers and discharges from seeps on the dam face measured several times per year beginning in 1990 supplemented measurements in 11 piezometers and less frequent seep discharge measurements made since 1983. Observations in one group of piezometers reveal heads above the land surface and head gradients favoring upward flow that correspond to factors of safety only slightly greater than 2. The steady-state ground-water flow system in the debris dam was simulated using a threedimensional finite difference computer program. A uniform, isotropic model having the same shape as the dam and a hydraulic conductivity of 1.55 ft/day simulates the correct water level at half the observation points, but is in error by 10 ft or more at other points. Spatial variations of hydraulic conductivity were required to calibrate the model. The model analysis suggests that ground water flows in both directions between the debris dam and Castle Lake. Factors of safety against heave and internal erosion were calculated where the model simulated upward flow of ground water. A critical gradient analysis yields factors of safety as low as 2 near the piezometers where water level observations indicate low factors of safety. Low safety factors are also computed near Castle Creek where slumping was caused by a storm in January, 1990. If hydraulic property contrasts are present in areas of the

19. Characterization of Jupiter's Deep Circulation and Static Stability through Wide Channel Numerical Simulations of the Dynamics and Interactions of Southern Midlatitudes Vortices

Morales-Juberias, Raul; Dowling, Timothy

2012-10-01

Previous studies have shown that the observed features and dynamics of Jovian vortices are sensitive to the underlying environmental structure of Jupiter's atmosphere, in particular to the vertical wind shear and the static stability, and that forward modeling techniques can be successfully used to eliminate a large range of possibilities in a self-consistent manner and hence constrain the atmospheric structure below the cloud regions (Youseff and Marcus 2003, Morales-Juberias et al. 2005). However, these studies have generally been done on a narrow latitude-band basis ( 15^o). Here we present wide channel simulations ( 40 ^o) of two major meteorological events observed in the southern atmosphere of Jupiter involving the interaction of the Great Red Spot (GRS) with other nearby vortices. By studying these two events using wide channel simulations, not unlike the strategy used in terrestrial synoptic meteorology, we show that we can gain new insights into the patterns governing Jupiter's global circulations, drawing a coherent picture of the vertical structure of the atmosphere for the whole southern mid-latitudinal regions of Jupiter over time.

20. Achieving Extreme Utilization of Excitons by an Efficient Sandwich-Type Emissive Layer Architecture for Reduced Efficiency Roll-Off and Improved Operational Stability in Organic Light-Emitting Diodes.

PubMed

Wu, Zhongbin; Sun, Ning; Zhu, Liping; Sun, Hengda; Wang, Jiaxiu; Yang, Dezhi; Qiao, Xianfeng; Chen, Jiangshan; Alshehri, Saad M; Ahamad, Tansir; Ma, Dongge

2016-02-10

It has been demonstrated that the efficiency roll-off is generally caused by the accumulation of excitons or charge carriers, which is intimately related to the emissive layer (EML) architecture in organic light-emitting diodes (OLEDs). In this article, an efficient sandwich-type EML structure with a mixed-host EML sandwiched between two single-host EMLs was designed to eliminate this accumulation, thus simultaneously achieving high efficiency, low efficiency roll-off and good operational stability in the resulting OLEDs. The devices show excellent electroluminescence performances, realizing a maximum external quantum efficiency (EQE) of 24.6% with a maximum power efficiency of 105.6 lm W(-1) and a maximum current efficiency of 93.5 cd A(-1). At the high brightness of 5,000 cd m(-2), they still remain as high as 23.3%, 71.1 lm W(-1), and 88.3 cd A(-1), respectively. And, the device lifetime is up to 2000 h at initial luminance of 1000 cd m(-2), which is significantly higher than that of compared devices with conventional EML structures. The improvement mechanism is systematically studied by the dependence of the exciton distribution in EML and the exciton quenching processes. It can be seen that the utilization of the efficient sandwich-type EML broadens the recombination zone width, thus greatly reducing the exciton quenching and increasing the probability of the exciton recombination. It is believed that the design concept provides a new avenue for us to achieve high-performance OLEDs. PMID:26828128

1. Characterization of Jupiter's Deep Circulation and Static Stability through Wide Channel Numerical Simulations of the Dynamics and Interactions of Southern Midlatitudes Vortices

Morales-Juberias, Raul; Dowling, T. E.

2012-10-01

Previous studies have shown that the observed features and dynamics of Jovian vortices are sensitive to the underlying environmental structure of Jupiter’s atmosphere, in particular to the vertical wind shear and the static stability, and that forward modeling techniques can be successfully used to eliminate a large range of possibilities in a self-consistent manner and hence constrain the atmospheric structure below the cloud regions (Youseff and Marcus 2003, Morales-Juberias et al. 2005). However, these studies have generally been done on a narrow latitude-band basis (˜15°). Here we present wide channel simulations (˜40°) of two major meteorological events observed in the southern atmosphere of Jupiter involving the interaction of the Great Red Spot (GRS) with other nearby vortices. Namely, the spots associated with the recirculation of the South Tropical Disturbance of 1979 (Smith et al. 1979) and the White Ovals (WOS) in 2000 when ovals BE and FA merged to form BA (Sanchez-Lavega et al. 2001). By studying these two events using wide channel simulations, not unlike the strategy used in terrestrial synoptic meteorology, we show that we can gain new insights into the patterns governing Jupiter's global circulations, drawing a coherent picture of the vertical structure of the atmosphere for the whole southern mid-latitudinal regions of Jupiter over time. In particular, we find that the model output best captures the dynamics of the individual vortices and the morphology of their interactions when the deformation length in this region is like that derived by Read et al. 2006 and the deep winds vary following a dependence like that derived by Dowling 1995 in which the westward jets remain constant with depth but the eastward jets increase with depth. Computational resources were provided by the New Mexico Computing Applications Center and NMT. This work was supported by PATM grants NNX08AE91G and NNX08AE64G.

2. Limits to Stability

ERIC Educational Resources Information Center

Cottey, Alan

2012-01-01

The author reflects briefly on what limited degree of global ecological stability and human cultural stability may be achieved, provided that humanity retains hope and does not give way to despair or hide in denial. These thoughts were triggered by a recent conference on International Stability and Systems Engineering. (Contains 5 notes.)

3. Stability of underwater periodic locomotion

Jing, Fangxu; Kanso, Eva

2013-07-01

Most aquatic vertebrates swim by lateral flapping of their bodies and caudal fins. While much effort has been devoted to understanding the flapping kinematics and its influence on the swimming efficiency, little is known about the stability (or lack of) of periodic swimming. It is believed that stability limits maneuverability and body designs/flapping motions that are adapted for stable swimming are not suitable for high maneuverability and vice versa. In this paper, we consider a simplified model of a planar elliptic body undergoing prescribed periodic heaving and pitching in potential flow. We show that periodic locomotion can be achieved due to the resulting hydrodynamic forces, and its value depends on several parameters including the aspect ratio of the body, the amplitudes and phases of the prescribed flapping.We obtain closedform solutions for the locomotion and efficiency for small flapping amplitudes, and numerical results for finite flapping amplitudes. This efficiency analysis results in optimal parameter values that are in agreement with values reported for some carangiform fish. We then study the stability of the (finite amplitude flapping) periodic locomotion using Floquet theory. We find that stability depends nonlinearly on all parameters. Interesting trends of switching between stable and unstable motions emerge and evolve as we continuously vary the parameter values. This suggests that, for live organisms that control their flapping motion, maneuverability and stability need not be thought of as disjoint properties, rather the organism may manipulate its motion in favor of one or the other depending on the task at hand.

4. Influence of Hydraulic Design on Stability and on Pressure Pulsations in Francis Turbines at Overload, Part Load and Deep Part Load based on Numerical Simulations and Experimental Model Test Results

Magnoli, M. V.; Maiwald, M.

2014-03-01

5. MHD stability of tokamak plasmas

SciTech Connect

Chance, M.S. Sun, Y.C.; Jardin, S.C.; Kessel, C.E.; Okabayashi, M.

1992-08-01

This paper will give an overview of the some of the methods which are used to simulate the ideal MHD properties of tokamak plasmas. A great deal of the research in this field is necessarily numerical and the substantial progress made during the past several years has roughly paralleled the continuing availability of more advanced supercomputers. These have become essential to accurately model the complex configurations necessary for achieving MHD stable reactor grade conditions. Appropriate tokamak MHD equilibria will be described. Then the stability properties is discussed in some detail, emphasizing the difficulties of obtaining stable high {beta} discharges in plasmas in which the current is mainly ohmically driven and thus demonstrating the need for tailoring the current and pressure profiles of the plasma away from the ohmic state. The outline of this paper will roughly follow the physics development to attain the second region of stability in the PBX-M device at The Princeton Plasmas Physics Laboratory.

6. Higher order accurate partial implicitization: An unconditionally stable fourth-order-accurate explicit numerical technique

NASA Technical Reports Server (NTRS)

Graves, R. A., Jr.

1975-01-01

The previously obtained second-order-accurate partial implicitization numerical technique used in the solution of fluid dynamic problems was modified with little complication to achieve fourth-order accuracy. The Von Neumann stability analysis demonstrated the unconditional linear stability of the technique. The order of the truncation error was deduced from the Taylor series expansions of the linearized difference equations and was verified by numerical solutions to Burger's equation. For comparison, results were also obtained for Burger's equation using a second-order-accurate partial-implicitization scheme, as well as the fourth-order scheme of Kreiss.

7. Hyperbolic heat conduction problems involving non-Fourier effects - Numerical simulations via explicit Lax-Wendroff/Taylor-Galerkin finite element formulations

NASA Technical Reports Server (NTRS)

Tamma, Kumar K.; Namburu, Raju R.

1989-01-01

Numerical simulations are presented for hyperbolic heat-conduction problems that involve non-Fourier effects, using explicit, Lax-Wendroff/Taylor-Galerkin FEM formulations as the principal computational tool. Also employed are smoothing techniques which stabilize the numerical noise and accurately predict the propagating thermal disturbances. The accurate capture of propagating thermal disturbances at characteristic time-step values is achieved; numerical test cases are presented which validate the proposed hyperbolic heat-conduction problem concepts.

8. Numerical Development

ERIC Educational Resources Information Center

Siegler, Robert S.; Braithwaite, David W.

2016-01-01

In this review, we attempt to integrate two crucial aspects of numerical development: learning the magnitudes of individual numbers and learning arithmetic. Numerical magnitude development involves gaining increasingly precise knowledge of increasing ranges and types of numbers: from non-symbolic to small symbolic numbers, from smaller to larger…

9. Numerical simulations in combustion

NASA Technical Reports Server (NTRS)

Chung, T. J.

1989-01-01

This paper reviews numerical simulations in reacting flows in general and combustion phenomena in particular. It is shown that use of implicit schemes and/or adaptive mesh strategies can improve convergence, stability, and accuracy of the solution. Difficulties increase as turbulence and multidimensions are considered, particularly when finite-rate chemistry governs the given combustion problem. Particular attention is given to the areas of solid-propellant combustion dynamics, turbulent diffusion flames, and spray droplet vaporization.

10. Feedback stabilization initiative

SciTech Connect

1997-06-01

Much progress has been made in attaining high confinement regimes in magnetic confinement devices. These operating modes tend to be transient, however, due to the onset of MHD instabilities, and their stabilization is critical for improved performance at steady state. This report describes the Feedback Stabilization Initiative (FSI), a broad-based, multi-institutional effort to develop and implement methods for raising the achievable plasma betas through active MHD feedback stabilization. A key element in this proposed effort is the Feedback Stabilization Experiment (FSX), a medium-sized, national facility that would be specifically dedicated to demonstrating beta improvement in reactor relevant plasmas by using a variety of MHD feedback stabilization schemes.

11. Stability analysis of the Euler discretization for SIR epidemic model

SciTech Connect

Suryanto, Agus

2014-06-19

In this paper we consider a discrete SIR epidemic model obtained by the Euler method. For that discrete model, existence of disease free equilibrium and endemic equilibrium is established. Sufficient conditions on the local asymptotical stability of both disease free equilibrium and endemic equilibrium are also derived. It is found that the local asymptotical stability of the existing equilibrium is achieved only for a small time step size h. If h is further increased and passes the critical value, then both equilibriums will lose their stability. Our numerical simulations show that a complex dynamical behavior such as bifurcation or chaos phenomenon will appear for relatively large h. Both analytical and numerical results show that the discrete SIR model has a richer dynamical behavior than its continuous counterpart.

12. Trapped particle stability for the kinetic stabilizer

Berk, H. L.; Pratt, J.

2011-08-01

A kinetically stabilized axially symmetric tandem mirror (KSTM) uses the momentum flux of low-energy, unconfined particles that sample only the outer end-regions of the mirror plugs, where large favourable field-line curvature exists. The window of operation is determined for achieving magnetohydrodynamic (MHD) stability with tolerable energy drain from the kinetic stabilizer. Then MHD stable systems are analysed for stability of the trapped particle mode. This mode is characterized by the detachment of the central-cell plasma from the kinetic-stabilizer region without inducing field-line bending. Stability of the trapped particle mode is sensitive to the electron connection between the stabilizer and the end plug. It is found that the stability condition for the trapped particle mode is more constraining than the stability condition for the MHD mode, and it is challenging to satisfy the required power constraint. Furthermore, a severe power drain may arise from the necessary connection of low-energy electrons in the kinetic stabilizer to the central region.

13. Numerical Integration

ERIC Educational Resources Information Center

Sozio, Gerry

2009-01-01

Senior secondary students cover numerical integration techniques in their mathematics courses. In particular, students would be familiar with the "midpoint rule," the elementary "trapezoidal rule" and "Simpson's rule." This article derives these techniques by methods which secondary students may not be familiar with and an approach that…

14. Numerical Relativity

NASA Technical Reports Server (NTRS)

Baker, John G.

2009-01-01

Recent advances in numerical relativity have fueled an explosion of progress in understanding the predictions of Einstein's theory of gravity, General Relativity, for the strong field dynamics, the gravitational radiation wave forms, and consequently the state of the remnant produced from the merger of compact binary objects. I will review recent results from the field, focusing on mergers of two black holes.

15. Graded Achievement, Tested Achievement, and Validity

ERIC Educational Resources Information Center

Brookhart, Susan M.

2015-01-01

Twenty-eight studies of grades, over a century, were reviewed using the argument-based approach to validity suggested by Kane as a theoretical framework. The review draws conclusions about the meaning of graded achievement, its relation to tested achievement, and changes in the construct of graded achievement over time. "Graded…

16. Adapting livestock behaviour to achieve management goals

Technology Transfer Automated Retrieval System (TEKTRAN)

Using livestock to efficiently achieve management goals requires melding animal behavior with mechanical and electronic equipment. Practices such as autonomously obtaining individual animal liveweight when combined with individual animal electronic identification can produce numerous cost saving ad...

17. Coupled trajectory and attitude stability of displaced orbits

Baoyin, Hexi; Li, Junfeng

2010-03-01

Coupled trajectory and attitude stability of displaced solar orbits is studied by using sailcraft with a kind of two-folding construction with two unequal rectangular plates forming a right angle. Three-dimensional coupled trajectory and attitude equations are developed for the coupled dynamical system, and the results show that all three types of displaced solar orbits widely referenced can be achieved through selecting an appropriate size of the two-folding sail. An analysis of the corresponding linear stability of the trajectory and attitude coupled system is carried out, and both trajectory and attitude linearly stable orbits are found to exist in a small range of parameters, whose non-linear stability is then examined via numerical simulations. Finally, passively stable orbits are found to have weak stability, and such passive means of station-keeping are attractive and useful in practice because of its simplicity.

18. Numerical simulation and modeling of combustion in scramjets

Clark, Ryan James

In the last fifteen years the development of a viable scramjet has quickly approached the following long term goals: responsive sub-orbital space access; long-range, prompt global strike; and high-speed transportation. Nonetheless, there are significant challenges that need to be resolved. These challenges include high skin friction drag and high heat transfer rates, inherent to vehicles in sustained, hypersonic flight. Another challenge is sustaining combustion. Numerical simulation and modeling was performed to provide insight into reducing skin friction drag and sustaining combustion. Numerical simulation was used to investigate boundary layer combustion, which has been shown to reduce skin friction drag. The objective of the numerical simulations was to quantify the effect of fuel injection parameters on boundary layer combustion and ultimately on the change in the skin friction coefficient and heat transfer rate. A qualitative analysis of the results suggest that the reduction in the skin friction coefficient depends on multiple parameters and potentially an interaction between parameters. Sustained combustion can be achieved through a stabilized detonation wave. Additionally, stabilizing a detonation wave will yield rapid combustion. This will allow for a shorter and lighter-weight engine system, resulting in less required combustor cooling. A stabilized detonation wave was numerically modeled for various inlet and geometric cases. The effect of fuel concentration, inlet Mach number, and geometric configuration on the stability of a detonation wave was quantified. Correlations were established between fuel concentration, inlet speed, geometric configuration and parameters characterizing the detonation wave. A linear relationship was quantified between the fuel concentration and the parameters characterizing the detonation wave.

19. Dynamics of Numerics and CFD

NASA Technical Reports Server (NTRS)

Yee, H. C.; Rai, Man Mohan (Technical Monitor)

1994-01-01

This lecture attempts to illustrate the basic ideas of how the recent advances in nonlinear dynamical systems theory (dynamics) can provide new insights into the understanding of numerical algorithms used in solving nonlinear differential equations (DEs). Examples will be given of the use of dynamics to explain unusual phenomena that occur in numerics. The inadequacy of the use of linearized analysis for the understanding of long time behavior of nonlinear problems will be illustrated, and the role of dynamics in studying the nonlinear stability, accuracy, convergence property and efficiency of using time- dependent approaches to obtaining steady-state numerical solutions in computational fluid dynamics (CFD) will briefly be explained.

20. Temperature-Stabilized Phase Detector

NASA Technical Reports Server (NTRS)

Yeeman, L.

1985-01-01

Precise temperature stabilized phase detector for clock signal distribution maintains 100-MHz signal with stability of 5 parts in 10 to the 16th power even for step changes of 20 degrees C in ambient temperature. Stabilization achieved by heating unit to 45 degrees C and maintaining temperature through bridge circuit.

1. School Achievement Strongly Predicts Midlife IQ

ERIC Educational Resources Information Center

Spinks, Ruth; Arndt, Stephan; Caspers, Kristin; Yucuis, Rebecca; McKirgan, L. William; Pfalzgraf, Christopher; Waterman, Elijah

2007-01-01

The relationship between measures of IQ and standardized school achievement tests is well established at around r~0.5 when the two are measured in close proximity. The current paper examined the stability of this correlation when comparing elementary school achievement (grades 3-8) and midlife IQ. Iowa Adoption Study participants who had…

2. Thermal-noise limit in the frequency stabilization of lasers with rigid cavities.

PubMed

Numata, Kenji; Kemery, Amy; Camp, Jordan

2004-12-17

We evaluate thermal noise (Brownian motion) in a rigid reference cavity used for frequency stabilization of lasers, based on the mechanical loss of cavity materials and the numerical analysis of the mirror-spacer mechanics with the direct application of the fluctuation dissipation theorem. This noise sets a fundamental limit for the frequency stability achieved with a rigid frequency-reference cavity of order 1 Hz/ square root Hz (0.01 Hz/ square root Hz) at 10 mHz (100 Hz) at room temperature. This level coincides with the world-highest level stabilization results. PMID:15697887

3. Investigation of the stability of AC repulsive force levitation systems for low-speed maglev.

SciTech Connect

He, J. L.; Wang, Z.; Rote, D. M.; Winkelman, S.; Energy Systems

1992-09-01

Discusses the stability of an AC induction levitation system, focusing on the analysis and optimum design of the secondary conductor. Several improved secondary conductor geometries are considered. A theoretical model with numerical results, as well as experimental observations and data are presented. Theoretical and experimental results indicate that only marginal stability can be achieved with a single-plate secondary conductor. Modifications of the single plate can enhance its stability at rest, but this design suffers from longitudinal instabilities when propelled. It is concluded that a double-plate secondary conductor is stable in all six degrees of freedom.

4. Rotational stability of a posterior stabilized total knee arthroplasty.

PubMed

Whiteside, L A; Amador, D D

1989-05-01

The effect of the posterior stabilizing mechanism on rotational stability in total knee arthroplasty (TKA) was investigated in six cadaver knees using a special knee-testing device. The device evaluated varus-valgus, rotational, and anteroposterior (AP) stability in the normal knee compared to a posterior stabilized TKA with either a rotationally constrained or an unconstrained articular surface. None of the stability parameters was significantly different from normal in either configuration of the tibial surface, but the constrained surface did decrease rotational deflection compared to the rotationally unconstrained surfaces. These findings show that rotational constraint in a posterior stabilized TKA is not necessary to achieve rotational stability as long as varus-valgus stability is achieved by appropriately tensioning the collateral ligaments. PMID:2706852

5. Elliptic systems and numerical transformations

NASA Technical Reports Server (NTRS)

Mastin, C. W.; Thompson, J. F.

1976-01-01

Properties of a transformation method, which was developed for solving fluid dynamic problems on general two dimensional regions, are discussed. These include construction error of the transformation and applications to mesh generation. An error and stability analysis for the numerical solution of a model parabolic problem is also presented.

6. Approaches to Numerical Relativity

d'Inverno, Ray

2005-07-01

Introduction Ray d'Inverno; Preface C. J. S. Clarke; Part I. Theoretical Approaches: 1. Numerical relativity on a transputer array Ray d'Inverno; 2. Some aspects of the characteristic initial value problem in numerical relativity Nigel Bishop; 3. The characteristic initial value problem in general relativity J. M. Stewart; 4. Algebraic approachs to the characteristic initial value problem in general relativity Jõrg Frauendiener; 5. On hyperboidal hypersurfaces Helmut Friedrich; 6. The initial value problem on null cones J. A. Vickers; 7. Introduction to dual-null dynamics S. A. Hayward; 8. On colliding plane wave space-times J. B. Griffiths; 9. Boundary conditions for the momentum constraint Niall O Murchadha; 10. On the choice of matter model in general relativity A. D. Rendall; 11. A mathematical approach to numerical relativity J. W. Barrett; 12. Making sense of the effects of rotation in general relativity J. C. Miller; 13. Stability of charged boson stars and catastrophe theory Franz E. Schunck, Fjodor V. Kusmartsev and Eckehard W. Mielke; Part II. Practical Approaches: 14. Numerical asymptotics R. Gómez and J. Winicour; 15. Instabilities in rapidly rotating polytropes Scott C. Smith and Joan M. Centrella; 16. Gravitational radiation from coalescing binary neutron stars Ken-Ichi Oohara and Takashi Nakamura; 17. 'Critical' behaviour in massless scalar field collapse M. W. Choptuik; 18. Goudunov-type methods applied to general relativistic gravitational collapse José Ma. Ibánez, José Ma. Martí, Juan A. Miralles and J. V. Romero; 19. Astrophysical sources of gravitational waves and neutrinos Silvano Bonazzola, Eric Gourgoulhon, Pawel Haensel and Jean-Alain Marck; 20. Gravitational radiation from triaxial core collapse Jean-Alain Marck and Silvano Bonazzola; 21. A vacuum fully relativistic 3D numerical code C. Bona and J. Massó; 22. Solution of elliptic equations in numerical relativity using multiquadrics M. R. Dubal, S. R. Oliveira and R. A. Matzner; 23

7. Unsteady boundary layers with an intelligent numerical scheme

Cebeci, T.

1986-02-01

A numerical method has been developed to represent unsteady boundary layers with large flow reversal. It makes use of the characteristic box scheme which examines the finite-difference grid in relation to the magnitude and direction of local velocity and reaches and implements a decision to ensure that the Courant, Friedricks and Lewey stability criterion is not violated. The method has been applied to the problem of an impulsively started circular cylinder and the results, though generally consistent with those of van Dommelen and Shen obtained with a Lagrangian method, show some differences. The time step is identified as very important and, with the present intelligent numerical scheme, the results were readily extended to times far beyond those previously achieved with Eulerian methods. Extrapolation of the results suggests that the much-discussed singularity for this unsteady flow is the same as that of the corresponding steady flow.

8. 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.

9. Ballooning Stability of the Compact Quasiaxially Symmetric Stellarator

SciTech Connect

M.H. Redi; J. Canik; R.L. Dewar; J.L. Johnson; S. Klasky; W.A. Cooper; W. Kerbichler

2001-09-19

The magnetohydrodynamic (MHD) ballooning stability of a compact, quasiaxially symmetric stellarator (QAS), expected to achieve good stability and particle confinement is examined with a method that can lead to estimates of global stability. Making use of fully 3D, ideal-MHD stability codes, the QAS beta is predicted to be limited above 4% by ballooning and high-n kink modes. Here MHD stability is analyzed through the calculation and examination of the ballooning mode eigenvalue isosurfaces in the 3-space [s, alpha, theta(subscript ''k'')]; s is the edge normalized toroidal flux, alpha is the field line variable, and theta(subscript ''k'') is the perpendicular wave vector or ballooning parameter. Broken symmetry, i.e., deviations from axisymmetry, in the stellarator magnetic field geometry causes localization of the ballooning mode eigenfunction, with new types of nonsymmetric, eigenvalue isosurfaces in both the stable and unstable spectrum. The isosurfaces around the most unstable points i n parameter space (well above marginal) are topologically spherical. In such cases, attempts to use ray tracing to construct global ballooning modes lead to a k-space runaway. Introduction of a reflecting cutoff in k(perpendicular) to model numerical truncation or finite Larmor radius (FLR) yields chaotic ray paths ergodically filling the allowed phase space, indicating that the global spectrum must be described using the language of quantum chaos theory. However, the isosurface for marginal stability in the cases studied are found to have a more complex topology, making estimation of FLR stabilization more difficult.

10. Design and testing of magnetic controllers for Satellite stabilization

Guelman, M.; Waller, R.; Shiryaev, A.; Psiaki, M.

2005-01-01

A study was carried out of attitude control algorithms that are able to provide 3-axis stabilization of a satellite equipped with a magnetometer as the only sensor, and magnetic torquers as the only actuators. Two different solutions to the problem were developed, namely Linear Quadratic Regulator and No Wheel controllers. Their aptitude to achieve the required performance was confirmed by multiple numerical simulations under different initial conditions and various scenarios. The new algorithms were tested onboard the Israeli Gurwin-TechSAT micro-satellite, nominally momentum-biased, stabilized within 2- 2.5∘ precision by the proportion-plus-derivative magnetic controller. In the flight tests of the new controllers, some valuable results were obtained, such as revealing the possibility to effectively maintain the satellite 3-axis stabilization even with a very small momentum bias, and the implementation and efficient performance of the properly modified extended and linear Kalman filters in the onboard computer.

11. Comparing Science Achievement Constructs: Targeted and Achieved

ERIC Educational Resources Information Center

Ferrara, Steve; Duncan, Teresa

2011-01-01

This article illustrates how test specifications based solely on academic content standards, without attention to other cognitive skills and item response demands, can fall short of their targeted constructs. First, the authors inductively describe the science achievement construct represented by a statewide sixth-grade science proficiency test.…

12. Varieties of Achievement Motivation.

ERIC Educational Resources Information Center

Kukla, Andre; Scher, Hal

1986-01-01

A recent article by Nicholls on achievement motivation is criticized on three points: (1) definitions of achievement motives are ambiguous; (2) behavioral consequences predicted do not follow from explicit theoretical assumptions; and (3) Nicholls's account of the relation between his theory and other achievement theories is factually incorrect.…

13. Motivation and School Achievement.

ERIC Educational Resources Information Center

Maehr, Martin L.; Archer, Jennifer

Addressing the question, "What can be done to promote school achievement?", this paper summarizes the literature on motivation relating to classroom achievement and school effectiveness. Particular attention is given to how values, ideology, and various cultural patterns impinge on classroom performance and serve to enhance motivation to achieve.…

ERIC Educational Resources Information Center

Waters, Theresa Z.

A study examined the effect of geographic mobility on elementary school students' achievement. Although such mobility, which requires students to make multiple moves among schools, can have a negative impact on academic achievement, the hypothesis for the study was that it was not a determining factor in reading achievement test scores. Subjects…

ERIC Educational Resources Information Center

Kirby, John R.

Two studies examined the effectiveness of the PASS (Planning, Attention, Simultaneous, and Successive cognitive processes) theory of intelligence in predicting reading achievement scores of normally achieving children and distinguishing children with reading disabilities from normally achieving children. The first study dealt with predicting…

16. Numerical methods in structural mechanics

Obraztsov, I. F.

The papers contained in this volume focus on numerical, numerical-analytical, and theoretical methods for dealing with strength, stability, and dynamics problems in the design of the structural elements of flight vehicles. Topics discussed include the solution of homogeneous boundary value problems for systems of ordinary differential equations modified by a difference factorization method, a study of the rupture strength of a welded joint between plates, singular solutions in mixed problems for a wedge and a half-strip, and a thermoelasticity problem for an open-profile cylindrical shell with a localized temperature field.

17. Numerical analysis of engine instability

Habiballah, M.; Dubois, I.

Following a literature review on numerical analyses of combustion instability, to give the state of the art in the area, the paper describes the ONERA methodology used to analyze the combustion instability in liquid propellant engines. Attention is also given to a model (named Phedre) which describes the unsteady turbulent two-phase reacting flow in a liquid rocket engine combustion chamber. The model formulation includes axial or radial propellant injection, baffles, and acoustic resonators modeling, and makes it possible to treat different engine types. A numerical analysis of a cryogenic engine stability is presented, and the results of the analysis are compared with results of tests of the Viking engine and the gas generator of the Vulcain engine, showing good qualitative agreement and some general trends between experiments and numerical analysis.

18. Heritability of Creative Achievement

ERIC Educational Resources Information Center

Piffer, Davide; Hur, Yoon-Mi

2014-01-01

Although creative achievement is a subject of much attention to lay people, the origin of individual differences in creative accomplishments remain poorly understood. This study examined genetic and environmental influences on creative achievement in an adult sample of 338 twins (mean age = 26.3 years; SD = 6.6 years). Twins completed the Creative…

19. Confronting the Achievement Gap

ERIC Educational Resources Information Center

Gardner, David

2007-01-01

This article talks about the large achievement gap between children of color and their white peers. The reasons for the achievement gap are varied. First, many urban minorities come from a background of poverty. One of the detrimental effects of growing up in poverty is receiving inadequate nourishment at a time when bodies and brains are rapidly…

ERIC Educational Resources Information Center

Christie, Kathy

2002-01-01

Summarizes 2 state initiatives to address the achievement gap: North Carolina's report by the Advisory Commission on Raising Achievement and Closing Gaps, containing an 11-point strategy, and Kentucky's legislation putting in place 10 specific processes. The North Carolina report is available at www.dpi.state.nc.us.closingthegap; Kentucky's…

1. Inverting the Achievement Pyramid

ERIC Educational Resources Information Center

White-Hood, Marian; Shindel, Melissa

2006-01-01

Attempting to invert the pyramid to improve student achievement and increase all students' chances for success is not a new endeavor. For decades, educators have strategized, formed think tanks, and developed school improvement teams to find better ways to improve the achievement of all students. Currently, the No Child Left Behind Act (NCLB) is…

2. Achievement Test Program.

ERIC Educational Resources Information Center

Ohio State Dept. of Education, Columbus. Trade and Industrial Education Service.

The Ohio Trade and Industrial Education Achievement Test battery is comprised of seven basic achievement tests: Machine Trades, Automotive Mechanics, Basic Electricity, Basic Electronics, Mechanical Drafting, Printing, and Sheet Metal. The tests were developed by subject matter committees and specialists in testing and research. The Ohio Trade and…

3. General Achievement Trends: Maryland

ERIC Educational Resources Information Center

Center on Education Policy, 2009

2009-01-01

This general achievement trends profile includes information that the Center on Education Policy (CEP) and the Human Resources Research Organization (HumRRO) obtained from states from fall 2008 through April 2009. Included herein are: (1) Bullet points summarizing key findings about achievement trends in that state at three performance…

4. General Achievement Trends: Arkansas

ERIC Educational Resources Information Center

Center on Education Policy, 2009

2009-01-01

This general achievement trends profile includes information that the Center on Education Policy (CEP) and the Human Resources Research Organization (HumRRO) obtained from states from fall 2008 through April 2009. Included herein are: (1) Bullet points summarizing key findings about achievement trends in that state at three performance…

5. General Achievement Trends: Idaho

ERIC Educational Resources Information Center

Center on Education Policy, 2009

2009-01-01

This general achievement trends profile includes information that the Center on Education Policy (CEP) and the Human Resources Research Organization (HumRRO) obtained from states from fall 2008 through April 2009. Included herein are: (1) Bullet points summarizing key findings about achievement trends in that state at three performance…

ERIC Educational Resources Information Center

Center on Education Policy, 2009

2009-01-01

This general achievement trends profile includes information that the Center on Education Policy (CEP) and the Human Resources Research Organization (HumRRO) obtained from states from fall 2008 through April 2009. Included herein are: (1) Bullet points summarizing key findings about achievement trends in that state at three performance…

ERIC Educational Resources Information Center

Center on Education Policy, 2009

2009-01-01

This general achievement trends profile includes information that the Center on Education Policy (CEP) and the Human Resources Research Organization (HumRRO) obtained from states from fall 2008 through April 2009. Included herein are: (1) Bullet points summarizing key findings about achievement trends in that state at three performance…

8. General Achievement Trends: Iowa

ERIC Educational Resources Information Center

Center on Education Policy, 2009

2009-01-01

This general achievement trends profile includes information that the Center on Education Policy (CEP) and the Human Resources Research Organization (HumRRO) obtained from states from fall 2008 through April 2009. Included herein are: (1) Bullet points summarizing key findings about achievement trends in that state at three performance…

9. General Achievement Trends: Hawaii

ERIC Educational Resources Information Center

Center on Education Policy, 2009

2009-01-01

This general achievement trends profile includes information that the Center on Education Policy (CEP) and the Human Resources Research Organization (HumRRO) obtained from states from fall 2008 through April 2009. Included herein are: (1) Bullet points summarizing key findings about achievement trends in that state at three performance…

10. General Achievement Trends: Kentucky

ERIC Educational Resources Information Center

Center on Education Policy, 2009

2009-01-01

This general achievement trends profile includes information that the Center on Education Policy (CEP) and the Human Resources Research Organization (HumRRO) obtained from states from fall 2008 through April 2009. Included herein are: (1) Bullet points summarizing key findings about achievement trends in that state at three performance…