62 Matching Results

Search Results

Advanced search parameters have been applied.

Asymptotic analysis of spatial discretizations in implicit Monte Carlo

Description: We perform an asymptotic analysis of spatial discretizations in Implicit Monte Carlo (IMC). We consider two asymptotic scalings: one that represents a time step that resolves the mean-free time, and one that corresponds to a fixed, optically large time step. We show that only the latter scaling results in a valid spatial discretization of the proper diffusion equation, and thus we conclude that IMC only yields accurate solutions when using optically large spatial cells if time steps are also optically large. We demonstrate the validity of our analysis with a set of numerical examples.
Date: January 1, 2009
Creator: Densmore, Jeffery D
Partner: UNT Libraries Government Documents Department

Modeling broadband poroelastic propagation using an asymptotic approach

Description: An asymptotic method, valid in the presence of smoothly-varying heterogeneity, is used to derive a semi-analytic solution to the equations for fluid and solid displacements in a poroelastic medium. The solution is defined along trajectories through the porous medium model, in the manner of ray theory. The lowest order expression in the asymptotic expansion provides an eikonal equation for the phase. There are three modes of propagation, two modes of longitudinal displacement and a single mode of transverse displacement. The two longitudinal modes define the Biot fast and slow waves which have very different propagation characteristics. In the limit of low frequency, the Biot slow wave propagates as a diffusive disturbance, in essence a transient pressure pulse. Conversely, at low frequencies the Biot fast wave and the transverse mode are modified elastic waves. At intermediate frequencies the wave characteristics of the longitudinal modes are mixed. A comparison of the asymptotic solution with analytic and numerical solutions shows reasonably good agreement for both homogeneous and heterogeneous Earth models.
Date: May 1, 2009
Creator: Vasco, Donald W.
Partner: UNT Libraries Government Documents Department

Zeroth-order inversion of transient head observations

Description: A high-frequency, asymptotic solution for transient head,appropriate for a medium containing smoothly varying heterogeneity,provides a basis for efficient inverse modeling. The semi analyticsolution is trajectory based, akin to ray methods used in modeling wavepropagation, and may be constructed by post processing the output of anumerical simulator. For high frequencies, the amplitude sensitivities,the relationship between changes in flow properties and changes in headampliude, are dominated by the phase term which may be computed directlyfrom the output of the simulator. Thus, transient head waveforms may beinverted with little more computation than is required to invert arrivaltimes. An applicatino to synthetic head values indicates that thetechnique can be used to improve the fit to waveforms. An application totransient head data from the Migration experiment in Switzerland revealsa narrow, high conductivity pathway within a 0.5 m thick zone offracturing.
Date: August 15, 2007
Creator: Vasco, D.W.
Partner: UNT Libraries Government Documents Department

A numerical study of bunched beam transverse e-p instability based on the centroid model

Description: In a recent theoretical study of the transverse electron-proton (e-p) instability, an asymptotic solution has been found for the equations describing the centroid motion of the traversing proton bunch and the stationary background electrons. It was shown that the combination of finite proton bunch length, non-uniform proton line density, and the single-pass e-p interaction cause the instability to evolve intricately in space and time even in the linear regime. This paper reports a numerical study of the e-p instability based on the same centroid equations. The purpose of the work is to compare the numerical solution with the analytic solution and to use the numerical approach to investigate the early development of the instability not covered by the asymptotic solution. In particular, the instability threshold and the initial growth of the instability are studied for various proton-beam conditions, fraction of charge neutralization, and initial perturbations.
Date: January 1, 2003
Creator: Wang, T. F. (Tai-Sen F.)
Partner: UNT Libraries Government Documents Department

Dynamics of the Ginzburg-Landau equations of superconductivity

