Volume 39, Issue 9, Pages 763-864(30 July 2002)
Research Article
Embedded turbulence model in numerical methods for hyperbolic conservation laws
NASA Astrophysics Data System (ADS)
Drikakis, D.
2002-07-01
The paper describes the use of numerical methods for hyperbolic conservation laws as an embedded turbulence modelling approach. Different Godunov-type schemes are utilized in computations of Burgers' turbulence and a two-dimensional mixing layer. The schemes include a total variation diminishing, characteristic-based scheme which is developed in this paper using the flux limiter approach. The embedded turbulence modelling property of the above methods is demonstrated through coarsely resolved large eddy simulations with and without subgrid scale models. Copyright
NASA Astrophysics Data System (ADS)
Wang, Gaili; Wong, Wai-Kin; Hong, Yang; Liu, Liping; Dong, Jili; Xue, Ming
2015-03-01
The primary objective of this study is to improve the performance of deterministic high resolution rainfall forecasts caused by severe storms by merging an extrapolation radar-based scheme with a storm-scale Numerical Weather Prediction (NWP) model. Effectiveness of Multi-scale Tracking and Forecasting Radar Echoes (MTaRE) model was compared with that of a storm-scale NWP model named Advanced Regional Prediction System (ARPS) for forecasting a violent tornado event that developed over parts of western and much of central Oklahoma on May 24, 2011. Then the bias corrections were performed to improve the forecast accuracy of ARPS forecasts. Finally, the corrected ARPS forecast and radar-based extrapolation were optimally merged by using a hyperbolic tangent weight scheme. The comparison of forecast skill between MTaRE and ARPS in high spatial resolution of 0.01° × 0.01° and high temporal resolution of 5 min showed that MTaRE outperformed ARPS in terms of index of agreement and mean absolute error (MAE). MTaRE had a better Critical Success Index (CSI) for less than 20-min lead times and was comparable to ARPS for 20- to 50-min lead times, while ARPS had a better CSI for more than 50-min lead times. Bias correction significantly improved ARPS forecasts in terms of MAE and index of agreement, although the CSI of corrected ARPS forecasts was similar to that of the uncorrected ARPS forecasts. Moreover, optimally merging results using hyperbolic tangent weight scheme further improved the forecast accuracy and became more stable.
Critical load: a novel approach to determining a sustainable intensity during resistance exercise.
Arakelian, Vivian M; Mendes, Renata G; Trimer, Renata; Rossi Caruso, Flavia C; de Sousa, Nuno M; Borges, Vanessa C; do Valle Gomes Gatto, Camila; Baldissera, Vilmar; Arena, Ross; Borghi-Silva, Audrey
2017-05-01
A hyperbolic function as well as a linear relationship between power output and time to exhaustion (Tlim) has been consistently observed during dynamic non-resistive exercises. However, little is known about its concept to resistance exercises (RE), which could be defined as critical load (CL). This study aimed to verify the existence of CL during dynamic RE and to verify the number of workbouts necessary to determine the optimal modeling to achieve it. Fifteen healthy men (23±2.5 yrs) completed 1 repetition maximum test (1RM) on a leg press and 3 (60%, 75% and 90% of 1RM) or 4 (+ 30% of 1RM) workbouts protocols to obtain the CL by hyperbolic and linear regression models between Tlim and load performed. Blood lactate and leg fatigue were also measured. CL was obtained during RE and 3 workbouts protocol estimate it at 53% while 4 tests at 38% of 1 RM. However, based on coefficients of determination, 3 protocols provided a better fit than the 4-parameter model, respectively (R2>0.95 vs. >0.77). Moreover, all intensities increased blood lactate and leg fatigue, however, when corrected by Tlim, were significantly lower at CL. It was possible to determinate CL during dynamic lower limbs RE and that 3 exhaustive workbouts can be used to better estimate the CL, constituting a new concept of determining this threshold during dynamic RE and reducing the physically demanding nature of the protocol. These findings may have important applications for functional performance evaluation and prescription of RE programs.
López Molina, Juan A; Rivera, María J; Trujillo, Macarena; Berjano, Enrique J
2009-04-01
The objectives of this study were to model the temperature progress of a pulsed radiofrequency (RF) power during RF heating of biological tissue, and to employ the hyperbolic heat transfer equation (HHTE), which takes the thermal wave behavior into account, and compare the results to those obtained using the heat transfer equation based on Fourier theory (FHTE). A theoretical model was built based on an active spherical electrode completely embedded in the biological tissue, after which HHTE and FHTE were analytically solved. We found three typical waveforms for the temperature progress depending on the relations between the dimensionless duration of the RF pulse delta(a) and the expression square root of lambda(rho-1), with lambda as the dimensionless thermal relaxation time of the tissue and rho as the dimensionless position. In the case of a unique RF pulse, the temperature at any location was the result of the overlapping of two different heat sources delayed for a duration delta(a) (each heat source being produced by a RF pulse of limitless duration). The most remarkable feature in the HHTE analytical solution was the presence of temperature peaks traveling through the medium at a finite speed. These peaks not only occurred during the RF power switch-on period but also during switch off. Finally, a physical explanation for these temperature peaks is proposed based on the interaction of forward and reverse thermal waves. All-purpose analytical solutions for FHTE and HHTE were obtained during pulsed RF heating of biological tissues, which could be used for any value of pulsing frequency and duty cycle.
An instability of hyperbolic space under the Yang-Mills flow
DOE Office of Scientific and Technical Information (OSTI.GOV)
Gegenberg, Jack; Day, Andrew C.; Liu, Haitao
2014-04-15
We consider the Yang-Mills flow on hyperbolic 3-space. The gauge connection is constructed from the frame-field and (not necessarily compatible) spin connection components. The fixed points of this flow include zero Yang-Mills curvature configurations, for which the spin connection has zero torsion and the associated Riemannian geometry is one of constant curvature. We analytically solve the linearized flow equations for a large class of perturbations to the fixed point corresponding to hyperbolic 3-space. These can be expressed as a linear superposition of distinct modes, some of which are exponentially growing along the flow. The growing modes imply the divergence ofmore » the (gauge invariant) perturbative torsion for a wide class of initial data, indicating an instability of the background geometry that we confirm with numeric simulations in the partially compactified case. There are stable modes with zero torsion, but all the unstable modes are torsion-full. This leads us to speculate that the instability is induced by the torsion degrees of freedom present in the Yang-Mills flow.« less
Hyperbolic umbilic caustics from oblate water drops with tilted illumination: Observations
NASA Astrophysics Data System (ADS)
Jobe, Oli; Thiessen, David B.; Marston, Philip L.
2017-11-01
Various groups have reported observations of hyperbolic umbilic diffraction catastrophe patterns in the far-field scattering by oblate acoustically levitated drops with symmetric illumination. In observations of that type the drop's symmetry axis is vertical and the illuminating light beam (typically an expanded laser beam) travels horizontally. In the research summarized here, scattering patterns in the primary rainbow region and drop measurements were recorded with vertically tilted laser beam illumination having a grazing angle as large as 4 degrees. The findings from these observations may be summarized as follows: (a) It remains possible to adjust the drop aspect ratio (diameter/height) = D/H so as to produce a V-shaped hyperbolic umbilic focal section (HUFS) in the far-field scattering. (b) The shift in the required D/H was typically an increase of less than 1% and was quadratic in the tilt. (c) The apex of the V-shaped HUFS was shifted vertically by an amount proportional to the tilt with a coefficient close to unity. The levitated drops had negligible up-down asymmetry. Our method of investigation should be useful for other generalized rainbows with tilted illumination.
Hyperbolic geometrical optics: Hyperbolic glass
NASA Astrophysics Data System (ADS)
De Micheli, Enrico; Scorza, Irene; Viano, Giovanni Alberto
2006-02-01
We study the geometrical optics generated by a refractive index of the form n (x,y)=1/y (y>0), where y is the coordinate of the vertical axis in an orthogonal reference frame in R2. We thus obtain what we call "hyperbolic geometrical optics" since the ray trajectories are geodesics in the Poincaré-Lobachevsky half-plane H2. Then we prove that the constant phase surface are horocycles and obtain the horocyclic waves, which are closely related to the classical Poisson kernel and are the analogs of the Euclidean plane waves. By studying the transport equation in the Beltrami pseudosphere, we prove (i) the conservation of the flow in the entire strip 0
NASA Technical Reports Server (NTRS)
Murad, P. A.
1993-01-01
Tsien's method is extended to treat the orbital motion of a body undergoing accelerations and decelerations. A generalized solution is discussed for the generalized case where a body undergoes azimuthal and radial thrust and the problem is further simplified for azimuthal thrust alone. Judicious selection of thrust could generate either an elliptic or hyperbolic trajectory. This is unexpected especially when the body has only enough energy for a lower state trajectory. The methodology is extended treating the problem of vehicle thrust for orbiting a sphere and vehicle thrust within the classical restricted three-body problem. Results for the latter situation can produce hyperbolic trajectories through eigen value decomposition. Since eigen values for no-thrust can be imaginary, thrust can generate real eigen values to describe hyperbolic trajectories. Keplerian dynamics appears to represent but a small subset of a much larger non-Keplerian domain especially when thrust effects are considered. The need for high thrust long duration space-based propulsion systems for changing a trajectory's canonical form is clearly demonstrated.
Ion beam figuring of Φ520mm convex hyperbolic secondary mirror
NASA Astrophysics Data System (ADS)
Meng, Xiaohui; Wang, Yonggang; Li, Ang; Li, Wenqing
2016-10-01
The convex hyperbolic secondary mirror is a Φ520-mm Zerodur lightweight hyperbolic convex mirror. Typically conventional methods like CCOS, stressed-lap polishing are used to manufacture this secondary mirror. Nevertheless, the required surface accuracy cannot be achieved through the use of conventional polishing methods because of the unpredictable behavior of the polishing tools, which leads to an unstable removal rate. Ion beam figuring is an optical fabrication method that provides highly controlled error of previously polished surfaces using a directed, inert and neutralized ion beam to physically sputter material from the optic surface. Several iterations with different ion beam size are selected and optimized to fit different stages of surface figure error and spatial frequency components. Before ion beam figuring, surface figure error of the secondary mirror is 2.5λ p-v, 0.23λ rms, and is improved to 0.12λ p-v, 0.014λ rms in several process iterations. The demonstration clearly shows that ion beam figuring can not only be used to the final correction of aspheric, but also be suitable for polishing the coarse surface of large, complex mirror.
A model for the plastic flow of landslides
Savage, William Z.; Smith, William K.
1986-01-01
To further the understanding of the mechanics of landslide flow, we present a model that predicts many of the observed attributes of landslides. The model is based on an integration of the hyperbolic differential equations for stress and velocity fields in a two-dimensional, inclined, semi-infinite half-space of Coulomb plastic material under elevated pore pressure and gravity. Our landslide model predicts commonly observed features. For example, compressive (passive), plug, or extending (active) flow will occur under appropriate longitudinal strain rates. Also, the model predicts that longitudinal stresses increase elliptically with depth to the basal slide plane, and that stress and velocity characteristics, surfaces along which discontinuities in stress and velocity are propagated, are coincident. Finally, the model shows how thrust and normal faults develop at the landslide surface in compressive and extending flow.
Dynamic Analysis of the Melanoma Model: From Cancer Persistence to Its Eradication
NASA Astrophysics Data System (ADS)
Starkov, Konstantin E.; Jimenez Beristain, Laura
In this paper, we study the global dynamics of the five-dimensional melanoma model developed by Kronik et al. This model describes interactions of tumor cells with cytotoxic T cells and respective cytokines under cellular immunotherapy. We get the ultimate upper and lower bounds for variables of this model, provide formulas for equilibrium points and present local asymptotic stability/hyperbolic instability conditions. Next, we prove the existence of the attracting set. Based on these results we come to global asymptotic melanoma eradication conditions via global stability analysis. Finally, we provide bounds for a locus of the melanoma persistence equilibrium point, study the case of melanoma persistence and describe conditions under which we observe global attractivity to the unique melanoma persistence equilibrium point.
NASA Technical Reports Server (NTRS)
Abarbanel, Saul; Gottlieb, David; Carpenter, Mark H.
1994-01-01
It has been previously shown that the temporal integration of hyperbolic partial differential equations (PDE's) may, because of boundary conditions, lead to deterioration of accuracy of the solution. A procedure for removal of this error in the linear case has been established previously. In the present paper we consider hyperbolic (PDE's) (linear and non-linear) whose boundary treatment is done via the SAT-procedure. A methodology is present for recovery of the full order of accuracy, and has been applied to the case of a 4th order explicit finite difference scheme.
Selective radiative heating of nanostructures using hyperbolic metamaterials
Ding, Ding; Minnich, Austin J
2015-01-01
Hyperbolic metamaterials (HMM) are of great interest due to their ability to break the diffraction limit for imaging and enhance near-field radiative heat transfer. Here we demonstrate that an annular, transparent HMM enables selective heating of a sub-wavelength plasmonic nanowire by controlling the angular mode number of a plasmonic resonance. A nanowire emitter, surrounded by an HMM, appears dark to incoming radiation from an adjacent nanowire emitter unless the second emitter is surrounded by an identical lens such that the wavelength and angular mode of the plasmonic resonance match. Our result can find applications in radiative thermal management.
A pseudospectral Legendre method for hyperbolic equations with an improved stability condition
NASA Technical Reports Server (NTRS)
Tal-Ezer, Hillel
1986-01-01
A new pseudospectral method is introduced for solving hyperbolic partial differential equations. This method uses different grid points than previously used pseudospectral methods: in fact the grid points are related to the zeroes of the Legendre polynomials. The main advantage of this method is that the allowable time step is proportional to the inverse of the number of grid points 1/N rather than to 1/n(2) (as in the case of other pseudospectral methods applied to mixed initial boundary value problems). A highly accurate time discretization suitable for these spectral methods is discussed.
A pseudospectral Legendre method for hyperbolic equations with an improved stability condition
NASA Technical Reports Server (NTRS)
Tal-Ezer, H.
1984-01-01
A new pseudospectral method is introduced for solving hyperbolic partial differential equations. This method uses different grid points than previously used pseudospectral methods: in fact the grid are related to the zeroes of the Legendre polynomials. The main advantage of this method is that the allowable time step is proportional to the inverse of the number of grid points 1/N rather than to 1/n(2) (as in the case of other pseudospectral methods applied to mixed initial boundary value problems). A highly accurate time discretization suitable for these spectral methods is discussed.
Error analysis of finite difference schemes applied to hyperbolic initial boundary value problems
NASA Technical Reports Server (NTRS)
Skollermo, G.
1979-01-01
Finite difference methods for the numerical solution of mixed initial boundary value problems for hyperbolic equations are studied. The reported investigation has the objective to develop a technique for the total error analysis of a finite difference scheme, taking into account initial approximations, boundary conditions, and interior approximation. Attention is given to the Cauchy problem and the initial approximation, the homogeneous problem in an infinite strip with inhomogeneous boundary data, the reflection of errors in the boundaries, and two different boundary approximations for the leapfrog scheme with a fourth order accurate difference operator in space.
Stability of hyperbolic-parabolic mixed type equations with partial boundary condition
NASA Astrophysics Data System (ADS)
Zhan, Huashui; Feng, Zhaosheng
2018-06-01
In this paper, we are concerned with the hyperbolic-parabolic mixed type equations with the non-homogeneous boundary condition. If it is degenerate on the boundary, the part of the boundary whose boundary value should be imposed, is determined by the entropy condition from the convection term. If there is no convection term in the equation, we show that the stability of solutions can be proved without any boundary condition. If the equation is completely degenerate, we show that the stability of solutions can be established just based on the partial boundary condition.
NASA Astrophysics Data System (ADS)
Zanraea, D. D. L.; Needham, D. J.
The depth-averaged hydraulic equations augmented with a suitable bed-load sediment transport function form a closed system which governs the one-dimensional flow in an alluvial river or channel. In this paper, it is shown that this system is hyperbolic and yields three families of shock-wave solutions. These are determined to be temporally stable in restricted regions of the (H, F0)-plane, via the Lax shock inequalities. Further, it is demonstrated that this criterion is equivalent to the energy dissipation criterion developed by Needham and Hey (1991).
Divergent conservation laws in hyperbolic thermoelasticity
NASA Astrophysics Data System (ADS)
Murashkin, E. V.; Radayev, Y. N.
2018-05-01
The present study is devoted to the problem of formulation of conservation laws in divergent form for hyperbolic thermoelastic continua. The field formalism is applied to study the problem. A natural density of thermoelastic action and the corresponding variational least action principle are formulated. A special form of the first variation of the action is employed to obtain 4-covariant divergent conservation laws. Differential field equations and constitutive laws are derived from a special form of the first variation of the action integral. The objectivity of constitutive equations is provided by the rotationally invariant forms of the Lagrangian employed.
NASA Technical Reports Server (NTRS)
Carpenter, Mark H.; Gottlieb, David; Abarbanel, Saul
1993-01-01
We present a systematic method for constructing boundary conditions (numerical and physical) of the required accuracy, for compact (Pade-like) high-order finite-difference schemes for hyperbolic systems. First, a roper summation-by-parts formula is found for the approximate derivative. A 'simultaneous approximation term' (SAT) is then introduced to treat the boundary conditions. This procedure leads to time-stable schemes even in the system case. An explicit construction of the fourth-order compact case is given. Numerical studies are presented to verify the efficacy of the approach.
