Numerical analysis of free vibrations of damped rotating structures
NASA Technical Reports Server (NTRS)
Gupta, K. K.
1977-01-01
This paper is concerned with the efficient numerical solution of damped and undamped free vibration problems of rotating structures. While structural discretization is achieved by the finite element method, the associated eigenproblem solution is effected by a combined Sturm sequence and inverse iteration technique that enables the computation of a few required roots only without having to compute any other. For structures of complex configurations, a modal synthesis technique is also presented, which is based on appropriate combinations of eigenproblem solution of various structural components. Such numerical procedures are general in nature, which fully exploit matrix sparsity inherent in finite element discretizations, and prove to be most efficient for the vibration analysis of any damped rotating structure, such as rotating machineries, helicopter and turbine blades, spinning space stations, among others.
Experimental validation of a numerical model for subway induced vibrations
NASA Astrophysics Data System (ADS)
Gupta, S.; Degrande, G.; Lombaert, G.
2009-04-01
This paper presents the experimental validation of a coupled periodic finite element-boundary element model for the prediction of subway induced vibrations. The model fully accounts for the dynamic interaction between the train, the track, the tunnel and the soil. The periodicity or invariance of the tunnel and the soil in the longitudinal direction is exploited using the Floquet transformation, which allows for an efficient formulation in the frequency-wavenumber domain. A general analytical formulation is used to compute the response of three-dimensional invariant or periodic media that are excited by moving loads. The numerical model is validated by means of several experiments that have been performed at a site in Regent's Park on the Bakerloo line of London Underground. Vibration measurements have been performed on the axle boxes of the train, on the rail, the tunnel invert and the tunnel wall, and in the free field, both at the surface and at a depth of 15 m. Prior to these vibration measurements, the dynamic soil characteristics and the track characteristics have been determined. The Bakerloo line tunnel of London Underground has been modelled using the coupled periodic finite element-boundary element approach and free field vibrations due to the passage of a train at different speeds have been predicted and compared to the measurements. The correspondence between the predicted and measured response in the tunnel is reasonably good, although some differences are observed in the free field. The discrepancies are explained on the basis of various uncertainties involved in the problem. The variation in the response with train speed is similar for the measurements as well as the predictions. This study demonstrates the applicability of the coupled periodic finite element-boundary element model to make realistic predictions of the vibrations from underground railways.
Results from Numerical General Relativity
NASA Technical Reports Server (NTRS)
Baker, John G.
2011-01-01
For several years numerical simulations have been revealing the details of general relativity's predictions for the dynamical interactions of merging black holes. I will review what has been learned of the rich phenomenology of these mergers and the resulting gravitational wave signatures. These wave forms provide a potentially observable record of the powerful astronomical events, a central target of gravitational wave astronomy. Asymmetric radiation can produce a thrust on the system which may accelerate the single black hole resulting from the merger to high relative velocity.
Numerical Simulation of Airfoil Vibrations Induced by Compressible Flow
NASA Astrophysics Data System (ADS)
Feistauer, Miloslav; Kučera, Václav; Šimánek, Petr
2010-09-01
The paper is concerned with the numerical solution of interaction of compressible flow and a vibrating airfoil with two degrees of freedom, which can rotate around an elastic axis and oscillate in the vertical direction. Compressible flow is described by the Euler or Navier-Stokes equations written in the ALE form. This system is discretized by the semi-implicit discontinuous Galerkin finite element method (DGFEM) and coupled with the solution of ordinary differential equations describing the airfoil motion. Computational results showing the flow induced airfoil vibrations are presented.
Numerical Approximations of Flow Induced Vibrations of Vocal Folds
NASA Astrophysics Data System (ADS)
Sváček, P.; Horáček, J.
2010-09-01
The paper is interested in numerical modelling of incompressible channel flow interacting with elastic part of its walls simulating vocal fold oscillations. The flow in moving domain is described with the aid of the Arbitrary Lagrangian-Eulerian method, see e.g. [1], and governed by the 2D incompressible Navier-Stokes equations. The flow model is coupled with the structural motion modelled by an aeroelastic two degrees of freedom model of the oscillating vocal folds, cf. [2], [9]. The described fluid-structure interaction problem is discretized in time and space, see also [1]. The numerical results of a channel flow modelling the glottal region of the human vocal tract including the vibrating vocal folds are shown. The vibrations of the channel walls are either prescribed (1st case) or induced by the aerodynamical forces (2nd case).
Experimental and numerical study of the effect of mold vibration on aluminum castings alloys
NASA Astrophysics Data System (ADS)
Abu-Dheir, Numan
2005-07-01
The recent advances in scientific and engineering tools have allowed researchers to integrate more science into manufacturing, leading to improved and new innovative processes. As a result, important accomplishments have been reached in the area of designing and engineering new materials for various industrial applications. This subject is of critical significance because of the impact it could have on the manufacturing industry. In the casting industry, obtaining the desired microstructure and properties during solidification may reduce or eliminate the need for costly thermo-mechanical processing prior to secondary manufacturing processes. Several techniques have been developed to alter and control the microstructure of castings during solidification including semi-solid processing, electromagnetic stirring, electromagnetic vibration, and mechanical vibration. Although it is established that mold vibration can significantly influence the structure and properties of castings, however, most of the studies are generally qualitative, limited to a small range of conditions and no attempts have been made to simulate the effect of vibration on casting microstructure. In this work, a detailed experimental and numerical investigation is carried out to advance the utilization of mold vibration as an effective tool for controlling and modifying the casting microstructure. The effects of a wide range of vibration amplitudes and frequencies on the solidification kinetics, microstructure formation and mechanical properties of Al-Si alloys are examined. Results show strong influence of mold vibration on the resulting casting. The presence of porosity was significantly reduced as a result of mold vibration. In addition, the changes in microstructure and mechanical properties can be successfully represented by the changes in solidification characteristics. Increasing the vibration amplitude tends to reduce the lamellar spacing and change the silicon morphology to become more
NASA Astrophysics Data System (ADS)
Ansari, R.; Hasrati, E.; Faghih Shojaei, M.; Gholami, R.; Shahabodini, A.
2015-05-01
In this paper, the nonlinear forced vibration behavior of composite plates reinforced by carbon nanotubes is investigated by a numerical approach. The reinforcement is considered to be functionally graded (FG) in the thickness direction according to a micromechanical model. The first-order shear deformation theory and von Kármán-type kinematic relations are employed. The governing equations and the corresponding boundary conditions are derived with the use of Hamilton's principle. The generalized differential quadrature (GDQ) method is utilized to achieve a discretized set of nonlinear governing equations. A Galerkin-based scheme is then applied to obtain a time-varying set of ordinary differential equations of Duffing-type. Subsequently, a time periodic discretization is done and the frequency response of plates is determined via the pseudo-arc length continuation method. Selected numerical results are given for the effects of different parameters on the nonlinear forced vibration characteristics of uniformly distributed carbon nanotube- and FG carbon nanotube-reinforced composite plates. It is found that with the increase of CNT volume fraction, the flexural stiffness of plate increases; and hence its natural frequency gets larger. Moreover, it is observed that the distribution type of CNTs significantly affects the vibrational behavior of plate. The results also show that when the mid-plane of plate is CNT-rich, the natural frequency takes its minimum value and the hardening-type response of plate is intensified.
Numerical Manifold Method for the Forced Vibration of Thin Plates during Bending
Jun, Ding; Song, Chen; Wei-Bin, Wen; Shao-Ming, Luo; Xia, Huang
2014-01-01
A novel numerical manifold method was derived from the cubic B-spline basis function. The new interpolation function is characterized by high-order coordination at the boundary of a manifold element. The linear elastic-dynamic equation used to solve the bending vibration of thin plates was derived according to the principle of minimum instantaneous potential energy. The method for the initialization of the dynamic equation and its solution process were provided. Moreover, the analysis showed that the calculated stiffness matrix exhibited favorable performance. Numerical results showed that the generalized degrees of freedom were significantly fewer and that the calculation accuracy was higher for the manifold method than for the conventional finite element method. PMID:24883403
NASA Astrophysics Data System (ADS)
Su, Zhu; Jin, Guoyong; Ye, Tiangui
2016-06-01
The paper presents a unified solution for free and transient vibration analyses of a functionally graded piezoelectric curved beam with general boundary conditions within the framework of Timoshenko beam theory. The formulation is derived by means of the variational principle in conjunction with a modified Fourier series which consists of standard Fourier cosine series and supplemented functions. The mechanical and electrical properties of functionally graded piezoelectric materials (FGPMs) are assumed to vary continuously in the thickness direction and are estimated by Voigt’s rule of mixture. The convergence, accuracy and reliability of the present formulation are demonstrated by comparing the present solutions with those from the literature and finite element analysis. Numerous results for FGPM beams with different boundary conditions, geometrical parameters as well as material distributions are given. Moreover, forced vibration of the FGPM beams subjected to dynamic loads and general boundary conditions are also investigated.
Static and Vibration Analyses of General Wing Structures Using Equivalent Plate Models
NASA Technical Reports Server (NTRS)
Kapania, Rakesh K.; Liu, Youhua
1999-01-01
An efficient method, using equivalent plate model, is developed for studying the static and vibration analyses of general built-up wing structures composed of skins, spars, and ribs. The model includes the transverse shear effects by treating the built-up wing as a plate following the Reissner-Mindlin theory, the so-called First-order Shear Deformation Theory (FSDT). The Ritz method is used with the Legendre polynomials being employed as the trial functions. This is in contrast to previous equivalent plate model methods which have used simple polynomials, known to be prone to numerical ill-conditioning, as the trial functions. The present developments are evaluated by comparing the results with those obtained using MSC/NASTRAN, for a set of examples. These examples are: (i) free-vibration analysis of a clamped trapezoidal plate with (a) uniform thickness, and (b) non-uniform thickness varying as an airfoil, (ii) free-vibration and static analyses (including skin stress distribution) of a general built-up wing, and (iii) free-vibration and static analyses of a swept-back box wing. The results obtained by the present equivalent plate model are in good agreement with those obtained by the finite element method.
Numerical simulation of vortex-induced vibration of a square cylinder at a low Reynolds number
NASA Astrophysics Data System (ADS)
Zhao, Ming; Cheng, Liang; Zhou, Tongming
2013-02-01
Vortex-induced vibrations (VIV) of a square cylinder at a Reynolds number of 100 and a low mass ratio of 3 are studied numerically by solving the Navier-Stokes equations using the finite element method. The equation of motion of the square cylinder is solved to simulate the vibration and the Arbitrary Lagrangian Eulerian scheme is employed to model the interaction between the vibrating cylinder and the fluid flow. The numerical model is validated against the published results of flow past a stationary square cylinder and the results of VIV of a circular cylinder at low Reynolds numbers. The effect of flow approaching angle (α) on the response of the square cylinder is investigated. It is found that α affects not only the vibration amplitude but also the lock-in regime. Among the three values of α (α = 0°, 45°, and 22.5°) that are studied, the smallest vibration amplitude and the narrowest lock-in regime occur at α = 0°. It is discovered that the vibration locks in with the natural frequency in two regimes of reduced velocity for α = 22.5°. Single loop vibration trajectories are observed in the lock-in regime at α = 22.5° and 45°, which is distinctively different from VIV of a circular cylinder. As a result, the vibration frequency in the in-line direction is the same as that in the cross-flow direction.
Numerical study of human vocal folds vibration using Immersed Finite Element Method
NASA Astrophysics Data System (ADS)
Wang, Xingshi; Zhang, Lucy; Krane, Michael
2011-11-01
The voice production procedure is a self-oscillating, fluid-structure interaction problem. In this study, the vocal folds vibration during phonation will be simulated by self-oscillated layered-structure vocal folds model, using Immersed Finite Element Method. With the numerical results, we will find out the vocal folds vibration pattern, and also show how the lung pressure, stiffness and geometry of vocal folds will affect the vocal folds vibration. With further analysis, we shall get better understanding of the dynamics of voice production. National Institute on Deafness and Other Communication Disorders.
Railway cuttings and embankments: Experimental and numerical studies of ground vibration.
Kouroussis, Georges; Connolly, David P; Olivier, Bryan; Laghrouche, Omar; Costa, Pedro Alves
2016-07-01
Railway track support conditions affect ground-borne vibration generation and propagation. Therefore this paper presents a combined experimental and numerical study into high speed rail vibrations for tracks on three types of support: a cutting, an embankment and an at grade section. Firstly, an experimental campaign is undertaken where vibrations and in-situ soil properties are measured at three Belgian rail sites. A finite element model is then developed to recreate the complex ground topology at each site. A validation is performed and it is found that although the at-grade and embankment cases show a correlation with the experimental results, the cutting case is more challenging to replicate. Despite this, each site is then analysed to determine the effect of earthworks profile on ground vibrations, with both the near and far fields being investigated. It is found that different earthwork profiles generate strongly differing ground-borne vibration characteristics, with the embankment profile generating lower vibration levels in comparison to the cutting and at-grade cases. Therefore it is concluded that it is important to consider earthwork profiles when undertaking vibration assessments. PMID:26994799
Numerical Investigation of Flapwise-Torsional Vibration Model of a Smart Section Blade with Microtab
Li, Nailu; Balas, Mark J.; Yang, Hua; Jiang, Wei; Magar, Kaman T.
2015-01-01
This study presents a method to develop an aeroelastic model of a smart section blade equipped with microtab. The model is suitable for potential passive vibration control study of the blade section in classic flutter. Equations of the model are described by the nondimensional flapwise and torsional vibration modes coupled with the aerodynamic model based on the Theodorsen theory and aerodynamic effects of the microtab based on the wind tunnel experimental data. The aeroelastic model is validated using numerical data available in the literature and then utilized to analyze the microtab control capability on flutter instability case and divergence instabilitymore » case. The effectiveness of the microtab is investigated with the scenarios of different output controllers and actuation deployments for both instability cases. The numerical results show that the microtab can effectively suppress both vibration modes with the appropriate choice of the output feedback controller.« less
Research on algorithm of blade vibration for general wind turbine
NASA Astrophysics Data System (ADS)
Wang, Long; Sun, Lun-ye; Wu, Guang; Li, Xue-bin; Lai, Yong-bin; Zhou, Yi-jun
2016-01-01
Evaluation of vibration characteristics for wind turbine blades is one of the important contents in the wind turbine research. This paper uses the compressible flow equations with the preconditioning technique, based on the finite volume method and combined with the LU-SGS algorithm for solving the flow area; meanwhile adopts the two degree of freedom of vibration equation with the vertical and torsional vibration for blades to simulate the vibration trajectory of blade under the aerodynamic force, uses the motion grid algorithm for changes in grid computing domain. Calculation program was developed autonomous in the C ++ platform, and the development of software correctness was verified by contrast the results of the classic cylindrical examples. Finally, the vibration characteristics of a wind turbine blade was given, and the software developed in this paper can provide technical support for wind turbine blade vibration study.
Methods for numerical study of tube bundle vibrations in cross-flows
NASA Astrophysics Data System (ADS)
Longatte, E.; Bendjeddou, Z.; Souli, M.
2003-11-01
In many industrial applications, mechanical structures like heat exchanger tube bundles are subjected to complex flows causing possible vibrations and damage. Part of fluid forces are coupled with tube motion and the so-called fluid-elastic forces can affect the structure dynamic behaviour generating possible instabilities and leading to possible short term failures through high amplitude vibrations. Most classical fluid force identification methods rely on structure response experimental measurements associated with convenient data processes. Owing to recent improvements in Computational Fluid Dynamics, numerical simulation of flow-induced vibrations is now practicable for industrial purposes. The present paper is devoted to the numerical identification of fluid-elastic effects affecting tube bundle motion in presence of fluid at rest and one-phase cross-flows. What is the numerical process? When fluid-elastic effects are not significant and are restricted to added mass effects, there is no strong coupling between structure and fluid motions. The structure displacement is not supposed to affect flow patterns. Thus it is possible to solve flow and structure problems separately by using a fixed nonmoving mesh for the fluid dynamic computation. Power spectral density and time record of lift and drag forces acting on tube bundles can be computed numerically by using an unsteady fluid computation involving for example a large Eddy simulation. Fluid force spectra or time record can then be introduced as inlet conditions into the structure code providing the tube dynamic response generated by flow. Such a computation is not possible in presence of strong flow structure coupling. When fluid-elastic effects cannot be neglected, in presence of tube bundles subjected to cross-flows for example, a coupling between flow and structure computations is required. Appropriate numerical methods are investigated in the present work. The purpose is to be able to provide a numerical
NASA Astrophysics Data System (ADS)
Dong, Xing-Jian; Meng, Guang; Peng, Juan-Chun
2006-11-01
The aim of this study is to investigate the efficiency of a system identification technique known as observer/Kalman filter identification (OKID) technique in the numerical simulation and experimental study of active vibration control of piezoelectric smart structures. Based on the structure responses determined by finite element method, an explicit state space model of the equivalent linear system is developed by employing OKID approach. The linear quadratic Gaussian (LQG) algorithm is employed for controller design. The control law is then incorporated into the ANSYS finite element model to perform closed loop simulations. Therefore, the control law performance can be evaluated in the context of a finite element environment. Furthermore, a complete active vibration control system comprising the cantilever plate, the piezoelectric actuators, the accelerometers and the digital signal processor (DSP) board is set up to conduct the experimental investigation. A state space model characterizing the dynamics of the physical system is developed from experimental results using OKID approach for the purpose of control law design. The controller is then implemented by using a floating point TMS320VC33 DSP. Numerical examples by employing the proposed numerical simulation method, together with the experimental results obtained by using the active vibration control system, have demonstrated the validity and efficiency of OKID method in application of active vibration control of piezoelectric smart structures.
Sharma, S.P.; Huo, W.M.; Park, C.
1988-01-01
A theoretical study of vibrational excitations and dissociations of nitrogen undergoing a nonequilibrium relaxation process upon heating and cooling is reported. The rate coefficients for collisional induced vibrational transitions and transitions from a bound vibrational state into a dissociative state have been calculated using an extension of the theory originally proposed by Schwarz (SSH) et al. (1952). High-lying vibrational states and dissociative states were explicitly included but rotational energy transfer was neglected. The transition probabilities calculated from the SSH theory were fed into the master equation, which was integrated numerically to determine the population distribution of the vibrational states as well as bulk thermodynamic properties. The results show that: (1) the transition rates have a minimum near the middle of the bound vibrational levels, causing a bottleneck in the vibrational relaxation and dissociation rates; (2) high vibrational states are always in equilibrium with the dissociative state; (3) for the heating case, only the low vibrational states relax according to the Landau-Teller theory; (4) for the cooling case, vibrational relaxation cannot be described by a rate equation; (5) Park's (1985, 1988) two-temperature model is approximately valid; and (6) the average vibrational energy removed in dissociation is about 30 percent of the dissociation energy. 29 references.
Miliordos, Evangelos; Xantheas, Sotiris S.
2013-08-15
We propose a general procedure for the numerical calculation of the harmonic vibrational frequencies that is based on internal coordinates and Wilson’s GF methodology via double differentiation of the energy. The internal coordinates are defined as the geometrical parameters of a Z-matrix structure, thus avoiding issues related to their redundancy. Linear arrangements of atoms are described using a dummy atom of infinite mass. The procedure has been automated in FORTRAN90 and its main advantage lies in the nontrivial reduction of the number of single-point energy calculations needed for the construction of the Hessian matrix when compared to the corresponding number using double differentiation in Cartesian coordinates. For molecules of C_{1} symmetry the computational savings in the energy calculations amount to 36N – 30, where N is the number of atoms, with additional savings when symmetry is present. Typical applications for small and medium size molecules in their minimum and transition state geometries as well as hydrogen bonded clusters (water dimer and trimer) are presented. Finally, in all cases the frequencies based on internal coordinates differ on average by <1 cm^{–1} from those obtained from Cartesian coordinates.
NASA Astrophysics Data System (ADS)
Wan, Qun
The purpose of this research is to numerically and analytically investigate the acoustic streaming and the associated heat transfer, which are induced by a beam vibrating in either standing or traveling waveforms. Analytical results show that the beam vibrating in standing waveforms scatters the acoustic waves into the free space, which have a larger attenuation coefficient and longer propagating traveling wavelength than those of the plane wave. In contrast to a constant Reynolds stress in the plane wave, the Reynolds stress generated by such acoustic wave is expected to drive the free space streaming away from the anti-nodes and towards nodes of the standing wave vibration. The sonic and ultrasonic streamings within the channel between the vibrating beam and a parallel stationary beam are also investigated. The acoustic streaming is utilized to cool the stationary beam, which has either a heat source attached to it or subjected to a uniform heat flux. The sonic streaming is found to be mainly the boundary layer streaming dominating the whole channel while the ultrasonic streaming is clearly composed of two boundary layer streamings near both beams and a core region streaming, which is driven by the streaming velocity at the edge of the boundary layer near the vibrating beam. The standing wave vibration of the beam induces acoustic streaming in a series of counterclockwise eddies, which is directed away from the anti-nodes and towards the nodes. The magnitude of the sonic streaming is proportional to o2A 2 while that of the ultrasonic streaming is proportional to o 3/2A2. Numerical results show that the acoustic streaming induced by the beam vibrating in either standing or traveling waveforms has almost the same cooling efficiency for the heat source and the heat flux cases although the flow and temperature fields within the channel are different. The hysteresis of the ultrasonic streaming flow patterns associated with the change of the aspect ratio of the channel
NASA Astrophysics Data System (ADS)
Kuo, K. A.; Verbraken, H.; Degrande, G.; Lombaert, G.
2016-07-01
Along with the rapid expansion of urban rail networks comes the need for accurate predictions of railway induced vibration levels at grade and in buildings. Current computational methods for making predictions of railway induced ground vibration rely on simplifying modelling assumptions and require detailed parameter inputs, which lead to high levels of uncertainty. It is possible to mitigate against these issues using a combination of field measurements and state-of-the-art numerical methods, known as a hybrid model. In this paper, two hybrid models are developed, based on the use of separate source and propagation terms that are quantified using in situ measurements or modelling results. These models are implemented using term definitions proposed by the Federal Railroad Administration and assessed using the specific illustration of a surface railway. It is shown that the limitations of numerical and empirical methods can be addressed in a hybrid procedure without compromising prediction accuracy.
Numerical Simulation of Planar and Axisymmetric Unsteady Flows Over Vibrating Bodies
NASA Astrophysics Data System (ADS)
Venkat, N. Kolluru
1991-02-01
Engineering applications of unsteady boundary layers are numerous and of great importance in hydro and aero dynamics. Relatively little research however is focussed on understanding the flow structure if the perturbations in the flow are caused by the motion of flexible boundaries rather than by time dependent variations in the flow itself. A model is therefore developed to examine the hydrodynamic characteristics of external laminar flow over an arbitrarily shaped body, a portion of which is subjected to harmonic motion. A vorticity-stream function formulation of the Navier-Stokes equations is used. A boundary fitted coordinate system is adopted to allow accurate modeling in the presence of the time dependent motion of the body surface. The flow equations are solved using a Thompson Tri Diagonal Finite Difference Algorithm. Inviscid-viscous interaction theory is used to split the model domain and save computational time. The model is tested by comparison to selected numerical, experimental or analytical results for flow over a cavity, boundary layer flow along a flat plate and time dependent flow in a channel with wavy walls. The model is then applied to predict the flow over a flat plate and circular cylinder with a section forced in simple harmonic motion. The nonlinear response of the flow is investigated for various Reynolds numbers, Re, Strouhal number, St, (ratio of the flow advective time scale to the plate or cylinder oscillation period), vibration amplitude ration, bf{H _ 0} (vibration amplitude divided by plate or cylinder vibrating section length) and the vibration mode number, n. The time varying friction (bf{C_ f}) and pressure ( bf{C_ p}) coefficients along the vibrating plate and cylinder are analyzed using Fast Fourier Transform techniques (FFT). Model results show that for low St and bf {H_ 0} the response of the external flow to the plate vibration amplitude is linear and there is little up or downstream influence. For high St, Re and bf{H_ 0
NASA Astrophysics Data System (ADS)
Silva, Gustavo H. C.; Paupitz Gonçalves, Paulo J.
2013-09-01
This paper develops a novel full analytic model for vibration analysis of solid-state electronic components. The model is just as accurate as finite element models and numerically light enough to permit for quick design trade-offs and statistical analysis. The paper shows the development of the model, comparison to finite elements and an application to a common engineering problem. A gull-wing flat pack component was selected as the benchmark test case, although the presented methodology is applicable to a wide range of component packages. Results showed very good agreement between the presented method and finite elements and demonstrated the usefulness of the method in how to use standard test data for a general application. The properties E, G, A, I, J and κ need not be constants; they may all be functions of s.
NASA Astrophysics Data System (ADS)
Haque, Md. Rejaul; Chowdhury, M. Arshad Zahangir; Goswami, Anjan
2016-07-01
A two-dimensional numerical study of flow induced vibration is reported in this paper to investigate flow over a semi-cricular D-shaped bluff body oriented at different angles-of-attack to determine an optimized design for energy harvesting. Bluff body structure governs fluid streamlines; therefore obtaining a suitable range of "lock in frequency" for energy harvesting purpose is dependent on refining and optimizing bluff body's shape and structure. A cantilever based novel energy harvester design incorporates the suitable angle-of-attack for optimized performance. This optimization was done by performing computations for 30°, 60° and 90° angles-of-attack. The frequency of vibration of the body was calculated at different Reynolds Number. A Fast Fourier Transformation yielded frequency of vortex shedding. From the wake velocity profile, lift oscillation and frequency of vortex shedding is estimated. Strouhal numbers of the body were analyzed at different angles-of-attack. A higher synchronized bandwidth of shedding frequencies is an indication of an optimized harvester design at different Reynolds number. The `D' shaped bluff bodies (with angle of attack of 30°,60° and 90°) are more suitable than that of cylindrical shaped bluff bodies. The research clearly stated that, bluff bodies shape has a prominent influence on vortex induced vibration and semicircular bluff body gives the highest vibration or energy under stated conditions.
NASA Astrophysics Data System (ADS)
Coulier, P.; Lombaert, G.; Degrande, G.
2014-06-01
The numerical prediction of vibrations in buildings due to railway traffic is a complicated problem where wave propagation in the soil couples the source (railway tunnel or track) and the receiver (building). This through-soil coupling is often neglected in state-of-the-art numerical models in order to reduce the computational cost. In this paper, the effect of this simplifying assumption on the accuracy of numerical predictions is investigated. A coupled finite element-boundary element methodology is employed to analyze the interaction between a building and a railway tunnel at depth or a ballasted track at the surface of a homogeneous halfspace, respectively. Three different soil types are considered. It is demonstrated that the dynamic axle loads can be calculated with reasonable accuracy using an uncoupled strategy in which through-soil coupling is disregarded. If the transfer functions from source to receiver are considered, however, large local variations in terms of vibration insertion gain are induced by source-receiver interaction, reaching up to 10 dB and higher, although the overall wave field is only moderately affected. A global quantification of the significance of through-soil coupling is made, based on the mean vibrational energy entering a building. This approach allows assessing the common assumption in seismic engineering that source-receiver interaction can be neglected if the distance between source and receiver is sufficiently large compared to the wavelength of waves in the soil. It is observed that the interaction between a source at depth and a receiver mainly affects the power flow distribution if the distance between source and receiver is smaller than the dilatational wavelength in the soil. Interaction effects for a railway track at grade are observed if the source-receiver distance is smaller than six Rayleigh wavelengths. A similar trend is revealed if the passage of a freight train is considered. The overall influence of dynamic
Free vibrations of beam-mass-spring systems: analytical analysis with numerical confirmation
NASA Astrophysics Data System (ADS)
Darabi, Mohammad A.; Kazemirad, Siavash; Ghayesh, Mergen H.
2012-04-01
Free vibrations of a beam-mass-spring system with different boundary conditions are analyzed both analytically and numerically. In the analytical analysis, the system is divided into three subsystems and the effects of the spring and the point mass are considered as internal boundary conditions between any two neighboring subsystems. The partial differential equations governing the motion of the subsystems and internal boundary conditions are then solved using the method of separation of variables. In the numerical analysis, the whole system is considered as a single system and the effects of the spring and point mass are introduced using the Dirac delta function. The Galerkin method is then employed to discretize the equation of motion and the resulting set of ordinary differential equations are solved via eigenvalue analysis. Analytical and numerical results are shown to be in very good agreement.
Numerical study of liquid-hydrogen droplet generation from a vibrating orifice
NASA Astrophysics Data System (ADS)
Xu, J.; Celik, D.; Hussaini, M. Y.; Van Sciver, S. W.
2005-08-01
Atomic hydrogen propellant feed systems for far-future spacecraft may utilize solid-hydrogen particle carriers for atomic species that undergo recombination to create hot rocket exhaust. Such technology will require the development of particle generation techniques. One such technique could involve the production of hydrogen droplets from a vibrating orifice that would then freeze in cryogenic helium vapor. Among other quantities, the shape and size of the droplet are of particular interest. The present paper addresses this problem within the framework of the incompressible Navier-Stokes equations for multiphase flows, in order to unravel the basic mechanisms of droplet formation with a view to control them. Surface tension, one of the most important mechanisms to determine droplet shape, is modeled as the source term in the momentum equation. Droplet shape is tracked using a volume-of-fluid approach. A dynamic meshing technique is employed to accommodate the vibration of the generator orifice. Numerically predicted droplet shapes show satisfactory agreement with photographs of droplets generated in experiments. A parametric study is carried out to understand the influence of injection velocity, nozzle vibrational frequency, and amplitude on the droplet shape and size. The computational model provides a definitive qualitative picture of the evolution of droplet shape as a function of the operating parameters. It is observed that, primarily, the orifice vibrational frequency affects the shape, the vibrational amplitude affects the time until droplet detachment from the orifice, and the injection velocity affects the size. However, it does not mean that, for example, there is no secondary effect of amplitude on shape or size.
Numerical and experimental investigation of natural flow-induced vibrations of flexible hydrofoils
NASA Astrophysics Data System (ADS)
Chae, Eun Jung; Akcabay, Deniz Tolga; Lelong, Alexandra; Astolfi, Jacques Andre; Young, Yin Lu
2016-07-01
The objective of this work is to present combined numerical and experimental studies of natural flow-induced vibrations of flexible hydrofoils. The focus is on identifying the dependence of the foil's vibration frequencies and damping characteristics on the inflow velocity, angle of attack, and solid-to-fluid added mass ratio. Experimental results are shown for a cantilevered polyacetate (POM) hydrofoil tested in the cavitation tunnel at the French Naval Academy Research Institute (IRENav). The foil is observed to primarily behave as a chordwise rigid body and undergoes spanwise bending and twisting deformations, and the flow is observed to be effectively two-dimensional (2D) because of the strong lift retention at the free tip caused by a small gap with a thickness less than the wall boundary layer. Hence, the viscous fluid-structure interaction (FSI) model is formulated by coupling a 2D unsteady Reynolds-averaged Navier-Stokes (URANS) model with a two degree-of-freedom (2-DOF) model representing the spanwise tip bending and twisting deformations. Good agreements were observed between viscous FSI predictions and experimental measurements of natural flow-induced vibrations in fully turbulent and attached flow conditions. The foil vibrations were found to be dominated by the natural frequencies in absence of large scale vortex shedding due to flow separation. The natural frequencies and fluid damping coefficients were found to vary with velocity, angle of attack, and solid-to-fluid added mass ratio. In addition, the numerical results showed that the in-water to in-air natural frequency ratios decreased rapidly, and the fluid damping coefficients increased rapidly, as the solid-to-fluid added mass ratio decreases. Uncoupled mode (UM) linear potential theory was found to significantly over-predict the fluid damping for cases of lightweight flexible hydrofoils, and this over-prediction increased with higher velocity and lower solid-to-fluid added mass ratio.
NASA Astrophysics Data System (ADS)
Ostasevicius, V.; Ubartas, M.; Gaidys, R.; Jurenas, V.; Samper, S.; Dauksevicius, R.
2012-11-01
This study is concerned with application of numerical-experimental approach for characterizing dynamic behavior of the developed piezoelectrically excited vibration drilling tool with the aim to identify the most effective conditions of tool vibration mode control for improved cutting efficiency. 3D finite element model of the tool was created on the basis of an elastically fixed pre-twisted cantilever (standard twist drill). The model was experimentally verified and used together with tool vibration measurements in order to reveal rich dynamic behavior of the pre-twisted structure, representing a case of parametric vibrations with axial, torsional and transverse natural vibrations accompanied by the additional dynamic effects arising due to the coupling of axial and torsional deflections ((un)twisting). Numerical results combined with extensive data from interferometric, accelerometric, dynamometric and surface roughness measurements allowed to determine critical excitation frequencies and the corresponding vibration modes, which have the largest influence on the performance metrics of the vibration drilling process. The most favorable tool excitation conditions were established: inducing the axial mode of the vibration tool itself through tailoring of driving frequency enables to minimize magnitudes of surface roughness, cutting force and torque. Research results confirm the importance of the tool mode control in enhancing the effectiveness of vibration cutting tools from the viewpoint of structural dynamics.
Numerical Dimension-Reduction Methods for Non-Linear Shell Vibrations
NASA Astrophysics Data System (ADS)
Foale, S.; Thompson, J. M. T.; McRobie, F. A.
1998-08-01
A number of methods are investigated for obtaining a low-dimensional dynamical system from a set of partial differential equations describing the non-linear vibrations of a shallow cylindrical panel under periodic axial forcing. In these approaches an initial (high-dimensional) spatial discretization of a (possibly irregular) domain is performed and a subsequent procedure is used to further reduce the resulting set of ordinary differential equations. In particular the results suggest that a numerical method based upon inertial manifold approximation is possible, but for the specific case studied, no advantage could be discerned over more direct dimension-reduction techniques.
Wang Lei; Babikov, Dmitri
2011-02-15
Anharmonicity of the quantized motional states of ions in a Paul trap can be utilized to address the state-to-state transitions selectively and control the motional modes of trapped ions coherently and adiabatically [Zhao and Babikov, Phys. Rev. A 77, 012338 (2008)]. In this paper we study two sources of the vibrational anharmonicity in the ion traps: the intrinsic Coulomb anharmonicity due to ion-ion interactions and the external anharmonicity of the trapping potential. An accurate numerical approach is used to compute energies and wave functions of vibrational eigenstates. The magnitude of the Coulomb anharmonicity is determined and shown to be insufficient for successful control. In contrast, anharmonicity of the trapping potential allows one to control the motion of ions very efficiently using the time-varying electric fields. Optimal control theory is used to derive the control pulses. One ion in a slightly anharmonic trap can be easily controlled. In the two- and three-ion systems the symmetric stretching mode is dark and cannot be controlled at all. The other two normal modes of the three-ion system can be controlled and used, for example, to encode a two-qubit system into the motional states of ions. A trap architecture that allows the necessary amount of vibrational anharmonicity to be achieved is proposed.
Numerical stability of pseudo-spectral PIC code generalizations
NASA Astrophysics Data System (ADS)
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.
NASA Technical Reports Server (NTRS)
Sutton, M. A.; Davis, P. K.
1976-01-01
Numerical solutions of the governing equations of motion of a liquid squeeze film damped forced vibration system were carried out to examine the feasibility of using a liquid squeeze film to cushion and protect large structures, such as buildings, located in areas of high seismic activity. The mathematical model used was that for a single degree of freedom squeeze film damped spring mass system. The input disturbance was simulated by curve fitting actual seismic data with an eleventh order Lagranging polynomial technique. Only the normal component of the seismic input was considered. The nonlinear, nonhomogeneous governing differential equation of motion was solved numerically to determine the transmissibility over a wide range of physical parameters using a fourth-order Runge-Kutta technique. It is determined that a liquid squeeze film used as a damping agent in a spring-mass system can significantly reduce the response amplitude for a seismic input disturbance.
An Asymptotic-Numerical Method for Large-Amplitude Free Vibrations of Thin Elastic Plates
NASA Astrophysics Data System (ADS)
Azrar, L.; Benamar, R.; Potier-Ferry, M.
1999-03-01
An Asymptotic-Numerical Method has been developed for large amplitude free vibrations of thin elastic plates. It is based on the perturbation method and the finite element method. This method eliminates the major difficulties of the classical perturbation methods, namely the complexity of the right hand sides and the limitation of the validity of the solution obtained. The applicability of this method to non-linear vibrations of plates is clearly presented. Based on the Von Karman theory and the harmonic balance method, a cubic non-linear operational formulation has been obtained. By using the mixed stress-displacement Hellinger-Reissner principle, a quadratic formulation is given. The displacement and frequency are expanded into power series with respect to a control parameter. The non-linear governing equation is then transformed into a sequence of linear problems having the same stiffness matrix, which can be solved by a classical FEM. Needing one matrix inversion, a large number of terms of the series can be easily computed with a small computation time. The non-linear mode and frequency are then obtained up to the radius of convergence. Taking the starting point in the zone of validity, the method is reapplied in order to determine a further part of the non-linear solution. Iteration of this method leads to a powerful incremental method. In order to increase the validity of the perturbed solution, another technique, called Padé approximants, is shrewdly incorporated. The solutions obtained by these two concepts coincide perfectly in a very large part of the backbone curve. Comprehensive numerical tests for non-linear free vibrations of circular, square, rectangular and annular plates with various boundary conditions are reported and discussed.
PERTURB: A program for calculating vibrational energies by generalized algebraic quantization
NASA Astrophysics Data System (ADS)
Fried, Laurence E.; Ezra, Gregory S.
1988-09-01
We describe PERTURB, a special purpose algebraic manipulation program which calculates vibrational eigenvalues in coupled oscillator systems. PERTURB implements the method of generalized algebraic quantization (AQ), in which Van Vleck perturbation theory is formulated in a mock phase space. The phase space formulation enables quantum and classical perturbation theory to be treated on the same footing, and allows the systematic calculation of corrections to classical perturbation results in powers of h̷. Generalized AQ is a powerful and efficient technique for calculating semiclassical vibrational energy levels. In many cases, including just the first correction to classical perturbation theory yields highly accurate energies.
Numerical simulation of vibration of horizontal cylinder induced by progressive waves
NASA Astrophysics Data System (ADS)
Chern, Ming-Jyh; Odhiambo, E. A.; Horng, Tzyy-Leng; Borthwick, A. G. L.
2016-02-01
Maritime structures often comprise cylinders of small diameter relative to the prevailing wave length. This paper describes the direct forcing immersed boundary simulation of the hydroelastic behaviour of a rigid, horizontal circular cylinder in regular progressive waves. Fluid motions are numerically solved by the full Navier-Stokes equations, and the free surface by the volume-of-fluid method. The Reynolds number Re = 110, Keulegan-Carpenter number KC = 10, Froude number Fr = 0.69 and Ursell number U rs ≈ 12. A single-degree-of-freedom model is used for the elastically mounted cylinder. Velocity profiles for the stationary cylinder case have been successfully validated using experimental results. The frequency response for reduced velocities 4.5\\lt {U}R*\\lt 5.3 have been compared with theoretical data. Three transverse vibration regimes are identified: lower beating (4\\lt {U}R*\\lt 4.5); lock-in (4.7\\lt {U}R*\\lt 4.8); and upper beating (5\\lt {U}R*\\lt 10) modes. The lower and upper beating regimes exhibit varying amplitude response. The lock-in mode represents the region of fixed and maximum response. The lower beating and lock-in modes have peaks at a common vibration to wave frequency ratio {f}{{w}}* = 2. For the upper beating mode, {f}{{w}}* = 1, except for {U}R*=10 when {f}{{w}}* = 2.
Wave interpretation of numerical results for the vibration in thin conical shells
NASA Astrophysics Data System (ADS)
Ni, Guangjian; Elliott, Stephen J.
2014-05-01
The dynamic behaviour of thin conical shells can be analysed using a number of numerical methods. Although the overall vibration response of shells has been thoroughly studied using such methods, their physical insight is limited. The purpose of this paper is to interpret some of these numerical results in terms of waves, using the wave finite element, WFE, method. The forced response of a thin conical shell at different frequencies is first calculated using the dynamic stiffness matrix method. Then, a wave finite element analysis is used to calculate the wave properties of the shell, in terms of wave type and wavenumber, as a function of position along it. By decomposing the overall results from the dynamic stiffness matrix analysis, the responses of the shell can then be interpreted in terms of wave propagation. A simplified theoretical analysis of the waves in the thin conical shell is also presented in terms of the spatially-varying ring frequency, which provides a straightforward interpretation of the wave approach. The WFE method provides a way to study the types of wave that travel in thin conical shell structures and to decompose the response of the numerical models into the components due to each of these waves. In this way the insight provided by the wave approach allows us to analyse the significance of different waves in the overall response and study how they interact, in particular illustrating the conversion of one wave type into another along the length of the conical shell.
Numerical solution of quadratic matrix equations for free vibration analysis of structures
NASA Technical Reports Server (NTRS)
Gupta, K. K.
1975-01-01
This paper is concerned with the efficient and accurate solution of the eigenvalue problem represented by quadratic matrix equations. Such matrix forms are obtained in connection with the free vibration analysis of structures, discretized by finite 'dynamic' elements, resulting in frequency-dependent stiffness and inertia matrices. The paper presents a new numerical solution procedure of the quadratic matrix equations, based on a combined Sturm sequence and inverse iteration technique enabling economical and accurate determination of a few required eigenvalues and associated vectors. An alternative procedure based on a simultaneous iteration procedure is also described when only the first few modes are the usual requirement. The employment of finite dynamic elements in conjunction with the presently developed eigenvalue routines results in a most significant economy in the dynamic analysis of structures.
Vibrational Spectra of Molecular Crystals with the Generalized Energy-Based Fragmentation Approach.
Fang, Tao; Jia, Junteng; Li, Shuhua
2016-05-01
The generalized energy-based fragmentation (GEBF) approach for molecular crystals with periodic boundary condition (PBC) (denoted as PBC-GEBF) is extended to allow vibrational spectra of molecular crystals to be easily computed at various theory levels. Within the PBC-GEBF approach, the vibrational frequencies of a molecular crystal can be directly evaluated from molecular quantum chemistry calculations on a series of nonperiodic molecular systems. With this approach, the vibrational spectra of molecular crystals can be calculated with much reduced computational costs at various theory levels, as compared to those required by the methods based on periodic electronic structure theory. By testing the performance of the PBC-GEBF method for two molecular crystals (CO2 and imidazole), we demonstrate that the PBC-GEBF approach can reproduce the results of the methods based on periodic electronic structure theory in predicting vibrational spectra of molecular crystals. We apply the PBC-GEBF method at second-order Møller-Plesset perturbation theory (PBC-GEBF-MP2 in short) to investigate the vibrational spectra of the urea and ammonia borane crystals. Our results show that the PBC-GEBF-MP2 method can provide quite accurate descriptions for the observed vibrational spectra of the two systems under study. PMID:27076120
NASA Astrophysics Data System (ADS)
Lin, Jian-Zhong; Jiang, Ren-Jie; Ku, Xiao-Ke
2014-03-01
Dynamics of vortex-induced vibrations (VIV) of one and two elastically mounted circular cylinders in tandem is numerically studied via a lattice Boltzmann method. The structural damping is set to zero and the cylinders are of one degree of freedom. The Reynolds number is in the range of Re = [5, 100], covering the sub- and super-critical regions. The effects of reduced velocity UR = [4, 10], center-to-center spacing between the two cylinders S/D = [1.1, 4] and mass ratio M = [0.5, 4] on the cylinder motion and flow structures are investigated. The VIV at subcritical Reynolds number exists extensively for both one- and two-cylinder systems. Besides the traditional symmetrical vibration regime, an anomalous biased vibration regime is observed for the case of two tandem cylinders. It is the first time that the existence of such an asymmetric vibration regime has been confirmed. The biased oscillation regime mainly exists in the region of subcritical Reynolds number, large reduced velocity, small spacing, and low mass ratio. As the reduced velocity is more than 7, such an asymmetric periodic vibration regime even occurs for Reynolds number as low as 10. Based on the positions of the oscillation axes and the flow structures, the biased vibration regime is classified into two subregimes and the mechanisms of the biased vibration subregimes are analysed separately.
Efficiency improvement in a vibration power generator for a linear MR damper: numerical study
NASA Astrophysics Data System (ADS)
Sapiński, Bogdan; Krupa, Stanisław
2013-04-01
This paper summarizes a numerical analysis of the electromagnetic field, voltage and circuit properties and the cogging force in a vibration power generator comprising permanent magnets and a coil with a foil winding. The device converts the energy harvested from vibrations into electrical energy which is next used to vary the damping characteristics of a linear MR damper attached to the generator. The objective of the study is to propose a sufficiently efficient generator whose finally developed (target) version could be integrated with a small-scale MR damper to build a single device. Two design options for the device are numerically studied, the previously engineered generator 1 and the newly devised generator 2. Generator 1 incorporates two magnet systems having four magnets each and a single-section coil, while generator 2 comprises three magnet systems with four magnets each and a two-section coil. Calculations were performed to determine the electromagnetic field, voltage and current properties and the cogging force in the generators. The electromagnetic field parameters include the distribution of the magnetic field, the electrical potential field and the current density in the open turn and closed turn of the generators’ coils. The voltage and current properties include electromotive force (emf) in the generators and the voltage, current, instantaneous power and energy of the magnetic field in the MR damper control coil which is represented by resistance parameter R and inductance parameter L. The cogging force expresses the magnetic interactions between the permanent magnet systems and ferromagnetic structural components of the generators. The occurrence of this force is very unfavourable and attempts should be made to reduce it through control of the parameters of the magnetic circuit components. On one hand, comparison of the numerical results for the electromagnetic field parameters and voltage and current properties revealed that for the
[Changes in specific sensation in pilots exposed to systematic general vibration].
Podshivalov, A A; Krylov, Iu V; Zaritskiĭ, V V
1995-01-01
Helicopter pilots exposed to excessive general vibration demonstrate changes of specific sensation (vibrotactile, vestibular, auditory), that could be signs of occupational disorder. Those changes are increased thresholds of vibrotactile sensation, lower vestibulovegetative stability, changed vestibulospinal reflexes, more common occurrence of consistent deafness for voice frequencies in comparison with jet aircraft pilots. Experimental vibration (50-1,800 (m/s2) 2 hour) caused no changes of the vibrotactile and auditory sensation, and the modified vestibular function could prove the increased vestibular reactivity. PMID:7620795
NASA Astrophysics Data System (ADS)
Biffi, Carlo Alberto; Bassani, P.; Tuissi, A.; Carnevale, M.; Lecis, N.; LoConte, A.; Previtali, B.
2012-12-01
Shape memory alloys (SMAs) are very interesting smart materials not only for their shape memory and superelastic effects but also because of their significant intrinsic damping capacity. The latter is exhibited upon martensitic transformations and especially in martensitic state. The combination of these SMA properties with the mechanical and the lightweight of fiberglass-reinforced polymer (FGRP) is a promising solution for manufacturing of innovative composites for vibration suppression in structural applications. CuZnAl sheets, after laser patterning, were embedded in a laminated composite between a thick FGRP core and two thin outer layers with the aim of maximizing the damping capacity of the beam for passive vibration suppression. The selected SMA Cu66Zn24Al10 at.% was prepared by vacuum induction melting; the ingot was subsequently hot-and-cold rolled down to 0.2 mm thickness tape. The choice of a copper alloy is related to some advantages in comparison with NiTiCu SMA alloys, which was tested for the similar presented application in a previous study: lower cost, higher storage modulus and consequently higher damping properties in martensitic state. The patterning of the SMA sheets was performed by means of a pulsed fiber laser. After the laser processing, the SMA sheets were heat treated to obtain the desired martensitic state at room temperature. The transformation temperatures were measured by differential scanning calorimetry (DSC). The damping properties were determined, at room temperature, on full-scale sheet, using a universal testing machine (MTS), with cyclic tensile tests at different deformation amplitudes. Damping properties were also determined as a function of the temperature on miniature samples with a dynamical mechanical analyzer (DMA). Numerical modeling of the laminated composite, done with finite element method analysis and modal strain energy approaches, was performed to estimate the corresponding total damping capacity and then
NASA Technical Reports Server (NTRS)
Sharma, Surendra P.; Huo, Winifred M.; Park, Chul
1988-01-01
A theoretical study of vibrational excitations and dissociations of nitrogen undergoing a nonequilibrium relaxation process upon heating and cooling is reported. The rate coefficients for collisional induced vibrational transitions and transitions from a bound vibrational state into a dissociative state have been calculated using an extension of the theory originally proposed by Schwarz (SSH) et al. (1952). High-lying vibrational states and dissociative states were explicitly included but rotational energy transfer was neglected. The transition probabilities calculated from the SSH theory were fed into the master equation, which was integrated numerically to determine the population distribution of the vibrational states as well as bulk thermodynamic properties. The results show that: (1) the transition rates have a minimum near the middle of the bound vibrational levels, causing a bottleneck in the vibrational relaxation and dissociation rates; (2) high vibrational states are always in equilibrium with the dissociative state; (3) for the heating case, only the low vibrational states relax according to the Landau-Teller theory; (4) for the cooling case, vibrational relaxation cannot be described by a rate equation; (5) Park's (1985, 1988) two-temperature model is approximately valid; and (6) the average vibrational energy removed in dissociation is about 30 percent of the dissociation energy.
NASA Astrophysics Data System (ADS)
Blažević, D.; Zelenika, S.
2015-05-01
Scavenging of low-level ambient vibrations i.e. the conversion of kinetic into electric energy, is proven as effective means of powering low consumption electronic devices such as wireless sensor nodes. Cantilever based scavengers are characterised by several advantages and thus thoroughly investigated; analytical models based on a distributed parameter approach, Euler-Bernoulli beam theory and eigenvalue analysis have thus been developed and experimentally verified. Finite element models (FEM) have also been proposed employing different modelling approaches and commercial software packages with coupled analysis capabilities. An approach of using a FEM analysis of a piezoelectric cantilever bimorph under harmonic excitation is used in this work. Modal, harmonic and linear and nonlinear transient analyses are performed. Different complex dynamic effects are observed and compared to the results obtained by using a distributed parameter model. The influence of two types of finite elements and three mesh densities is also investigated. A complex bimorph cantilever, based on commercially available Midé Technology® Volture energy scavengers, is then considered. These scavengers are characterised by an intricate multilayer structure not investigated so far in literature. An experimental set-up is developed to evaluate the behaviour of the considered class of devices. The results of the modal and the harmonic FEM analyses of the behaviour of the multilayer scavengers are verified experimentally for three different tip masses and 12 different electrical load values. A satisfying agreement between numerical and experimental results is achieved.
A Numerical Simulator for Three-Dimensional Flows Through Vibrating Blade Rows
NASA Technical Reports Server (NTRS)
Chuang, H. Andrew; Verdon, Joseph M.
1998-01-01
The three-dimensional, multi-stage, unsteady, turbomachinery analysis, TURBO, has been extended to predict the aeroelastic and aeroacoustic response behaviors of a single blade row operating within a cylindrical annular duct. In particular, a blade vibration capability has been incorporated so that the TURBO analysis can be applied over a solution domain that deforms with a vibratory blade motion. Also, unsteady far-field conditions have been implemented to render the computational boundaries at inlet and exit transparent to outgoing unsteady disturbances. The modified TURBO analysis is applied herein to predict unsteady subsonic and transonic flows. The intent is to partially validate this nonlinear analysis for blade flutter applications, via numerical results for benchmark unsteady flows, and to demonstrate the analysis for a realistic fan rotor. For these purposes, we have considered unsteady subsonic flows through a 3D version of the 10th Standard Cascade, and unsteady transonic flows through the first stage rotor of the NASA Lewis, Rotor 67, two-stage fan.
Vortex-induced vibration of two parallel risers: Experimental test and numerical simulation
NASA Astrophysics Data System (ADS)
Huang, Weiping; Zhou, Yang; Chen, Haiming
2016-04-01
The vortex-induced vibration of two identical rigidly mounted risers in a parallel arrangement was studied using Ansys- CFX and model tests. The vortex shedding and force were recorded to determine the effect of spacing on the two-degree-of-freedom oscillation of the risers. CFX was used to study the single riser and two parallel risers in 2-8 D spacing considering the coupling effect. Because of the limited width of water channel, only three different riser spacings, 2 D, 3 D, and 4 D, were tested to validate the characteristics of the two parallel risers by comparing to the numerical simulation. The results indicate that the lift force changes significantly with the increase in spacing, and in the case of 3 D spacing, the lift force of the two parallel risers reaches the maximum. The vortex shedding of the risers in 3 D spacing shows that a variable velocity field with the same frequency as the vortex shedding is generated in the overlapped area, thus equalizing the period of drag force to that of lift force. It can be concluded that the interaction between the two parallel risers is significant when the risers are brought to a small distance between them because the trajectory of riser changes from oval to curve 8 as the spacing is increased. The phase difference of lift force between the two risers is also different as the spacing changes.
A General Framework for Multiphysics Modeling Based on Numerical Averaging
NASA Astrophysics Data System (ADS)
Lunati, I.; Tomin, P.
2014-12-01
In the last years, multiphysics (hybrid) modeling has attracted increasing attention as a tool to bridge the gap between pore-scale processes and a continuum description at the meter-scale (laboratory scale). This approach is particularly appealing for complex nonlinear processes, such as multiphase flow, reactive transport, density-driven instabilities, and geomechanical coupling. We present a general framework that can be applied to all these classes of problems. The method is based on ideas from the Multiscale Finite-Volume method (MsFV), which has been originally developed for Darcy-scale application. Recently, we have reformulated MsFV starting with a local-global splitting, which allows us to retain the original degree of coupling for the local problems and to use spatiotemporal adaptive strategies. The new framework is based on the simple idea that different characteristic temporal scales are inherited from different spatial scales, and the global and the local problems are solved with different temporal resolutions. The global (coarse-scale) problem is constructed based on a numerical volume-averaging paradigm and a continuum (Darcy-scale) description is obtained by introducing additional simplifications (e.g., by assuming that pressure is the only independent variable at the coarse scale, we recover an extended Darcy's law). We demonstrate that it is possible to adaptively and dynamically couple the Darcy-scale and the pore-scale descriptions of multiphase flow in a single conceptual and computational framework. Pore-scale problems are solved only in the active front region where fluid distribution changes with time. In the rest of the domain, only a coarse description is employed. This framework can be applied to other important problems such as reactive transport and crack propagation. As it is based on a numerical upscaling paradigm, our method can be used to explore the limits of validity of macroscopic models and to illuminate the meaning of
Balunov, V D; Barsukov, A F; Artamonova, V G
1998-01-01
The article covers complex evaluation of health state in building industry workers engaged into ferro-concrete production in St. Petersburg. The health state was considered under combined action of infrasound, noise and general vibration. Clinical and functional evaluation included medical examination by doctors, blood biochemistry and CBC, ECG, computer integral rheography, voice audiometry. Data for 62 moulders helped to assess acoustic environment at workplace and to reveal the morbidity structure. PMID:9662931
NASA Astrophysics Data System (ADS)
Lumentut, Mikail F.; Howard, Ian M.
2016-02-01
This paper focuses on the primary development of novel numerical and analytical techniques of the modal damped vibration energy harvesters with arbitrary proof mass offset. The key equations of electromechanical finite element discretisation using the extended Lagrangian principle are revealed and simplified to give matrix and scalar forms of the coupled system equations, indicating the most relevant numerical technique for the power harvester research. To evaluate the performance of the numerical study, the analytical closed-form boundary value equations have been developed using the extended Hamiltonian principle. The results from the electromechanical frequency response functions (EFRFs) derived from two theoretical studies show excellent agreement with experimental studies. The benefit of the numerical technique is in providing effective and quick predictions for analysing parametric designs and physical properties of piezoelectric materials. Although analytical technique provides a challenging process for analysing the complex smart structure, it shows complementary study for validating the numerical technique.
NASA Astrophysics Data System (ADS)
Nghiem, H. T. M.; Costi, T. A.
2014-02-01
The time-dependent numerical renormalization group (TDNRG) method [Anders et al., Phys. Rev. Lett. 95, 196801 (2005), 10.1103/PhysRevLett.95.196801] offers the prospect of investigating in a nonperturbative manner the time dependence of local observables of interacting quantum impurity models at all time scales following a quantum quench. Here, we present a generalization of this method to arbitrary finite temperature by making use of the full density matrix approach [Weichselbaum et al., Phys. Rev. Lett. 99, 076402 (2007), 10.1103/PhysRevLett.99.076402]. We show that all terms in the projected full density matrix ρi →f=ρ+++ρ--+ρ+-+ρ-+ appearing in the time evolution of a local observable may be evaluated in closed form at finite temperature, with ρ+-=ρ-+=0. The expression for ρ-- is shown to be finite at finite temperature, becoming negligible only in the limit of vanishing temperatures. We prove that this approach recovers the short-time limit for the expectation value of a local observable exactly at arbitrary temperatures. In contrast, the corresponding long-time limit is recovered exactly only for a continuous bath, i.e., when the logarithmic discretization parameter Λ →1+. Since the numerical renormalization group approach breaks down in this limit, and calculations have to be carried out at Λ >1, the long-time behavior following an arbitrary quantum quench has a finite error, which poses an obstacle for the method, e.g., in its application to the scattering-states numerical renormalization group method for describing steady-state nonequilibrium transport through correlated impurities [Anders, Phys. Rev. Lett. 101, 066804 (2008), 10.1103/PhysRevLett.101.066804]. We suggest a way to overcome this problem by noting that the time dependence, in general, and the long-time limit, in particular, become increasingly more accurate on reducing the size of the quantum quench. This suggests an improved generalized TDNRG approach in which the system is time
Thermally induced vibrations in a generalized thermoelastic solid with a cavity
Erbay, H.A.; Erbay, S.; Dost, S. )
1991-06-01
The present work deals with thermally induced vibrations in an infinite solid with a cavity. The medium is assumed to be linear, isotropic, temperature-rate-dependent thermoelastic. The problem is solved for the cases of cylindrical and spherical cavities. The surface of the cavity is assumed to be subjected to a temperature varying harmonically with time, and free of stress. For the cases considered, the coupled field equations admit exact solutions in terms of Hankel and the spherical Hankel functions, respectively. Numerical results are compared with those of classical thermoelasticity. The contribution of the second sound parameters in these problems becomes more significant as the frequency of applied temperature increases. 8 refs.
NASA Astrophysics Data System (ADS)
Tomioka, Takahiro; Takigami, Tadao
2015-05-01
The effects of passengers on vertical flexural vibrations of railway vehicle carbodies have been investigated experimentally and numerically. The primary focus is the damping effect caused by passengers. Vibration measurement tests, including stationary excitation tests and a running test on a commercial line, were conducted on three different types of actual railway vehicles for varying the numbers, postures and distributions of passengers. The measurement results showed that the peak frequencies in the measured FRF (frequency response function) and acceleration PSD (power spectral density) in response to excitation force or acceleration changed only a little; in contrast large reduction of the peak values was observed when passengers were aboard. These results show that passengers behave not as additional mass but as damping elements upon the carbody flexural vibrations. The damping effect by passengers appeared over several flexural modes and was significantly induced even by few passengers. Numerical studies were also carried out to develop analytical models for representing the measured results of the vibrations of carbody with passengers. It was shown that the change of the carbody FRF due to passengers could be simulated well by using a simple one-degree-of-freedom passenger model comprising a mass-spring-dashpot with a large damping ratio.
NASA Astrophysics Data System (ADS)
Zhu, Lihua; Yang, Qingxin; Yan, Rongge; Li, Yongjian; Zhang, Xian; Yan, Weili; Zhu, Jianguo
2013-05-01
Magnetostriction (MS) caused by the global magnetization of limbs and yokes and magnetic forces are the undisputed causes of the vibration and noise in power transformer cores. This paper presents a novel way to reduce the vibration and noise, in which nanocrystalline soft magnetic composite (NSMC) material with high permeability is used to fill the step-lap joint gaps of the power transformer magnetic cores. In order to numerically predict the effectiveness of the proposed method, a 3-D magneto-mechanical strong coupled model including MS and magnetic anisotropy of steel sheet was founded. Then, the numerical model was applied to analyze the step-lap joint region of the corner of magnetic cores. The analysis results illustrated that the deformation and noise of core with NSMC are lower than with the traditional epoxy damping material. Moreover, the validity of the proposed new way was verified by the simplified step-lap joint cores, which were achieved based on Epstein Frames.
Generalized Database Management System Support for Numeric Database Environments.
ERIC Educational Resources Information Center
Dominick, Wayne D.; Weathers, Peggy G.
1982-01-01
This overview of potential for utilizing database management systems (DBMS) within numeric database environments highlights: (1) major features, functions, and characteristics of DBMS; (2) applicability to numeric database environment needs and user needs; (3) current applications of DBMS technology; and (4) research-oriented and…
Natkaniec, I.; Holderna-Natkaniec, K.; Kalus, J.; Khavryutchenko, V. D.
1999-06-15
Vibrational densities of states of solid xylenes were determined from the inelastic neutron scattering spectra measured on the NERA spectrometer at the IBR-2 pulsed reactor. These spectra were used to test the semi-empirical quantum-chemistry calculations of internal vibrations of xylene molecules with differently deuterated sub-units. Rotations of methyl groups were found to be strongly affected by intermolecular interactions in the crystals and mixed with phenyl ring deformations.
Experimental and numerical study on vibration of the full-revolving propulsion ship stern
NASA Astrophysics Data System (ADS)
Liu, Chang-qing; Che, Chi-dong; Shen, Xiao-han
2015-03-01
In order to solve the severe vibration problems of an ocean engineering ship with a full-revolving propulsion system, the navigation tests, including forced vibration response test and modal test, are carried out in its stern. It is concluded from the comparison of the time-domain waveform and spectrum from different measurement points that three main factors lead to a high-level stern vibration. Firstly, the specific dynamic stiffness of a water tank is relatively small compared with its neighbor hold, which makes it act like a vibration isolator preventing vibrational energy transmitting to the main hold. Secondly, there exists high-density local modes in the working frequency range of the main engine and thus the local resonance occurs. Thirdly, the abnormal engagement of gears caused by the large deflection of the shaft bearing due to its low mounting rigidity leads to violent extra impulse excitations at high speeds. Then the modification against the dynamic defects is given by simply improving the specific stiffness of the water tanks. And the effect is validated by the FEM calculation. Some important experience is obtained with the problems being solved, which is useful in the design of ships with the same propulsion system. It is also believed that the dynamic consideration is as important as the static analysis for the ships, and that most of the vibration problems may be avoided with a proper acoustic design.
A numerical model for calculating vibration from a railway tunnel embedded in a full-space
NASA Astrophysics Data System (ADS)
Hussein, M. F. M.; Hunt, H. E. M.
2007-08-01
Vibration generated by underground railways transmits to nearby buildings causing annoyance to inhabitants and malfunctioning to sensitive equipment. Vibration can be isolated through countermeasures by reducing the stiffness of railpads, using floating-slab tracks and/or supporting buildings on springs. Modelling of vibration from underground railways has recently gained more importance on account of the need to evaluate accurately the performance of vibration countermeasures before these are implemented. This paper develops an existing model, reported by Forrest and Hunt, for calculating vibration from underground railways. The model, known as the Pipe-in-Pipe model, has been developed in this paper to account for anti-symmetrical inputs and therefore to model tangential forces at the tunnel wall. Moreover, three different arrangements of supports are considered for floating-slab tracks, one which can be used to model directly-fixed slabs. The paper also investigates the wave-guided solution of the track, the tunnel, the surrounding soil and the coupled system. It is shown that the dynamics of the track have significant effect on the results calculated in the wavenumber-frequency domain and therefore an important role on controlling vibration from underground railways.
NASA Astrophysics Data System (ADS)
Verbraken, H.; Lombaert, G.; Degrande, G.
2011-04-01
Vibrations induced by the passage of trains are a major environmental concern in urban areas. In practice, vibrations are often predicted using empirical methods such as the detailed vibration assessment procedure of the Federal Railroad Administration (FRA) of the U.S. Department of Transportation. This procedure allows predicting ground surface vibrations and re-radiated noise in buildings. Ground vibrations are calculated based on force densities, measured when a vehicle is running over a track, and line source transfer mobilities, measured on site to account for the effect of the local geology on wave propagation. Compared to parametric models, the advantage of this approach is that it inherently takes into account all important parameters. It can only be used, however, when an appropriate estimation of the force density is available. In this paper, analytical expressions are derived for the force density and the line source transfer mobility of the FRA procedure. The derivation of these expressions is verified using a coupled finite element-boundary element method.
A general numerical model for wave rotor analysis
NASA Technical Reports Server (NTRS)
Paxson, Daniel W.
1992-01-01
Wave rotors represent one of the promising technologies for achieving very high core temperatures and pressures in future gas turbine engines. Their operation depends upon unsteady gas dynamics and as such, their analysis is quite difficult. This report describes a numerical model which has been developed to perform such an analysis. Following a brief introduction, a summary of the wave rotor concept is given. The governing equations are then presented, along with a summary of the assumptions used to obtain them. Next, the numerical integration technique is described. This is an explicit finite volume technique based on the method of Roe. The discussion then focuses on the implementation of appropriate boundary conditions. Following this, some results are presented which first compare the numerical approximation to the governing differential equations and then compare the overall model to an actual wave rotor experiment. Finally, some concluding remarks are presented concerning the limitations of the simplifying assumptions and areas where the model may be improved.
NASA Astrophysics Data System (ADS)
Zhao, Ming; Yan, Guirong
2013-08-01
Two-degree-of-freedom Vortex-Induced Vibration (VIV) of two rigidly coupled circular cylinders of different diameters at a low Reynolds number of 250 is investigated numerically. While the diameter ratio and the mass ratio are kept constant, the study is focused on the effect of the position angle of the small cylinder on the lock-in regime of the VIV. Simulations are carried out for position angles α of the small cylinder ranging from 0° to 180° with an interval of 22.5° and the reduced velocities ranging from 1 to 15 with an increment of 1. In order to find the effect of the gap between the two cylinders on the vibration, two gap-to-diameter ratios (0 and 0.2) are considered. It is found that compared with a single cylinder case, the lock-in regime of the reduced velocity is widened significantly when the position angle of the small cylinder is α = 0°, 22.5°, 90°, or 112.5°. Pulsed beating phenomenon characterized by regular vibration with occasional high-amplitude disturbances at regular or irregular intervals is observed at G = 0 and α = 90°. At α = 135°, more than one lock-in regimes are observed in the computed range of reduced velocity for both gaps (G = 0 and 0.2). Setting a small gap (gap-to-diameter ratio of 0.2) between the two cylinders mitigates the vibration by narrowing the lock-in regime and reducing the vibration amplitude.
A GENERAL MASS-CONSERVATIVE NUMERICAL SOLUTION FOR THE UNSATURATED FLOW EQUATION
Numerical approximations based on different forms of the governing partial differential equation can lead to significantly different results for unsaturated flow problems. Numerical solution based on the standard h-based form of Richards equation generally yields poor results, ch...
Free Vibration of Simply Supported General Triangular Thin Plates: AN Accurate Simplified Solution
NASA Astrophysics Data System (ADS)
Saliba, H. T.
1996-09-01
In this paper, a highly accurate, simplified, and economical solution is provided for the free vibration problem of simply supported thin general triangular plates. The method is applicable to thin plates with linear boundaries regardless of there geometrical shapes. Results are compared with previously published reliable data for both the isosceles as well as the general triangles. Excellent agreements are reported. Eigenvalues are provided for a wide range of place aspect ratios. The first five mode shapes for some isosceles triangles are also provided for illustrative purposes. The advantages of the solution presented in the paper over previously published solutions are briefly discussed. Although the paper deals only with simple support conditions, it is mentioned that any combination of classical boundary conditions, with or without complicating factors, can easily be handled.
Numerical and semiclassical analysis of some generalized Casimir pistons
Schaden, M.
2009-05-15
The Casimir force due to a scalar field in a cylinder of radius r with a spherical cap of radius R>r is computed numerically in the world-line approach. A geometrical subtraction scheme gives the finite interaction energy that determines the Casimir force. The spectral function of convex domains is obtained from a probability measure on convex surfaces that is induced by the Wiener measure on Brownian bridges the convex surfaces are the hulls of. Due to reflection positivity, the vacuum force on the piston by a scalar field satisfying Dirichlet boundary conditions is attractive in these geometries, but the strength and short-distance behavior of the force depend strongly on the shape of the piston casing. For a cylindrical casing with a hemispherical head, the force on the piston does not depend on the dimension of the casing at small piston elevation a<
Numerical methods for a general class of porous medium equations
Rose, M. E.
1980-03-01
The partial differential equation par. deltau/par. deltat + par. delta(f(u))/par. deltax = par. delta(g(u)par. deltau/par. deltax)/par. deltax, where g(u) is a non-negative diffusion coefficient that may vanish for one or more values of u, was used to model fluid flow through a porous medium. Error estimates for a numerical procedure to approximate the solution are derived. A revised version of this report will appear in Computers and Mathematics with Applications.
Numerical implementation of generalized Coddington equations for ophthalmic lens design
NASA Astrophysics Data System (ADS)
Rojo, P.; Royo, S.; Ramírez, J.; Madariaga, I.
2014-02-01
A method for general implementation in any software platform of the generalized Coddington equations is presented, developed, and validated within a Matlab environment. The ophthalmic lens design strategy is presented thoroughly, and the basic concepts of generalized ray tracing are introduced. The methodology for ray tracing is shown to include two inter-related processes. Firstly, finite ray tracing is used to provide the main direction of propagation of the considered ray at the incidence point of interest. Afterwards, generalized ray tracing provides the principal curvatures of the local wavefront at that point, and its orientation after being refracted by the lens. The curvature values of the local wavefront are interpreted as the sagital and tangential powers of the lens at the point of interest. The proposed approach is validated using a double-check of the calculated lens performance in the spherical lens case: while finite ray tracing is validated using a commercial ray tracing software, generalized ray tracing is validated using a software application for ophthalmic lens design based on the classical version of Coddington equations. Equations of the complete tracing process are developed in detail for the case of generic astigmatic ophthalmic lenses as an example. Three-dimensional representation of the sagital and tangential powers of the ophthalmic lens at all directions of gaze then becomes possible, and results are presented for lenses with different geometries.
NASA Astrophysics Data System (ADS)
Kouroussis, G.; Bergeret, E.; Conti, C.; Verlinden, O.
With the development of continuous technological innovation, the railway transport is presented as an interesting alternative to the road traffic. Some drawbacks exist, one of the most problematic being certainly the vibrations induced by the railway traffic. The presented research wants to establish a reliable methodology in order to evaluate, from the design stage of a vehicle or of a track, the efforts transmitted by the vehicle to the track/soil system and consequently the level of vibrations in the surroundings. An analysis of the interaction between the track and the soil has been performed in order to show when the track/soil uncoupling can be assumed, with the aim of working in two stages. The first step is based on the vertical dynamic behaviour of the vehicle/track subsystem, taking into account any irregularity in the rail surface. For the soil subsystem (second step), recent publications showed that the finite/infinite element method can be an interesting alternative to boundary element method. The objective of this paper is to demonstrate the real benefit of the vehicle modelling in this kind of problem. Typical railway applications (Brussels tram, Thalys HST) are proposed, showing among others that significant reduction of ground vibration level can be obtained by modifying the dynamic characteristics of the vehicle.
NASA Astrophysics Data System (ADS)
Guilmineau, E.; Queutey, P.
2004-05-01
In this paper, we present some numerical results from a study of the dynamics and fluid forcing on an elastically mounted rigid cylinder with low mass-damping, constrained to oscillate transversely to a free stream. The vortex shedding around the cylinder is investigated numerically by the incompressible two-dimensional Reynolds-Averaged Navier-Stokes (RANS) equations. These equations are written in a primitive formulation in which the Cartesian velocity components and pressure share the same location at the center of the control volume. The numerical method uses a consistent physical reconstruction for the mass and momentum fluxes: the so-called consistent physical interpolation (CPI) approach in a conservative discretization using finite volumes on structured grids. The turbulence modeling is carried out by the SST K-ω model of Menter (AIAA 24th Fluid Dynamics Conference, Orlando, FL, USA). The numerical results are compared with the 1996 experimental results of Khalak and Williamson (J. Fluids Struct. 10 (1996) 455). The Reynolds number is in the range 900-15 000, the reduced velocity is including between 1.0 and 17.0. The mass ratio is 2.4 and the mass-damping is 0.013. Several initial conditions are used. According the initial condition used, the simulations predict correctly the maximum amplitude. On the other hand, the numerical results do not match the upper branch found experimentally. However, these results are encouraging, because no simulations have yet predicted such a high amplitude of vibration.
NASA Astrophysics Data System (ADS)
Bekhoucha, Ferhat; Rechak, Said; Duigou, Laëtitia; Cadou, Jean-Marc
2015-05-01
This paper deals with the computation of backbone curves bifurcated from a Hopf bifurcation point in the framework of nonlinear free vibrations of a rotating flexible beams. The intrinsic and geometrical equations of motion for anisotropic beams subjected to large displacements are used and transformed with Galerkin and harmonic balance methods to one quadratic algebraic equation involving one parameter, the pulsation. The latter is treated with the asymptotic numerical method using Padé approximants. An algorithm, equivalent to the Lyapunov-Schmidt reduction is proposed, to compute the bifurcated branches accurately from a Hopf bifurcation point, with singularity of co-rank 2, related to a conservative and gyroscopic dynamical system steady state, toward a nonlinear periodic state. Numerical tests dealing with clamped, isotropic and composite, rotating beams show the reliability of the proposed method reinforced by accurate results.
On a New Numerical Method for Solving General Variational Inequalities
NASA Astrophysics Data System (ADS)
Bnouhachem, Abdellah; Noor, Muhammad Aslam; Khalfaoui, Mohamed; Sheng, Zhaohan
In this paper, we suggest and analyze a new extragradient method for solving the general variational inequalities involving two operators. We also prove the global convergence of the proposed modified method under certain mild conditions. We used a self-adaptive technique to adjust parameter ρ at each iteration. It is proved theoretically that the lower-bound of the progress obtained by the proposed method is greater than that by the extragradient method. An example is given to illustrate the efficiency and its comparison with the extragradient method. Since the general variational inequalities include the classical variational inequalities and complementarity problems as special cases, our results obtained in this paper continue to hold for these problems. Results obtained in this paper may be viewed as an improvement and refinement of the previously known results in this field.
A general numerical analysis program for the superconducting quasiparticle mixer
NASA Technical Reports Server (NTRS)
Hicks, R. G.; Feldman, M. J.; Kerr, A. R.
1986-01-01
A user-oriented computer program SISCAP (SIS Computer Analysis Program) for analyzing SIS mixers is described. The program allows arbitrary impedance terminations to be specified at all LO harmonics and sideband frequencies. It is therefore able to treat a much more general class of SIS mixers than the widely used three-frequency analysis, for which the harmonics are assumed to be short-circuited. An additional program, GETCHI, provides the necessary input data to program SISCAP. The SISCAP program performs a nonlinear analysis to determine the SIS junction voltage waveform produced by the local oscillator. The quantum theory of mixing is used in its most general form, treating the large signal properties of the mixer in the time domain. A small signal linear analysis is then used to find the conversion loss and port impedances. The noise analysis includes thermal noise from the termination resistances and shot noise from the periodic LO current. Quantum noise is not considered. Many aspects of the program have been adequately verified and found accurate.
NASA Technical Reports Server (NTRS)
White, W. F., Jr.; Malatino, R. E.
1975-01-01
A method is presented for determining the free vibration characteristics of a rotating blade having nonuniform spanwise properties and cantilever boundary conditions. The equations which govern the coupled flapwise, chordwise, and torsional motion of such a blade are solved using an integrating matrix method. By expressing the equations of motion and matrix notation, utilizing the integrating matrix as an operator, and applying the boundary conditions, the equations are formulated into an eigenvalue problem whose solutions may be determined by conventional methods. Computer results are compared with experimental data.
WhiskyMHD: Numerical Code for General Relativistic Magnetohydrodynamics
NASA Astrophysics Data System (ADS)
Baiotti, Luca; Giacomazzo, Bruno; Hawke, Ian; et al.
2010-10-01
Whisky is a code to evolve the equations of general relativistic hydrodynamics (GRHD) and magnetohydrodynamics (GRMHD) in 3D Cartesian coordinates on a curved dynamical background. It was originally developed by and for members of the EU Network on Sources of Gravitational Radiation and is based on the Cactus Computational Toolkit. Whisky can also implement adaptive mesh refinement (AMR) if compiled together with Carpet. Whisky has grown from earlier codes such as GR3D and GRAstro_Hydro, but has been rewritten to take advantage of some of the latest research performed here in the EU. The motivation behind Whisky is to compute gravitational radiation waveforms for systems that involve matter. Examples would include the merger of a binary system containing a neutron star, which are expected to be reasonably common in the universe and expected to produce substantial amounts of radiation. Other possible sources are given in the projects list.
A numerical investigation of flow induced vibrations in a rocket engine manifold
NASA Astrophysics Data System (ADS)
Peugeot, John W.
2011-12-01
Flow induced vibrations are common in liquid rocket engine components and have been the subject of several recent studies within the Space Shuttle and Delta launch vehicle programs. Understanding how unsteady flow phenomena develop is important when investigating failures in existing hardware and in the design of new propulsion systems. In this study, a subsonic turbulent flow in a rocket engine manifold is analyzed using a compressible form of the viscous flow equations coupled with a hybrid RANS-DES turbulence model. It is found that vortex shedding and pressure perturbations within a manifold significantly influence the stability of shear layers and flow through exit cooling tubes. By adding a chamfer to the inlet of the cooling tubes, it was demonstrated that greater shear layer stability can be obtaIned at a given pressure ratio.
Transverse vibrations of shear-deformable beams using a general higher order theory
NASA Technical Reports Server (NTRS)
Kosmatka, J. B.
1993-01-01
A general higher order theory is developed to study the static and vibrational behavior of beam structures having an arbitrary cross section that utilizes both out-of-plane shear-dependent warping and in-plane (anticlastic) deformations. The equations of motion are derived via Hamilton's principle, where the full 3D constitutive relations are used. A simplified version of the general higher-order theory is also presented for beams having an arbitrary cross section that includes out-of-plane shear deformation but assumes that stresses within the cross section and in-plane deformations are negligible. This simplified model, which is accurate for long to moderately short wavelengths, offers substantial improvements over existing higher order theories that are limited to beams with thin rectangular cross sections. The current approach will be very useful in the study of thin-wall closed-cell beams such as airfoil-type sections where the magnitude of shear-related cross-sectional warping is significant.
NASA Technical Reports Server (NTRS)
Rosen, I. G.
1984-01-01
A cubic spline based Galerkin-like method is developed for the identification of a class of hybrid systems which describe the transverse vibration to flexible beams with attached tip bodies. The identification problem is formulated as a least squares fit to data subject to the system dynamics given by a coupled system of ordnary and partial differential equations recast as an abstract evolution equation (AEE) in an appropriate infinite dimensional Hilbert space. Projecting the AEE into spline-based subspaces leads naturally to a sequence of approximating finite dimensional identification problems. The solutions to these problems are shown to exist, are relatively easily computed, and are shown to, in some sense, converge to solutions to the original identification problem. Numerical results for a variety of examples are discussed.
Numerical investigation of nonlinear vibration for rotor-seal system of centrifugal pump
NASA Astrophysics Data System (ADS)
Zhou, W. J.; Yang, Y. C.; Xing, G. K.; Wang, L. Q.
2013-12-01
The exciting force in the seal is an important factor for the stability of a multiple stage centrifugal pump. With the speed increasing, the rotor system of multiple stage centrifugal pump presents some nonlinear characters. In order to provide supports for the research of nonlinear characters of multiple stage centrifugal pump, a rotor-seal system model of centrifugal pump is presented and the Muszynska nonlinear seal model is used to express the seal exciting force with multiple parameters in the paper. The fourth-order Runge-Kutta method is also used to determine the vibration response at the impeller place and obtain bifurcation diagram, axis orbit, phase diagram as well as Poincaré Map. The bifurcation results show that the rotor-seal system would be stable under a lower speed and change to be unstable as the rotor speed increases. Various multi-periodic motions and quasi-periodic motions are found showing the complicated motions in the rotor-seal system under nonlinear seal forces.
Experimental evidence, numerics, and theory of vibrational resonance in bistable systems.
Baltanás, J P; López, L; Blechman, I I; Landa, P S; Zaikin, A; Kurths, J; Sanjuán, M A F
2003-06-01
We consider an overdamped bistable oscillator subject to the action of a biharmonic force with very different frequencies, and study the response of the system when the parameters of the high-frequency force are varied. A resonantlike behavior is obtained when the amplitude or the frequency of this force is modified in an experiment performed by means of an analog circuit. This behavior, confirmed by numerical simulations, is explained on the basis of a theoretical approach. PMID:16241316
NASA Astrophysics Data System (ADS)
Jiang, Xiaohai; Lee, Taehun; Andreopoulos, Yiannis; Wang, Zhexuan
2013-11-01
Vortex-induced vibrations (VIV) phenomena related to self-excited energy harvesters consisting of circular or square cylinders have been investigated numerically by using the BGK or MRT Lattice Boltzmann Method. In the present work such a harvester is placed inside a channel flow and is allowed to oscillate without a structural restoring force in a direction normal to the flow. Currently the half-way bounce-back boundary scheme and interpolations are being used to model the moving boundary. The numerical results were compared to the ones by classical CFD methods and experiments. A good agreement was obtained. The vortex dynamics and the development of the flow patterns for different flow parameters such as Reynolds number, blockage and aspect ratios will be presented. Particular emphasis is given to the dynamics of vortex pairing observed in several of the simulations. The present approach will be extended to simulate the flexible beam with the Immersed Boundary Method. Sponsored by the National Science Foundation (CBET #1033117) and a fellowship support from China Scholarship Council.
Thomson, S L; Tack, J W; Verkerke, G J
2007-01-01
A computational model for exploring the design of a voice-producing voice prosthesis, or voice-producing element (VPE), is presented. The VPE is intended for use by laryngectomized patients who cannot benefit from current speech rehabilitation techniques. Previous experiments have focused on the design of a double-membrane voice generator as a VPE. For optimization studies, a numerical model has been developed. The numerical model introduced incorporates the finite element (FE) method to solve for the flow-induced vibrations of the VPE system, including airflow coupled with a mass-loaded membrane. The FE model includes distinct but coupled fluid and solid domains. The flow solver is governed by the incompressible, laminar, unsteady Navier-Stokes equations. The solid solver allows for large deformation, large strain, and collision. It is first shown that the model satisfactorily represents previously published experimental results in terms of frequency and flow rate, enabling the model for use as a design tool. The model is then used to study the influence of geometric scaling, membrane thickness, membrane stiffness, and slightly convergent or divergent channel geometry on the model response. It is shown that physiological allowable changes in the latter three device parameters alone will not be sufficient to generate the desired reduction in fundamental frequency. However, their effects are quantified and it is shown that membrane stiffness and included angle should be considered in future designs. PMID:17662296
Chapman, Craig T; Cheng, Xiaolu; Cina, Jeffrey A
2011-04-28
A recently framed quantum/semiclassical treatment for the internal nuclear dynamics of a small molecule and the induced small-amplitude coherent motion of a low-temperature host medium (Chapman, C. T.; Cina, J. A. J. Chem. Phys.2007,127, 114502) is further analyzed and subjected to initial tests of its numerical implementation. In the illustrative context of a 1D system interacting with a 1D medium, we rederive the fixed vibrational basis/gaussian bath (FVB/GB) equations of motion for the parameters defining the gaussian bath wave packet accompanying each of the energy eigenkets of the quantum mechanical system. The conditions of validity for the gaussian-bath approximation are shown to coincide with those supporting approximate population conservation. We perform initial numerical tests of the FVB/GB scheme and illustrate the semiclassical description it provides of coherent motion in the medium by comparing its predictions with the exact results for a high-frequency system harmonic oscillator bilinearly coupled to a lower-frequency bath oscillator. Linear vibronic absorption spectra or, equivalently, ultrafast wave packet interferometry signals are shown to be readily and accurately calculable within the FVB/GB framework. PMID:21462985
NASA Technical Reports Server (NTRS)
Fowlis, W. W. (Editor); Davis, M. H. (Editor)
1981-01-01
The atmospheric general circulation experiment (AGCE) numerical design for Spacelab flights was studied. A spherical baroclinic flow experiment which models the large scale circulations of the Earth's atmosphere was proposed. Gravity is simulated by a radial dielectric body force. The major objective of the AGCE is to study nonlinear baroclinic wave flows in spherical geometry. Numerical models must be developed which accurately predict the basic axisymmetric states and the stability of nonlinear baroclinic wave flows. A three dimensional, fully nonlinear, numerical model and the AGCE based on the complete set of equations is required. Progress in the AGCE numerical design studies program is reported.
NASA Astrophysics Data System (ADS)
Zhao, Ming
2013-05-01
Two-degree-of-freedom (2dof) vortex-induced vibration (VIV) of a circular cylinder in oscillatory flow is investigated numerically. The direction of the oscillatory flow is perpendicular to the spanwise direction of the circular cylinder. Simulations are carried out for the Keulegan-Carpenter (KC) numbers of 10, 20 and 40 and the Reynolds numbers ranging from 308 to 9240. The ratio of the Reynolds number to the reduced velocity is 308. At KC=10, the amplitude of the primary frequency component is much larger than those of other frequency components. Most vibrations for KC=20 and 40 have multiple frequencies. The primary frequency of the response in the cross-flow direction decreases with the increasing reduced velocity, except when the reduced velocity is very small. Because the calculated primary frequencies of the response in the cross-flow direction are multiple of the oscillatory flow frequency in most of the calculated cases, the responses are classified into single-frequency mode, double-frequency mode, triple frequency mode, etc. If the reduced velocity is in the range where the VIV is transiting from one mode to another, the vibration is very irregular.For each KC number the range of the reduced velocity can be divided into a cross-flow-in-phase regime (low Vr), where the response and the hydrodynamic force in the cross-flow direction synchronize, and a cross-flow-anti-phase regime (high Vr), where the response and the hydrodynamic force in the cross-flow direction are in anti-phase with each other. The boundary values of Vr between the cross-flow-in-phase and the cross-flow-anti-phase regimes are 7, 9 and 11 for KC=10, 20 and 40, respectively. For KC=20, another cross-flow-anti-phase regime is found between 15≤Vr≤19. Similarly the in-line-in-phase and the in-line-anti-phase regimes are also identified for the response in the in-line direction. It is found that the boundary value of Vr between the in-line-in-phase and the in-line-anti-phase regimes is
NASA Astrophysics Data System (ADS)
Zhu, Dibin; Roberts, Stephen; Mouille, Thomas; Tudor, Michael J.; Beeby, Stephen P.
2012-10-01
This paper presents a general model and its experimental validation for electrically tunable electromagnetic energy harvesters. Electrical tuning relies on the adjustment of the electrical load so that the maximum output power of the energy harvester occurs at a frequency which is different from the mechanical resonant frequency of the energy harvester. Theoretical analysis shows that for this approach to be feasible the electromagnetic vibration energy harvester’s coupling factor must be maximized so that its resonant frequency can be tuned with the minimum decrease of output power. Two different-sized electromagnetic energy harvesters were built and tested to validate the model. Experimentally, the micro-scale energy harvester has a coupling factor of 0.0035 and an untuned resonant frequency of 70.05 Hz. When excited at 30 mg, it was tuned by 0.23 Hz by changing its capacitive load from 0 to 4000 nF its effective tuning range is 0.15 Hz for a capacitive load variation from 0 to 1500 nF. The macro-scale energy harvester has a coupling factor of 552.25 and an untuned resonant frequency of 95.1 Hz and 95.5 Hz when excited at 10 mg and 25 mg, respectively. When excited at 10 mg, it was tuned by 3.8 Hz by changing its capacitive load from 0 to 1400 nF it has an effective tuning range of 3.5 Hz for a capacitive load variation from 0 to 1200 nF. When excited at 25 mg, its resonant frequency was tuned by 4.2 Hz by changing its capacitive load from 0 to 1400 nF it has an effective tuning range of about 5 Hz. Experimental results were found to agree with the theoretical analysis to within 10%.
Modal wavefront reconstruction over general shaped aperture by numerical orthogonal polynomials
NASA Astrophysics Data System (ADS)
Ye, Jingfei; Li, Xinhua; Gao, Zhishan; Wang, Shuai; Sun, Wenqing; Wang, Wei; Yuan, Qun
2015-03-01
In practical optical measurements, the wavefront data are recorded by pixelated imaging sensors. The closed-form analytical base polynomial will lose its orthogonality in the discrete wavefront database. For a wavefront with an irregularly shaped aperture, the corresponding analytical base polynomials are laboriously derived. The use of numerical orthogonal polynomials for reconstructing a wavefront with a general shaped aperture over the discrete data points is presented. Numerical polynomials are orthogonal over the discrete data points regardless of the boundary shape of the aperture. The performance of numerical orthogonal polynomials is confirmed by theoretical analysis and experiments. The results demonstrate the adaptability, validity, and accuracy of numerical orthogonal polynomials for estimating the wavefront over a general shaped aperture from regular boundary to an irregular boundary.
NASA Astrophysics Data System (ADS)
Bian, Xuecheng; Chen, Yunmin; Hu, Ting
2008-06-01
An efficient 2.5D finite element numerical modeling approach was developed to simulate wave motions generated in ground by high-speed train passages. Fourier transform with respect to the coordinate in the track direction was applied to reducing the three-dimensional dynamic problem to a plane strain problem which has been solved in a section perpendicular to the track direction. In this study, the track structure and supporting ballast layer were simplified as a composite Euler beam resting on the ground surface, while the ground with complicated geometry and physical properties was modeled by 2.5D quadrilateral elements. Wave dissipation into the far field was dealt with the transmitting boundary constructed with frequency-dependent dashpots. Three-dimensional responses of track structure and ground were obtained from the wavenumber expansion in the track direction. The simulated wave motions in ground were interpreted for train moving loads traveling at speeds below or above the critical velocity of a specific track-ground system. It is found that, in the soft ground area, the high-speed train operations can enter the transonic range, which can lead to resonances of the track structure and the supporting ground. The strong vibration will endanger the safe operations of high-speed train and accelerate the deterioration of railway structure.
NASA Astrophysics Data System (ADS)
Wang, Yanrong; Liu, Bin; Tian, Aimei; Tang, Wei
2016-06-01
Particle damping (PD) has been well known for its simplicity and high efficiency in attenuating structure vibration. Recent studies on PD have focused mainly on new types of dampers and applications. Meanwhile, excitation applied to the primary structure is still limited to either horizontal or vertical direction, perpendicular or parallel to gravity. In this study, the characteristics of PD under horizontal-vertical excitations (HVE) are investigated numerically and experimentally. The particle damper, which is attached to the top free end of an L-shaped cantilever beam, is simultaneously excited in the horizontal and vertical directions in the context of free decay. An equivalent model capable of motion in both the horizontal and vertical directions is generated. Given an initial displacement disturbance, this model starts vibrating freely in the vertical plane. A code based on the 3D discrete element method is programmed, and the high coincidence between the numerical and experimental results shows that this equivalent model is capable of high-fidelity simulation for PD under HVE. Parametric studies have been implemented to characterize the basic nonlinear damping capacity of particle dampers under this new operating condition. The effects of seven dimensionless independent parameters on the specific damping capacity (SDC) are investigated, including dimensionless acceleration amplitude, particle mass ratio, dimensionless horizontal and vertical impact clearances, coefficients of friction and restitution, and amplitude ratio of the horizontal excitation to the vertical excitation. The results show that the basic damping properties of PD under HVE are similar to those of PD under only vertical excitation. However, PD under HVE signifies its own characteristics because of the existence of horizontal excitation: (1) The impact clearances in both the horizontal and vertical directions have significant effects on the SDC because of the significant increase in oblique
Generalized thick strip modelling for vortex-induced vibration of long flexible cylinders
NASA Astrophysics Data System (ADS)
Bao, Y.; Palacios, R.; Graham, M.; Sherwin, S.
2016-09-01
We propose a generalized strip modelling method that is computationally efficient for the VIV prediction of long flexible cylinders in three-dimensional incompressible flow. In order to overcome the shortcomings of conventional strip-theory-based 2D models, the fluid domain is divided into "thick" strips, which are sufficiently thick to locally resolve the small scale turbulence effects and three dimensionality of the flow around the cylinder. An attractive feature of the model is that we independently construct a three-dimensional scale resolving model for individual strips, which have local spanwise scale along the cylinder's axial direction and are only coupled through the structural model of the cylinder. Therefore, this approach is able to cover the full spectrum for fully resolved 3D modelling to 2D strip theory. The connection between these strips is achieved through the calculation of a tensioned beam equation, which is used to represent the dynamics of the flexible body. In the limit, however, a single "thick" strip would fill the full 3D domain. A parallel Fourier spectral/hp element method is employed to solve the 3D flow dynamics in the strip-domain, and then the VIV response prediction is achieved through the strip-structure interactions. Numerical tests on both laminar and turbulent flows as well as the comparison against the fully resolved DNS are presented to demonstrate the applicability of this approach.
NASA Astrophysics Data System (ADS)
Wildman, R. D.; Jenkins, J. T.; Krouskop, P. E.; Talbot, J.
2006-07-01
A comparison of the predictions of a simple kinetic theory with experimental and numerical results for a vibrated granular bed consisting of nearly elastic particles of two sizes has been performed. The results show good agreement between the data sets for a range of numbers of each size of particle, and are particularly good for particle beds containing similar proportions of each species. The agreement suggests that such a model may be a good starting point for describing polydisperse systems of granular flows.
Drolshagen, G.; Mayne, H.R.; Toennies, J.P.
1981-07-01
We extend the theory of inelastic rainbows to include vibrationally inelastic scattering, showing how the existence of vibrational rainbows can be deduced from collinear classical scattering theory. Exact close-coupling calculations are carried out for a breathing sphere potential, and rainbow structures are, in fact, observed. The location of the rainbows generally agrees well with the classical prediction. In addition, the sensitivity of the location of the rainbow to changes in the vibrational coupling has been investigated. It is shown that vibrational rainbows persist in the presence of anisotropy. Experimental results (R. David, M. Faubel, and J. P. Toennies, Chem. Phys. Lett. 18, 87 (1973)) are examined for evidence of vibrational rainbow structure, and it is shown that vibrational rainbow theory is not inconsistent with these results.
NASA Astrophysics Data System (ADS)
Takada, Shoji; Nakamura, Hiroki
1995-03-01
Tunneling energy splittings of vibrationally excited states are calculated quantum mechanically using several models of two-dimensional symmetric double well potentials. Various effects of vibrational excitation on tunneling are found to appear, depending on the topography of potential energy surface; the symmetry of the mode coupling plays an essential role. Especially, oscillation of tunneling splitting with respect to vibrational quantum number can occur and is interpreted by a clear physical picture based on the semiclassical theory formulated recently [Takada and Nakamura, J. Chem. Phys. 100, 98 (1994)]. The mixed tunneling in the C region found there allows the wave functions to have nodal lines in classically inaccessible region and can cause the suppression of the tunneling. The above analysis is followed by the interpretation of recent experiments of proton tunneling in tropolone. Ab initio molecular orbital calculations are carried out for the electronically ground state. A simple three-dimensional model potential is constructed and employed to analyze the proton tunneling dynamics. Some of the experimentally observed intriguing features can be explained by the typical mechanisms discussed above.
Numerically pricing American options under the generalized mixed fractional Brownian motion model
NASA Astrophysics Data System (ADS)
Chen, Wenting; Yan, Bowen; Lian, Guanghua; Zhang, Ying
2016-06-01
In this paper, we introduce a robust numerical method, based on the upwind scheme, for the pricing of American puts under the generalized mixed fractional Brownian motion (GMFBM) model. By using portfolio analysis and applying the Wick-Itô formula, a partial differential equation (PDE) governing the prices of vanilla options under the GMFBM is successfully derived for the first time. Based on this, we formulate the pricing of American puts under the current model as a linear complementarity problem (LCP). Unlike the classical Black-Scholes (B-S) model or the generalized B-S model discussed in Cen and Le (2011), the newly obtained LCP under the GMFBM model is difficult to be solved accurately because of the numerical instability which results from the degeneration of the governing PDE as time approaches zero. To overcome this difficulty, a numerical approach based on the upwind scheme is adopted. It is shown that the coefficient matrix of the current method is an M-matrix, which ensures its stability in the maximum-norm sense. Remarkably, we have managed to provide a sharp theoretic error estimate for the current method, which is further verified numerically. The results of various numerical experiments also suggest that this new approach is quite accurate, and can be easily extended to price other types of financial derivatives with an American-style exercise feature under the GMFBM model.
Harleston, H. Instituto de Ciencias Nucleares, Universidad Nacional Autonoma de Mexico, Distrito Federal ); Vishniac, E.T. )
1992-06-15
The Arnowitt-Deser-Misner formalism is used to write the Einstein-Boltzmann coupled system of equations. The sources of gravitational field are represented by ordinary matter described by a perfect-fluid approximation together with a particle gas described by a phase-space distribution function obeying the general-relativistic Boltzmann transport equation. Through the use of the Liouville operator in phase space, we obtain a form of the Boltzmann equation that makes it very amenable for numerical treatment. The resulting system of equations can be used for the numerical study of either massless or massive particles interacting with ordinary matter.
Steady and Unsteady Numerical Solution of Generalized Newtonian Fluids Flow by Runge-Kutta method
NASA Astrophysics Data System (ADS)
Keslerová, R.; Kozel, K.; Prokop, V.
2010-09-01
In this paper the laminar viscous incompressible flow for generalized Newtonian (Newtonian and non-Newtonian) fluids is considered. The governing system of equations is the system of Navier-Stokes equations and the continuity equation. The steady and unsteady numerical solution for this system is computed by finite volume method combined with an artificial compressibility method. For time discretization the explicit multistage Runge-Kutta numerical scheme is considered. Steady state solution is achieved for t→∞ using steady boundary conditions and followed by steady residual behavior. The dual time-stepping method is considered for unsteady computation. The high artificial compressibility coefficient is used in the artificial compressibility method applied in the dual time τ. The steady and unsteady numerical results of Newtonian and non-Newtonian (shear thickening and shear thinning) fluids flow in the branching channel are presented.
Numerical studies of constraints and gravitational wave extraction in general relativity
NASA Astrophysics Data System (ADS)
Fiske, David Robert
Within classical physics, general relativity is the theory of gravity. Its equations are non-linear partial differential equations for which relatively few closed form solutions are known. Because of the growing observational need for solutions representing gravitational waves from astrophysically plausible sources, a subfield of general relativity; numerical relativity, has a emerged with the goal of generating numerical solutions to the Einstein equations. This dissertation focuses on two fundamental problems in modern numerical relativity: (1)Creating a theoretical treatment of the constraints in the presence of constraint-violating numerical errors, and (2)Designing and implementing an algorithm to compute the spherical harmonic decomposition of radiation quantities for comparison with observation. On the issue of the constraints, I present a novel and generic procedure for incorporating the constraints into the equations of motion of the theory in a way designed to make the constraint hypersurface an attractor of the evolution. In principle, the prescription generates non- linear corrections for the Einstein equations. The dissertation presents numerical evidence that the correction terms do work in the case of two formulations of the Maxwell equations and two formulations of the linearized Einstein equations. On the issue of radiation extraction, I provide the first in-depth analysis of a novel algorithm, due originally to Misner, for computing spherical harmonic components on a cubic grid. I compute explicitly how the truncation error in the algorithm depends on its various parameters, and I also provide a detailed analysis showing how to implement the method on grids in which explicit symmetries are enforced via boundary conditions. Finally, I verify these error estimates and symmetry arguments with a numerical study using a solution of the linearized Einstein equations known as a Teukolsky wave. The algorithm performs well and the estimates prove true both
Numerical solution of shock and ramp compression for general material properties
Swift, D C
2009-01-28
A general formulation was developed to represent material models for applications in dynamic loading. Numerical methods were devised to calculate response to shock and ramp compression, and ramp decompression, generalizing previous solutions for scalar equations of state. The numerical methods were found to be flexible and robust, and matched analytic results to a high accuracy. The basic ramp and shock solution methods were coupled to solve for composite deformation paths, such as shock-induced impacts, and shock interactions with a planar interface between different materials. These calculations capture much of the physics of typical material dynamics experiments, without requiring spatially-resolving simulations. Example calculations were made of loading histories in metals, illustrating the effects of plastic work on the temperatures induced in quasi-isentropic and shock-release experiments, and the effect of a phase transition.
Chinyoka, T.; Makinde, O. D.
2013-01-01
The thermodynamic second law analysis is utilized to investigate the inherent irreversibility in an unsteady hydromagnetic generalized Couette flow with variable electrical conductivity in the presence of induced electric field. Based on some simplified assumption, the model nonlinear governing equations are obtained and solved numerically using semidiscretization finite difference techniques. Effects of various thermophysical parameters on the fluid velocity, temperature, current density, skin friction, the Nusselt number, entropy generation number, and the Bejan number are presented graphically and discussed quantitatively. PMID:23956691
NASA Astrophysics Data System (ADS)
Sakaris, Christos S.; Sakellariou, John S.; Fassois, Spilios D.
2015-07-01
A Generalized Functional Model Based Method for vibration-based damage precise localization on structures consisting of 1D, 2D, or 3D elements is introduced. The method generalizes previous versions applicable to structures consisting of 1D elements, thus allowing for 2D and 3D elements as well. It is based on scalar (single sensor) or vector (multiple sensor) Functional Models which - in the inspection phase - incorporate the mathematical form of the specific structural topology. Precise localization is then based on coordinate estimation within this model structure, and confidence bounds are also obtained. The effectiveness of the method is demonstrated through experiments on a 3D truss structure where damage corresponds to single bolt loosening. Both the scalar and vector versions of the method are shown to be effective even within a very limited, low frequency, bandwidth of 3-59 Hz. The improvement achieved through the use of multiple sensors is also demonstrated.
Jendrzejczyk, J.A.; Wambsganss, M.W.; Smith, R.K.
1993-03-01
The primary objective of this study was to measure the coupling between the experimental hall structure and the storage ring tunnel. The study was motivated in part by concern over the potential for wind excitation of the building to couple through the ground and/or piping/ductwork to the storage ring basemat. Unfortunately, there was very little wind on the measurement date. However, it can be assumed that the wind turbulence would excite the experimental hall structure at one or more of its resonant frequencies. Similar coupling information can be obtained by using ambient vibration as the excitation source and measuring the resultant structure/basemat response and associated coupling. This was the approach taken in this study.
Dong, S.
2015-02-15
We present a family of physical formulations, and a numerical algorithm, based on a class of general order parameters for simulating the motion of a mixture of N (N⩾2) immiscible incompressible fluids with given densities, dynamic viscosities, and pairwise surface tensions. The N-phase formulations stem from a phase field model we developed in a recent work based on the conservations of mass/momentum, and the second law of thermodynamics. The introduction of general order parameters leads to an extremely strongly-coupled system of (N−1) phase field equations. On the other hand, the general form enables one to compute the N-phase mixing energy density coefficients in an explicit fashion in terms of the pairwise surface tensions. We show that the increased complexity in the form of the phase field equations associated with general order parameters in actuality does not cause essential computational difficulties. Our numerical algorithm reformulates the (N−1) strongly-coupled phase field equations for general order parameters into 2(N−1) Helmholtz-type equations that are completely de-coupled from one another. This leads to a computational complexity comparable to that for the simplified phase field equations associated with certain special choice of the order parameters. We demonstrate the capabilities of the method developed herein using several test problems involving multiple fluid phases and large contrasts in densities and viscosities among the multitude of fluids. In particular, by comparing simulation results with the Langmuir–de Gennes theory of floating liquid lenses we show that the method using general order parameters produces physically accurate results for multiple fluid phases.
ERIC Educational Resources Information Center
Fuchs, Lynn S.; Geary, David C.; Compton, Donald L.; Fuchs, Douglas; Hamlett, Carol L.; Seethaler, Pamela M.; Bryant, Joan D.; Schatschneider, Christopher
2010-01-01
The purpose of this study was to examine the interplay between basic numerical cognition and domain-general abilities (such as working memory) in explaining school mathematics learning. First graders (N = 280; mean age = 5.77 years) were assessed on 2 types of basic numerical cognition, 8 domain-general abilities, procedural calculations, and word…
A fully general relativistic numerical simulation code for spherically symmetric matter
NASA Astrophysics Data System (ADS)
Park, Dong-Ho; Cho, Inyong; Kang, Gungwon; Lee, Hyung Mok
2013-02-01
We present a fully general relativistic open-source code that can be used for simulating a system of spherically symmetric perfect fluid matter. It is based on the Arnowitt-Deser-Misner 3+1 formalism with maximal slicing and isotropic spatial coordinates. For hydrodynamic matter High Resolution Shock Capturing (HRSC) schemes with a monotonized central-difference limiter and approximated Riemann solvers are used in the Eulerian viewpoint. The accuracy and the convergence of our numerical code are verified by performing several test problems. These include a relativistic blast wave, relativistic spherical accretion of matter into a black hole, Tolman-Oppenheimer-Volkoff (TOV) stars and Oppenheimer-Snyder (OS) dust collapses. In particular, a dynamical code test is done for the OS collapse by explicitly performing numerical coordinate transformations between our coordinate 8system and the one used for the analytic solution. Finally, some TOV star solutions are presented for the Eddington-inspired Born-Infeld gravity theory.
Generalization Evaluation of Machine Learning Numerical Observers for Image Quality Assessment.
Kalayeh, Mahdi M; Marin, Thibault; Brankov, Jovan G
2013-06-01
In this paper, we present two new numerical observers (NO) based on machine learning for image quality assessment. The proposed NOs aim to predict human observer performance in a cardiac perfusion-defect detection task for single-photon emission computed tomography (SPECT) images. Human observer (HumO) studies are now considered to be the gold standard for task-based evaluation of medical images. However such studies are impractical for use in early stages of development for imaging devices and algorithms, because they require extensive involvement of trained human observers who must evaluate a large number of images. To address this problem, numerical observers (also called model observers) have been developed as a surrogate for human observers. The channelized Hotelling observer (CHO), with or without internal noise model, is currently the most widely used NO of this kind. In our previous work we argued that development of a NO model to predict human observers' performance can be viewed as a machine learning (or system identification) problem. This consideration led us to develop a channelized support vector machine (CSVM) observer, a kernel-based regression model that greatly outperformed the popular and widely used CHO. This was especially evident when the numerical observers were evaluated in terms of generalization performance. To evaluate generalization we used a typical situation for the practical use of a numerical observer: after optimizing the NO (which for a CHO might consist of adjusting the internal noise model) based upon a broad set of reconstructed images, we tested it on a broad (but different) set of images obtained by a different reconstruction method. In this manuscript we aim to evaluate two new regression models that achieve accuracy higher than the CHO and comparable to our earlier CSVM method, while dramatically reducing model complexity and computation time. The new models are defined in a Bayesian machine-learning framework: a channelized
NASA Astrophysics Data System (ADS)
Wang, Qingshan; Shi, Dongyan; Pang, Fuzhen; Liang, Qian
2016-04-01
A Fourier-Ritz method for predicting the free vibration of composite laminated circular panels and shells of revolution subjected to various combinations of classical and non-classical boundary conditions is presented in this paper. A modified Fourier series approach in conjunction with a Ritz technique is employed to derive the formulation based on the first-order shear deformation theory. The general boundary condition can be achieved by the boundary spring technique in which three types of liner and two types of rotation springs along the edges of the composite laminated circular panels and shells of revolution are set to imitate the boundary force. Besides, the complete shells of revolution can be achieved by using the coupling spring technique to imitate the kinematic compatibility and physical compatibility conditions of composite laminated circular panels at the common meridian with θ = 0 and 2π. The comparisons established in a sufficiently conclusive manner show that the present formulation is capable of yielding highly accurate solutions with little computational effort. The influence of boundary and coupling restraint parameters, circumference angles, stiffness ratios, numbers of layer and fiber orientations on the vibration behavior of the composite laminated circular panels and shells of revolution are also discussed.
NASA Astrophysics Data System (ADS)
Davoodi, H.; Noori, M.
1990-07-01
The work presented in this paper constitutes the second phase of on-going research aimed at developing mathematical models for representing general hysteretic behavior of structures and approximation techniques for the computation and analysis of the response of hysteretic systems to random excitations. In this second part, the technique previously developed by the authors for the Gaussian response analysis of non-linear systems with general hysteretic behavior is extended for the non-Gaussian analysis of these systems. This approximation technique is based on the approach proposed independently by Ibrahim and Wu-Lin. In this work up to fourth order moments of the response co-ordinates are obtained for the Bouc-Baber-Wen smooth hysteresis model. These higher order statistics previously have not been made available for general hysteresis models by using existing approximation methods. Second order moments obtained for the model by this non-Gaussian closure scheme are compared with equivalent linearization and Gaussian closure results via Monte Carlo simulation (MCS). Higher order moments are compared with the simulation results. The study performed for a wide range of degradation parameters and input power spectral density ( PSD) levels shows that the non-Gaussian responses obtained by this approach are in better agreement with the MCS results than the linearized and Gaussian ones. This approximation technique can provide information on higher order moments for general hysteretic systems. This information is valuable in random vibration and the reliability analysis of hysteretically yielding structures.
NASA Astrophysics Data System (ADS)
Xiong, Xingyu; Oyadiji, S. Olutunde
2014-10-01
Piezoelectric vibration energy harvesters with multi-layer stacked structures have been developed. They consist of multi-layer beams, of zigzag configurations, with rigid masses attached between the beams. The rigid masses, which also serve as spacers, are attached to each layer to tune the frequencies of the harvester. Close resonance frequencies and considerable power output can be achieved in multiple modes by varying the positions of the masses. A modal approach is introduced to determine the modal performance conveniently using the mass ratio and the modal electromechanical coupling coefficient, and the required modal parameters are derived using the finite element method. Mass ratio represents the influence of modal mechanical behaviour on the power density. Since the modes with larger mass ratios cause the remaining modes to have smaller mass ratios and lower power densities, a screening process using the modal approach is developed to determine the optimal or near-optimal performance of the harvesters when altering mass positions. This procedure obviates the need for full analysis by pre-selecting the harvester configurations with close resonances and favourable values of mass ratio initially. Furthermore, the multi-layer stacked designs using the modal approach can be used to develop harvesters with different sizes with the power ranging from microwatts to milliwatts.
NASA Astrophysics Data System (ADS)
Nagy, Péter R.; Surján, Péter R.; Szabados, Ágnes
2014-01-01
Cross sections of inelastic light scattering accompanied by vibronic excitation in large conjugated carbon structures is assessed at the π-electron level. Intensities of Raman and vibrational Raman optical activity (VROA) spectra of fullerenes are computed, relying on a single electron per atom. When considering only first neighbor terms in the Hamiltonian (a tight-binding (TB) type or Hückel-model), Raman intensities are captured remarkably well, based on comparison with frequency-dependent linear response of the self-consistent field (SCF) method. Resorting to π-electron levels when computing spectral intensities brings a beneficial reduction in computational cost as compared to linear response SCF. At difference with total intensities, the first neighbor TB model is found inadequate for giving the left and right circularly polarized components of the scattered light, especially when the molecular surface is highly curved. To step beyond first neighbor approximation, an effective π-electron Hamiltonian, including interaction of all sites is derived from the all-electron Fockian, in the spirit of the Bloch-equation. Chiroptical cross-sections computed by this novel π-electron method improve upon first-neighbor TB considerably, with no increase in computational cost. Computed VROA spectra of chiral fullerenes, such as C76 and C28, are reported for the first time, both by conventional linear response SCF and effective π-electron models.
NASA Technical Reports Server (NTRS)
Zhang, Y. C.; Zhang, J. Z. H.; Kouri, D. J.; Haug, K.; Schwenke, D. W.
1988-01-01
Numerically exact, fully three-dimensional quantum mechanicl reactive scattering calculations are reported for the H2Br system. Both the exchange (H + H-prime Br to H-prime + HBr) and abstraction (H + HBR to H2 + Br) reaction channels are included in the calculations. The present results are the first completely converged three-dimensional quantum calculations for a system involving a highly exoergic reaction channel (the abstraction process). It is found that the production of vibrationally hot H2 in the abstraction reaction, and hence the extent of population inversion in the products, is a sensitive function of initial HBr rotational state and collision energy.
Code of Federal Regulations, 2013 CFR
2013-04-01
... the structure to withstand transportation shock and vibration. 3280.903 Section 3280.903 Housing and... structure to withstand transportation shock and vibration. (a) The cumulative effect of highway transportation shock and vibration upon a manufactured home structure may result in incremental degradation...
Code of Federal Regulations, 2010 CFR
2010-04-01
... the structure to withstand transportation shock and vibration. 3280.903 Section 3280.903 Housing and... structure to withstand transportation shock and vibration. (a) The cumulative effect of highway transportation shock and vibration upon a manufactured home structure may result in incremental degradation...
Code of Federal Regulations, 2011 CFR
2011-04-01
... the structure to withstand transportation shock and vibration. 3280.903 Section 3280.903 Housing and... structure to withstand transportation shock and vibration. (a) The cumulative effect of highway transportation shock and vibration upon a manufactured home structure may result in incremental degradation...
Code of Federal Regulations, 2014 CFR
2014-04-01
... the structure to withstand transportation shock and vibration. 3280.903 Section 3280.903 Housing and... structure to withstand transportation shock and vibration. (a) The cumulative effect of highway transportation shock and vibration upon a manufactured home structure may result in incremental degradation...
An efficient numerical method for general L(p) regularization in fluorescence molecular tomography.
Baritaux, Jean-Charles; Hassler, Kat; Unser, Michael
2010-04-01
Reconstruction algorithms for fluorescence tomography have to address two crucial issues: 1) the ill-posedness of the reconstruction problem, 2) the large scale of numerical problems arising from imaging of 3-D samples. Our contribution is the design and implementation of a reconstruction algorithm that incorporates general Lp regularization (p ¿ 1). The originality of this work lies in the application of general Lp constraints to fluorescence tomography, combined with an efficient matrix-free strategy that enables the algorithm to deal with large reconstruction problems at reduced memory and computational costs. In the experimental part, we specialize the application of the algorithm to the case of sparsity promoting constraints (L (1)). We validate the adequacy of L (1) regularization for the investigation of phenomena that are well described by a sparse model, using data acquired during phantom experiments. PMID:20236875
Analytical and Numerical Solutions of a Generalized Hyperbolic Non-Newtonian Fluid Flow
NASA Astrophysics Data System (ADS)
Pakdemirli, Mehmet; Sarı, Pınar; Solmaz, Bekir
2010-03-01
The generalized hyperbolic non-Newtonian fluid model first proposed by Al-Zahrani [J. Petroleum Sci. Eng. 17, 211 (1997)] is considered. This model was successfully applied to some drilling fluids with a better performance in relating shear stress and velocity gradient compared to power-law and the Hershel-Bulkley model. Special flow geometries namely pipe flow, parallel plate flow, and flow between two rotating cylinders are treated. For the first two cases, analytical solutions of velocity profiles and discharges in the form of integrals are presented. These quantities are calculated by numerically evaluating the integrals. For the flow between two rotating cylinders, the differential equation is solved by the Runge-Kutta method combined with shooting. For all problems, the power-law approximation of the model is compared with the generalized hyperbolic model, too.
NASA Astrophysics Data System (ADS)
Polome, J.
2012-07-01
This paper presents the performance methodology and the achieved results of a vibration test campaign. The test performed by shaker simulates an unusual acoustic excitation containing a harmonic serial superimposed on random noise seen by electronic flight equipment. The paper is focused on main experimental aspects resulting of “helicopter simulation” capability applied on (representative) dummy electronic equipment. Wide internal instrumentation shows that the equipment is effectively answering to the stimuli by resonances excitation.
Numerical methods for the weakly compressible Generalized Langevin Model in Eulerian reference frame
NASA Astrophysics Data System (ADS)
Azarnykh, Dmitrii; Litvinov, Sergey; Adams, Nikolaus A.
2016-06-01
A well established approach for the computation of turbulent flow without resolving all turbulent flow scales is to solve a filtered or averaged set of equations, and to model non-resolved scales by closures derived from transported probability density functions (PDF) for velocity fluctuations. Effective numerical methods for PDF transport employ the equivalence between the Fokker-Planck equation for the PDF and a Generalized Langevin Model (GLM), and compute the PDF by transporting a set of sampling particles by GLM (Pope (1985) [1]). The natural representation of GLM is a system of stochastic differential equations in a Lagrangian reference frame, typically solved by particle methods. A representation in a Eulerian reference frame, however, has the potential to significantly reduce computational effort and to allow for the seamless integration into a Eulerian-frame numerical flow solver. GLM in a Eulerian frame (GLMEF) formally corresponds to the nonlinear fluctuating hydrodynamic equations derived by Nakamura and Yoshimori (2009) [12]. Unlike the more common Landau-Lifshitz Navier-Stokes (LLNS) equations these equations are derived from the underdamped Langevin equation and are not based on a local equilibrium assumption. Similarly to LLNS equations the numerical solution of GLMEF requires special considerations. In this paper we investigate different numerical approaches to solving GLMEF with respect to the correct representation of stochastic properties of the solution. We find that a discretely conservative staggered finite-difference scheme, adapted from a scheme originally proposed for turbulent incompressible flow, in conjunction with a strongly stable (for non-stochastic PDE) Runge-Kutta method performs better for GLMEF than schemes adopted from those proposed previously for the LLNS. We show that equilibrium stochastic fluctuations are correctly reproduced.
The generalized Dirichlet-Neumann map for linear elliptic PDEs and its numerical implementation
NASA Astrophysics Data System (ADS)
Sifalakis, A. G.; Fokas, A. S.; Fulton, S. R.; Saridakis, Y. G.
2008-09-01
A new approach for analyzing boundary value problems for linear and for integrable nonlinear PDEs was introduced in Fokas [A unified transform method for solving linear and certain nonlinear PDEs, Proc. Roy. Soc. London Ser. A 53 (1997) 1411-1443]. For linear elliptic PDEs, an important aspect of this approach is the characterization of a generalized Dirichlet to Neumann map: given the derivative of the solution along a direction of an arbitrary angle to the boundary, the derivative of the solution perpendicularly to this direction is computed without solving on the interior of the domain. This is based on the analysis of the so-called global relation, an equation which couples known and unknown components of the derivative on the boundary and which is valid for all values of a complex parameter k. A collocation-type numerical method for solving the global relation for the Laplace equation in an arbitrary bounded convex polygon was introduced in Fulton et al. [An analytical method for linear elliptic PDEs and its numerical implementation, J. Comput. Appl. Math. 167 (2004) 465-483]. Here, by choosing a different set of the "collocation points" (values for k), we present a significant improvement of the results in Fulton et al. [An analytical method for linear elliptic PDEs and its numerical implementation, J. Comput. Appl. Math. 167 (2004) 465-483]. The new collocation points lead to well-conditioned collocation methods. Their combination with sine basis functions leads to a collocation matrix whose diagonal blocks are point diagonal matrices yielding efficient implementation of iterative methods; numerical experimentation suggests quadratic convergence. The choice of Chebyshev basis functions leads to higher order convergence, which for regular polygons appear to be exponential.
NASA Astrophysics Data System (ADS)
Shintaku, Hirofumi; Yonemura, Tsubasa; Tsuru, Kazuaki; Isoyama, Takashi; Yambe, Tomoyuki; Kawano, Satoyuki
In this study, we construct an experimental apparatus for a prototype artificial heart and lung (AHL) by installing hollow fibers into the cylindrical tube of the vibrating flow pump (VFP). The oxygenation characteristics are investigated both by experiments using bovine blood and by numerical analyses based on the computational fluid dynamics. The analyses are carried out at the Reynolds numbers Re ranged from O(1) to O(103), which are determined based on the experimental conditions. The blood flow and the diffusion of oxygen gas are analyzed based on the Newtonian/non-Newtonian, unsteady, incompressible and axisymmetric Navier-Stokes equations, and the advection-diffusion equation. The results show that the oxygenation rate increases in proportion to Re1/3, where the phenomenon corresponds to the decreasing thickness of the concentration boundary layer with Re. Although the effects of the vibrating flow and the rheology of the blood are clearly appeared on the velocity field, their effects on the gas exchange are relatively small at the ranges of prescribed Reynolds numbers. Furthermore, the numerical results in terms of the oxygenation rate are compared with the experimental ones. The basic design data of VFP were accumulated for the development of AHL in the clinical applications.
NASA Astrophysics Data System (ADS)
Monsalve-Cano, J. F.; Darío Aristizábal-Ochoa, J.
2009-12-01
The stability and free vibration analyses (i.e., buckling, natural frequencies and modal shapes) of an orthotropic singly symmetric 3D Timoshenko beam-column with generalized boundary conditions (i.e., with bending and torsional semirigid restraints and lateral bracings as well as lumped masses at both ends) subjected to an eccentric end axial load are presented in a classical manner. The five governing equations of dynamic equilibrium (i.e., two transverse shear equations, two bending moment equations and pure torsional moment equation) are sufficient to determine the natural frequencies and the corresponding modal shapes of the beam-column in the two principal planes of bending and torsion about its longitudinal axis. The proposed model includes the coupling effects among: (1) the deformations due to bending, shear and pure torsion; (2) inertias (translational, rotational and torsional) of all masses considered; (3) eccentric axial loads applied at the ends, and (4) restraints at the supports (bending, torsional and lateral bracings at both ends of the member). However, the effects of axial deformations and warping torsion produced by the axial load are not included; consequently the proposed model is not capable of capturing the phenomena of torsional buckling or combined lateral bending-torsional buckling. The proposed analytical model indicates that the stability and dynamic response of beam-columns are highly sensitive to the coupling effects, particularly in members with both ends free to rotate. The natural frequencies and modal shapes can be determined from the eigenvalues of a full 4×4 matrix for vibration in the plane of symmetry (using the uncoupled equations of transverse force and moment equilibrium at both ends) and from a full 6×6 matrix for the coupled shear-bending-torsional vibration (using the coupled equations of transverse shear, bending and torsional moment equilibrium at both ends). Also, it is shown that the proposed method reproduces the
NASA Astrophysics Data System (ADS)
Fleury, Gérard; Mistrot, Pierre
2006-12-01
While driving off-road vehicles, operators are exposed to whole-body vibration acting in the fore-and-aft direction. Seat manufacturers supply products equipped with fore-and-aft suspension but only a few studies report on their performance. This work proposes a computational approach to design fore-and-aft suspensions for wheel loader seats. Field tests were conducted in a quarry to analyse the nature of vibration to which the driver was exposed. Typical input signals were recorded to be reproduced in the laboratory. Technical specifications are defined for the suspension. In order to evaluate the suspension vibration attenuation performance, a model of a sitting human body was developed and coupled to a seat model. The seat model combines the models of each suspension component. A linear two-degree-of-freedom model is used to describe the dynamic behaviour of the sitting driver. Model parameters are identified by fitting the computed apparent mass frequency response functions to the measured values. Model extensions are proposed to investigate postural effects involving variations in hands and feet positions and interaction of the driver's back with the backrest. Suspension design parameters are firstly optimized by computing the seat/man model response to sinusoidal acceleration. Four criteria including transmissibility, interaction force between the driver's back and the backrest and relative maximal displacement of the suspension are computed. A new suspension design with optimized features is proposed. Its performance is checked from calculations of the response of the seat/man model subjected to acceleration measured on the wheel loader during real work conditions. On the basis of the computed values of the SEAT factors, it is found possible to design a suspension that would increase the attenuation provided by the seat by a factor of two.
Fixtureless geometric inspection of nonrigid parts using "generalized numerical inspection fixture"
NASA Astrophysics Data System (ADS)
Radvar Esfahlan, Hassan
Free-form nonrigid parts form the substance of today's automotive and aerospace industries. These parts have different shapes in free state due to their dimensional and geometric variations, gravity and residual strains. For the geometric inspection of such compliant parts, special inspection fixtures, in combination with coordinate measuring systems (CMM) and/or optical data acquisition devices (scanners) are used. This inevitably causes additional costs and delays that result in a lack of competitiveness in the industry. The goal of this thesis is to facilitate the dimensional and geometrical inspection of flexible components from a point cloud without using a jig or secondary conformation operation. More specifically, we aim to develop a methodology to localize and quantify the profile defects in the case of thin shells which are typical to the aerospace and automotive industries. The presented methodology is based on the fact that the interpoint geodesic distance between any two points of a shape remains unchangeable during an isometric deformation. This study elaborates on the theory and general methods for the metrology of nonrigid parts. We have developed a Generalized Numerical Inspection Fixture (GNIF), a robust methodology which merges existing technologies in metric and computational geometry, nonlinear dimensionality reduction techniques, and finite element methods to introduce a general approach to the fixtureless geometrical inspection of nonrigid parts.
Shurtz, Timothy E.; Thomson, Scott L.
2012-01-01
Computational vocal fold models are often used to study the physics of voice production. In this paper the sensitivity of predicted vocal fold flow-induced vibration and resulting airflow patterns to several modeling selections is explored. The location of contact lines used to prevent mesh collapse and assumptions of symmetry were found to influence airflow patterns. However, these variables had relatively little effect on the vibratory response of the vocal fold model itself. Model motion was very sensitive to Poisson’s ratio. The importance of these parameter sensitivities in the context of vocal fold modeling is discussed. PMID:23794762
Shurtz, Timothy E; Thomson, Scott L
2013-06-01
Computational vocal fold models are often used to study the physics of voice production. In this paper the sensitivity of predicted vocal fold flow-induced vibration and resulting airflow patterns to several modeling selections is explored. The location of contact lines used to prevent mesh collapse and assumptions of symmetry were found to influence airflow patterns. However, these variables had relatively little effect on the vibratory response of the vocal fold model itself. Model motion was very sensitive to Poisson's ratio. The importance of these parameter sensitivities in the context of vocal fold modeling is discussed. PMID:23794762
Analytical and Numerical Solutions of Generalized Fokker-Planck Equations - Final Report
Prinja, Anil K.
2000-12-31
The overall goal of this project was to develop advanced theoretical and numerical techniques to quantitatively describe the spreading of a collimated beam of charged particles in space, in angle, and in energy, as a result of small deflection, small energy transfer Coulomb collisions with the target nuclei and electrons. Such beams arise in several applications of great interest in nuclear engineering, and include electron and ion radiotherapy, ion beam modification of materials, accelerator transmutation of waste, and accelerator production of tritium, to name some important candidates. These applications present unique and difficult modeling challenges, but from the outset are amenable to the language of ''transport theory'', which is very familiar to nuclear engineers and considerably less-so to physicists and material scientists. Thus, our approach has been to adopt a fundamental description based on transport equations, but the forward peakedness associated with charged particle interactions precludes a direct application of solution methods developed for neutral particle transport. Unique problem formulations and solution techniques are necessary to describe the transport and interaction of charged particles. In particular, we have developed the Generalized Fokker-Planck (GFP) approach to describe the angular and radial spreading of a collimated beam and a renormalized transport model to describe the energy-loss straggling of an initially monoenergetic distribution. Both analytic and numerical solutions have been investigated and in particular novel finite element numerical methods have been developed. In the first phase of the project, asymptotic methods were used to develop closed form solutions to the GFP equation for different orders of expansion, and was described in a previous progress report. In this final report we present a detailed description of (i) a novel energy straggling model based on a Fokker-Planck approximation but which is adapted for a
The sensitivity of the general circulation to Arctic Sea ice boundaries - A numerical experiment
NASA Technical Reports Server (NTRS)
Herman, G. F.; Johnson, W. T.
1978-01-01
Results are presented for a set of numerical experiments conducted with the Goddard (formerly GISS) general circulation model. The experiments were designed to test the model atmospheric response to a single fixed and specified parameter, the total ice cover in the Davis Strait, Barents Sea, East Greenland Sea, Sea of Okhotsk and Bering Sea. Margin variations are considered that are substantially smaller than those involved in ice age or ice-free Arctic simulations. Anomaly is defined as the mean of two runs corresponding to climatological maximum sea ice conditions. Model results indicate that the ice margin anomalies are capable of altering local climates in certain regions of high and middle latitudes. Possible interactions between high latitudes and subtropical regions are suggested.
Numerical model of Zeeman splitting of ro-vibrational lines in the fundamental band of NO molecule
NASA Astrophysics Data System (ADS)
Borkov, Yu. G.; Sulakshina, O. N.; Klimachev, Yu. M.
2016-07-01
This paper presents the results of calculation the LMR spectrograms of NO molecule in a variable magnetic field with maximum induction up to 6 T for probed CO laser lines. For the simulation of the LMR spectrum a numerical model was developed. This model is based on the numerical diagonalization the matrix of the effective molecular Hamiltonian, which includes Zeeman operator corresponding to interaction an external magnetic field with NO molecule. The comparison of calculated and experimental spectrograms has shown that the numerical model is very reliable and can reproduce the location of absorption peaks measured in a damped oscillating magnetic field.
NASA Astrophysics Data System (ADS)
Sourie, Aurélien; Oertel, Micaela; Novak, Jérôme
2016-04-01
We present a numerical model for uniformly rotating superfluid neutron stars in a fully general relativistic framework with, for the first time, realistic microphysics including entrainment. We compute stationary and axisymmetric configurations of neutron stars composed of two fluids, namely superfluid neutrons and charged particles (protons and electrons), rotating with different rates around a common axis. Both fluids are coupled by entrainment, a nondissipative interaction which in the case of a nonvanishing relative velocity between the fluids causes the fluid momenta to be not aligned with the respective fluid velocities. We extend the formalism put forth by Comer and Joynt in order to calculate the equation of state (EOS) and entrainment parameters for an arbitrary relative velocity as far as superfluidity is maintained. The resulting entrainment matrix fulfills all necessary sum rules, and in the limit of small relative velocity our results agree with Fermi liquid theory ones derived to lowest order in the velocity. This formalism is applied to two new nuclear equations of state which are implemented in the numerical model, which enables us to obtain precise equilibrium configurations. The resulting density profiles and moments of inertia are discussed employing both EOSs, showing the impact of entrainment and the dependence on the EOS.
NASA Astrophysics Data System (ADS)
Capron, Alexandre; Valentin, Johann; Jongmans, Denis; Baillet, Laurent; Larose, Eric; Bottelin, Pierre; Donze, Frédéric; Mangeney, Anne
2015-04-01
During the last two decades, seismic noise measurements have been increasingly used in gravitational hazard assessment for both investigation and monitoring purposes. The wide frequency range allows ambient vibrations to be applied for investigating geological and civil engineering structures in a great variety of sizes, from the lithospheric or crust scale to a few m-thick landslide and rock column or buildings. On unstable slopes, ambient vibrations have been applied in very different ways for reconnaissance, depending on the investigation purpose and the landslide type. The simplest way to extract information from ambient vibrations on a given site is to perform single-station measurements with a 3-C sensor and to process the records computing Fourier spectra of the three components or the spectral ratio between the horizontal and vertical components (the so-called H/V method). On landslide sites, several studies revealed significant spectral amplification at given frequency and polarization of the wave-field in the direction of maximum slope displacement. They show that different characteristics of the seismic noise (resonant frequencies, polarization, and spectral amplification) could be used from the spectral analysis of the motion or of spectral ratios for characterizing the landslides. For cliff-like sites, this study aims to identify the pertinent and applicable parameters that could be extracted from ambient vibrations and used to gain information on the prone-to-fall column geometry. We first use 2D numerical modeling for better understanding the influence of the rear fracture characteristics (wideness w and depth L) on the horizontal motion H(f), as well as on the spectral ratios H(f)/V(f) and H(f)/Hr(f), where Hr(f) is the horizontal motion measured at a reference site. We then identify the seismic parameters able to characterize the column decoupling and we compare numerical results to data acquired at two rocky sites exhibiting cliff-like geometry
NASA Astrophysics Data System (ADS)
El Mouhayar, Rabih; Jurdak, Murad
2016-02-01
This paper explored variation of student numerical and figural reasoning approaches across different pattern generalization types and across grade level. An instrument was designed for this purpose. The instrument was given to a sample of 1232 students from grades 4 to 11 from five schools in Lebanon. Analysis of data showed that the numerical reasoning approach seems to be more dominant than the figural reasoning approach for the near and far pattern generalization types but not for the immediate generalization type. The findings showed that for the recursive strategy, the numerical reasoning approach seems to be more dominant than the figural reasoning approach for each of the three pattern generalization types. However, the figural reasoning approach seems to be more dominant than the numerical reasoning approach for the functional strategy, for each generalization type. The findings also showed that the numerical reasoning was more dominant than the figural reasoning in lower grade levels (grades 4 and 5) for each generalization type. In contrast, the figural reasoning became more dominant than the numerical reasoning in the upper grade levels (grades 10 and 11).
Generalized ``thick'' strip modelling for vortex-induced vibration of long flexible cylinders
NASA Astrophysics Data System (ADS)
Bao, Yan; Palocios, Rafael; Sherwin, Spencer; Nektar++ Collaboration
2015-11-01
We propose a generalized strip modelling method that is computationally efficient for the VIV prediction of long flexible cylinders in three-dimensional incompressible flow. In order to overcome the shortcomings of conventional strip theory-based 2D models, the fluid domain is divided into ``thick'' strips, which are sufficiently thick to locally resolve the small scale turbulence effects and three dimensionality of the flow around the cylinder. An attractive feature of the model is that we independently construct a three-dimensional scale resolving model for individual strips, which have local spanwise scale along the cylinder's axial direction and are only coupled through the structural model of the cylinder. Therefore, this model is able to cover the situations of fully resolved 3D model and 2D strip theory model. The connection between these strips is achieved through the calculation of a tensioned beam equation, which is used to represent the dynamics of the flexible body. In the limit, however, a single ``thick'' strip would request the full 3D domain. A parallel Fourier spectral/hp element method is employed to solve the 3D flow dynamics in the strip-domain, and then the VIV response prediction is achieved through the strip-structure interactions. This work is supported by EPSRC grant EP/K037536/1. Acknowledge UK Turbulence Consortium (UKTC) for ARCHER time under EPSRC grant EP/L000261/1.
NASA Astrophysics Data System (ADS)
Li, Zhaorui; Livescu, Daniel
2014-11-01
By using the second-law of thermodynamics and the Onsager reciprocal method for irreversible processes, we have developed a set of physically consistent multicomponent compressible generalized Cahn-Hilliard Navier-Stokes (CGCHNS) equations from basic thermodynamics. The new equations can describe not only flows with pure miscible and pure immiscible materials but also complex flows in which mass diffusion and surface tension or Korteweg stresses effects may coexist. Furthermore, for the first time, the incompressible generalized Cahn-Hilliard Navier-Stokes (IGCHNS) equations are rigorously derived from the incompressible limit of the CGCHNS equations (as the infinite sound speed limit) and applied to the immiscible Rayleigh-Taylor instability problem. Extensive good agreements between numerical results and the linear stability theory (LST) predictions for the Rayleigh-Taylor instability are achieved for a wide range of wavenumber, surface tension, and viscosity values. The late-time results indicate that the IGCHNS equations can naturally capture complex interface topological changes including merging and breaking-up and are free of singularity problems.
NASA Technical Reports Server (NTRS)
Wie, Yong-Sun
1990-01-01
A procedure for calculating 3-D, compressible laminar boundary layer flow on general fuselage shapes is described. The boundary layer solutions can be obtained in either nonorthogonal 'body oriented' coordinates or orthogonal streamline coordinates. The numerical procedure is 'second order' accurate, efficient and independent of the cross flow velocity direction. Numerical results are presented for several test cases, including a sharp cone, an ellipsoid of revolution, and a general aircraft fuselage at angle of attack. Comparisons are made between numerical results obtained using nonorthogonal curvilinear 'body oriented' coordinates and streamline coordinates.
NASA Astrophysics Data System (ADS)
Zhang, Han; Zhang, Zhao-Hui; Zhao, Xiao-Yan; Zhang, Tian-Yao; Yan, Fang; Shen, Jiang
2015-07-01
The crystal structure of L-glutamine is stabilized by a three-dimensional network of intermolecular hydrogen bonds. We utilize plane-wave density functional theory lattice-dynamics calculations within the generalized-gradient approximation (GGA), Perdew-Burke-Ernzerhof (PBE), PBE for solids (PBEsol), PBE with Wu-Cohen exchange (WC), and dispersion-corrected PBE, to investigate the effect of these intermolecular contacts on the absorption spectra of glutamine in the terahertz frequency range. Among these calculations, the solid-state simulated results obtained using the WC method exhibit a good agreement with the measured absorption spectra, and the absorption features are assigned with the help of WC. This indicates that the vibrational modes of glutamine were related to the combination of intramolecular and intermolecular motions, the intramolecular modes were dominated by rocking or torsion involving functional groups; the intermolecular modes mainly result from the translational motions of individual molecules, and the rocking of the hydrogen-bonded functional groups. Project supported by the National Natural Science Foundation of China (Grant Nos. 61302007 and 60977065), the Fundamental Research Funds for the Central Universities of China (Grant No. FRF-SD-12-016A), and the Engineering Research Center of Industrial Spectrum Imaging of Beijing, China.
NASA Astrophysics Data System (ADS)
Degrande, G.; Clouteau, D.; Othman, R.; Arnst, M.; Chebli, H.; Klein, R.; Chatterjee, P.; Janssens, B.
2006-06-01
A numerical model is presented to predict vibrations in the free field from excitation due to metro trains in tunnels. The three-dimensional dynamic tunnel-soil interaction problem is solved with a subdomain formulation, using a finite element formulation for the tunnel and a boundary element method for the soil. The periodicity of the geometry in the longitudinal direction of the tunnel is exploited using the Floquet transform, limiting the discretization to a single-bounded reference cell. The responses of two different types of tunnel due to a harmonic load on the tunnel invert are compared, both in the frequency-wavenumber and spatial domains. The first tunnel is a shallow cut-and-cover masonry tunnel on the Paris metro network, embedded in layers of sand, while the second tunnel is a deep bored tunnel of London Underground, with a cast iron lining and embedded in the London clay.
Aref's chaotic orbits tracked by a general ellipsoid using 3D numerical simulations
NASA Astrophysics Data System (ADS)
Shui, Pei; Popinet, Stéphane; Govindarajan, Rama; Valluri, Prashant
2015-11-01
The motion of an ellipsoidal solid in an ideal fluid has been shown to be chaotic (Aref, 1993) under the limit of non-integrability of Kirchhoff's equations (Kozlov & Oniscenko, 1982). On the other hand, the particle could stop moving when the damping viscous force is strong enough. We present numerical evidence using our in-house immersed solid solver for 3D chaotic motion of a general ellipsoidal solid and suggest criteria for triggering such motion. Our immersed solid solver functions under the framework of the Gerris flow package of Popinet et al. (2003). This solver, the Gerris Immersed Solid Solver (GISS), resolves 6 degree-of-freedom motion of immersed solids with arbitrary geometry and number. We validate our results against the solution of Kirchhoff's equations. The study also shows that the translational/ rotational energy ratio plays the key role on the motion pattern, while the particle geometry and density ratio between the solid and fluid also have some influence on the chaotic behaviour. Along with several other benchmark cases for viscous flows, we propose prediction of chaotic Aref's orbits as a key benchmark test case for immersed boundary/solid solvers.
Numerical study of fourth-order linearized compact schemes for generalized NLS equations
NASA Astrophysics Data System (ADS)
Liao, Hong-lin; Shi, Han-sheng; Zhao, Ying
2014-08-01
The fourth-order compact approximation for the spatial second-derivative and several linearized approaches, including the time-lagging method of Zhang et al. (1995), the local-extrapolation technique of Chang et al. (1999) and the recent scheme of Dahlby et al. (2009), are considered in constructing fourth-order linearized compact difference (FLCD) schemes for generalized NLS equations. By applying a new time-lagging linearized approach, we propose a symmetric fourth-order linearized compact difference (SFLCD) scheme, which is shown to be more robust in long-time simulations of plane wave, breather, periodic traveling-wave and solitary wave solutions. Numerical experiments suggest that the SFLCD scheme is a little more accurate than some other FLCD schemes and the split-step compact difference scheme of Dehghan and Taleei (2010). Compared with the time-splitting pseudospectral method of Bao et al. (2003), our SFLCD method is more suitable for oscillating solutions or the problems with a rapidly varying potential.
NASA Technical Reports Server (NTRS)
Padovan, J.; Adams, M.; Fertis, J.; Zeid, I.; Lam, P.
1982-01-01
Finite element codes are used in modelling rotor-bearing-stator structure common to the turbine industry. Engine dynamic simulation is used by developing strategies which enable the use of available finite element codes. benchmarking the elements developed are benchmarked by incorporation into a general purpose code (ADINA); the numerical characteristics of finite element type rotor-bearing-stator simulations are evaluated through the use of various types of explicit/implicit numerical integration operators. Improving the overall numerical efficiency of the procedure is improved.
A general spectral method for the numerical simulation of one-dimensional interacting fermions
NASA Astrophysics Data System (ADS)
Clason, Christian; von Winckel, Gregory
2012-08-01
This software implements a general framework for the direct numerical simulation of systems of interacting fermions in one spatial dimension. The approach is based on a specially adapted nodal spectral Galerkin method, where the basis functions are constructed to obey the antisymmetry relations of fermionic wave functions. An efficient Matlab program for the assembly of the stiffness and potential matrices is presented, which exploits the combinatorial structure of the sparsity pattern arising from this discretization to achieve optimal run-time complexity. This program allows the accurate discretization of systems with multiple fermions subject to arbitrary potentials, e.g., for verifying the accuracy of multi-particle approximations such as Hartree-Fock in the few-particle limit. It can be used for eigenvalue computations or numerical solutions of the time-dependent Schrödinger equation. The new version includes a Python implementation of the presented approach. New version program summaryProgram title: assembleFermiMatrix Catalogue identifier: AEKO_v1_1 Program summary URL:http://cpc.cs.qub.ac.uk/summaries/AEKO_v1_1.html Program obtainable from: CPC Program Library, Queen's University, Belfast, N. Ireland Licensing provisions: Standard CPC licence, http://cpc.cs.qub.ac.uk/licence/licence.html No. of lines in distributed program, including test data, etc.: 332 No. of bytes in distributed program, including test data, etc.: 5418 Distribution format: tar.gz Programming language: MATLAB/GNU Octave, Python Computer: Any architecture supported by MATLAB, GNU Octave or Python Operating system: Any supported by MATLAB, GNU Octave or Python RAM: Depends on the data Classification: 4.3, 2.2. External routines: Python 2.7+, NumPy 1.3+, SciPy 0.10+ Catalogue identifier of previous version: AEKO_v1_0 Journal reference of previous version: Comput. Phys. Commun. 183 (2012) 405 Does the new version supersede the previous version?: Yes Nature of problem: The direct numerical
NASA Astrophysics Data System (ADS)
Blake, W. K.
Physical and mathematical analyses of the vibration and sound induced by different types of fluid flow are presented in a comprehensive introduction intended primarily for practicing engineers. The elementary concepts are explained, and chapters are devoted to the theory of sound and its generation by flow; shear-layer instabilities, flow tones, and jet noise; dipole sound from cylinders; the fundamentals of flow-induced vibration and noise; bubble dynamics and cavitation; hydrodynamically induced cavitation and bubble noise; turbulent wall-pressure fluctuations; structural response to turbulent wall flow and random sound; noise radiation from pipe and duct systems; noncavitating lifting sections; and noise from rotating machinery. Graphs, diagrams, drawings, and tables of numerical data are provided.
Iwasa, Takeshi; Takenaka, Masato; Taketsugu, Tetsuya
2016-03-28
A theoretical method to compute infrared absorption spectra when a molecule is interacting with an arbitrary nonuniform electric field such as near-fields is developed and numerically applied to simple model systems. The method is based on the multipolar Hamiltonian where the light-matter interaction is described by a spatial integral of the inner product of the molecular polarization and applied electric field. The computation scheme is developed under the harmonic approximation for the molecular vibrations and the framework of modern electronic structure calculations such as the density functional theory. Infrared reflection absorption and near-field infrared absorption are considered as model systems. The obtained IR spectra successfully reflect the spatial structure of the applied electric field and corresponding vibrational modes, demonstrating applicability of the present method to analyze modern nanovibrational spectroscopy using near-fields. The present method can use arbitral electric fields and thus can integrate two fields such as computational chemistry and electromagnetics. PMID:27036436
TOUGH2: A general-purpose numerical simulator for multiphase nonisothermal flows
Pruess, K.
1991-06-01
Numerical simulators for multiphase fluid and heat flows in permeable media have been under development at Lawrence Berkeley Laboratory for more than 10 yr. Real geofluids contain noncondensible gases and dissolved solids in addition to water, and the desire to model such `compositional` systems led to the development of a flexible multicomponent, multiphase simulation architecture known as MULKOM. The design of MULKOM was based on the recognition that the mass-and energy-balance equations for multiphase fluid and heat flows in multicomponent systems have the same mathematical form, regardless of the number and nature of fluid components and phases present. Application of MULKOM to different fluid mixtures, such as water and air, or water, oil, and gas, is possible by means of appropriate `equation-of-state` (EOS) modules, which provide all thermophysical and transport parameters of the fluid mixture and the permeable medium as a function of a suitable set of primary thermodynamic variables. Investigations of thermal and hydrologic effects from emplacement of heat-generating nuclear wastes into partially water-saturated formations prompted the development and release of a specialized version of MULKOM for nonisothermal flow of water and air, named TOUGH. TOUGH is an acronym for `transport of unsaturated groundwater and heat` and is also an allusion to the tuff formations at Yucca Mountain, Nevada. The TOUGH2 code is intended to supersede TOUGH. It offers all the capabilities of TOUGH and includes a considerably more general subset of MULKOM modules with added capabilities. The paper briefly describes the simulation methodology and user features.
ERIC Educational Resources Information Center
Cowan, Richard; Powell, Daisy
2014-01-01
Explanations of the marked individual differences in elementary school mathematical achievement and mathematical learning disability (MLD or dyscalculia) have involved domain-general factors (working memory, reasoning, processing speed, and oral language) and numerical factors that include single-digit processing efficiency and multidigit skills…
ERIC Educational Resources Information Center
El Mouhayar, Rabih; Jurdak, Murad
2016-01-01
This paper explored variation of student numerical and figural reasoning approaches across different pattern generalization types and across grade level. An instrument was designed for this purpose. The instrument was given to a sample of 1232 students from grades 4 to 11 from five schools in Lebanon. Analysis of data showed that the numerical…
Cowan, Richard; Powell, Daisy
2014-02-01
Explanations of the marked individual differences in elementary school mathematical achievement and mathematical learning disability (MLD or dyscalculia) have involved domain-general factors (working memory, reasoning, processing speed, and oral language) and numerical factors that include single-digit processing efficiency and multidigit skills such as number system knowledge and estimation. This study of 3rd graders (N = 258) finds both domain-general and numerical factors contribute independently to explaining variation in 3 significant arithmetic skills: basic calculation fluency, written multidigit computation, and arithmetic word problems. Estimation accuracy and number system knowledge show the strongest associations with every skill, and their contributions are independent of both each other and other factors. Different domain-general factors independently account for variation in each skill. Numeral comparison, a single digit processing skill, uniquely accounts for variation in basic calculation. Subsamples of children with MLD (at or below 10th percentile, n = 29) are compared with low achievement (LA, 11th to 25th percentiles, n = 42) and typical achievement (above 25th percentile, n = 187). Examination of these and subsets with persistent difficulties supports a multiple deficits view of number difficulties: Most children with number difficulties exhibit deficits in both domain-general and numerical factors. The only factor deficit common to all persistent MLD children is in multidigit skills. These findings indicate that many factors matter but multidigit skills matter most in 3rd grade mathematical achievement. PMID:24532854
Numeric and symbolic evaluation of the pfaffian of general skew-symmetric matrices
NASA Astrophysics Data System (ADS)
González-Ballestero, C.; Robledo, L. M.; Bertsch, G. F.
2011-10-01
Evaluation of pfaffians arises in a number of physics applications, and for some of them a direct method is preferable to using the determinantal formula. We discuss two methods for the numerical evaluation of pfaffians. The first is tridiagonalization based on Householder transformations. The main advantage of this method is its numerical stability that makes unnecessary the implementation of a pivoting strategy. The second method considered is based on Aitken's block diagonalization formula. It yields to a kind of LU (similar to Cholesky's factorization) decomposition (under congruence) of arbitrary skew-symmetric matrices that is well suited both for the numeric and symbolic evaluations of the pfaffian. Fortran subroutines (FORTRAN 77 and 90) implementing both methods are given. We also provide simple implementations in Python and Mathematica for purpose of testing, or for exploratory studies of methods that make use of pfaffians.
NASA Astrophysics Data System (ADS)
Sakaris, C. S.; Sakellariou, J. S.; Fassois, S. D.
2016-06-01
This study focuses on the problem of vibration-based damage precise localization via data-based, time series type, methods for structures consisting of 1D, 2D, or 3D elements. A Generalized Functional Model Based method is postulated based on an expanded Vector-dependent Functionally Pooled ARX (VFP-ARX) model form, capable of accounting for an arbitrary structural topology. The FP model's operating parameter vector elements are properly constrained to reflect any given topology. Damage localization is based on operating parameter vector estimation within the specified topology, so that the location estimate and its uncertainty bounds are statistically optimal. The method's effectiveness is experimentally demonstrated through damage precise localization on a laboratory spatial truss structure using various damage scenarios and a single pair of random excitation - vibration response signals in a low and limited frequency bandwidth.
Numerical Convergence of the Block-Maxima Approach to the Generalized Extreme Value Distribution
NASA Astrophysics Data System (ADS)
Faranda, Davide; Lucarini, Valerio; Turchetti, Giorgio; Vaienti, Sandro
2011-12-01
In this paper we perform an analytical and numerical study of Extreme Value distributions in discrete dynamical systems. In this setting, recent works have shown how to get a statistics of extremes in agreement with the classical Extreme Value Theory. We pursue these investigations by giving analytical expressions of Extreme Value distribution parameters for maps that have an absolutely continuous invariant measure. We compare these analytical results with numerical experiments in which we study the convergence to limiting distributions using the so called block-maxima approach, pointing out in which cases we obtain robust estimation of parameters. In regular maps for which mixing properties do not hold, we show that the fitting procedure to the classical Extreme Value Distribution fails, as expected. However, we obtain an empirical distribution that can be explained starting from a different observable function for which Nicolis et al. (Phys. Rev. Lett. 97(21): 210602, 2006) have found analytical results.
Analytical and experimental study of vibrations in a gear transmission
NASA Technical Reports Server (NTRS)
Choy, F. K.; Ruan, Y. F.; Zakrajsek, J. J.; Oswald, F. B.; Coy, J. J.
1991-01-01
An analytical simulation of the dynamics of a gear transmission system is presented and compared to experimental results from a gear noise test rig at NASA Lewis. The analytical procedure developed couples the dynamic behaviors of the rotor-bearing-gear system with the response of the gearbox structure. Transient and steady-state vibrations of the gearbox system are presented in the time and frequency domains. The vibration characteristics of a simple single-mesh-gear noise test rig are modeled. The numerical simulations are compared to experimental data measured under typical operating conditions. The system natural frequencies, peak vibration amplitudes, and gear mesh frequencies are generally in good agreement.
NASA Technical Reports Server (NTRS)
Padovan, J.; Adams, M.; Lam, P.; Fertis, D.; Zeid, I.
1982-01-01
Second-year efforts within a three-year study to develop and extend finite element (FE) methodology to efficiently handle the transient/steady state response of rotor-bearing-stator structure associated with gas turbine engines are outlined. The two main areas aim at (1) implanting the squeeze film damper element into a general purpose FE code for testing and evaluation; and (2) determining the numerical characteristics of the FE-generated rotor-bearing-stator simulation scheme. The governing FE field equations are set out and the solution methodology is presented. The choice of ADINA as the general-purpose FE code is explained, and the numerical operational characteristics of the direct integration approach of FE-generated rotor-bearing-stator simulations is determined, including benchmarking, comparison of explicit vs. implicit methodologies of direct integration, and demonstration problems.
NASA Astrophysics Data System (ADS)
Zarepour, Misagh; Amirhosein Hosseini, Seyed
2016-08-01
This study presents an examination of nonlinear free vibration of a nanobeam under electro-thermo-mechanical loading with elastic medium and various boundary conditions, especially the elastic boundary condition. The nanobeam is modeled as an Euler–Bernoulli beam. The von Kármán strain-displacement relationship together with Hamilton’s principle and Eringen’s theory are employed to derive equations of motion. The nonlinear free vibration frequency is obtained for simply supported (S-S) and elastic supported (E-E) boundary conditions. E-E boundary condition is a general and actual form of boundary conditions and it is chosen because of more realistic behavior. By applying the differential transform method (DTM), the nanobeam’s natural frequencies can be easily obtained for the two different boundary conditions mentioned above. Performing a precise study led to investigation of the influences of nonlocal parameter, temperature change, spring constants (either for elastic medium or boundary condition) and imposed electric potential on the nonlinear free vibration characteristics of nanobeam. The results for S-S and E-E nanobeams are compared with each other. In order to validate the results, some comparisons are presented between DTM results and open literature to show the accuracy of this new approach. It has been discovered that DTM solves the equations with minimum calculation cost.
NASA Astrophysics Data System (ADS)
Zarepour, Misagh; Amirhosein Hosseini, Seyed
2016-08-01
This study presents an examination of nonlinear free vibration of a nanobeam under electro-thermo-mechanical loading with elastic medium and various boundary conditions, especially the elastic boundary condition. The nanobeam is modeled as an Euler–Bernoulli beam. The von Kármán strain-displacement relationship together with Hamilton’s principle and Eringen’s theory are employed to derive equations of motion. The nonlinear free vibration frequency is obtained for simply supported (S-S) and elastic supported (E-E) boundary conditions. E-E boundary condition is a general and actual form of boundary conditions and it is chosen because of more realistic behavior. By applying the differential transform method (DTM), the nanobeam’s natural frequencies can be easily obtained for the two different boundary conditions mentioned above. Performing a precise study led to investigation of the influences of nonlocal parameter, temperature change, spring constants (either for elastic medium or boundary condition) and imposed electric potential on the nonlinear free vibration characteristics of nanobeam. The results for S-S and E-E nanobeams are compared with each other. In order to validate the results, some comparisons are presented between DTM results and open literature to show the accuracy of this new approach. It has been discovered that DTM solves the equations with minimum calculation cost.
Chau, W.Y.; Lake, K.; Stone, J.
1984-06-15
The Oppenheimer-Volkoff equations for isothermal, partially degenerate neutral lepton configurations have been solved numerically for arbitrary temperature and degree of degeneracy. For the specific case where the leptons are massive neutrions (approx.30 eV) with properties quite tightly constrained by cosmological considerations, we find that it possible to have a neutrino halo surrounding a normal galaxy with the right values of mass and radius required of the ''invisible halo'' in the missing mass problem. The density distribution, however, yields a rotational curve which, with the present chosen values of parameters, does not fit the observations for spiral galaxies.
A novel numerical flux for the 3D Euler equations with general equation of state
NASA Astrophysics Data System (ADS)
Toro, Eleuterio F.; Castro, Cristóbal E.; Lee, Bok Jik
2015-12-01
Here we extend the flux vector splitting approach recently proposed in E.F. Toro and M.E. Vázquez-Cendón (2012) [42]. The scheme was originally presented for the 1D Euler equations for ideal gases and its extension presented in this paper is threefold: (i) we solve the three-dimensional Euler equations on general meshes; (ii) we use a general equation of state; and (iii) we achieve high order of accuracy in both space and time through application of the semi-discrete ADER methodology on general meshes. The resulting methods are systematically assessed for accuracy, robustness and efficiency on a carefully selected suite of test problems. Formal high accuracy is assessed through convergence rates studies for schemes of up to 4th order of accuracy in both space and time on unstructured meshes.
A Numerical Treatment of the Rf SQUID: I. General Properties andNoise Energy
Kleiner, Reinhold; Koelle, Dieter; Clarke, John
2007-01-15
We investigate the characteristics and noise performance of rf Superconducting Quantum Interference Devices (SQUIDs) by solving the corresponding Langevin equations numerically and optimizing the model parameters with respect to noise energy. After introducing the basic concepts of the numerical simulations, we give a detailed discussion of the performance of the SQUID as a function of all relevant parameters. The best performance is obtained in the crossover region between the dispersive and dissipative regimes, characterized by an inductance parameter {beta}{prime}{sub L} {triple_bond} 2{pi}LI{sub 0}/{Phi}{sub 0} {approx} 1; L is the loop inductance, I{sub 0} the critical current of the Josephson junction, and {phi}{sub 0} the flux quantum. In this regime, which is not well explored by previous analytical approaches, the lowest (intrinsic) values of noise energy are a factor of about 2 above previous estimates based on analytical approaches. However, several other analytical predictions, such as the inverse proportionality of the noise energy on the tank circuit quality factor and the square of the coupling coefficient between the tank circuit and the SQUID loop, could not be well reproduced. The optimized intrinsic noise energy of the rf SQUID is superior to that of the dc SQUID at all temperatures. Although for technologically achievable parameters this advantage shrinks, particularly at low thermal fluctuation levels, we give an example for realistic parameters that leads to a noise energy comparable to that of the dc SQUID even in this regime.
NASA Technical Reports Server (NTRS)
Green, C.
1971-01-01
Guidelines of the methods and applications used in vibration technology at the MSFC are presented. The purpose of the guidelines is to provide a practical tool for coordination and understanding between industry and government groups concerned with vibration of systems and equipments. Topics covered include measuring, reducing, analyzing, and methods for obtaining simulated environments and formulating vibration specifications. Methods for vibration and shock testing, theoretical aspects of data processing, vibration response analysis, and techniques of designing for vibration are also presented.
TOUGH2: A general-purpose numerical simulator for multiphase fluid and heat flow
Pruess, K.
1991-05-01
TOUGH2 is a numerical simulation program for nonisothermal flows of multicomponent, multiphase fluids in porous and fractured media. The chief applications for which TOUGH2 is designed are in geothermal reservoir engineering, nuclear waste disposal, and unsaturated zone hydrology. A successor to the TOUGH program, TOUGH2 offers added capabilities and user features, including the flexibility to handle different fluid mixtures, facilities for processing of geometric data (computational grids), and an internal version control system to ensure referenceability of code applications. This report includes a detailed description of governing equations, program architecture, and user features. Enhancements in data inputs relative to TOUGH are described, and a number of sample problems are given to illustrate code applications. 46 refs., 29 figs., 12 tabs.
Numeric promoter description - A comparative view on concepts and general application.
Beier, Rico; Labudde, Dirk
2016-01-01
Nucleic acid molecules play a key role in a variety of biological processes. Starting from storage and transfer tasks, this also comprises the triggering of biological processes, regulatory effects and the active influence gained by target binding. Based on the experimental output (in this case promoter sequences), further in silico analyses aid in gaining new insights into these processes and interactions. The numerical description of nucleic acids thereby constitutes a bridge between the concrete biological issues and the analytical methods. Hence, this study compares 26 descriptor sets obtained by applying well-known numerical description concepts to an established dataset of 38 DNA promoter sequences. The suitability of the description sets was evaluated by computing partial least squares regression models and assessing the model accuracy. We conclude that the major importance regarding the descriptive power is attached to positional information rather than to explicitly incorporated physico-chemical information, since a sufficient amount of implicit physico-chemical information is already encoded in the nucleobase classification. The regression models especially benefited from employing the information that is encoded in the sequential and structural neighborhood of the nucleobases. Thus, the analyses of n-grams (short fragments of length n) suggested that they are valuable descriptors for DNA target interactions. A mixed n-gram descriptor set thereby yielded the best description of the promoter sequences. The corresponding regression model was checked and found to be plausible as it was able to reproduce the characteristic binding motifs of promoter sequences in a reasonable degree. As most functional nucleic acids are based on the principle of molecular recognition, the findings are not restricted to promoter sequences, but can rather be transferred to other kinds of functional nucleic acids. Thus, the concepts presented in this study could provide
NASA Astrophysics Data System (ADS)
Brown, Sandra E.; Georgescu, IonuÅ£; Mandelshtam, Vladimir A.
2013-01-01
The self-consistent phonons (SCP) method provides a consistent way to include anharmonic effects when treating a many-body quantum system at thermal equilibrium. The system is then described by an effective temperature-dependent harmonic Hamiltonian, which can be used to estimate the system's properties, such as its free energy or its vibrational spectrum. The numerical bottleneck of the method is the evaluation of Gaussian averages of the potential energy and its derivatives. Several algorithmic ideas/tricks are introduced to reduce the cost of such integration by orders of magnitude, e.g., relative to that of the previous implementation of the SCP approach by Calvo et al. [J. Chem. Phys. 133, 074303 (2010), 10.1063/1.3465554]. One such algorithmic improvement is the replacement of standard Monte Carlo integration by quasi-Monte Carlo integration utilizing low-discrepancy sequences. The performance of the method is demonstrated on the calculation of vibrational frequencies of pyrene. It is then applied to compute the free energies of five isomers of water hexamer using the WHBB potential of Bowman and co-workers [J. Chem. Phys. 134, 094509 (2011), 10.1063/1.3554905]. The present results predict the hexamer prism being thermodynamically most stable, with the free energy of the hexamer cage being about 0.2 kcal mol-1 higher at all temperatures below T = 200 K.
Piccardo, Matteo; Bloino, Julien; Barone, Vincenzo
2015-01-01
Models going beyond the rigid-rotor and the harmonic oscillator levels are mandatory for providing accurate theoretical predictions for several spectroscopic properties. Different strategies have been devised for this purpose. Among them, the treatment by perturbation theory of the molecular Hamiltonian after its expansion in power series of products of vibrational and rotational operators, also referred to as vibrational perturbation theory (VPT), is particularly appealing for its computational efficiency to treat medium-to-large systems. Moreover, generalized (GVPT) strategies combining the use of perturbative and variational formalisms can be adopted to further improve the accuracy of the results, with the first approach used for weakly coupled terms, and the second one to handle tightly coupled ones. In this context, the GVPT formulation for asymmetric, symmetric, and linear tops is revisited and fully generalized to both minima and first-order saddle points of the molecular potential energy surface. The computational strategies and approximations that can be adopted in dealing with GVPT computations are pointed out, with a particular attention devoted to the treatment of symmetry and degeneracies. A number of tests and applications are discussed, to show the possibilities of the developments, as regards both the variety of treatable systems and eligible methods. © 2015 Wiley Periodicals, Inc. PMID:26345131
NASA Technical Reports Server (NTRS)
Rogers, S. E.; Kwak, D.; Chang, J. L. C.
1986-01-01
Numerically solving the incompressible Navier-Stokes equations is known to be time consuming and expensive. Testing of the INS3D computers code, which solves these equations with the use of the pseudocompressibility method, shows this method to be an efficient way to obtain the steady state solution. The effects of the waves introduced by the pseudocompressibility method are analyzed and criteria are set and tested for the choice of the pseudocompressibility parameter which governs the artificial sound speed. The code is tested using laminar flow over a two dimensional backward-facing step, and laminar flow over a two dimensional circular cylinder. The results of the computations over the backward-facing step are in excellent agreement with experimental results. The transient solution of the flow over the cylinder impulsively started from rest is in good agreement with experimental results. However, the computed frequency of periodic shedding of vortices behind the cylinder is not in agreement with the experimental value. For a three dimensional test case, computations were conducted for a cylinder end wall junction. The saddle point separation and horseshoe vortex system appear in the computed field. The solution also shows secondary vortex filaments which wrap around the cylinder and spiral up in the wake.
An efficient and general numerical method to compute steady uniform vortices
NASA Astrophysics Data System (ADS)
Luzzatto-Fegiz, Paolo; Williamson, Charles H. K.
2011-07-01
Steady uniform vortices are widely used to represent high Reynolds number flows, yet their efficient computation still presents some challenges. Existing Newton iteration methods become inefficient as the vortices develop fine-scale features; in addition, these methods cannot, in general, find solutions with specified Casimir invariants. On the other hand, available relaxation approaches are computationally inexpensive, but can fail to converge to a solution. In this paper, we overcome these limitations by introducing a new discretization, based on an inverse-velocity map, which radically increases the efficiency of Newton iteration methods. In addition, we introduce a procedure to prescribe Casimirs and remove the degeneracies in the steady vorticity equation, thus ensuring convergence for general vortex configurations. We illustrate our methodology by considering several unbounded flows involving one or two vortices. Our method enables the computation, for the first time, of steady vortices that do not exhibit any geometric symmetry. In addition, we discover that, as the limiting vortex state for each flow is approached, each family of solutions traces a clockwise spiral in a bifurcation plot consisting of a velocity-impulse diagram. By the recently introduced "IVI diagram" stability approach [Phys. Rev. Lett. 104 (2010) 044504], each turn of this spiral is associated with a loss of stability for the steady flows. Such spiral structure is suggested to be a universal feature of steady, uniform-vorticity flows.
Using a Vibration Device to Ease Pain During Facial Needling and Injection
Ogawa, Rei
2016-01-01
Objective: In general, needling and injection are painful procedures, especially when the face is the target. Although local anesthetics (cream or tape) can be used to reduce the pain, they are not sufficiently effective. It has been suggested that vibration can reduce pain. The aim of this case study was to determine whether application of a vibration device to an area adjacent to the facial target area to be injected/needled would relieve pain. Methods: Consecutive women scheduled to undergo facial injection with hyaluronic acid or botulinum toxin were recruited. Half of the face was injected with concomitant vibration, whereas the other half was injected without vibration. The pain experienced by the women during both procedures was assessed using the Numeric Rating Scale. The safety of injection with vibration was also assessed. Results: Of the 32 patients, 28 indicated that vibration relieved the pain, 3 stated that it had no effect, and 1 (who received deep botulinum toxin injections to the masseter muscle) complained that it made the pain worse. Vibration did not affect the safety of the injections. The average Numeric Rating Scale scores for the no-vibration and vibration injections were 4.5 ± 1.5 and 2.3 ± 0.9, respectively (P < .001). Conclusions: The Gate Control Theory of Pain explains why vibration reduces pain. PMID:26933468
NASA Technical Reports Server (NTRS)
Zeng, S.; Wesseling, P.
1993-01-01
The performance of a linear multigrid method using four smoothing methods, called SCGS (Symmetrical Coupled GauBeta-Seidel), CLGS (Collective Line GauBeta-Seidel), SILU (Scalar ILU), and CILU (Collective ILU), is investigated for the incompressible Navier-Stokes equations in general coordinates, in association with Galerkin coarse grid approximation. Robustness and efficiency are measured and compared by application to test problems. The numerical results show that CILU is the most robust, SILU the least, with CLGS and SCGS in between. CLGS is the best in efficiency, SCGS and CILU follow, and SILU is the worst.
Cryns, Jackson W.; Hatchell, Brian K.; Santiago-Rojas, Emiliano; Silvers, Kurt L.
2013-07-01
Formal journal article Experimental analysis of a piezoelectric energy harvesting system for harmonic, random, and sine on random vibration Abstract: Harvesting power with a piezoelectric vibration powered generator using a full-wave rectifier conditioning circuit is experimentally compared for varying sinusoidal, random and sine on random (SOR) input vibration scenarios. Additionally, the implications of source vibration characteristics on harvester design are discussed. Studies in vibration harvesting have yielded numerous alternatives for harvesting electrical energy from vibrations but piezoceramics arose as the most compact, energy dense means of energy transduction. The rise in popularity of harvesting energy from ambient vibrations has made piezoelectric generators commercially available. Much of the available literature focuses on maximizing harvested power through nonlinear processing circuits that require accurate knowledge of generator internal mechanical and electrical characteristics and idealization of the input vibration source, which cannot be assumed in general application. In this manuscript, variations in source vibration and load resistance are explored for a commercially available piezoelectric generator. We characterize the source vibration by its acceleration response for repeatability and transcription to general application. The results agree with numerical and theoretical predictions for in previous literature that load optimal resistance varies with transducer natural frequency and source type, and the findings demonstrate that significant gains are seen with lower tuned transducer natural frequencies for similar source amplitudes. Going beyond idealized steady state sinusoidal and simplified random vibration input, SOR testing allows for more accurate representation of real world ambient vibration. It is shown that characteristic interactions from more complex vibrational sources significantly alter power generation and power processing
Numerical simulation of the general circulation of the atmosphere of Titan.
Hourdin, F; Talagrand, O; Sadourny, R; Courtin, R; Gautier, D; McKay, C P
1995-10-01
The atmospheric circulation of Titan is investigated with a general circulation model. The representation of the large-scale dynamics is based on a grid point model developed and used at Laboratoire de Météorologie Dynamique for climate studies. The code also includes an accurate representation of radiative heating and cooling by molecular gases and haze as well as a parametrization of the vertical turbulent mixing of momentum and potential temperature. Long-term simulations of the atmospheric circulation are presented. Starting from a state of rest, the model spontaneously produces a strong superrotation with prograde equatorial winds (i.e., in the same sense as the assumed rotation of the solid body) increasing from the surface to reach 100 m sec-1 near the 1-mbar pressure level. Those equatorial winds are in very good agreement with some indirect observations, especially those of the 1989 occultation of Star 28-Sgr by Titan. On the other hand, the model simulates latitudinal temperature contrasts in the stratosphere that are significantly weaker than those observed by Voyager 1 which, we suggest, may be partly due to the nonrepresentation of the spatial and temporal variations of the abundances of molecular species and haze. We present diagnostics of the simulated atmospheric circulation underlying the importance of the seasonal cycle and a tentative explanation for the creation and maintenance of the atmospheric superrotation based on a careful angular momentum budget. PMID:11538593
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
NASA Astrophysics Data System (ADS)
Klein, C.; Peter, R.
2015-06-01
We present a detailed numerical study of solutions to general Korteweg-de Vries equations with critical and supercritical nonlinearity, both in the context of dispersive shocks and blow-up. We study the stability of solitons and show that they are unstable against being radiated away and blow-up. In the L2 critical case, the blow-up mechanism by Martel, Merle and Raphaël can be numerically identified. In the limit of small dispersion, it is shown that a dispersive shock always appears before an eventual blow-up. In the latter case, always the first soliton to appear will blow up. It is shown that the same type of blow-up as for the perturbations of the soliton can be observed which indicates that the theory by Martel, Merle and Raphaël is also applicable to initial data with a mass much larger than the soliton mass. We study the scaling of the blow-up time t∗ in dependence of the small dispersion parameter ɛ and find an exponential dependence t∗(ɛ) and that there is a minimal blow-up time t0∗ greater than the critical time of the corresponding Hopf solution for ɛ → 0. To study the cases with blow-up in detail, we apply the first dynamic rescaling for generalized Korteweg-de Vries equations. This allows to identify the type of the singularity.
49 CFR 178.819 - Vibration test.
Code of Federal Regulations, 2014 CFR
2014-10-01
... 49 Transportation 3 2014-10-01 2014-10-01 false Vibration test. 178.819 Section 178.819... Vibration test. (a) General. The vibration test must be conducted for the qualification of all rigid IBC design types. Flexible IBC design types must be capable of withstanding the vibration test. (b)...
49 CFR 178.819 - Vibration test.
Code of Federal Regulations, 2013 CFR
2013-10-01
... 49 Transportation 3 2013-10-01 2013-10-01 false Vibration test. 178.819 Section 178.819... Vibration test. (a) General. The vibration test must be conducted for the qualification of all rigid IBC design types. Flexible IBC design types must be capable of withstanding the vibration test. (b)...
NASA Astrophysics Data System (ADS)
Holm, D. D.; Ivanov, R. I.
2010-12-01
The Lax pair formulation of the two-component Camassa-Holm equation (CH2) is generalized to produce an integrable multi-component family, CH(n, k), of equations with n components and 1 <= |k| <= n velocities. All of the members of the CH(n, k) family show fluid-dynamics properties with coherent solitons following particle characteristics. We determine their Lie-Poisson Hamiltonian structures and give numerical examples of their soliton solution behaviour. We concentrate on the CH(2, k) family with one or two velocities, including the CH(2, -1) equation in the Dym position of the CH2 hierarchy. A brief discussion of the CH(3, 1) system reveals the underlying graded Lie-algebraic structure of the Hamiltonian formulation for CH(n, k) when n >= 3. Fondly recalling our late friend Jerry Marsden.
NASA Astrophysics Data System (ADS)
Holmedal, Bjørn; Osmundsen, Elisa; Du, Qiang
2016-01-01
Particles precipitated during aging treatments often have non-spherical shapes, e.g., needles or plates, while in the classical Kampmann-Wagner Numerical (KWN) precipitation model, it is assumed that the particles are of spherical shape. This model is here generalized resulting in two correction factors accounting for the effects induced by the particles' non-spherical shape on their growth kinetics. The first one is for the correction of the growth rate. It is derived from the approximate solution of the diffusion problem on spheroidal coordinate and verified by the three-dimensional numerical solutions for cuboid particles. The second factor is for the energetic correction due to the particle surface curvature. It is derived from chemical potential equality (or Gibbs energy minimization principle) at equilibrium for non-spherical particles and provides a correction factor for the Gibbs-Thomson effect. In the accompanying paper, the two correction factors are implemented into a multi-component KWN modeling framework, and the resulting improvements on the model's predictive power are demonstrated.
Vertical Distribution of Vibrationally Excited Hydroxyl
NASA Astrophysics Data System (ADS)
Grygalashvyly, Mykhaylo; Becker, Erich; Sonnemann, Gerd
2016-04-01
Knowledge about the vertical distribution of the vibrationally excited states of hydroxyl (OH*) is important for the interpretation of airglow measurements with respect to dynamical processes in the mesopause region. We derive an approximate analytical expression for the distribution of OH* that highlights the dependence on atomic oxygen and temperature. In addition, we use an advanced numerical model for the formation and relaxation of OH* and investigate the distributions of the different vibrationally exited states of OH*. For the production of OH*, the model includes the reaction of atomic hydrogen with ozone, as well as the reaction of atomic oxygen with hydroperoxy radicals. As loss processes we include 1) deactivation by atomic oxygen, molecular oxygen, and molecular nitrogen, 2) spontaneous emission, and 3) loss due to chemical reaction with atomic oxygen. All these processes take the dependence on the vibrational number into account. The quenching by molecular and atomic oxygen is parameterized by a multi-quantum relaxation scheme. This diagnostic model for OH* has been implemented as part of a chemistry-transport model that is driven by the dynamics simulated with the KMCM (Kühlungsborn Mechanistic general Circulation Model). Numerical results confirm that emission from excited states with higher vibrational number is weaker and emanates from higher altitudes. In addition we find that the OH*-peak altitudes depend significantly on season and latitude. This behavior is mainly controlled by the corresponding variations of atomic oxygen and temperature, as is also confirmed by the aforementioned approximate theory.
NASA Astrophysics Data System (ADS)
Sun, Liangming; Xie, Weiping; He, Xingwen; Hayashikawa, Toshiro
2016-03-01
In this study a 3D numerical analysis approach is developed to predict the ground vibration around rigid-frame viaducts induced by running high-speed trains. The train-bridge-ground interaction system is divided into two subsystems: the train-bridge interaction and the soil-structure interaction. First, the analytical program to simulate bridge vibration with consideration of train-bridge interaction is developed to obtain the vibration reaction forces at the pier bottoms. The highspeed train is described by a multi-DOFs vibration system and the rigid-frame viaduct is modeled with 3D beam elements. Second, applying these vibration reaction forces as input external excitations, the ground vibration is simulated by using a general-purpose program that includes soil-structure interaction effects. The validity of the analytical procedure is confirmed by comparing analytical and experimental results. The characteristics of high-speed train-induced vibrations, including the location of predominant vibration, are clarified. Based on this information a proposed vibration countermeasure using steel strut and new barrier is found effective in reducing train-induced vibrations and it satisfies environmental vibration requirements. The vibration screening efficiency is evaluated by reduction VAL based on 1/3 octave band spectral analysis.
NASA Astrophysics Data System (ADS)
Pessoa, Renato; Castro, Marcos A.; Amaral, Orlando A. V.; Fonseca, Tertius L.
2005-08-01
In this work we report results of CPHF calculations, including vibrational corrections, for the dynamic polarizability and first hyperpolarizability of the LiF molecule. Vibrational corrections were computed through the Numerov-Cooley and perturbation-theoretic methods. Comparison between the results obtained using both methods shows that first-order perturbation-theoretic provides a good approximation for the zpva contribution. The double-harmonic-oscillator approximation is reliable for the pv contributions of αzz, βxxz and βzxx, while [ μα] 0 + [ μ3] I is a good approximation for βzzz. CCSD(T) results obtained for the pv contribution show that the electron correlation effects are small for the polarizability but significant for the first hyperpolarizability.
Malik, Suheel Abdullah; Qureshi, Ijaz Mansoor; Amir, Muhammad; Malik, Aqdas Naveed; Haq, Ihsanul
2015-01-01
In this paper, a new heuristic scheme for the approximate solution of the generalized Burgers'-Fisher equation is proposed. The scheme is based on the hybridization of Exp-function method with nature inspired algorithm. The given nonlinear partial differential equation (NPDE) through substitution is converted into a nonlinear ordinary differential equation (NODE). The travelling wave solution is approximated by the Exp-function method with unknown parameters. The unknown parameters are estimated by transforming the NODE into an equivalent global error minimization problem by using a fitness function. The popular genetic algorithm (GA) is used to solve the minimization problem, and to achieve the unknown parameters. The proposed scheme is successfully implemented to solve the generalized Burgers'-Fisher equation. The comparison of numerical results with the exact solutions, and the solutions obtained using some traditional methods, including adomian decomposition method (ADM), homotopy perturbation method (HPM), and optimal homotopy asymptotic method (OHAM), show that the suggested scheme is fairly accurate and viable for solving such problems. PMID:25811858
NASA Astrophysics Data System (ADS)
Casavant, D.; Brodsky, I.; MacDougall, G. J.
Many important details regarding magnetism in a material can be inferred from the magnetic excitation spectrum, and in this context, general calculations of the classical spinwave spectrum are often necessary. Beyond the simplest of lattices, however, it is difficult to numerically determine the full spinwave spectrum, due primarily to the non-linearity of the problem. In this talk, I will present MATLAB code, developed over the last few years at the University of Illinois, that calculates the dispersions of spinwave excitations out of an arbitrarily defined ordered spin system. The calculation assumes a standard Heisenberg exchange Hamiltonian with the incorporation of a single-ion anisotropy term which can be varied site-by-site and can also simulate the application of an applied field. An overview of the calculation method and the structure of the code will be given, with emphasis on its general applicability. Extensions to the code enable the simulation of both single-crystal and powder-averaged neutron scattering intensity patterns. As a specfic example, I will present the calculated neutron scattering spectrum for powders of CoV2O4, where good agreement between the simulated and experimental data suggests a self-consistent picture of the low-temperature magnetism.
Free vibrations of delaminated beams
NASA Technical Reports Server (NTRS)
Shen, M.-H. H.; Grady, J. E.
1992-01-01
Free vibration of laminated composite beams is studied. The effect of interply delaminations on natural frequencies and mode shapes is evaluated both analytically and experimentally. A generalized vibrational principle is used to formulate the equation of motion and associated boundary conditions for the free vibration of a composite beam with a delamination of arbitrary size and location. The effect of coupling between longitudinal vibration and bending vibration is considered. This coupling effect is shown to significantly affect the calculated natural frequencies and mode shapes of the delaminated beam.
NASA Astrophysics Data System (ADS)
Noori, M.; Davoodi, H.; Saffar, A.
1988-12-01
The cumulant-neglect closure scheme independently developed by Ibrahim and Lin is extended for determining the stationary and non-stationary response of non-linear systems with hysteric restoring force characteristics. The method is applied to the analysis of a hysteresis and model with strength and/or stiffness degradation capabilities. This model has been studied in the past by Baber and Wen for the analysis of hysterically degrading systems using equivalent linearization. The same model has also been used for stochastic seismic performance evaluation of reinforced concrete buildings. Response statistics obtained for the model by using this closure scheme are compared with results of equivalent linearization via Monte Carlo simulation. The study performed, for a wide range of degradation parameters and input power spectral density levels, shows that the Gaussian responses obtained by this approach are identical with the linearized results. This general approximation technique, however, can provide information on higher order statistics for hysteretic systems. These non-Gaussian statistics have not been made available so far by the existing approximation techniques. In this paper the Gaussian statistics are presented.
NASA Astrophysics Data System (ADS)
Borromeo, M.; Marchesoni, F.
2006-01-01
Transport in one-dimensional symmetric devices can be activated by the combination of thermal noise and a biharmonic drive. For the study case of an overdamped Brownian particle diffusing on a periodic one-dimensional substrate, we distinguish two apparently different biharmonic regimes: (i) Harmonic mixing, where the two drive frequencies are commensurate and of the order of some intrinsic relaxation rate. Earlier predictions based on perturbation expansions seem inadequate to interpret our simulation results; (ii) Vibrational mixing, where one harmonic drive component is characterized by high frequency but finite amplitude-to-frequency ratio. Its effect on the device response to either a static or a low-frequency additional input signal is accurately reproduced by rescaling each spatial Fourier component of the substrate potential, separately. Contrary to common wisdom, based on the linear response theory, we show that extremely high-frequency modulations can indeed influence the response of slowly (or dc) operated devices, with potential applications in sensor technology and cellular physiology. Finally, the mixing of two high-frequency beating signal is also investigated both numerically and analytically.
NASA Astrophysics Data System (ADS)
Xu, Feng; Chen, Wen-Li; Xiao, Yi-Qing; Li, Hui; Ou, Jin-Ping
2014-01-01
In the present paper, the commercial CFD code "Fluent" was employed to perform 2-D simulations of an entire process that included the flow around a fixed circular cylinder, the oscillating cylinder (vortex-induced vibration, VIV) and the oscillating cylinder subjected to shape control by a traveling wave wall (TWW) method. The study mainly focused on using the TWW control method to suppress the VIV of an elastically supported circular cylinder with two degrees of freedom at a low Reynolds number of 200. The cross flow (CF) and the inline flow (IL) displacements, the centroid motion trajectories and the lift and drag forces of the cylinder that changed with the frequency ratios were analyzed in detail. The results indicate that a series of small-scale vortices will be formed in the troughs of the traveling wave located on the rear part of the circular cylinder; these vortices can effectively control the flow separation from the cylinder surface, eliminate the oscillating wake and suppress the VIV of the cylinder. A TWW starting at the initial time or at some time halfway through the time interval can significantly suppress the CF and IL vibrations of the cylinder and can remarkably decrease the fluctuations of the lift coefficients and the average values of the drag coefficients; however, it will simultaneously dramatically increase the fluctuations of the drag coefficients.
Formulation of numerical procedures for dynamic analysis of spinning structures
NASA Technical Reports Server (NTRS)
Gupta, K. K.
1986-01-01
The paper presents the descriptions of recently developed numerical algorithms that prove to be useful for the solution of the free vibration problem of spinning structures. First, a generalized procedure for the computation of nodal centrifugal forces in a finite element owing to any specified spin rate is derived in detail. This is followed by a description of an improved eigenproblem solution procedure that proves to be economical for the free vibration analysis of spinning structures. Numerical results are also presented which indicate the efficacy of the currently developed procedures.
Numerical study of the structural and vibrational properties of amorphous Ta2O5 and TiO2-doped Ta2O5
NASA Astrophysics Data System (ADS)
Damart, T.; Coillet, E.; Tanguy, A.; Rodney, D.
2016-05-01
Using classical molecular dynamics simulations, we synthesized amorphous Ta2O5 and amorphous TiO2-doped Ta2O5. We show that Ta2O5 is composed primarily of six-folded Ta atoms forming octahedra that are either organized in chain-like structures or share edges or faces. When Ta2O5 is doped with TiO2, Ti atoms form equally five- and six-folded polyhedra that perturb but do not break the network structure of the glass. Performing a vibrational eigenmode analysis and projecting the eigenmodes on the rocking, stretching, and bending motions of the Ta-2O and Ta-3O bonds, we provide an atomic-scale analysis that substantiates the interpretations of Raman spectra of amorphous Ta2O5. This eigenmode analysis also reveals the key role played by Ti atoms in the 5 to 12 THz range.
Torsional vibration of aircraft engines
NASA Technical Reports Server (NTRS)
Lurenbaum, Karl
1932-01-01
Exhaustive torsional-vibration investigations are required to determine the reliability of aircraft engines. A general outline of the methods used for such investigations and of the theoretical and mechanical means now available for this purpose is given, illustrated by example. True vibration diagrams are usually obtained from vibration measurements on the completed engine. Two devices for this purpose and supplementing each other, the D.V.L. torsiograph and the D.V.L. torsion recorder, are described in this report.
Numerical simulation of 137Cs and (239,240)Pu concentrations by an ocean general circulation model.
Tsumune, Daisuke; Aoyama, Michio; Hirose, Katsumi
2003-01-01
We simulated the spatial distributions and the temporal variations of 137Cs and (239,240)Pu concentrations in the ocean by using the ocean general circulation model which was developed by National Center of Atmospheric Research. These nuclides are introduced into seawaters from global fallout due to atmospheric nuclear weapons tests. The distribution of radioactive deposition on the world ocean is estimated from global precipitation data and observed values of annual deposition of radionuclides at the Meteorological Research Institute in Japan and several observed points in New Zealand. Radionuclides from global fallout have been transported by advection, diffusion and scavenging, and this concentration reduces by radioactive decay in the ocean. We verified the results of the model calculations by comparing simulated values of 137Cs and (239,240)Pu in seawater with the observed values included in the Historical Artificial Radionuclides in the HAM database, which has been constructed by the Meteorological Research Institute. The vertical distributions of the calculated 137Cs concentrations were in good agreement and are in good agreement with the observed profiles in the 1960s up to 250 m, in the 1970s up to 500 m, in the 1980s up to 750 m and in the 1990s up to 750 m. However, the calculated 137Cs concentrations were underestimated compared with the observed 137Cs at the deeper layer. This may suggest other transport processes of 137Cs to deep waters. The horizontal distributions of 137Cs concentrations in surface water could be simulated. A numerical tracer release experiment was performed to explain the horizontal distribution pattern. A maximum (239,240)Pu concentration layer occurs at an intermediate depth for both observed and calculated values, which is formed by particle scavenging. The horizontal distributions of the calculated (239,240)Pu concentrations in surface water could be simulated by considering the scavenging effect. PMID:12860090
NASA Technical Reports Server (NTRS)
Bastin, Paul
1990-01-01
Viewgraphs on vibration isolation are presented. Techniques to control and isolate centrifuge disturbances were identified. Topics covered include: disturbance sources in the microgravity environment; microgravity assessment criteria; life sciences centrifuge; flight support equipment for launch; active vibration isolation system; active balancing system; and fuzzy logic control.
Hauth, J.J.
1962-07-01
A method of compacting a powder in a metal container is described including the steps of vibrating the container at above and below the resonant frequency and also sweeping the frequency of vibration across the resonant frequency several times thereby following the change in resonant frequency caused by compaction of the powder. (AEC)
2011-01-01
By homing in on the distribution patterns of electrons around an atom, a team of scientists team with Berkeley Lab's Molecular Foundry showed how certain vibrations from benzene thiol cause electrical charge to "slosh" onto a gold surface (left), while others do not (right). The vibrations that cause this "sloshing" behavior yield a stronger SERS signal.
A shell model for tyre belt vibrations
NASA Astrophysics Data System (ADS)
Lecomte, C.; Graham, W. R.; Dale, M.
2010-05-01
We present a new formulation for the prediction of tyre belt vibrations in the frequency range 0-500 Hz. Our representation includes the effects of belt width, curvature and anisotropy, and also explicitly models the tyre sidewalls. Many of the associated numerical parameters are fixed by physical considerations; the remainder require empirical input. A systematic and general approach to this problem is developed, and illustrated for the specific example of a Goodyear Wrangler tyre. The resulting predictions for the radial response to radial forcing show good correspondence with experiment up to 300 Hz, and satisfactory agreement up to 1 kHz.
NASA Astrophysics Data System (ADS)
Malagnini, Luca; Dreger, Douglas S.
2016-03-01
Although optimal, computing the moment tensor solution is not always a viable option for the calculation of the size of an earthquake, especially for small events (say, below MW 2.0). Here we show an alternative approach to the calculation of the moment-rate spectra of small earthquakes, and thus of their scalar moments, that uses a network-based calibration of crustal wave propagation. The method works best when applied to a relatively small crustal volume containing both the seismic sources and the recording sites. In this study we present the calibration of the crustal volume monitored by the High Resolution Seismic Network (HRSN), along the San Andreas Fault (SAF) at Parkfield. After the quantification of the attenuation parameters within the crustal volume under investigation, we proceed to the spectral correction of the observed Fourier amplitude spectra for the 100 largest events in our data set. Multiple estimates of seismic moment for the all events (1811 events total) are obtained by calculating the ratio of rms-averaged spectral quantities based on the peak values of the ground velocity in the time domain, as they are observed in narrowband-filtered time series. The mathematical operations allowing the described spectral ratios are obtained from Random Vibration Theory (RVT). Due to the optimal conditions of the HRSN, in terms of signal-to-noise ratios, our network-based calibration allows the accurate calculation of seismic moments down to MW< 0. However, because the HRSN is equipped only with borehole instruments, we define a frequency-dependent Generalized Free-Surface Effect (GFSE), to be used instead of the usual free-surface constant F = 2. Our spectral corrections at Parkfield need a different GFSE for each side of the SAF, which can be quantified by means of the analysis of synthetic seismograms. The importance of the GFSE of borehole instruments increases for decreasing earthquake's size, because for smaller earthquakes the bandwidth available
Mobile high frequency vibrator system
Fair, D.W.; Buller, P.L.
1985-01-08
A carrier mounted seismic vibrator system that is primarily adapted for generation of high force, high frequency seismic energy into an earth medium. The apparatus includes first and second vibrators as supported by first and second lift systems disposed in tandem juxtaposition generally centrally in said vehicle, and the lift systems are designed to maintain equal hold-down force on the vibrator coupling baseplates without exceeding the weight of the carrier vehicle. The juxtaposed vibrators are then energized in synchronized relationship to propagate increased amounts of higher frequency seismic energy into an earth medium.
NASA Astrophysics Data System (ADS)
Hussein, M. F. M.; Hunt, H. E. M.
2009-03-01
This paper presents a new method for modelling floating-slab tracks with discontinuous slabs in underground railway tunnels. The track is subjected to a harmonic load moving with a constant velocity. The model consists of two sub-models. The first is an infinite track with periodic double-beam unit formulated as a periodic infinite structure. The second is modelled with a new version of the Pipe-in-Pipe (PiP) model that accounts for a tunnel wall embedded in a half-space. The two sub-models are coupled by writing the force transmitted from the track to the tunnel as a continuous function using Fourier series representation and satisfying the compatibility condition. The displacements at the free surface are calculated for a track with discontinuous slab and compared with those of a track with continuous slab. The results show that the far-field vibration can be significantly increased due to resonance frequencies of slabs for tracks with discontinuous slabs.
Lerwill, W.E.
1980-09-16
Apparatus for generating vibrations in a medium, such as the ground, comprises a first member which contacts the medium, means , preferably electromagnetic, which includes two relatively movable members for generating vibrations in the apparatus and means operatively connecting the said two members to said first member such that the relatively amplitudes of the movements of said three members can be adjusted to match the impedances of the apparatus and the medium.
14 CFR 29.907 - Engine vibration.
Code of Federal Regulations, 2011 CFR
2011-01-01
... 14 Aeronautics and Space 1 2011-01-01 2011-01-01 false Engine vibration. 29.907 Section 29.907... STANDARDS: TRANSPORT CATEGORY ROTORCRAFT Powerplant General § 29.907 Engine vibration. (a) Each engine must be installed to prevent the harmful vibration of any part of the engine or rotorcraft. (b)...
14 CFR 27.907 - Engine vibration.
Code of Federal Regulations, 2014 CFR
2014-01-01
... 14 Aeronautics and Space 1 2014-01-01 2014-01-01 false Engine vibration. 27.907 Section 27.907... STANDARDS: NORMAL CATEGORY ROTORCRAFT Powerplant General § 27.907 Engine vibration. (a) Each engine must be installed to prevent the harmful vibration of any part of the engine or rotorcraft. (b) The addition of...
14 CFR 29.907 - Engine vibration.
Code of Federal Regulations, 2013 CFR
2013-01-01
... 14 Aeronautics and Space 1 2013-01-01 2013-01-01 false Engine vibration. 29.907 Section 29.907... STANDARDS: TRANSPORT CATEGORY ROTORCRAFT Powerplant General § 29.907 Engine vibration. (a) Each engine must be installed to prevent the harmful vibration of any part of the engine or rotorcraft. (b)...
14 CFR 27.907 - Engine vibration.
Code of Federal Regulations, 2013 CFR
2013-01-01
... 14 Aeronautics and Space 1 2013-01-01 2013-01-01 false Engine vibration. 27.907 Section 27.907... STANDARDS: NORMAL CATEGORY ROTORCRAFT Powerplant General § 27.907 Engine vibration. (a) Each engine must be installed to prevent the harmful vibration of any part of the engine or rotorcraft. (b) The addition of...
49 CFR 178.985 - Vibration test.
Code of Federal Regulations, 2014 CFR
2014-10-01
... 49 Transportation 3 2014-10-01 2014-10-01 false Vibration test. 178.985 Section 178.985... Packagings § 178.985 Vibration test. (a) General. All rigid Large Packaging and flexible Large Packaging design types must be capable of withstanding the vibration test. (b) Test method. (1) A sample...
14 CFR 29.907 - Engine vibration.
Code of Federal Regulations, 2014 CFR
2014-01-01
... 14 Aeronautics and Space 1 2014-01-01 2014-01-01 false Engine vibration. 29.907 Section 29.907... STANDARDS: TRANSPORT CATEGORY ROTORCRAFT Powerplant General § 29.907 Engine vibration. (a) Each engine must be installed to prevent the harmful vibration of any part of the engine or rotorcraft. (b)...
49 CFR 178.985 - Vibration test.
Code of Federal Regulations, 2011 CFR
2011-10-01
... 49 Transportation 3 2011-10-01 2011-10-01 false Vibration test. 178.985 Section 178.985... Packagings § 178.985 Vibration test. (a) General. All rigid Large Packaging and flexible Large Packaging design types must be capable of withstanding the vibration test. (b) Test method. (1) A sample...
49 CFR 178.985 - Vibration test.
Code of Federal Regulations, 2012 CFR
2012-10-01
... 49 Transportation 3 2012-10-01 2012-10-01 false Vibration test. 178.985 Section 178.985... Packagings § 178.985 Vibration test. (a) General. All rigid Large Packaging and flexible Large Packaging design types must be capable of withstanding the vibration test. (b) Test method. (1) A sample...
14 CFR 27.907 - Engine vibration.
Code of Federal Regulations, 2012 CFR
2012-01-01
... 14 Aeronautics and Space 1 2012-01-01 2012-01-01 false Engine vibration. 27.907 Section 27.907... STANDARDS: NORMAL CATEGORY ROTORCRAFT Powerplant General § 27.907 Engine vibration. (a) Each engine must be installed to prevent the harmful vibration of any part of the engine or rotorcraft. (b) The addition of...
14 CFR 29.907 - Engine vibration.
Code of Federal Regulations, 2012 CFR
2012-01-01
... 14 Aeronautics and Space 1 2012-01-01 2012-01-01 false Engine vibration. 29.907 Section 29.907... STANDARDS: TRANSPORT CATEGORY ROTORCRAFT Powerplant General § 29.907 Engine vibration. (a) Each engine must be installed to prevent the harmful vibration of any part of the engine or rotorcraft. (b)...
49 CFR 178.985 - Vibration test.
Code of Federal Regulations, 2013 CFR
2013-10-01
... 49 Transportation 3 2013-10-01 2013-10-01 false Vibration test. 178.985 Section 178.985... Packagings § 178.985 Vibration test. (a) General. All rigid Large Packaging and flexible Large Packaging design types must be capable of withstanding the vibration test. (b) Test method. (1) A sample...
14 CFR 29.907 - Engine vibration.
Code of Federal Regulations, 2010 CFR
2010-01-01
... 14 Aeronautics and Space 1 2010-01-01 2010-01-01 false Engine vibration. 29.907 Section 29.907... STANDARDS: TRANSPORT CATEGORY ROTORCRAFT Powerplant General § 29.907 Engine vibration. (a) Each engine must be installed to prevent the harmful vibration of any part of the engine or rotorcraft. (b)...
NASA Astrophysics Data System (ADS)
Gupta, Amita; Singh, Ranvir; Ahmad, Amir; Kumar, Mahesh
2003-10-01
Today, vibration sensors with low and medium sensitivities are in great demand. Their applications include robotics, navigation, machine vibration monitoring, isolation of precision equipment & activation of safety systems e.g. airbags in automobiles. Vibration sensors have been developed at SSPL, using silicon micromachining to sense vibrations in a system in the 30 - 200 Hz frequency band. The sensing element in the silicon vibration sensor is a seismic mass suspended by thin silicon hinges mounted on a metallized glass plate forming a parallel plate capacitor. The movement of the seismic mass along the vertical axis is monitored to sense vibrations. This is obtained by measuring the change in capacitance. The movable plate of the parallel plate capacitor is formed by a block connected to a surrounding frame by four cantilever beams located on sides or corners of the seismic mass. This element is fabricated by silicon micromachining. Several sensors in the chip sizes 1.6 cm x 1.6 cm, 1 cm x 1 cm and 0.7 cm x 0.7 cm have been fabricated. Work done on these sensors, techniques used in processing and silicon to glass bonding are presented in the paper. Performance evaluation of these sensors is also discussed.
NASA Astrophysics Data System (ADS)
Matthews, Devin A.; Gong, Justin Z.; Stanton, John F.
2014-06-01
The derivation of analytic expressions for vibrational and rovibrational constants, for example the anharmonicity constants χij and the vibration-rotation interaction constants α^B_r, from second-order vibrational perturbation theory (VPT2) can be accomplished with pen and paper and some practice. However, the corresponding quantities from fourth-order perturbation theory (VPT4) are considerably more complex, with the only known derivations by hand extensively using many layers of complicated intermediates and for rotational quantities requiring specialization to orthorhombic cases or the form of Watson's reduced Hamiltonian. We present an automatic computer program for generating these expressions with full generality based on the adaptation of an existing numerical program based on the sum-over-states representation of the energy to a computer algebra context. The measures taken to produce well-simplified and factored expressions in an efficient manner are discussed, as well as the framework for automatically checking the correctness of the generated equations.
NASA Technical Reports Server (NTRS)
Cao, Yiding; Faghri, Amir; Chang, Won Soon
1989-01-01
An enthalpy transforming scheme is proposed to convert the energy equation into a nonlinear equation with the enthalpy, E, being the single dependent variable. The existing control-volume finite-difference approach is modified so it can be applied to the numerical performance of Stefan problems. The model is tested by applying it to a three-dimensional freezing problem. The numerical results are in agreement with those existing in the literature. The model and its algorithm are further applied to a three-dimensional moving heat source problem showing that the methodology is capable of handling complicated phase-change problems with fixed grids.
Force Limited Vibration Testing
NASA Technical Reports Server (NTRS)
Scharton, Terry; Chang, Kurng Y.
2005-01-01
This slide presentation reviews the concept and applications of Force Limited Vibration Testing. The goal of vibration testing of aerospace hardware is to identify problems that would result in flight failures. The commonly used aerospace vibration tests uses artificially high shaker forces and responses at the resonance frequencies of the test item. It has become common to limit the acceleration responses in the test to those predicted for the flight. This requires an analysis of the acceleration response, and requires placing accelerometers on the test item. With the advent of piezoelectric gages it has become possible to improve vibration testing. The basic equations have are reviewed. Force limits are analogous and complementary to the acceleration specifications used in conventional vibration testing. Just as the acceleration specification is the frequency spectrum envelope of the in-flight acceleration at the interface between the test item and flight mounting structure, the force limit is the envelope of the in-flight force at the interface . In force limited vibration tests, both the acceleration and force specifications are needed, and the force specification is generally based on and proportional to the acceleration specification. Therefore, force limiting does not compensate for errors in the development of the acceleration specification, e.g., too much conservatism or the lack thereof. These errors will carry over into the force specification. Since in-flight vibratory force data are scarce, force limits are often derived from coupled system analyses and impedance information obtained from measurements or finite element models (FEM). Fortunately, data on the interface forces between systems and components are now available from system acoustic and vibration tests of development test models and from a few flight experiments. Semi-empirical methods of predicting force limits are currently being developed on the basis of the limited flight and system test
NASA Technical Reports Server (NTRS)
Messaro. Semma; Harrison, Phillip
2010-01-01
Ares I Zonal Random vibration environments due to acoustic impingement and combustion processes are develop for liftoff, ascent and reentry. Random Vibration test criteria for Ares I Upper Stage pyrotechnic components are developed by enveloping the applicable zonal environments where each component is located. Random vibration tests will be conducted to assure that these components will survive and function appropriately after exposure to the expected vibration environments. Methodology: Random Vibration test criteria for Ares I Upper Stage pyrotechnic components were desired that would envelope all the applicable environments where each component was located. Applicable Ares I Vehicle drawings and design information needed to be assessed to determine the location(s) for each component on the Ares I Upper Stage. Design and test criteria needed to be developed by plotting and enveloping the applicable environments using Microsoft Excel Spreadsheet Software and documenting them in a report Using Microsoft Word Processing Software. Conclusion: Random vibration liftoff, ascent, and green run design & test criteria for the Upper Stage Pyrotechnic Components were developed by using Microsoft Excel to envelope zonal environments applicable to each component. Results were transferred from Excel into a report using Microsoft Word. After the report is reviewed and edited by my mentor it will be submitted for publication as an attachment to a memorandum. Pyrotechnic component designers will extract criteria from my report for incorporation into the design and test specifications for components. Eventually the hardware will be tested to the environments I developed to assure that the components will survive and function appropriately after exposure to the expected vibration environments.
NASA Astrophysics Data System (ADS)
Malagnini, Luca; Dreger, Douglas S.
2016-07-01
Although optimal, computing the moment tensor solution is not always a viable option for the calculation of the size of an earthquake, especially for small events (say, below Mw 2.0). Here we show an alternative approach to the calculation of the moment-rate spectra of small earthquakes, and thus of their scalar moments, that uses a network-based calibration of crustal wave propagation. The method works best when applied to a relatively small crustal volume containing both the seismic sources and the recording sites. In this study we present the calibration of the crustal volume monitored by the High-Resolution Seismic Network (HRSN), along the San Andreas Fault (SAF) at Parkfield. After the quantification of the attenuation parameters within the crustal volume under investigation, we proceed to the spectral correction of the observed Fourier amplitude spectra for the 100 largest events in our data set. Multiple estimates of seismic moment for the all events (1811 events total) are obtained by calculating the ratio of rms-averaged spectral quantities based on the peak values of the ground velocity in the time domain, as they are observed in narrowband-filtered time-series. The mathematical operations allowing the described spectral ratios are obtained from Random Vibration Theory (RVT). Due to the optimal conditions of the HRSN, in terms of signal-to-noise ratios, our network-based calibration allows the accurate calculation of seismic moments down to Mw < 0. However, because the HRSN is equipped only with borehole instruments, we define a frequency-dependent Generalized Free-Surface Effect (GFSE), to be used instead of the usual free-surface constant F = 2. Our spectral corrections at Parkfield need a different GFSE for each side of the SAF, which can be quantified by means of the analysis of synthetic seismograms. The importance of the GFSE of borehole instruments increases for decreasing earthquake's size because for smaller earthquakes the bandwidth available
Digital vibration control techniques
NASA Technical Reports Server (NTRS)
Chapman, P.; Kim, B. K.; Boctor, W.
1974-01-01
Analog vibration control techniques are reviewed and are compared with digital techniques. The advantages of the digital methods over the analog methods are demonstrated. The following topics are covered: (1) methods of computer-controlled random vibration and reverberation acoustic testing; (2) methods of computer-controlled sinewave vibration testing; and (3) methods of computer-controlled shock testing. Basic concepts are stressed rather than specific techniques or equipment. General algorithms are described in the form of block diagrams and flow diagrams. Specific problems and potential problems are discussed. The material is computer sciences oriented but is kept at a level that facilitates an understanding of the basic concepts of computer-controlled induced environmental test systems.
Electron-vibration relaxation in oxygen plasmas
NASA Astrophysics Data System (ADS)
Laporta, V.; Heritier, K. L.; Panesi, M.
2016-06-01
An ideal chemical reactor model is used to study the vibrational relaxation of oxygen molecules in their ground electronic state, X3 Σg-, in presence of free electrons. The model accounts for vibrational non-equilibrium between the translational energy mode of the gas and the vibrational energy mode of individual molecules. The vibrational levels of the molecules are treated as separate species, allowing for non-Boltzmann distributions of their population. The electron and vibrational temperatures are varied in the range [0-20,000] K. Numerical results show a fast energy transfer between oxygen molecules and free electron, which causes strong deviation of the vibrational distribution function from Boltzmann distribution, both in heating and cooling conditions. Comparison with Landau-Teller model is considered showing a good agreement for electron temperature range [2000-12,000] K. Finally analytical fit of the vibrational relaxation time is given.
Optimum vibration control of flexible beams by piezo-electric actuators
NASA Technical Reports Server (NTRS)
Baz, A.; Poh, S.
1987-01-01
The utilization of piezoelectric actuators in controlling the structural vibrations of flexible beams is examined. A Modified Independent Modal Space Control (MIMSC) method is devised to enable the selection of the optimal location, control gains and excitation voltage of the piezoelectric actuators in a way that would minimize the amplitudes of vibrations of beams to which these actuators are bonded, as well as the input control energy necessary to suppress these vibrations. The developed method accounts for the effects that the piezoelectric actuators have on changing the elastic and inertial properties of the flexible beams. Numerical examples are presented to illustrate the application of the developed MIMSC method in minimizing the structural vibrations of beams of different materials when subjected to different loading and end conditions using ceramic or polymeric piezoelectric actuators. The obtained results emphasize the importance of the devised method in designing more realistic active control systems for flexible beams, in particular, and large flexible structures in general.
Liu, J.; Strzalka, J; Tronin, A; Johansson, J; Blasie, J
2009-01-01
We demonstrate that cyano-phenylalanine (PheCN) can be utilized to probe the binding of the inhalational anesthetic halothane to an anesthetic-binding, model ion channel protein hbAP-PheCN. The Trp to PheCN mutation alters neither the a-helical conformation nor the 4-helix bundle structure. The halothane binding properties of this PheCN mutant hbAP-PheCN, based on fluorescence quenching, are consistent with those of the prototype, hbAP1. The dependence of fluorescence lifetime as a function of halothane concentration implies that the diffusion of halothane in the nonpolar core of the protein bundle is one-dimensional. As a consequence, at low halothane concentrations, the quenching of the fluorescence is dynamic, whereas at high concentrations the quenching becomes static. The 4-helix bundle structure present in aqueous detergent solution and at the air-water interface, is preserved in multilayer films of hbAP-PheCN, enabling vibrational spectroscopy of both the protein and its nitrile label (-CN). The nitrile groups' stretching vibration band shifts to higher frequency in the presence of halothane, and this blue-shift is largely reversible. Due to the complexity of this amphiphilic 4-helix bundle model membrane protein, where four PheCN probes are present adjacent to the designed cavity forming the binding site within each bundle, all contributing to the infrared absorption, molecular dynamics (MD) simulation is required to interpret the infrared results. The MD simulations indicate that the blue-shift of -CN stretching vibration induced by halothane arises from an indirect effect, namely an induced change in the electrostatic protein environment averaged over the four probe oscillators, rather than a direct interaction with the oscillators. hbAP-PheCN therefore provides a successful template for extending these investigations of the interactions of halothane with the model membrane protein via vibrational spectroscopy, using cyano-alanine residues to form the
Liu, Jing; Strzalka, Joseph; Tronin, Andrey; Johansson, Jonas S.; Blasie, J. Kent
2009-01-01
We demonstrate that cyano-phenylalanine (PheCN) can be utilized to probe the binding of the inhalational anesthetic halothane to an anesthetic-binding, model ion channel protein hbAP-PheCN. The Trp to PheCN mutation alters neither the α-helical conformation nor the 4-helix bundle structure. The halothane binding properties of this PheCN mutant hbAP-PheCN, based on fluorescence quenching, are consistent with those of the prototype, hbAP1. The dependence of fluorescence lifetime as a function of halothane concentration implies that the diffusion of halothane in the nonpolar core of the protein bundle is one-dimensional. As a consequence, at low halothane concentrations, the quenching of the fluorescence is dynamic, whereas at high concentrations the quenching becomes static. The 4-helix bundle structure present in aqueous detergent solution and at the air-water interface, is preserved in multilayer films of hbAP-PheCN, enabling vibrational spectroscopy of both the protein and its nitrile label (-CN). The nitrile groups' stretching vibration band shifts to higher frequency in the presence of halothane, and this blue-shift is largely reversible. Due to the complexity of this amphiphilic 4-helix bundle model membrane protein, where four PheCN probes are present adjacent to the designed cavity forming the binding site within each bundle, all contributing to the infrared absorption, molecular dynamics (MD) simulation is required to interpret the infrared results. The MD simulations indicate that the blue-shift of -CN stretching vibration induced by halothane arises from an indirect effect, namely an induced change in the electrostatic protein environment averaged over the four probe oscillators, rather than a direct interaction with the oscillators. hbAP-PheCN therefore provides a successful template for extending these investigations of the interactions of halothane with the model membrane protein via vibrational spectroscopy, using cyano-alanine residues to form the
Parameter Reconstruction of Vibration Systems from Partial Eigeninformation
Lin, Matthew M.
2013-01-01
Quadratic matrix polynomials are fundamental to vibration analysis. Because of the predetermined interconnectivity among the constituent elements and the mandatory nonnegativity of the physical parameters, most given vibration systems will impose some inherent structure on the coefficients of the corresponding quadratic matrix polynomials. In the inverse problem of reconstructing a vibration system from its observed or desirable dynamical behavior, respecting the intrinsic structure becomes important and challenging both theoretically and practically. The issue of whether a structured inverse eigenvalue problem is solvable is problem dependent and has to be addressed structure by structure. In an earlier work, physical systems that can be modeled under the paradigm of a serially linked mass-spring system have been considered via specifically formulated inequality systems. In this paper, the framework is generalized to arbitrary generally linked systems. In particular, given any configuration of interconnectivity in a mass-spring system, this paper presents a mechanism that systematically and automatically generates a corresponding inequality system. A numerical approach is proposed to determine whether the inverse problem is solvable and, if it is so, computes the coefficient matrices while providing an estimate of the residual error. The most important feature of this approach is that it is problem independent, that is, the approach is general and robust for any kind of physical configuration. The ideas discussed in this paper have been implemented into a software package by which some numerical experiments are reported. PMID:23966750
Development of a numerical procedure to map a general 3-d body onto a near-circle
NASA Technical Reports Server (NTRS)
Hommel, M. J.
1986-01-01
Conformal mapping is a classical technique utilized for solving problems in aerodynamics and hydrodynamics. Conformal mapping is utilized in the construction of grids around airfoils, engine inlets and other aircraft configurations. These shapes are transformed onto a near-circle image for which the equations of fluid motion are discretized on the mapped plane and solved numerically by utilizing the appropriate techniques. In comparison to other grid-generation techniques such as algerbraic or differential type, conformal mapping offers an analytical and accurate form even if the grid deformation is large. One of the most appealing features is that the grid can be constrained to remain orthogonal to the body after the transformation. Hence, the grid is suitable for analyzing the supersonic flow past a blunt object. The associated shock as a coordinate surface adjusts its position in the course of computation until convergence is reached. The present work applied conformal mapping to 3-D bodies with no axis of symmetry such as the Aerobraking Flight Experiment (AFE) vehicle, transforming the AFE shape onto a near-circle image. A numerical procedure and code are used to generate grids around the AFE body.
Natural vibration dynamics of Rainbow Bridge, Utah
NASA Astrophysics Data System (ADS)
Moore, J. R.; Thorne, M. S.; Wood, J. R.; Doyle, S.; Stanfield, E.; White, B.
2015-12-01
We measured resonant frequencies of Rainbow Bridge, Utah, one of the world's longest rock spans, during a field experiment recording ambient vibration data. Measurements were generated over 20 hours on March 23-24, 2015 using two broadband three-component seismometers placed on the bridge, and compared to concurrent data from nearby reference stations 20 and 220 m distant. We identified seven distinct modes of vibration for Rainbow Bridge between 1 and 6 Hz. Data for each resonant frequency was then analyzed to determine the frequency-dependent polarization vector in an attempt to clarify mode shapes; e.g. the fundamental mode represents out-of-plane horizontal flexure. We compared experimental data to results of 3D numerical modal analysis, using a new photogrammetric model of Rainbow Bridge generated in this study imported into COMSOL Multiphysics. Results compare well with measured data for seven of the first eight modeled modes, matching vibrational frequencies and polarization orientations generally within 10%. Only predicted mode 6 was not explicitly apparent in our experimental data. Large site-to-reference spectral ratios resolved from experimental data indicate high amplification on the bridge as compared to nearby bedrock.
NASA Technical Reports Server (NTRS)
Bozeman, Richard J., Jr. (Inventor)
1990-01-01
The invention relates to monitoring circuitry for the real time detection of vibrations of a predetermined frequency and which are greater than a predetermined magnitude. The circuitry produces an instability signal in response to such detection. The circuitry is particularly adapted for detecting instabilities in rocket thrusters, but may find application with other machines such as expensive rotating machinery, or turbines. The monitoring circuitry identifies when vibration signals are present having a predetermined frequency of a multi-frequency vibration signal which has an RMS energy level greater than a predetermined magnitude. It generates an instability signal only if such a vibration signal is identified. The circuitry includes a delay circuit which responds with an alarm signal only if the instability signal continues for a predetermined time period. When used with a rocket thruster, the alarm signal may be used to cut off the thruster if such thruster is being used in flight. If the circuitry is monitoring tests of the thruster, it generates signals to change the thruster operation, for example, from pulse mode to continuous firing to determine if the instability of the thruster is sustained once it is detected.
NASA Technical Reports Server (NTRS)
Moura, A. D.; Shukla, J.
1981-01-01
The establishment of a thermally direct local circulation which has its ascending branch at about 10 deg N and its descending branch over northeast Brazil and the adjoining oceanic region is proposed as a possible mechanism for the occurrence of severe droughts over this Brazilian region. The driving for this anomalous circulation is provided by enhanced moist convection due to the effect of warmer sea surface anomalies over the northern tropical Atlantic and cooling associated with colder sea surface temperature anomalies in the southern tropical Atlantic. A simple primitive equation model is used to calculate the frictionally-controlled and thermally-driven circulation due to a prescribed heating function in a resting atmosphere, and a series of numerical experiments are carried out to test the sensitivity of the Goddard Laboratory's model to prescribed sea surface temperature anomalies over the tropical Atlantic.
NASA Astrophysics Data System (ADS)
Makinde, O. D.; Onyejekwe, O. O.
2011-11-01
The steady flow and heat transfer of an electrically conducting fluid with variable viscosity and electrical conductivity between two parallel plates in the presence of a transverse magnetic field is investigated. It is assumed that the flow is driven by combined action of axial pressure gradient and uniform motion of the upper plate. The governing nonlinear equations of momentum and energy transport are solved numerically using a shooting iteration technique together with a sixth-order Runge-Kutta integration algorithm. Solutions are presented in graphical form and given in terms of fluid velocity, fluid temperature, skin friction and heat transfer rate for various parametric values. Our results reveal that the combined effect of magnetic field, viscosity, exponents of variable properties, various fluid and heat transfer dimensionless quantities and the electrical conductivity variation, have significant impact on the hydromagnetic and electrical properties of the fluid.
NASA Astrophysics Data System (ADS)
Eremin, M.; Makarov, P.
2016-04-01
The results of 2D modelling of rock mass elements fracture are shown in the article. The results of modelling are in a good agreement with empirical and theoretical estimations of roof caving steps for the flat-dipping coal seams when the horizons are not so deep (less than 600 m). The estimations of the general and set steps of roof caving are given for the lava conditions at different lengths of the main roof containing sandstone.
NASA Astrophysics Data System (ADS)
May, V.
2002-12-01
To fully account for electron-vibrational coupling and vibrational relaxation in the course of electron motion through a molecular wire a density operator approach is utilized. If combined with a particular projection operator technique a generalized master equation can be derived which governs the populations of the electronic wire states. The respective memory kernels are determined beyond any perturbation theory with respect to the electron-vibrational coupling and can be classified via so-called Liouville space pathways. An ordering of the different contributions to the current-voltage characteristics becomes possible by introducing an electron transmission coefficient which describes ballistic as well as inelastic electron transport through the wire. The general derivations are illustrated by numerical calculations which demonstrate the drastic influence of the electron-vibrational coupling on the wire transmission coefficient as well as on the current-voltage characteristics.
NASA Astrophysics Data System (ADS)
Geroyannis, Vassilis; Tzelati, Eleftheria
In this paper we compute general-relativistic polytropic models simulating rigidly rotating, pulsating neutron stars. These relativistic compact objects, with a radius of $\\sim 10 \\, \\mathrm{km}$ and mass between $\\sim 1.4$ and $3.2$ solar masses, are closely related to pulsars. We emphasize on computing the change in the pulsation eigenfrequencies owing to a rigid rotation, which, in turn, is a decisive issue for studying stability of such objects. In our computations, we keep rotational perturbation terms of up to second order in the angular velocity.
NASA Technical Reports Server (NTRS)
2001-01-01
A Small Business Innovation Research (SBIR) sponsorship from NASA's Dryden Flight Research Center, assisted MetroLaser, of Irvine, California, in the development of a self-aligned laser vibrometer system. VibroMet, capable of measuring surface vibrations in a variety of industries, provides information on the structural integrity and acoustical characteristics of manufactured products. This low-cost, easy-to-use sensor performs vibration measurement from distances of up to three meters without the need for adjustment. The laser beam is simply pointed at the target and the system then uses a compact laser diode to illuminate the surface and to subsequently analyze the reflected light. The motion of the surface results in a Doppler shift that is measured with very high precision. VibroMet is considered one of the many behind-the-scenes tools that can be relied on to assure the quality, reliability and safety of everything from airplane panels to disk brakes
Benedetti, Andrea; Platt, Robert; Atherton, Juli
2014-01-01
Background Over time, adaptive Gaussian Hermite quadrature (QUAD) has become the preferred method for estimating generalized linear mixed models with binary outcomes. However, penalized quasi-likelihood (PQL) is still used frequently. In this work, we systematically evaluated whether matching results from PQL and QUAD indicate less bias in estimated regression coefficients and variance parameters via simulation. Methods We performed a simulation study in which we varied the size of the data set, probability of the outcome, variance of the random effect, number of clusters and number of subjects per cluster, etc. We estimated bias in the regression coefficients, odds ratios and variance parameters as estimated via PQL and QUAD. We ascertained if similarity of estimated regression coefficients, odds ratios and variance parameters predicted less bias. Results Overall, we found that the absolute percent bias of the odds ratio estimated via PQL or QUAD increased as the PQL- and QUAD-estimated odds ratios became more discrepant, though results varied markedly depending on the characteristics of the dataset Conclusions Given how markedly results varied depending on data set characteristics, specifying a rule above which indicated biased results proved impossible. This work suggests that comparing results from generalized linear mixed models estimated via PQL and QUAD is a worthwhile exercise for regression coefficients and variance components obtained via QUAD, in situations where PQL is known to give reasonable results. PMID:24416249
Gakh, G. I.; Rekalo, A. P.; Tomasi-Gustafsson, E.; Boucher, J.; Gakh, A. G.
2011-02-15
A general formalism is developed to calculate the cross section and the polarization observables for the reaction N-bar+N{yields}{pi}+l{sup +}+l{sup -}. The matrix element and the observables are expressed in terms of six scalar amplitudes (complex functions of three kinematical variables) that determine the reaction dynamics. The numerical predictions are given in the frame of a particular model in the kinematical range accessible in the antiproton annihilation at Darmstadt (PANDA) experiment at the Facility for Antiproton and Ion Research (FAIR).
Optimal active vibration absorber - Design and experimental results
NASA Technical Reports Server (NTRS)
Lee-Glauser, Gina; Juang, Jer-Nan; Sulla, Jeffrey L.
1993-01-01
An optimal active vibration absorber can provide guaranteed closed-loop stability and control for large flexible space structures with collocated sensors/actuators. The active vibration absorber is a second-order dynamic system which is designed to suppress any unwanted structural vibration. This can be designed with minimum knowledge of the controlled system. Two methods for optimizing the active vibration absorber parameters are illustrated: minimum resonant amplitude and frequency matched active controllers. The Controls-Structures Interaction Phase-1 Evolutionary Model at NASA LaRC is used to demonstrate the effectiveness of the active vibration absorber for vibration suppression. Performance is compared numerically and experimentally using acceleration feedback.
Optimal active vibration absorber: Design and experimental results
NASA Technical Reports Server (NTRS)
Lee-Glauser, Gina; Juang, Jer-Nan; Sulla, Jeffrey L.
1992-01-01
An optimal active vibration absorber can provide guaranteed closed-loop stability and control for large flexible space structures with collocated sensors/actuators. The active vibration absorber is a second-order dynamic system which is designed to suppress any unwanted structural vibration. This can be designed with minimum knowledge of the controlled system. Two methods for optimizing the active vibration absorber parameters are illustrated: minimum resonant amplitude and frequency matched active controllers. The Controls-Structures Interaction Phase-1 Evolutionary Model at NASA LaRC is used to demonstrate the effectiveness of the active vibration absorber for vibration suppression. Performance is compared numerically and experimentally using acceleration feedback.
Coupling between plate vibration and acoustic radiation
NASA Technical Reports Server (NTRS)
Frendi, Abdelkader; Maestrello, Lucio; Bayliss, Alvin
1992-01-01
A detailed numerical investigation of the coupling between the vibration of a flexible plate and the acoustic radiation is performed. The nonlinear Euler equations are used to describe the acoustic fluid while the nonlinear plate equation is used to describe the plate vibration. Linear, nonlinear, and quasi-periodic or chaotic vibrations and the resultant acoustic radiation are analyzed. We find that for the linear plate response, acoustic coupling is negligible. However, for the nonlinear and chaotic responses, acoustic coupling has a significant effect on the vibration level as the loading increases. The radiated pressure from a plate undergoing nonlinear or chaotic vibrations is found to propagate nonlinearly into the far-field. However, the nonlinearity due to wave propagation is much weaker than that due to the plate vibrations. As the acoustic wave propagates into the far-field, the relative difference in level between the fundamental and its harmonics and subharmonics decreases with distance.
Sekiguchi, Yu-ichirou; Shibata, Masaru
2005-04-15
We perform axisymmetric simulations for gravitational collapse of a massive iron core to a black hole in full general relativity. The iron cores are modeled by {gamma}=4/3 equilibrium polytrope for simplicity. The hydrodynamic equations are solved using a high-resolution shock-capturing scheme with a parametric equation of state. The Cartoon method is adopted for solving the Einstein equations. Simulations are performed for a wide variety of initial conditions changing the mass ({approx_equal}2.0-3.0M{sub {center_dot}}), the angular momentum, the rotational velocity profile of the core, and the parameters of the equations of state which are chosen so that the maximum mass of the cold spherical polytrope is {approx_equal}1.6M{sub {center_dot}}. Then, the criterion for the prompt black hole formation is clarified in terms of the mass and the angular momentum for several rotational velocity profile of the core and equations of state. It is found that (i) with the increase of the thermal energy generated by shocks, the threshold mass for the prompt black hole formation is increased by 20-40%, (ii) the rotational centrifugal force increases the threshold mass by < or approx. 25%, (iii) with the increase of the degree of differential rotation, the threshold mass is also increased, and (iv) the amplification factors shown in the results (i)-(iii) depend sensitively on the equation of state. We also find that the collapse dynamics and the structure of the shock formed at the bounce depend strongly on the stiffness of the adopted equation of state. In particular, as a new feature, a strong bipolar explosion is observed for the collapse of rapidly rotating iron cores with an equation of state which is stiff in subnuclear density and soft in supranuclear density. Gravitational waves are computed in terms of a quadrupole formula. It is also found that the waveform depends sensitively on the equations of state.
Reduced elbow extension torque during vibrations.
Friesenbichler, Bernd; Coza, Aurel; Nigg, Benno M
2012-08-31
Impact sports and vibration platforms trigger vibrations within soft tissues and the skeleton. Although the long-term effects of vibrations on the body have been studied extensively, the acute effects of vibrations are little understood. This study determined the influence of acute vibrations at different frequencies and elbow angles on maximal isometric elbow extension torque and muscle activity. Vibrations were generated by a pneumatic vibrator attached to the lever of a dynamometer, and were applied on the forearm of 15 healthy female subjects. The subjects were instructed to push maximally against the lever at three different elbow angles, while extension torque and muscle activity were quantified and compared between vibration and non-vibration (control) conditions. A change in vibration frequency had no significant effects on torque and muscle activity although vibrations in general decreased the maximal extension torque relative to the control by 1.8% (±5.7%, p>0.05), 7.4% (±7.9%, p<0.01), and 5.0% (±8.2%, p<0.01) at elbow angles of 60°, 90°, and 120°, respectively. Electromyographic activity increased significantly between ∼30% and 40% in both triceps and biceps with vibrations. It is speculated that a similar increase in muscle activity between agonist and antagonist, in combination with an unequal increase in muscle moment arms about the elbow joint, limit the maximal extension torque during exposure to vibrations. This study showed that maximal extension torque decreased during vibration exposure while muscle activity increased and suggests that vibrations may be counterproductive during activities requiring maximal strength but potentially beneficial for strength training. PMID:22771229
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.
Eckart frame vibration-rotation Hamiltonians: Contravariant metric tensor
Pesonen, Janne
2014-02-21
Eckart frame is a unique embedding in the theory of molecular vibrations and rotations. It is defined by the condition that the Coriolis coupling of the reference structure of the molecule is zero for every choice of the shape coordinates. It is far from trivial to set up Eckart kinetic energy operators (KEOs), when the shape of the molecule is described by curvilinear coordinates. In order to obtain the KEO, one needs to set up the corresponding contravariant metric tensor. Here, I derive explicitly the Eckart frame rotational measuring vectors. Their inner products with themselves give the rotational elements, and their inner products with the vibrational measuring vectors (which, in the absence of constraints, are the mass-weighted gradients of the shape coordinates) give the Coriolis elements of the contravariant metric tensor. The vibrational elements are given as the inner products of the vibrational measuring vectors with themselves, and these elements do not depend on the choice of the body-frame. The present approach has the advantage that it does not depend on any particular choice of the shape coordinates, but it can be used in conjunction with all shape coordinates. Furthermore, it does not involve evaluation of covariant metric tensors, chain rules of derivation, or numerical differentiation, and it can be easily modified if there are constraints on the shape of the molecule. Both the planar and non-planar reference structures are accounted for. The present method is particular suitable for numerical work. Its computational implementation is outlined in an example, where I discuss how to evaluate vibration-rotation energies and eigenfunctions of a general N-atomic molecule, the shape of which is described by a set of local polyspherical coordinates.
Eckart frame vibration-rotation Hamiltonians: contravariant metric tensor.
Pesonen, Janne
2014-02-21
Eckart frame is a unique embedding in the theory of molecular vibrations and rotations. It is defined by the condition that the Coriolis coupling of the reference structure of the molecule is zero for every choice of the shape coordinates. It is far from trivial to set up Eckart kinetic energy operators (KEOs), when the shape of the molecule is described by curvilinear coordinates. In order to obtain the KEO, one needs to set up the corresponding contravariant metric tensor. Here, I derive explicitly the Eckart frame rotational measuring vectors. Their inner products with themselves give the rotational elements, and their inner products with the vibrational measuring vectors (which, in the absence of constraints, are the mass-weighted gradients of the shape coordinates) give the Coriolis elements of the contravariant metric tensor. The vibrational elements are given as the inner products of the vibrational measuring vectors with themselves, and these elements do not depend on the choice of the body-frame. The present approach has the advantage that it does not depend on any particular choice of the shape coordinates, but it can be used in conjunction with all shape coordinates. Furthermore, it does not involve evaluation of covariant metric tensors, chain rules of derivation, or numerical differentiation, and it can be easily modified if there are constraints on the shape of the molecule. Both the planar and non-planar reference structures are accounted for. The present method is particular suitable for numerical work. Its computational implementation is outlined in an example, where I discuss how to evaluate vibration-rotation energies and eigenfunctions of a general N-atomic molecule, the shape of which is described by a set of local polyspherical coordinates. PMID:24559332
Vibrational Control of a Nonlinear Elastic Panel
NASA Technical Reports Server (NTRS)
Chow, P. L.; Maestrello, L.
1998-01-01
The paper is concerned with the stabilization of the nonlinear panel oscillation by an active control. The control is actuated by a combination of additive and parametric vibrational forces. A general method of vibrational control is presented for stabilizing panel vibration satisfying a nonlinear beam equation. To obtain analytical results, a perturbation technique is used in the case of weak nonlinearity. Possible application to other types of problems is briefly discussed.
Optimization of a vacuum chamber for vibration measurements
NASA Astrophysics Data System (ADS)
Danyluk, Mike; Dhingra, Anoop
2011-10-01
A 200 °C high vacuum chamber has been built to improve vibration measurement sensitivity. The optimized design addresses two significant issues: (i) vibration measurements under high vacuum conditions and (ii) use of design optimization tools to reduce operating costs. A test rig consisting of a cylindrical vessel with one access port has been constructed with a welded-bellows assembly used to seal the vessel and enable vibration measurements in high vacuum that are comparable with measurements in air. The welded-bellows assembly provides a force transmissibility of 0.1 or better at 15 Hz excitation under high vacuum conditions. Numerical results based on design optimization of a larger diameter chamber are presented. The general constraints on the new design include material yield stress, chamber first natural frequency, vibration isolation performance, and forced convection heat transfer capabilities over the exterior of the vessel access ports. Operating costs of the new chamber are reduced by 50% compared to a preexisting chamber of similar size and function.
Microgravity vibration isolation: Optimal preview and feedback control
NASA Technical Reports Server (NTRS)
Hampton, R. D.; Knospe, C. R.; Grodsinsky, C. M.; Allaire, P. E.; Lewis, D. W.
1992-01-01
In order to achieve adequate low-frequency vibration isolation for certain space experiments an active control is needed, due to inherent passive-isolator limitations. Proposed here are five possible state-space models for a one-dimensional vibration isolation system with a quadratic performance index. The five models are subsets of a general set of nonhomogeneous state space equations which includes disturbance terms. An optimal control is determined, using a differential equations approach, for this class of problems. This control is expressed in terms of constant, Linear Quadratic Regulator (LQR) feedback gains and constant feedforward (preview) gains. The gains can be easily determined numerically. They result in a robust controller and offers substantial improvements over a control that uses standard LQR feedback alone.
Damping in micro-scale generalized thermoelastic circular plate resonators.
Sharma, J N; Sharma, R
2011-04-01
The out-of-plane vibrations of a generalized thermoelastic circular plate are studied under different environmental temperature, plate dimensions and boundary conditions. The analytical expressions for thermoelastic damping of vibration and phase velocity of circumferential surface wave modes are obtained. It is noticed that the damping of vibrations and phase velocities of circumferential surface wave modes significantly depend on thermal relaxation time in addition to thermoelastic coupling in circular plates under resonance conditions. The surface conditions also impose significant effects on the vibrations of such resonators. The expressions for displacement and temperature fields in the plate resonator are also derived and obtained. Some numerical results have also been presented for illustration purpose in case of silicon material plate. PMID:21168892
Diagrammatic Vibrational Coupled-Cluster
NASA Astrophysics Data System (ADS)
Faucheaux, Jacob A.; Hirata, So
2015-06-01
A diagrammatic vibrational coupled-cluster method for calculation of zero-point energies and an equation-of-motion coupled-cluster method for calculation of anharmonic vibrational frequencies are developed. The methods, which we refer to as XVCC and EOM-XVCC respectively, rely on the size-extensive vibrational self-consistient field (XVSCF) method for reference wave functions. The methods retain the efficiency advantages of XVSCF making them suitable for applications to large molecules and solids, while they are numerically shown to accurately predict zero-point energies and frequencies of small molecules as well. In particular, EOM-XVCC is shown to perform well for modes which undergo Fermi resonance where traditional perturbative methods fail. Rules for the systematic generation and interpretation of the XVCC and EOM-XVCC diagrams to any order are presented.
Nieh, C.D.; Zengyan, H.
1986-04-01
Based on the classical beam theory, a simple method for calculating the natural frequency of unequally spanned tubes is presented. The method is suitable for various boundary conditions. Accuracy of the calculations is sufficient for practical applications. This method will help designers and operators estimate the vibration of tubular exchangers. In general, there are three reasons why a tube vibrates in cross flow: vortex shedding, fluid elasticity and turbulent buffeting. No matter which is the cause, the basic reason is that the frequency of exciting force is approximately the same as or equal to the natural frequency of the tube. To prevent the heat exchanger from vibrating, it is necessary to select correctly the shell-side fluid velocity so that the frequency of exciting force is different from the natural frequency of the tube, or to vary the natural frequency of the heat exchanger tube. So precisely determining the natural frequency of the heat exchanger, especially its foundational frequency under various supporting conditions, is of significance.
Vibrational Dynamics of Tricyanomethanide
NASA Astrophysics Data System (ADS)
Weidinger, Daniel; Houchins, Cassidy; Owrutsky, Jeffrey C.
2011-06-01
Time-resolved and steady-state IR spectroscopy have been used to characterize vibrational spectra and energy relaxation dynamics of the CN stretching band of the tricyanomethanide (TCM, C(CN)3-) anion near 2170 Cm-1 in solutions of water, heavy water, methanol, formamide, dimethyl sulfoxide (DMSO) and the ionic liquid 1-butyl methyl imidazolium tetrafluoroborate ([BMIM][BF4]). The band intensity is strong (˜1500 M-1Cm-1) and the vibrational energy relaxation times are relatively long (˜5 ps in water, 12 ps in heavy water, and ˜30 ps in DMSO and [BMIM][BF4]). They are longer than those previously reported for dicyanamide in the same solvents. Although the static TCM frequency generally shifts to higher frequency with more strongly interacting solvents, the shift does not follow the same trend as the vibrational dynamics. The results for the experimental frequencies and intensities agree well with results from ab initio calculations. Proton and electron affinities for TCM are also calculated because they are relevant to potential applications of this anion in low viscosity ionic liquids.
49 CFR 178.985 - Vibration test.
Code of Federal Regulations, 2010 CFR
2010-10-01
... Testing of Large Packagings § 178.985 Vibration test. (a) General. All rigid Large Packaging and flexible Large Packaging design types must be capable of withstanding the vibration test. (b) Test method. (1) A sample Large Packaging, selected at random, must be filled and closed as for shipment. Large...
Li, Junjie; Li, Xiaohu; Iyengar, Srinivasan S
2014-06-10
agreement with experiments. Spectroscopic features are computed using a unified velocity/flux autocorrelation function and include vibrational fundamentals and combination bands. These agree well with experiments and other theories. PMID:26580749
Chaotic vortex induced vibrations
Zhao, J.; Sheridan, J.; Leontini, J. S.; Lo Jacono, D.
2014-12-15
This study investigates the nature of the dynamic response of an elastically mounted cylinder immersed in a free stream. A novel method is utilized, where the motion of the body during a free vibration experiment is accurately recorded, and then a second experiment is conducted where the cylinder is externally forced to follow this recorded trajectory. Generally, the flow response during both experiments is identical. However, particular regimes exist where the flow response is significantly different. This is taken as evidence of chaos in these regimes.
Structural Acoustics and Vibrations
NASA Astrophysics Data System (ADS)
Chaigne, Antoine
This chapter is devoted to vibrations of structures and to their coupling with the acoustic field. Depending on the context, the radiated sound can be judged as desirable, as is mostly the case for musical instruments, or undesirable, like noise generated by machinery. In architectural acoustics, one main goal is to limit the transmission of sound through walls. In the automobile industry, the engineers have to control the noise generated inside and outside the passenger compartment. This can be achieved by means of passive or active damping. In general, there is a strong need for quieter products and better sound quality generated by the structures in our daily environment.
Vibration-Torsion Dynamics of Ethane-like Molecules in Degenerate Vibrational States
NASA Astrophysics Data System (ADS)
Dilauro, C.; Lattanzi, F.; Avellino, R.
1994-10-01
We examine the mechanisms of interaction between the deformations of the two identical halves of ethane-like molecules in degenerate vibrational states, in order to determine the most convenient symmetries of the degenerate vibrational wavefunctions, in the molecular symmetry group G36( EM). In moderate barrier molecules this is related to the problem of the vibration-torsion symmetry labeling of the tunneling split components, in the order of increasing energies, in given rotational states. Numerical calculations have been performed as a guide to explore several different situations. It has been proved that the torsional splittings in degenerate vibrational states of molecules with a moderate barrier depend dramatically on the value of the γ-Coriolis coefficient and, unlike the splittings in nondegenerate vibrational states, are not determined only by the effective torsional potential function. Theory and numerical predictions support the experimental result that all normal modes of ethane behave as E1 d, E2 d pairs.
NASA Astrophysics Data System (ADS)
Tornabene, Francesco; Liverani, Alfredo; Caligiana, Gianni
2012-10-01
The Generalized Differential Quadrature (GDQ) method is applied to study the dynamic behavior of anisotropic doubly-curved shells and panels of revolution with a free-form meridian. The First-order Shear Deformation Theory (FSDT) is used to analyze the above mentioned moderately thick structural elements. In order to include the effect of the initial curvature in the evaluation of the stress resultants three different approaches, specifically Qatu approach, Toorani-Lakis approach and Reissner-Mindlin approach, are considered and compared. An improvement of the Classical Reissner-Mindlin Theory (CRMT) using a different kinematical model is considered. By so doing a generalization of the theory of anisotropic doubly-curved shells and panels of revolution is proposed. Four different anisotropic shell theories, namely General First-order Shear Deformation Theory by Qatu (GFSDTQ), General First-order Shear Deformation Theory by Toorani-Lakis (GFSDTTL), General First-order Shear Deformation Theory by Reissner-Mindlin (GFSDTRM) and Classical Reissner-Mindlin Theory (CRMT), are compared in order to show the differences and the accuracy of these theories. The solution is given in terms of generalized displacement components of points lying on the middle surface of the shell. Simple Rational Bézier curves are used to define the meridian curve of the revolution structures. Results are obtained taking the meridian and circumferential coordinates into account, without using the Fourier modal expansion methodology. Furthermore, GDQ results are compared with those obtained by using commercial programs such as Abaqus, Ansys, Nastran, Straus and Pro/Mechanica. Very good agreement is observed.
Noncontact Electromagnetic Vibration Source
NASA Technical Reports Server (NTRS)
Namkung, Min; Fulton, James P.; Wincheski, Buzz A.
1994-01-01
Metal aircraft skins scanned rapidly in vibration tests. Relatively simple combination of permanent magnets and electromagnet serves as noncontact vibration source for nondestructive testing of metal aircraft skins. In test, source excites vibrations, and vibration waveforms measured, then analyzed for changes in resonances signifying cracks and other flaws.
Transverse and longitudinal vibrations in amorphous silicon
NASA Astrophysics Data System (ADS)
Beltukov, Y. M.; Fusco, C.; Tanguy, A.; Parshin, D. A.
2015-12-01
We show that harmonic vibrations in amorphous silicon can be decomposed to transverse and longitudinal components in all frequency range even in the absence of the well defined wave vector q. For this purpose we define the transverse component of the eigenvector with given ω as a component, which does not change the volumes of Voronoi cells around atoms. The longitudinal component is the remaining orthogonal component. We have found the longitudinal and transverse components of the vibrational density of states for numerical model of amorphous silicon. The vibrations are mostly transverse below 7 THz and above 15 THz. In the frequency interval in between the vibrations have a longitudinal nature. Just this sudden transformation of vibrations at 7 THz from almost transverse to almost longitudinal ones explains the prominent peak in the diffusivity of the amorphous silicon just above 7 THz.