Description: This article is concerned with the dynamical properties of solutions of the time-dependent Ginzburg-Landau (TDGL) equations of superconductivity. It is shown that the TDGL equations define a dynamical process when the applied magnetic field varies with time, and a dynamical system when the applied magnetic field is stationary. The dynamical system describes the large-time asymptotic behavior: Every solution of the TDGL equations is attracted to a set of stationary solutions, which are divergence free. These results are obtained in the {open_quotes}{phi} = -{omega}({gradient}{center_dot}A){close_quotes} gauge, which reduces to the standard {close_quotes}{phi} = -{gradient}{center_dot}A{close_quotes} gauge if {omega} = 1 and to the zero-electric potential gauge if {omega} = 0; the treatment captures both in a unified framework. This gauge forces the London gauge, {gradient}{center_dot}A = 0, for any stationary solution of the TDGL equations.
Date: August 1997
Creator: Fleckinger-Pelle, J.; Kaper, H. G. & Takac, P.
Partner: UNT Libraries Government Documents Department

Asymptotic persistence of collective modes in shear flows

Description: A new nonasymptotic method is presented that reveals an unexpected richness in the spectrum of fluctuations sustained by a shear flow with nontrivial arbitrary mean kinematics. The vigor of the method is illustrated by analyzing a two-dimensional, compressible hydrodynamic shear flow. The temporal evolution of perturbations spans a wide range of nonexponential behavior from growth-cum oscillations to monotonic growth. The principal characteristic of the revealed exotic collective modes in their asymptotic persistence. {open_quotes}Echoing{close_quotes} as well as unstable (including parametrically-driven) solutions are displayed. Further areas of application, for both the method and the new physics, are outlined.
Date: March 31, 1998
Creator: Mahajan, S.M. & Rogava, A.D.
Partner: UNT Libraries Government Documents Department

Detonation front theories: Using high-resolution DNS to define extended asymptotic scalings and models

Description: When the detonation reaction-zone length, {eta}{sub r}, is short in comparison to the dimensions of the explosive piece being burnt, the detonation can be viewed as a propagating surface (or front) separating burnt from unburnt material. If the product of the shock curvature, {kappa} and {eta}{sub r} is small (i.e., the scaled shock curvature satisfies the {vert_bar}{kappa}{eta}{sub r}{vert_bar} {much_lt} 1), then to leading order the speed of this surface, D{sub n}({kappa}) is a function only of {kappa}. It is in this limit that the original version of the asymptotic detonation front theory, called detonation shock dynamics (DSD), derives the propagation law, D{sub n}({kappa}). In this lecture, the authors compare D{sub n}({kappa})-theory with the results obtained with high-resolution direct numerical simulations (DNS), and then use the DNS results to guide the development of extended asymptotic front theories with enhanced predictive capabilities.
Date: February 1, 1998
Creator: Aslam, T.D. & Bdzil, J.B.
Partner: UNT Libraries Government Documents Department

Sinuous oscillations and steady warps of polytropic disks

Description: In an asymptotic development of the equations governing the equilibria and linear stability of rapidly rotating polytropes we employed the slender aspect of these objects to reduce the three-dimensional partial differential equations to a somewhat simpler, ordinary integro-differential form. The earlier calculations dealt with isolated objects that were in centrifugal balance, that is the centrifugal acceleration of the configuration was balanced largely by self gravity with small contributions from the pressure gradient. Another interesting situation is that in which the polytrope rotates subject to externally imposed gravitational fields. In astrophysics, this is common in the theory of galactic dynamics because disks are unlikely to be isolated objects. The dark halos associated with disks also provide one possible explanation of the apparent warping of many galaxies. If the axis of the highly flattened disk is not aligned with that of the much less flattened halo, then the resultant torque of the halo gravity on the disk might provide a nonaxisymmetric distortion or disk warp. Motivated by these possibilities we shall here build models of polytropic disks of small but finite thickness which are subjected to prescribed, external gravitational fields. First we estimate how a symmetrical potential distorts the structure of the disk, then we examine its sinuous oscillations to confirm that they freely decay, hence suggesting that a warp must be externally forced. Finally, we consider steady warps of the disk plane when the axis of the disk does not coincide with that of the halo.
Date: May 1, 1995
Creator: Balmforth, N.J. & Spiegel, E.A.
Partner: UNT Libraries Government Documents Department

Application of the Two-Parameter J-A2 Description to Ductile Crack Growth

Description: Typical ASTM fracture testing determines J-integral resistance (J-R) curve or fracture toughness (JIC) based on specimens with high constraint geometry such as those specified in ASTM E1737-96. A three-term asymptotic solution with two parameters J and A2 (a constraint parameter) has been developed for characterizing the constraint effect of various geometries. The present paper extends the J-A2 characterization of a stationary crack tip to the regime of stable crack growth. Similar to the concept of J-controlled crack growth, the J-A2 description can be approximately used to characterize ductile crack growth under certain amount of crack extension. The region of J-A2 controlled crack growth is much larger than that controlled by J-integral alone. From the relationships between A2 and the test data, JIC and tearing modulus (TR), the coefficients used to define a J-R curve can be determined. For non-standard specimens or actual structures, once the constraint parameter A2 is determined, the J-R curves appropriate for these geometries can then be obtained. A procedure of transferring J-R curves determined from the standard ASTM procedure to non-standard specimens or flawed structures is outlined in the paper.
Date: July 8, 1999
Creator: Lam, P.S.
Partner: UNT Libraries Government Documents Department

Scattering of color dipoles: from low to high energies

Description: A dipole-dipole scattering amplitude is calculated exactly in the first two orders of perturbation theory. This amplitude is an analytic function of the relative energy and the dipoles' sizes. The cross section of the dipole-dipole scattering approached the high-energy BFKL asymptotics starting from relatively large rapidity {approx}4.
Date: November 1, 2002
Creator: Babansky, Alexander & Balitsky, Ian
Partner: UNT Libraries Government Documents Department

Axymptotic normality of X/sup 2/ in mxn tables with n large and small cell expectations

Description: Asymptotic normality for chi/sup 2/ used as a test for homogeneity is established under nonstandard conditions. The case of an mxn table with m fixed and the total number of observations proportional to n is studied for n large. Results are obtained under very mild assumptions on the marginal totals.
Date: January 1, 1977
Creator: Cuzick, J.
Partner: UNT Libraries Government Documents Department

Jet fragmentation and MLLA

Description: Recent CDF results in inclusive momentum distributions and multiplicities of particles in restricted cones around jets are compared to predictions using the Modified Leading Log Approximation. The authors found that MLLA gives a very reasonable description of jet fragmentation for a wide range of energies. Model parameters are extracted separately from the multiplicity and from the shape of the momentum distributions and are found to agree. The ratio of charged particle multiplicities in the gluon and quark jets measured in the context of MLLA is compared to the model-independent result and also found to agree.
Date: July 18, 2000
Creator: Safonov, Alexei N.
Partner: UNT Libraries Government Documents Department

Bounds for approximation in total variation distance by quantum circuits

Description: It was recently shown that for reasonable notions of approximation of states and functions by quantum circuits, almost all states and,functions are exponentially hard to approximate. The bounds obtained are asymptotically tight except for the one based on total variation distance (TVD). TVD is the most relevant metric for the performance of a quantum circuit. In this paper we obtain asymptotically tight bounds for TVD. We show that in a natural sense, almost all states are hard to approximate to within a TVD of 2/e -- {epsilon} even for exponentially small {epsilon}. The quantity 2/e -- {epsilon} is asymptotically the average distance to the uniform distribution. Almost all states with probability amplitudes concentrated in a small fraction of the space are hard to approximate to within a TVD of 2 -- {epsilon}. These results imply that non-uniform quantum circuit complexity is non-trivial in any reasonable model. They also reinforce the notion that the relative information distance between states (which is based on the difficulty of transforming one state to another) fully reflects the dimensionality of the space of qubits, not the number of qubits.
Date: September 1, 1995
Creator: Knill, E.
Partner: UNT Libraries Government Documents Department

Asymptotic distribution of a histogram density estimator

Description: Two theorems on the asymptotic distribution of a histogram density estimator based on randomly determined spacings introduced by Van Ryzin in 1973 are stated and proved. One theorem gives conditions for the pointwise asymptotic normality of the density estimator for points in the support of the density at which the density is continuously differentiable. A second theorem gives conditions for the pointwise asymptotic normality of the density estimator with a faster convergence rate for points in the support of the density at which the density is twice continuously differentiable. The results are used to compare the relative asymptotic efficiencies of the histogram estimator with the kernal method of density estimation.
Date: June 1, 1980
Creator: Kim, B K & Van Ryzin, J
Partner: UNT Libraries Government Documents Department

Investigation of an empirical probability measure based test for multivariate normality

Description: Foutz (1980) derived a goodness of fit test for a hypothesis specifying a continuous, p-variate distribution. The test statistic is both distribution-free and independent of p. In adapting the Foutz test for multivariate normality, we consider using chi/sup 2/ and rescaled beta variates in constructing statistically equivalent blocks. The Foutz test is compared to other multivariate normality tests developed by Hawkins (1981) and Malkovich and Afifi (1973). The set of alternative distributions tested include Pearson type II and type VII, Johnson translations, Plackett, and distributions arising from Khintchine's theorem. Univariate alternatives from the general class developed by Johnson et al. (1980) were also used. An empirical study confirms the independence of the test statistic on p even when parameters are estimated. In general, the Foutz test is less conservative under the null hypothesis but has poorer power under most alternatives than the other tests.
Date: January 1, 1984
Creator: Booker, J.M.; Johnson, M.E. & Beckman, R.J.
Partner: UNT Libraries Government Documents Department

CALCULATING POINT-CHARGE WAKEFIELDS FROM FINITE LENGTH BUNCH WAKE-POTENTIALS

Description: Starting from analytical properties of high frequency geometric impedance we show how one can accurately calculate short bunch wake-potentials (and even point-charge wakefields) from time domain calculations performed with a much longer bunch. In many practical instances this drastically reduces the need for computer resources, speeds up the calculations, and improves their accuracy. To illustrate this method we give examples for 2D accelerator structures of various complexities. We describe preliminary results of a new method that allows us to accurately obtain longitudinal wakefields of short bunches by adding a long-bunch result from an EM solver and a singular analytical wake model. In the future this work will be generalized to 3D geometries as well. Similarly, the method should be equally applicable to the calculations of transverse wakefields. Periodic structures with a significant number of periods (2 {ge} a{sup 2}/{sigma}L, where L is the period length) have not been considered so far. They have asymptotic wakefields that differ from the examples described above. We believe this method is applicable to such geometries as well, as long as correct asymptotic solutions are used.
Date: March 28, 2011
Creator: Podobedov, B. & Stupakov, G.
Partner: UNT Libraries Government Documents Department

An Order-of-Magnitude Estimation of Benzene Concentration in Saltstone Vault

Description: The contents of Tank 48H that include the tetraphenylborate (TPB) precipitates of potassium and cesium will be grouted and stored in the Saltstone vault. The grouting process is exothermic, which should accelerate the rate of decomposition of TPB precipitates eventually to benzene. Because the vault is not currently outfitted with an active ventilation system, there is a concern that a mixture of flammable gases may form in the vapor space of each cell filled with the curing grout. The purpose of this study was to determine if passive breathing induced by the diurnal fluctuations of barometric pressure would provide any mitigating measure against potential flammability in the cell vapor space. In Revision 0 of this document, a set of algorithms were presented that would predict the equilibrium concentration of benzene in the cell vapor space as a function of benzene generation rate, fill height, and passive breathing rate. The algorithms were derived based on several simplifying assumptions so that order of magnitude estimates could be made quickly for scoping purposes. In particular, it was assumed that passive breathing would occur solely due to barometric pressure fluctuations that were sinusoidal; the resulting algorithm for estimating the rate of passive breathing into or out of each cell is given in Eq. (10). Since Revision 0 was issued, the validity of this critical assumption on the mode of passive breathing was checked against available passive ventilation data for the Hanford waste tanks. It was found that the passive breathing rates estimated from Eq. (10) were on average 50 to 90% lower than those measured for 5 out of 6 Hanford tanks considered in this study (see Table 1); for Tank U-106, the estimated passive breathing rates were on average 20% lower than the measured data. These results indicate that Eq. (10) would most ...
Date: March 20, 2006
Creator: CHOI, A
Partner: UNT Libraries Government Documents Department

An Integrated Approach to Characterizing Bypassed Oil in Heterogeneous and Fractured Reservoirs Using Partitioning Tracers

Description: This report presents an efficient trajectory-based approach to integrate transient pressure data into high-resolution reservoir and aquifer models. The method involves alternating travel time and peak amplitude matching of pressure response using inverse modeling and is particularly well-suited for high resolution subsurface characterization using hydraulic tomography or pressure interference tests. Compared to travel time inversion only, our proposed approach results in a significantly improved match of the pressure response at the wells and also better estimates of subsurface properties. This is accomplished with very little increase in computational cost. Utilizing the concept of a ''diffusive'' time of flight derived from an asymptotic solution of the diffusivity equation, we develop analytical approaches to estimate the sensitivities for travel time and peak amplitude of pressure response to subsurface properties. The sensitivities are then used in an iterative least-squared minimization to match the pressure data. We illustrate our approach using synthetic and field examples. In the field application at a fractured limestone formation, the predominant fracture patterns emerging from the inversion are shown to be consistent with independent geophysical experiments and borehole data.
Date: August 1, 2006
Creator: Datta-Gupta, Akhil
Partner: UNT Libraries Government Documents Department

Modeling quasi-static poroelastic propagation using an asymptotic approach

Description: Since the formulation of poroelasticity (Biot(1941)) and its reformulation (Rice & Cleary(1976)), there have been many efforts to solve the coupled system of equations. Perhaps because of the complexity of the governing equations, most of the work has been directed towards finding numerical solutions. For example, Lewis and co-workers published early papers (Lewis & Schrefler(1978); Lewis et al.(1991)Lewis, Schrefler, & Simoni) concerned with finite-element methods for computing consolidation, subsidence, and examining the importance of coupling. Other early work dealt with flow in a deformable fractured medium (Narasimhan & Witherspoon 1976); Noorishad et al.(1984)Noorishad, Tsang, & Witherspoon. This effort eventually evolved into a general numerical approach for modeling fluid flow and deformation (Rutqvist et al.(2002)Rutqvist, Wu, Tsang, & Bodvarsson). As a result of this and other work, numerous coupled, computer-based algorithms have emerged, typically falling into one of three categories: one-way coupling, loose coupling, and full coupling (Minkoff et al.(2003)Minkoff, Stone, Bryant, Peszynska, & Wheeler). In one-way coupling the fluid flow is modeled using a conventional numerical simulator and the resulting change in fluid pressures simply drives the deformation. In loosely coupled modeling distinct geomechanical and fluid flow simulators are run for a sequence of time steps and at the conclusion of each step information is passed between the simulators. In full coupling, the fluid flow and geomechanics equations are solved simultaneously at each time step (Lewis & Sukirman(1993); Lewis & Ghafouri(1997); Gutierrez & Lewis(2002)). One disadvantage of a purely numerical approach to solving the governing equations of poroelasticity is that it is not clear how the various parameters interact and influence the solution. Analytic solutions have an advantage in that respect; the relationship between the medium and fluid properties is clear from the form of the solution. Unfortunately, analytic solutions are only available for highly idealized conditions, such as a ...
Date: November 1, 2007
Creator: Vasco, D.W.
Partner: UNT Libraries Government Documents Department

A variational principle for the Ackerberg-O'Malley resonance problem

Description: A new variational principle is proposed for determining the asymptotic expansion of the solution of the Ackerberg-O'Malley resonance problem to any order in epsilon. The method used differs from that originally proposed by Grasman and Matkowsky in that variations are allowed in the boundary values of a transformed problem, making the calculation applicable to any order in epsilon. Explicit results using the new approach are reported to O(epsilon) and confirmed with asymptotic expansions of the exact solution. The O(1) results agree with those reported in the literature, and the O(epsilon) results are new. It is not difficult to perform the higher order calculations using the new variational approach, in contrast to the current methods in use. In the case where the coefficient functions are analytic in the independent variable and epsilon, and the outer expansion satisfies a certain restriction, the exact solution presented provides an alternative and equally simple method for calculating the asymptotic expansion to any order in epsilon. 14 refs.
Date: August 1, 1987
Creator: Srinivasan, R.
Partner: UNT Libraries Government Documents Department

Numerical and asymptotic studies of complex flow dynamics. Progress report, January 1--December 31, 1997

Description: If one has computed a stationary state of an evolution equation, a major question is if the state is stable, and many different techniques have been developed to discuss this important question. Together with H.-O. Kreiss the author completed a survey article on nonlinear stability for time dependent PDEs. They give an account of the Lyapunov technique and the resolvent technique to study stability questions. It has been known for some time that the strength of the resolvent technique is the control of the small-wave-number projection of the solution, whereas Lyapunov`s technique is good for high wave numbers. In ref. 6 they succeeded in combining the two techniques. Together with T. Hagstrom, the author has studied somewhat related questions for flows at low Mach number. They show all time existence of classical solutions when the initial data are almost incompressible. The result is stronger than Hoff`s in that they can allow for a ball of slightly compressible data with a radius independent of the Mach number M. In contrast, the radius of Hoff`s ball shrinks to zero as M approaches 0. Together with H.-J. Schroll, RWTH Aachen, the author worked on conservation laws with stiff source term. From a mathematical point of view, these are singular perturbation problems for which the reduced problem is singular. In ref. 2 they introduce the notion of stiff well-posedness, which--in the linear constant coefficient case--characterizes all hyperbolic systems with stiff source term 1/{epsilon} Bu whose solutions converge for {epsilon} approaches 0. There are many applications. A rather new field, which has received considerable attention in the physics community, is self-organized criticality. Bak gives a popular account. The aim is to study the self-organization of systems towards a critical state and to understand the scaling properties at criticality. Together with the students the author has ...
Date: January 1, 1997
Creator: Lorenz, J.
Partner: UNT Libraries Government Documents Department

Heterogeneous domain decomposition for singularly perturbed elliptic boundary value problems

Description: A heterogeneous domain-decomposition method is presented for the numerical solution of singularly perturbed elliptic boundary value problem. The method, which is parallelizable at various levels, uses several ideas of asymptotic analysis. The sub-domains match the domains of validity of the local ({open_quotes}inner{close_quotes} and {open_quotes}outer{close_quotes}) asymptotic expansions, and cut-off functions are used to match solutions in neighboring subdomains. The positions of the interfaces, as well as the mesh widths, depend on the small parameter, {epsilon}. On the subdomains, iterative solution techniques are used, which may vary from one subdomain to another. The global convergence rate depends on {epsilon}; it generally increases like some power of (log({epsilon}{sup -1})){sup -1} as {epsilon} {down_arrow} 0. The method is illustrated on several two-dimensional singular perturbation problems.
Date: April 14, 1995
Creator: Garbey, M. & Kaper, H.G.
Partner: UNT Libraries Government Documents Department