Numerical and analytical simulation of the production process of ZrO2 hollow particles

NASA Astrophysics Data System (ADS)

Safaei, Hadi; Emami, Mohsen Davazdah

2017-12-01

In this paper, the production process of hollow particles from the agglomerated particles is addressed analytically and numerically. The important parameters affecting this process, in particular, the initial porosity level of particles and the plasma gun types are investigated. The analytical model adopts a combination of quasi-steady thermal equilibrium and mechanical balance. In the analytical model, the possibility of a solid core existing in agglomerated particles is examined. In this model, a range of particle diameters (50μm ≤ D_{p0} ≤ 160 μ m) and various initial porosities ( 0.2 ≤ p ≤ 0.7) are considered. The numerical model employs the VOF technique for two-phase compressible flows. The production process of hollow particles from the agglomerated particles is simulated, considering an initial diameter of D_{p0} = 60 μm and initial porosity of p = 0.3, p = 0.5, and p = 0.7. Simulation results of the analytical model indicate that the solid core diameter is independent of the initial porosity, whereas the thickness of the particle shell strongly depends on the initial porosity. In both models, a hollow particle may hardly develop at small initial porosity values ( p < 0.3), while the particle disintegrates at high initial porosity values ( p > 0.6.

Distribution of velocities and acceleration for a particle in Brownian correlated disorder: Inertial case

NASA Astrophysics Data System (ADS)

Le Doussal, Pierre; Petković, Aleksandra; Wiese, Kay Jörg

2012-06-01

We study the motion of an elastic object driven in a disordered environment in presence of both dissipation and inertia. We consider random forces with the statistics of random walks and reduce the problem to a single degree of freedom. It is the extension of the mean-field Alessandro-Beatrice- Bertotti-Montorsi (ABBM) model in presence of an inertial mass m. While the ABBM model can be solved exactly, its extension to inertia exhibits complicated history dependence due to oscillations and backward motion. The characteristic scales for avalanche motion are studied from numerics and qualitative arguments. To make analytical progress, we consider two variants which coincide with the original model whenever the particle moves only forward. Using a combination of analytical and numerical methods together with simulations, we characterize the distributions of instantaneous acceleration and velocity, and compare them in these three models. We show that for large driving velocity, all three models share the same large-deviation function for positive velocities, which is obtained analytically for small and large m, as well as for m=6/25. The effect of small additional thermal and quantum fluctuations can be treated within an approximate method.

Studies in nonlinear problems of energy. Final report

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

Matkowsky, B.J.

1998-12-01

The author completed a successful research program on Nonlinear Problems of Energy, with emphasis on combustion and flame propagation. A total of 183 papers associated with the grant has appeared in the literature, and the efforts have twice been recognized by DOE`s Basic Science Division for Top Accomplishment. In the research program the author concentrated on modeling, analysis and computation of combustion phenomena, with particular emphasis on the transition from laminar to turbulent combustion. Thus he investigated the nonlinear dynamics and pattern formation in the successive stages of transition. He described the stability of combustion waves, and transitions to wavesmore » exhibiting progressively higher degrees of spatio-temporal complexity. Combustion waves are characterized by large activation energies, so that chemical reactions are significant only in thin layers, termed reaction zones. In the limit of infinite activation energy, the zones shrink to moving surfaces, termed fronts, which must be found during the course of the analysis, so that the problems are moving free boundary problems. The analytical studies were carried out for the limiting case with fronts, while the numerical studies were carried out for the case of finite, though large, activation energy. Accurate resolution of the solution in the reaction zone(s) is essential, otherwise false predictions of dynamical behavior are possible. Since the reaction zones move, and their location is not known a-priori, the author has developed adaptive pseudo-spectral methods, which have proven to be very useful for the accurate, efficient computation of solutions of combustion, and other, problems. The approach is based on a combination of analytical and numerical methods. The numerical computations built on and extended the information obtained analytically. Furthermore, the solutions obtained analytically served as benchmarks for testing the accuracy of the solutions determined computationally. Finally, the computational results suggested new analysis to be considered. A cumulative list of publications citing the grant make up the contents of this report.« less

The development of a rail passenger coach car crush zone

DOT National Transportation Integrated Search

2002-03-21

There is currently an active program within the United States to improve the passive safety of train systems. In particular, the Volpe Center is supporting the Federal Railroad Administrations research that includes numerous analytical studies and...

Modeling and analysis to quantify MSE wall behavior and performance.

DOT National Transportation Integrated Search

2009-08-01

To better understand potential sources of adverse performance of mechanically stabilized earth (MSE) walls, a suite of analytical models was studied using the computer program FLAC, a numerical modeling computer program widely used in geotechnical en...

Ground-state magnetization of the Ising spin glass: A recursive numerical method and Chen-Ma scaling

NASA Astrophysics Data System (ADS)

Sepehrinia, Reza; Chalangari, Fartash

2018-03-01

The ground-state properties of quasi-one-dimensional (Q1D) Ising spin glass are investigated using an exact numerical approach and analytical arguments. A set of coupled recursive equations for the ground-state energy are introduced and solved numerically. For various types of coupling distribution, we obtain accurate results for magnetization, particularly in the presence of a weak external magnetic field. We show that in the weak magnetic field limit, similar to the 1D model, magnetization exhibits a singular power-law behavior with divergent susceptibility. Remarkably, the spectrum of magnetic exponents is markedly different from that of the 1D system even in the case of two coupled chains. The magnetic exponent makes a crossover from being dependent on a distribution function to a constant value independent of distribution. We provide an analytic theory for these observations by extending the Chen-Ma argument to the Q1D case. We derive an analytical formula for the exponent which is in perfect agreement with the numerical results.

A differential delay equation arising from the sieve of Eratosthenes

NASA Technical Reports Server (NTRS)

Cheer, A. Y.; Goldston, D. A.

1990-01-01

Consideration is given to the differential delay equation introduced by Buchstab (1937) in connection with an asymptotic formula for the uncanceled terms in the sieve of Eratosthenes. Maier (1985) used this result to show there is unexpected irreqularity in the distribution of primes in short intervals. The function omega(u) is studied in this paper using numerical and analytical techniques. The results are applied to give some numerical constants in Maier's theorem.

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

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

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

Symmetry-plane model of 3D Euler flows: Mapping to regular systems and numerical solutions of blowup

NASA Astrophysics Data System (ADS)

Mulungye, Rachel M.; Lucas, Dan; Bustamante, Miguel D.

2014-11-01

We introduce a family of 2D models describing the dynamics on the so-called symmetry plane of the full 3D Euler fluid equations. These models depend on a free real parameter and can be solved analytically. For selected representative values of the free parameter, we apply the method introduced in [M.D. Bustamante, Physica D: Nonlinear Phenom. 240, 1092 (2011)] to map the fluid equations bijectively to globally regular systems. By comparing the analytical solutions with the results of numerical simulations, we establish that the numerical simulations of the mapped regular systems are far more accurate than the numerical simulations of the original systems, at the same spatial resolution and CPU time. In particular, the numerical integrations of the mapped regular systems produce robust estimates for the growth exponent and singularity time of the main blowup quantity (vorticity stretching rate), converging well to the analytically-predicted values even beyond the time at which the flow becomes under-resolved (i.e. the reliability time). In contrast, direct numerical integrations of the original systems develop unstable oscillations near the reliability time. We discuss the reasons for this improvement in accuracy, and explain how to extend the analysis to the full 3D case. Supported under the programme for Research in Third Level Institutions (PRTLI) Cycle 5 and co-funded by the European Regional Development Fund.

Solar neutrino masses and mixing from bilinear R-parity broken supersymmetry: Analytical versus numerical results

NASA Astrophysics Data System (ADS)

Díaz, M.; Hirsch, M.; Porod, W.; Romão, J.; Valle, J.

2003-07-01

We give an analytical calculation of solar neutrino masses and mixing at one-loop order within bilinear R-parity breaking supersymmetry, and compare our results to the exact numerical calculation. Our method is based on a systematic perturbative expansion of R-parity violating vertices to leading order. We find in general quite good agreement between the approximate and full numerical calculations, but the approximate expressions are much simpler to implement. Our formalism works especially well for the case of the large mixing angle Mikheyev-Smirnov-Wolfenstein solution, now strongly favored by the recent KamLAND reactor neutrino data.

Fourth order Douglas implicit scheme for solving three dimension reaction diffusion equation with non-linear source term

NASA Astrophysics Data System (ADS)

Hasnain, Shahid; Saqib, Muhammad; Mashat, Daoud Suleiman

2017-07-01

This research paper represents a numerical approximation to non-linear three dimension reaction diffusion equation with non-linear source term from population genetics. Since various initial and boundary value problems exist in three dimension reaction diffusion phenomena, which are studied numerically by different numerical methods, here we use finite difference schemes (Alternating Direction Implicit and Fourth Order Douglas Implicit) to approximate the solution. Accuracy is studied in term of L2, L∞ and relative error norms by random selected grids along time levels for comparison with analytical results. The test example demonstrates the accuracy, efficiency and versatility of the proposed schemes. Numerical results showed that Fourth Order Douglas Implicit scheme is very efficient and reliable for solving 3-D non-linear reaction diffusion equation.

A Numerical Study of Coupled Non-Linear Equations of Thermo-Viscous Fluid Flow in Cylindrical Geometry

NASA Astrophysics Data System (ADS)

Pothanna, N.; Aparna, P.; Gorla, R. S. R.

2017-12-01

In this paper we present numerical solutions to coupled non-linear governing equations of thermo-viscous fluid flow in cylindrical geometry using MATHEMATICA software solver. The numerical results are presented in terms of velocity, temperature and pressure distribution for various values of the material parameters such as the thermo-mechanical stress coefficient, thermal conductivity coefficient, Reiner Rivlin cross viscosity coefficient and the Prandtl number in the form of tables and graphs. Also, the solutions to governing equations for slow steady motion of a fluid have been obtained numerically and compared with the existing analytical results and are found to be in excellent agreement. The results of the present study will hopefully enable a better understanding applications of the flow under consideration.

Analytical and numerical prediction of harmonic sound power in the inlet of aero-engines with emphasis on transonic rotation speeds

NASA Astrophysics Data System (ADS)

Lewy, Serge; Polacsek, Cyril; Barrier, Raphael

2014-12-01

Tone noise radiated through the inlet of a turbofan is mainly due to rotor-stator interactions at subsonic regimes (approach flight), and to the shock waves attached to each blade at supersonic helical tip speeds (takeoff). The axial compressor of a helicopter turboshaft engine is transonic as well and can be studied like turbofans at takeoff. The objective of the paper is to predict the sound power at the inlet radiating into the free field, with a focus on transonic conditions because sound levels are much higher. Direct numerical computation of tone acoustic power is based on a RANS (Reynolds averaged Navier-Stokes) solver followed by an integration of acoustic intensity over specified inlet cross-sections, derived from Cantrell and Hart equations (valid in irrotational flows). In transonic regimes, sound power decreases along the intake because of nonlinear propagation, which must be discriminated from numerical dissipation. This is one of the reasons why an analytical approach is also suggested. It is based on three steps: (i) appraisal of the initial pressure jump of the shock waves; (ii) 2D nonlinear propagation model of Morfey and Fisher; (iii) calculation of the sound power of the 3D ducted acoustic field. In this model, all the blades are assumed to be identical such that only the blade passing frequency and its harmonics are predicted (like in the present numerical simulations). However, transfer from blade passing frequency to multiple pure tones can be evaluated in a fourth step through a statistical analysis of irregularities between blades. Interest of the analytical method is to provide a good estimate of nonlinear acoustic propagation in the upstream duct while being easy and fast to compute. The various methods are applied to two turbofan models, respectively in approach (subsonic) and takeoff (transonic) conditions, and to a Turbomeca turboshaft engine (transonic case). The analytical method in transonic appears to be quite reliable by comparison with the numerical solution and with available experimental data.

Numerical evaluation of heating in the human head due to magnetic resonance imaging (MRI)

NASA Astrophysics Data System (ADS)

Nguyen, Uyen; Brown, Steve; Chang, Isaac; Krycia, Joe; Mirotznik, Mark S.

2003-06-01

In this paper we present a numerical model for evaluating tissue heating during magnetic resonance imaging (MRI). Our method, which included a detailed anatomical model of a human head, calculated both the electromagnetic power deposition and the associated temperature elevations during a MRI head examination. Numerical studies were conducted using a realistic birdcage coil excited at frequencies ranging from 63 MHz to 500 MHz. The model was validated both experimentally and analytically. The experimental validation was performed at the MR test facility located at the FDA's Center for Devices and Radiological Health (CDRH).

A Conceptual Analytics Model for an Outcome-Driven Quality Management Framework as Part of Professional Healthcare Education.

PubMed

Hervatis, Vasilis; Loe, Alan; Barman, Linda; O'Donoghue, John; Zary, Nabil

2015-10-06

Preparing the future health care professional workforce in a changing world is a significant undertaking. Educators and other decision makers look to evidence-based knowledge to improve quality of education. Analytics, the use of data to generate insights and support decisions, have been applied successfully across numerous application domains. Health care professional education is one area where great potential is yet to be realized. Previous research of Academic and Learning analytics has mainly focused on technical issues. The focus of this study relates to its practical implementation in the setting of health care education. The aim of this study is to create a conceptual model for a deeper understanding of the synthesizing process, and transforming data into information to support educators' decision making. A deductive case study approach was applied to develop the conceptual model. The analytics loop works both in theory and in practice. The conceptual model encompasses the underlying data, the quality indicators, and decision support for educators. The model illustrates how a theory can be applied to a traditional data-driven analytics approach, and alongside the context- or need-driven analytics approach.

A Conceptual Analytics Model for an Outcome-Driven Quality Management Framework as Part of Professional Healthcare Education

PubMed Central

Loe, Alan; Barman, Linda; O'Donoghue, John; Zary, Nabil

2015-01-01

Background Preparing the future health care professional workforce in a changing world is a significant undertaking. Educators and other decision makers look to evidence-based knowledge to improve quality of education. Analytics, the use of data to generate insights and support decisions, have been applied successfully across numerous application domains. Health care professional education is one area where great potential is yet to be realized. Previous research of Academic and Learning analytics has mainly focused on technical issues. The focus of this study relates to its practical implementation in the setting of health care education. Objective The aim of this study is to create a conceptual model for a deeper understanding of the synthesizing process, and transforming data into information to support educators’ decision making. Methods A deductive case study approach was applied to develop the conceptual model. Results The analytics loop works both in theory and in practice. The conceptual model encompasses the underlying data, the quality indicators, and decision support for educators. Conclusions The model illustrates how a theory can be applied to a traditional data-driven analytics approach, and alongside the context- or need-driven analytics approach. PMID:27731840

Chemical transport in a fissured rock: Verification of a numerical model

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

Rasmuson, A.; Narasimhan, T. N.; Neretnieks, I.

1982-10-01

Numerical models for simulating chemical transport in fissured rocks constitute powerful tools for evaluating the acceptability of geological nuclear waste repositories. Due to the very long-term, high toxicity of some nuclear waste products, the models are required to predict, in certain cases, the spatial and temporal distribution of chemical concentration less than 0.001% of the concentration released from the repository. Whether numerical models can provide such accuracies is a major question addressed in the present work. To this end, we have verified a numerical model, TRUMP, which solves the advective diffusion equation in general three dimensions with or without decaymore » and source terms. The method is based on an integrated finite-difference approach. The model was verified against known analytic solution of the one-dimensional advection-diffusion problem as well as the problem of advection-diffusion in a system of parallel fractures separated by spherical particles. The studies show that as long as the magnitude of advectance is equal to or less than that of conductance for the closed surface bounding any volume element in the region (that is, numerical Peclet number <2), the numerical method can indeed match the analytic solution within errors of ±10{sup -3} % or less. The realistic input parameters used in the sample calculations suggest that such a range of Peclet numbers is indeed likely to characterize deep groundwater systems in granitic and ancient argillaceous systems. Thus TRUMP in its present form does provide a viable tool for use in nuclear waste evaluation studies. A sensitivity analysis based on the analytic solution suggests that the errors in prediction introduced due to uncertainties in input parameters is likely to be larger than the computational inaccuracies introduced by the numerical model. Currently, a disadvantage in the TRUMP model is that the iterative method of solving the set of simultaneous equations is rather slow when time constants vary widely over the flow region. Although the iterative solution may be very desirable for large three-dimensional problems in order to minimize computer storage, it seems desirable to use a direct solver technique in conjunction with the mixed explicit-implicit approach whenever possible. work in this direction is in progress.« less

Testing density-dependent groundwater models: Two-dimensional steady state unstable convection in infinite, finite and inclined porous layers

USGS Publications Warehouse

Weatherill, D.; Simmons, C.T.; Voss, C.I.; Robinson, N.I.

2004-01-01

This study proposes the use of several problems of unstable steady state convection with variable fluid density in a porous layer of infinite horizontal extent as two-dimensional (2-D) test cases for density-dependent groundwater flow and solute transport simulators. Unlike existing density-dependent model benchmarks, these problems have well-defined stability criteria that are determined analytically. These analytical stability indicators can be compared with numerical model results to test the ability of a code to accurately simulate buoyancy driven flow and diffusion. The basic analytical solution is for a horizontally infinite fluid-filled porous layer in which fluid density decreases with depth. The proposed test problems include unstable convection in an infinite horizontal box, in a finite horizontal box, and in an infinite inclined box. A dimensionless Rayleigh number incorporating properties of the fluid and the porous media determines the stability of the layer in each case. Testing the ability of numerical codes to match both the critical Rayleigh number at which convection occurs and the wavelength of convection cells is an addition to the benchmark problems currently in use. The proposed test problems are modelled in 2-D using the SUTRA [SUTRA-A model for saturated-unsaturated variable-density ground-water flow with solute or energy transport. US Geological Survey Water-Resources Investigations Report, 02-4231, 2002. 250 p] density-dependent groundwater flow and solute transport code. For the case of an infinite horizontal box, SUTRA results show a distinct change from stable to unstable behaviour around the theoretical critical Rayleigh number of 4??2 and the simulated wavelength of unstable convection agrees with that predicted by the analytical solution. The effects of finite layer aspect ratio and inclination on stability indicators are also tested and numerical results are in excellent agreement with theoretical stability criteria and with numerical results previously reported in traditional fluid mechanics literature. ?? 2004 Elsevier Ltd. All rights reserved.

Regional scale impact of tidal forcing on groundwater flow in unconfined coastal aquifers

NASA Astrophysics Data System (ADS)

Pauw, P. S.; Oude Essink, G. H. P.; Leijnse, A.; Vandenbohede, A.; Groen, J.; van der Zee, S. E. A. T. M.

2014-09-01

This paper considers the impact of tidal forcing on regional groundwater flow in an unconfined coastal aquifer. Numerical models are used to quantify this impact for a wide range of hydrogeological conditions. Both a shallow and a deep aquifer are investigated with regard to three dimensionless parameter groups that determine the groundwater flow to a large extent. Analytical expressions are presented that allow for a quick estimate of the regional scale effect of tidal forcing under the same conditions as used in the numerical models. Quantitatively, the results in this paper are complementary to previous studies by taking into account variable density groundwater flow, dispersive salt transport and a seepage face in the intertidal area. Qualitatively, the results are in line with previous investigations. The time-averaged hydraulic head at the high tide mark increases upon a decrease of each of the three considered dimensionless parameter groups: R (including the ratio of the hydraulic conductivity and the precipitation excess), α (the slope of the intertidal area) and AL (the ratio of the width of the fresh water lens and the tidal amplitude). The relative change of the location and the hydraulic head of the groundwater divide, which together characterize regional groundwater flow, increase as α and AL decrease, but decrease as R decreases. The difference between the analytical solutions and numerical results is small. Therefore, the presented analytical solutions can be used to estimate the bias that is introduced in a numerical model if tidal forcing is neglected. The results should be used with caution in case of significant wave forcing, as this was not considered.

Modeling of the Global Water Cycle - Analytical Models

Treesearch

Yongqiang Liu; Roni Avissar

2005-01-01

Both numerical and analytical models of coupled atmosphere and its underlying ground components (land, ocean, ice) are useful tools for modeling the global and regional water cycle. Unlike complex three-dimensional climate models, which need very large computing resources and involve a large number of complicated interactions often difficult to interpret, analytical...

The contribution of Raman spectroscopy to the analytical quality control of cytotoxic drugs in a hospital environment: eliminating the exposure risks for staff members and their work environment.

PubMed

Bourget, Philippe; Amin, Alexandre; Vidal, Fabrice; Merlette, Christophe; Troude, Pénélope; Baillet-Guffroy, Arlette

2014-08-15

The purpose of the study was to perform a comparative analysis of the technical performance, respective costs and environmental effect of two invasive analytical methods (HPLC and UV/visible-FTIR) as compared to a new non-invasive analytical technique (Raman spectroscopy). Three pharmacotherapeutic models were used to compare the analytical performances of the three analytical techniques. Statistical inter-method correlation analysis was performed using non-parametric correlation rank tests. The study's economic component combined calculations relative to the depreciation of the equipment and the estimated cost of an AQC unit of work. In any case, analytical validation parameters of the three techniques were satisfactory, and strong correlations between the two spectroscopic techniques vs. HPLC were found. In addition, Raman spectroscopy was found to be superior as compared to the other techniques for numerous key criteria including a complete safety for operators and their occupational environment, a non-invasive procedure, no need for consumables, and a low operating cost. Finally, Raman spectroscopy appears superior for technical, economic and environmental objectives, as compared with the other invasive analytical methods. Copyright © 2014 Elsevier B.V. All rights reserved.

Analytical and Numerical solutions of a nonlinear alcoholism model via variable-order fractional differential equations

NASA Astrophysics Data System (ADS)

Gómez-Aguilar, J. F.

2018-03-01

In this paper, we analyze an alcoholism model which involves the impact of Twitter via Liouville-Caputo and Atangana-Baleanu-Caputo fractional derivatives with constant- and variable-order. Two fractional mathematical models are considered, with and without delay. Special solutions using an iterative scheme via Laplace and Sumudu transform were obtained. We studied the uniqueness and existence of the solutions employing the fixed point postulate. The generalized model with variable-order was solved numerically via the Adams method and the Adams-Bashforth-Moulton scheme. Stability and convergence of the numerical solutions were presented in details. Numerical examples of the approximate solutions are provided to show that the numerical methods are computationally efficient. Therefore, by including both the fractional derivatives and finite time delays in the alcoholism model studied, we believe that we have established a more complete and more realistic indicator of alcoholism model and affect the spread of the drinking.

Distribution of diameters for Erdős-Rényi random graphs.

PubMed

Hartmann, A K; Mézard, M

2018-03-01

We study the distribution of diameters d of Erdős-Rényi random graphs with average connectivity c. The diameter d is the maximum among all the shortest distances between pairs of nodes in a graph and an important quantity for all dynamic processes taking place on graphs. Here we study the distribution P(d) numerically for various values of c, in the nonpercolating and percolating regimes. Using large-deviation techniques, we are able to reach small probabilities like 10^{-100} which allow us to obtain the distribution over basically the full range of the support, for graphs up to N=1000 nodes. For values c<1, our results are in good agreement with analytical results, proving the reliability of our numerical approach. For c>1 the distribution is more complex and no complete analytical results are available. For this parameter range, P(d) exhibits an inflection point, which we found to be related to a structural change of the graphs. For all values of c, we determined the finite-size rate function Φ(d/N) and were able to extrapolate numerically to N→∞, indicating that the large-deviation principle holds.

Distribution of diameters for Erdős-Rényi random graphs

NASA Astrophysics Data System (ADS)

Hartmann, A. K.; Mézard, M.

2018-03-01

We study the distribution of diameters d of Erdős-Rényi random graphs with average connectivity c . The diameter d is the maximum among all the shortest distances between pairs of nodes in a graph and an important quantity for all dynamic processes taking place on graphs. Here we study the distribution P (d ) numerically for various values of c , in the nonpercolating and percolating regimes. Using large-deviation techniques, we are able to reach small probabilities like 10-100 which allow us to obtain the distribution over basically the full range of the support, for graphs up to N =1000 nodes. For values c <1 , our results are in good agreement with analytical results, proving the reliability of our numerical approach. For c >1 the distribution is more complex and no complete analytical results are available. For this parameter range, P (d ) exhibits an inflection point, which we found to be related to a structural change of the graphs. For all values of c , we determined the finite-size rate function Φ (d /N ) and were able to extrapolate numerically to N →∞ , indicating that the large-deviation principle holds.

On the superposition principle in interference experiments.

PubMed

Sinha, Aninda; H Vijay, Aravind; Sinha, Urbasi

2015-05-14

The superposition principle is usually incorrectly applied in interference experiments. This has recently been investigated through numerics based on Finite Difference Time Domain (FDTD) methods as well as the Feynman path integral formalism. In the current work, we have derived an analytic formula for the Sorkin parameter which can be used to determine the deviation from the application of the principle. We have found excellent agreement between the analytic distribution and those that have been earlier estimated by numerical integration as well as resource intensive FDTD simulations. The analytic handle would be useful for comparing theory with future experiments. It is applicable both to physics based on classical wave equations as well as the non-relativistic Schrödinger equation.

Low-density lipoprotein transport through an arterial wall under hyperthermia and hypertension conditions--An analytical solution.

PubMed

Iasiello, Marcello; Vafai, Kambiz; Andreozzi, Assunta; Bianco, Nicola

2016-01-25

An analytical solution for Low-Density Lipoprotein transport through an arterial wall under hyperthermia conditions is established in this work. A four-layer model is used to characterize the arterial wall. Transport governing equations are obtained as a combination between Staverman-Kedem-Katchalsky membrane equations and volume-averaged porous media equations. Temperature and solute transport fields are coupled by means of Ludwig-Soret effect. Results are in excellent agreement with numerical and analytical literature data under isothermal conditions, and with numerical literature data for the hyperthermia case. Effects of hypertension combined with hyperthermia, are also analyzed in this work. Copyright © 2015 Elsevier Ltd. All rights reserved.

Variations on Debris Disks. IV. An Improved Analytical Model for Collisional Cascades

NASA Astrophysics Data System (ADS)

Kenyon, Scott J.; Bromley, Benjamin C.

2017-04-01

We derive a new analytical model for the evolution of a collisional cascade in a thin annulus around a single central star. In this model, r max the size of the largest object changes with time, {r}\\max \\propto {t}-γ , with γ ≈ 0.1-0.2. Compared to standard models where r max is constant in time, this evolution results in a more rapid decline of M d , the total mass of solids in the annulus, and L d , the luminosity of small particles in the annulus: {M}d\\propto {t}-(γ +1) and {L}d\\propto {t}-(γ /2+1). We demonstrate that the analytical model provides an excellent match to a comprehensive suite of numerical coagulation simulations for annuli at 1 au and at 25 au. If the evolution of real debris disks follows the predictions of the analytical or numerical models, the observed luminosities for evolved stars require up to a factor of two more mass than predicted by previous analytical models.

Analytical approximation for the Einstein-dilaton-Gauss-Bonnet black hole metric

NASA Astrophysics Data System (ADS)

Kokkotas, K. D.; Konoplya, R. A.; Zhidenko, A.

2017-09-01

We construct an analytical approximation for the numerical black hole metric of P. Kanti et al. [Phys. Rev. D 54, 5049 (1996), 10.1103/PhysRevD.54.5049] in the four-dimensional Einstein-dilaton-Gauss-Bonnet (EdGB) theory. The continued fraction expansion in terms of a compactified radial coordinate, used here, converges slowly when the dilaton coupling approaches its extremal values, but for a black hole far from the extremal state, the analytical formula has a maximal relative error of a fraction of one percent already within the third order of the continued fraction expansion. The suggested analytical representation of the numerical black hole metric is relatively compact and a good approximation in the whole space outside the black hole event horizon. Therefore, it can serve in the same way as an exact solution when analyzing particles' motion, perturbations, quasinormal modes, Hawking radiation, accreting disks, and many other problems in the vicinity of a black hole. In addition, we construct the approximate analytical expression for the dilaton field.

On the analytic and numeric optimisation of airplane trajectories under real atmospheric conditions

NASA Astrophysics Data System (ADS)

Gonzalo, J.; Domínguez, D.; López, D.

2014-12-01

From the beginning of aviation era, economic constraints have forced operators to continuously improve the planning of the flights. The revenue is proportional to the cost per flight and the airspace occupancy. Many methods, the first started in the middle of last century, have explore analytical, numerical and artificial intelligence resources to reach the optimal flight planning. In parallel, advances in meteorology and communications allow an almost real-time knowledge of the atmospheric conditions and a reliable, error-bounded forecast for the near future. Thus, apart from weather risks to be avoided, airplanes can dynamically adapt their trajectories to minimise their costs. International regulators are aware about these capabilities, so it is reasonable to envisage some changes to allow this dynamic planning negotiation to soon become operational. Moreover, current unmanned airplanes, very popular and often small, suffer the impact of winds and other weather conditions in form of dramatic changes in their performance. The present paper reviews analytic and numeric solutions for typical trajectory planning problems. Analytic methods are those trying to solve the problem using the Pontryagin principle, where influence parameters are added to state variables to form a split condition differential equation problem. The system can be solved numerically -indirect optimisation- or using parameterised functions -direct optimisation-. On the other hand, numerical methods are based on Bellman's dynamic programming (or Dijkstra algorithms), where the fact that two optimal trajectories can be concatenated to form a new optimal one if the joint point is demonstrated to belong to the final optimal solution. There is no a-priori conditions for the best method. Traditionally, analytic has been more employed for continuous problems whereas numeric for discrete ones. In the current problem, airplane behaviour is defined by continuous equations, while wind fields are given in a discrete grid at certain time intervals. The research demonstrates advantages and disadvantages of each method as well as performance figures of the solutions found for typical flight conditions under static and dynamic atmospheres. This provides significant parameters to be used in the selection of solvers for optimal trajectories.

Microwave reflection, transmission, and absorption by human brain tissue

NASA Astrophysics Data System (ADS)

Ansari, M. A.; Akhlaghipour, N.; Zarei, M.; Niknam, A. R.

2018-04-01

These days, the biological effects of electromagnetic (EM) radiations on the brain, especially in the frequency range of mobile communications, have caught the attention of many scientists. Therefore, in this paper, the propagation of mobile phone electromagnetic waves in the brain tissues is investigated analytically and numerically. The brain is modeled by three layers consisting of skull, grey and white matter. First, we have analytically calculated the microwave reflection, transmission, and absorption coefficients using signal flow graph technique. The effect of microwave frequency and variations in the thickness of layers on the propagation of microwave through brain are studied. Then, the penetration of microwave in the layers is numerically investigated by Monte Carlo method. It is shown that the analytical results are in good agreement with those obtained by Monte Carlo method. Our results indicate the absorbed microwave energy depends on microwave frequency and thickness of brain layers, and the absorption coefficient is optimized at a number of frequencies. These findings can be used for comparing the microwave absorbed energy in a child's and adult's brain.

Phosphorescence Kinetics of Singlet Oxygen Produced by Photosensitization in Spherical Nanoparticles. Part I. Theory.

PubMed

Hovan, Andrej; Datta, Shubhashis; Kruglik, Sergei G; Jancura, Daniel; Miskovsky, Pavol; Bánó, Gregor

2018-05-24

The singlet oxygen produced by energy transfer between an excited photosensitizer (pts) and ground-state oxygen molecules plays a key role in photodynamic therapy. Different nanocarrier systems are extensively studied to promote targeted pts delivery in a host body. The phosphorescence kinetics of the singlet oxygen produced by the short laser pulse photosensitization of pts inside nanoparticles is influenced by singlet oxygen diffusion from the particles to the surrounding medium. Two theoretical models are presented in this work: a more complex numerical one and a simple analytical one. Both the models predict the time course of singlet oxygen concentration inside and outside of the spherical particles following short-pulse excitation of pts. On the basis of the comparison of the numerical and analytical results, a semiempirical analytical formula is derived to calculate the characteristic diffusion time of singlet oxygen from the nanoparticles to the surrounding solvent. The phosphorescence intensity of singlet oxygen produced in pts-loaded nanocarrier systems can be calculated as a linear combination of the two concentrations (inside and outside the particles), taking the different phosphorescence emission rate constants into account.

Vector solitons in coupled nonlinear Schrödinger equations with spatial stimulated scattering and inhomogeneous dispersion

NASA Astrophysics Data System (ADS)

Gromov, E. M.; Malomed, B. A.; Tyutin, V. V.

2018-01-01

The dynamics of two-component solitons is studied, analytically and numerically, in the framework of a system of coupled extended nonlinear Schrödinger equations, which incorporate the cross-phase modulation, pseudo-stimulated-Raman-scattering (pseudo-SRS), cross-pseudo-SRS, and spatially inhomogeneous second-order dispersion (SOD). The system models co-propagation of electromagnetic waves with orthogonal polarizations in plasmas. It is shown that the soliton's wavenumber downshift, caused by pseudo-SRS, may be compensated by an upshift, induced by the inhomogeneous SOD, to produce stable stationary two-component solitons. The corresponding approximate analytical solutions for stable solitons are found. Analytical results are well confirmed by their numerical counterparts. Further, the evolution of inputs composed of spatially even and odd components is investigated by means of systematic simulations, which reveal three different outcomes: formation of a breather which keeps opposite parities of the components; splitting into a pair of separating vector solitons; and spreading of the weak odd component into a small-amplitude pedestal with an embedded dark soliton.

Skin friction enhancement in a model problem of undulatory swimming

NASA Astrophysics Data System (ADS)

Ehrenstein, Uwe; Eloy, Christophe

2013-10-01

To calculate the energy costs of swimming, it is crucial to evaluate the drag force originating from skin friction. In this paper we examine the assumption, known as the 'Bone-Lighthill boundary-layer thinning hypothesis', that undulatory swimming motions induce a drag increase because of the compression of the boundary layer. Studying analytically an incoming flow along a flat plate moving at a normal velocity as a limit case of a yawed cylinder in uniform flow under the laminar boundary layer assumption, we demonstrate that the longitudinal drag scales as the square root of the normal velocity component. This analytical prediction is interpreted in the light of a three-dimensional numerical simulation result for a plate of finite length and width. An analogous two-dimensional Navier-Stokes problem by artificially accelerating the flow in a channel of finite height is proposed and solved numerically, showing the robustness of the analytical results. Solving the problem for an undulatory plate motion similar to fish swimming, we find a drag enhancement which can be estimated to be of the order of 20 %.

Coronal heating by the resonant absorption of Alfven waves - Importance of the global mode and scaling laws

NASA Technical Reports Server (NTRS)

Steinolfson, Richard S.; Davila, Joseph M.

1993-01-01

Numerical simulations of the MHD equations for a fully compressible, low-beta, resistive plasma are used to study the resonance absorption process for the heating of coronal active region loops. Comparisons with more approximate analytic models show that the major predictions of the analytic theories are, to a large extent, confirmed by the numerical computations. The simulations demonstrate that the dissipation occurs primarily in a thin resonance layer. Some of the analytically predicted features verified by the simulations are (a) the position of the resonance layer within the initial inhomogeneity; (b) the importance of the global mode for a large range of loop densities; (c) the dependence of the resonance layer thickness and the steady-state heating rate on the dissipation coefficient; and (d) the time required for the resonance layer to form. In contrast with some previous analytic and simulation results, the time for the loop to reach a steady state is found to be the phase-mixing time rather than a dissipation time. This disagreement is shown to result from neglect of the existence of the global mode in some of the earlier analyses. The resonant absorption process is also shown to behave similar to a classical driven harmonic oscillator.

Numerical investigation of band gaps in 3D printed cantilever-in-mass metamaterials

NASA Astrophysics Data System (ADS)

Qureshi, Awais; Li, Bing; Tan, K. T.

2016-06-01

In this research, the negative effective mass behavior of elastic/mechanical metamaterials is exhibited by a cantilever-in-mass structure as a proposed design for creating frequency stopping band gaps, based on local resonance of the internal structure. The mass-in-mass unit cell model is transformed into a cantilever-in-mass model using the Bernoulli-Euler beam theory. An analytical model of the cantilever-in-mass structure is derived and the effects of geometrical dimensions and material parameters to create frequency band gaps are examined. A two-dimensional finite element model is created to validate the analytical results, and excellent agreement is achieved. The analytical model establishes an easily tunable metamaterial design to realize wave attenuation based on locally resonant frequency. To demonstrate feasibility for 3D printing, the analytical model is employed to design and fabricate 3D printable mechanical metamaterial. A three-dimensional numerical experiment is performed using COMSOL Multiphysics to validate the wave attenuation performance. Results show that the cantilever-in-mass metamaterial is capable of mitigating stress waves at the desired resonance frequency. Our study successfully presents the use of one constituent material to create a 3D printed cantilever-in-mass metamaterial with negative effective mass density for stress wave mitigation purposes.

Long Term Evolution of Planetary Systems with a Terrestrial Planet and a Giant Planet

NASA Technical Reports Server (NTRS)

Georgakarakos, Nikolaos; Dobbs-Dixon, Ian; Way, Michael J.

2016-01-01

We study the long term orbital evolution of a terrestrial planet under the gravitational perturbations of a giant planet. In particular, we are interested in situations where the two planets are in the same plane and are relatively close. We examine both possible configurations: the giant planet orbit being either outside or inside the orbit of the smaller planet. The perturbing potential is expanded to high orders and an analytical solution of the terrestrial planetary orbit is derived. The analytical estimates are then compared against results from the numerical integration of the full equations of motion and we find that the analytical solution works reasonably well. An interesting finding is that the new analytical estimates improve greatly the predictions for the timescales of the orbital evolution of the terrestrial planet compared to an octupole order expansion. Finally, we briefly discuss possible applications of the analytical estimates in astrophysical problems.

TimeBench: a data model and software library for visual analytics of time-oriented data.

PubMed

Rind, Alexander; Lammarsch, Tim; Aigner, Wolfgang; Alsallakh, Bilal; Miksch, Silvia

2013-12-01

Time-oriented data play an essential role in many Visual Analytics scenarios such as extracting medical insights from collections of electronic health records or identifying emerging problems and vulnerabilities in network traffic. However, many software libraries for Visual Analytics treat time as a flat numerical data type and insufficiently tackle the complexity of the time domain such as calendar granularities and intervals. Therefore, developers of advanced Visual Analytics designs need to implement temporal foundations in their application code over and over again. We present TimeBench, a software library that provides foundational data structures and algorithms for time-oriented data in Visual Analytics. Its expressiveness and developer accessibility have been evaluated through application examples demonstrating a variety of challenges with time-oriented data and long-term developer studies conducted in the scope of research and student projects.

The possible equilibrium shapes of static pendant drops

NASA Astrophysics Data System (ADS)

Sumesh, P. T.; Govindarajan, Rama

2010-10-01

Analytical and numerical studies are carried out on the shapes of two-dimensional and axisymmetric pendant drops hanging under gravity from a solid surface. Drop shapes with both pinned and equilibrium contact angles are obtained naturally from a single boundary condition in the analytical energy optimization procedure. The numerical procedure also yields optimum energy shapes, satisfying Young's equation without the explicit imposition of a boundary condition at the plate. It is shown analytically that a static pendant two-dimensional drop can never be longer than 3.42 times the capillary length. A related finding is that a range of existing solutions for long two-dimensional drops correspond to unphysical drop shapes. Therefore, two-dimensional drops of small volume display only one static solution. In contrast, it is known that axisymmetric drops can display multiple solutions for a given volume. We demonstrate numerically that there is no limit to the height of multiple-lobed Kelvin drops, but the total volume is finite, with the volume of successive lobes forming a convergent series. The stability of such drops is in question, though. Drops of small volume can attain large heights. A bifurcation is found within the one-parameter space of Laplacian shapes, with a range of longer drops displaying a minimum in energy in the investigated space. Axisymmetric Kelvin drops exhibit an infinite number of bifurcations.

N -jettiness subtractions for g g →H at subleading power

NASA Astrophysics Data System (ADS)

Moult, Ian; Rothen, Lorena; Stewart, Iain W.; Tackmann, Frank J.; Zhu, Hua Xing

2018-01-01

N -jettiness subtractions provide a general approach for performing fully-differential next-to-next-to-leading order (NNLO) calculations. Since they are based on the physical resolution variable N -jettiness, TN , subleading power corrections in τ =TN/Q , with Q a hard interaction scale, can also be systematically computed. We study the structure of power corrections for 0-jettiness, T0, for the g g →H process. Using the soft-collinear effective theory we analytically compute the leading power corrections αsτ ln τ and αs2τ ln3τ (finding partial agreement with a previous result in the literature), and perform a detailed numerical study of the power corrections in the g g , g q , and q q ¯ channels. This includes a numerical extraction of the αsτ and αs2τ ln2τ corrections, and a study of the dependence on the T0 definition. Including such power suppressed logarithms significantly reduces the size of missing power corrections, and hence improves the numerical efficiency of the subtraction method. Having a more detailed understanding of the power corrections for both q q ¯ and g g initiated processes also provides insight into their universality, and hence their behavior in more complicated processes where they have not yet been analytically calculated.

Surface response of a viscoelastic medium to subsurface acoustic sources with application to medical diagnosis

NASA Astrophysics Data System (ADS)

Royston, Thomas J.; Yazicioglu, Yigit; Loth, Francis

2003-02-01

The response at the surface of an isotropic viscoelastic medium to buried fundamental acoustic sources is studied theoretically, computationally and experimentally. Finite and infinitesimal monopole and dipole sources within the low audible frequency range (40-400 Hz) are considered. Analytical and numerical integral solutions that account for compression, shear and surface wave response to the buried sources are formulated and compared with numerical finite element simulations and experimental studies on finite dimension phantom models. It is found that at low audible frequencies, compression and shear wave propagation from point sources can both be significant, with shear wave effects becoming less significant as frequency increases. Additionally, it is shown that simple closed-form analytical approximations based on an infinite medium model agree well with numerically obtained ``exact'' half-space solutions for the frequency range and material of interest in this study. The focus here is on developing a better understanding of how biological soft tissue affects the transmission of vibro-acoustic energy from biological acoustic sources below the skin surface, whose typical spectral content is in the low audible frequency range. Examples include sound radiated from pulmonary, gastro-intestinal and cardiovascular system functions, such as breath sounds, bowel sounds and vascular bruits, respectively.

Analytical approach on the stiffness of MR fluid filled spring

NASA Astrophysics Data System (ADS)

Sikulskyi, Stanislav; Kim, Daewon

2017-04-01

A solid mechanical spring generally exhibits uniform stiffness. This paper studies a mechanical spring filled with magnetorheological (MR) fluid to achieve controllable stiffness. The hollow spring filled with MR fluid is subjected to a controlled magnetic field in order to change the viscosity of the MR fluid and thereby to change the overall stiffness of the spring. MR fluid is considered as a Bingham viscoplastic linear material in the mathematical model. The goal of this research is to study the feasibility of such spring system by analytically computing the effects of MR fluid on overall spring stiffness. For this purpose, spring mechanics and MR fluid behavior are studied to increase the accuracy of the analysis. Numerical simulations are also performed to generate some assumptions, which simplify calculations in the analytical part. The accuracy of the present approach is validated by comparing the analytical results to previously known experimental results. Overall stiffness variations of the spring are also discussed for different spring designs.

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

Canavan, G.H.

Optimizations of missile allocation based on linearized exchange equations produce accurate allocations, but the limits of validity of the linearization are not known. These limits are explored in the context of the upload of weapons by one side to initially small, equal forces of vulnerable and survivable weapons. The analysis compares analytic and numerical optimizations and stability induces based on aggregated interactions of the two missile forces, the first and second strikes they could deliver, and they resulting costs. This note discusses the costs and stability indices induced by unilateral uploading of weapons to an initially symmetrical low force configuration.more » These limits are quantified for forces with a few hundred missiles by comparing analytic and numerical optimizations of first strike costs. For forces of 100 vulnerable and 100 survivable missiles on each side, the analytic optimization agrees closely with the numerical solution. For 200 vulnerable and 200 survivable missiles on each side, the analytic optimization agrees with the induces to within about 10%, but disagrees with the allocation of the side with more weapons by about 50%. The disagreement comes from the interaction of the possession of more weapons with the shift of allocation from missiles to value that they induce.« less

Analytical and numerical solutions for heat transfer and effective thermal conductivity of cracked media

NASA Astrophysics Data System (ADS)

Tran, A. B.; Vu, M. N.; Nguyen, S. T.; Dong, T. Q.; Le-Nguyen, K.

2018-02-01

This paper presents analytical solutions to heat transfer problems around a crack and derive an adaptive model for effective thermal conductivity of cracked materials based on singular integral equation approach. Potential solution of heat diffusion through two-dimensional cracked media, where crack filled by air behaves as insulator to heat flow, is obtained in a singular integral equation form. It is demonstrated that the temperature field can be described as a function of temperature and rate of heat flow on the boundary and the temperature jump across the cracks. Numerical resolution of this boundary integral equation allows determining heat conduction and effective thermal conductivity of cracked media. Moreover, writing this boundary integral equation for an infinite medium embedding a single crack under a far-field condition allows deriving the closed-form solution of temperature discontinuity on the crack and particularly the closed-form solution of temperature field around the crack. These formulas are then used to establish analytical effective medium estimates. Finally, the comparison between the developed numerical and analytical solutions allows developing an adaptive model for effective thermal conductivity of cracked media. This model takes into account both the interaction between cracks and the percolation threshold.

An innovative hybrid 3D analytic-numerical model for air breathing parallel channel counter-flow PEM fuel cells.

PubMed

Tavčar, Gregor; Katrašnik, Tomaž

2014-01-01

The parallel straight channel PEM fuel cell model presented in this paper extends the innovative hybrid 3D analytic-numerical (HAN) approach previously published by the authors with capabilities to address ternary diffusion systems and counter-flow configurations. The model's core principle is modelling species transport by obtaining a 2D analytic solution for species concentration distribution in the plane perpendicular to the cannel gas-flow and coupling consecutive 2D solutions by means of a 1D numerical pipe-flow model. Electrochemical and other nonlinear phenomena are coupled to the species transport by a routine that uses derivative approximation with prediction-iteration. The latter is also the core of the counter-flow computation algorithm. A HAN model of a laboratory test fuel cell is presented and evaluated against a professional 3D CFD simulation tool showing very good agreement between results of the presented model and those of the CFD simulation. Furthermore, high accuracy results are achieved at moderate computational times, which is owed to the semi-analytic nature and to the efficient computational coupling of electrochemical kinetics and species transport.

Preliminary numerical analysis of improved gas chromatograph model

NASA Technical Reports Server (NTRS)

Woodrow, P. T.

1973-01-01

A mathematical model for the gas chromatograph was developed which incorporates the heretofore neglected transport mechanisms of intraparticle diffusion and rates of adsorption. Because a closed-form analytical solution to the model does not appear realizable, techniques for the numerical solution of the model equations are being investigated. Criteria were developed for using a finite terminal boundary condition in place of an infinite boundary condition used in analytical solution techniques. The class of weighted residual methods known as orthogonal collocation is presently being investigated and appears promising.

Application of differential transformation method for solving dengue transmission mathematical model

NASA Astrophysics Data System (ADS)

Ndii, Meksianis Z.; Anggriani, Nursanti; Supriatna, Asep K.

2018-03-01

The differential transformation method (DTM) is a semi-analytical numerical technique which depends on Taylor series and has application in many areas including Biomathematics. The aim of this paper is to employ the differential transformation method (DTM) to solve system of non-linear differential equations for dengue transmission mathematical model. Analytical and numerical solutions are determined and the results are compared to that of Runge-Kutta method. We found a good agreement between DTM and Runge-Kutta method.

Non-degenerate two-photon absorption in silicon waveguides. Analytical and experimental study

DOE PAGES

Zhang, Yanbing; Husko, Chad; Lefrancois, Simon; ...

2015-06-22

We theoretically and experimentally investigate the nonlinear evolution of two optical pulses in a silicon waveguide. We provide an analytic solution for the weak probe wave undergoing non-degenerate two-photon absorption (TPA) from the strong pump. At larger pump intensities, we employ a numerical solution to study the interplay between TPA and photo-generated free carriers. We develop a simple and powerful approach to extract and separate out the distinct loss contributions of TPA and free-carrier absorption from readily available experimental data. Our analysis accounts accurately for experimental results in silicon photonic crystal waveguides.

Cogeneration Technology Alternatives Study (CTAS). Volume 2: Analytical approach

NASA Technical Reports Server (NTRS)

Gerlaugh, H. E.; Hall, E. W.; Brown, D. H.; Priestley, R. R.; Knightly, W. F.

1980-01-01

The use of various advanced energy conversion systems were compared with each other and with current technology systems for their savings in fuel energy, costs, and emissions in individual plants and on a national level. The ground rules established by NASA and assumptions made by the General Electric Company in performing this cogeneration technology alternatives study are presented. The analytical methodology employed is described in detail and is illustrated with numerical examples together with a description of the computer program used in calculating over 7000 energy conversion system-industrial process applications. For Vol. 1, see 80N24797.

Internal-external stimulus competition in a system of interacting moving particles: Persuasion versus propaganda

NASA Astrophysics Data System (ADS)

Clementi, N. C.; Revelli, J. A.; Sibona, G. J.

2015-07-01

We propose a general nonlinear analytical framework to study the effect of an external stimulus in the internal state of a population of moving particles. This novel scheme allows us to study a broad range of excitation transport phenomena. In particular, considering social systems, it gives insight of the spatial dynamics influence in the competition between propaganda (mass media) and convincement. By extending the framework presented by Terranova et al. [Europhys. Lett. 105, 30007 (2014), 10.1209/0295-5075/105/30007], we now allow changes in individual's opinions due to a reflection induced by mass media. The equations of the model could be solved numerically, and, for some special cases, it is possible to derive analytical solutions for the steady states. We implement computational simulations for different social and dynamical systems to check the accuracy of our scheme and to study a broader variety of scenarios. In particular, we compare the numerical outcome with the analytical results for two possible real cases, finding a good agreement. From the results, we observe that mass media dominates the opinion state in slow dynamics communities; whereas, for higher agent active speeds, the rate of interactions increases and the opinion state is determined by a competition between propaganda and persuasion. This difference suggests that kinetics can not be neglected in the study of transport of any excitation over a particle system.

Internal-external stimulus competition in a system of interacting moving particles: Persuasion versus propaganda.

PubMed

Clementi, N C; Revelli, J A; Sibona, G J

2015-07-01

We propose a general nonlinear analytical framework to study the effect of an external stimulus in the internal state of a population of moving particles. This novel scheme allows us to study a broad range of excitation transport phenomena. In particular, considering social systems, it gives insight of the spatial dynamics influence in the competition between propaganda (mass media) and convincement. By extending the framework presented by Terranova et al. [Europhys. Lett. 105, 30007 (2014)], we now allow changes in individual's opinions due to a reflection induced by mass media. The equations of the model could be solved numerically, and, for some special cases, it is possible to derive analytical solutions for the steady states. We implement computational simulations for different social and dynamical systems to check the accuracy of our scheme and to study a broader variety of scenarios. In particular, we compare the numerical outcome with the analytical results for two possible real cases, finding a good agreement. From the results, we observe that mass media dominates the opinion state in slow dynamics communities; whereas, for higher agent active speeds, the rate of interactions increases and the opinion state is determined by a competition between propaganda and persuasion. This difference suggests that kinetics can not be neglected in the study of transport of any excitation over a particle system.

Study of tethered satellite active attitude control

NASA Technical Reports Server (NTRS)

Colombo, G.

1982-01-01

Existing software was adapted for the study of tethered subsatellite rotational dynamics, an analytic solution for a stable configuration of a tethered subsatellite was developed, the analytic and numerical integrator (computer) solutions for this "test case' was compared in a two mass tether model program (DUMBEL), the existing multiple mass tether model (SKYHOOK) was modified to include subsatellite rotational dynamics, the analytic "test case,' was verified, and the use of the SKYHOOK rotational dynamics capability with a computer run showing the effect of a single off axis thruster on the behavior of the subsatellite was demonstrated. Subroutines for specific attitude control systems are developed and applied to the study of the behavior of the tethered subsatellite under realistic on orbit conditions. The effect of all tether "inputs,' including pendular oscillations, air drag, and electrodynamic interactions, on the dynamic behavior of the tether are included.

The thermoelectric properties of strongly correlated systems

NASA Astrophysics Data System (ADS)

Cai, Jianwei

Strongly correlated systems are among the most interesting and complicated systems in physics. Large Seebeck coefficients are found in some of these systems, which highlight the possibility for thermoelectric applications. In this thesis, we study the thermoelectric properties of these strongly correlated systems with various methods. We derived analytic formulas for the resistivity and Seebeck coefficient of the periodic Anderson model based on the dynamic mean field theory. These formulas were possible as the self energy of the single impurity Anderson model could be given by an analytic ansatz derived from experiments and numerical calculations instead of complicated numerical calculations. The results show good agreement with the experimental data of rare-earth compound in a restricted temperature range. These formulas help to understand the properties of periodic Anderson model. Based on the study of rare-earth compounds, we proposed a design for the thermoelectric meta-material. This manmade material is made of quantum dots linked by conducting linkers. The quantum dots act as the rare-earth atoms with heavier mass. We set up a model similar to the periodic Anderson model for this new material. The new model was studied with the perturbation theory for energy bands. The dynamic mean field theory with numerical renormalization group as the impurity solver was used to study the transport properties. With these studies, we confirmed the improved thermoelectric properties of the designed material.

Adaptive numerical algorithms to simulate the dynamical Casimir effect in a closed cavity with different boundary conditions

NASA Astrophysics Data System (ADS)

Villar, Paula I.; Soba, Alejandro

2017-07-01

We present an alternative numerical approach to compute the number of particles created inside a cavity due to time-dependent boundary conditions. The physical model consists of a rectangular cavity, where a wall always remains still while the other wall of the cavity presents a smooth movement in one direction. The method relies on the setting of the boundary conditions (Dirichlet and Neumann) and the following resolution of the corresponding equations of modes. By a further comparison between the ground state before and after the movement of the cavity wall, we finally compute the number of particles created. To demonstrate the method, we investigate the creation of particle production in vibrating cavities, confirming previously known results in the appropriate limits. Within this approach, the dynamical Casimir effect can be investigated, making it possible to study a variety of scenarios where no analytical results are known. Of special interest is, of course, the realistic case of the electromagnetic field in a three-dimensional cavity, with transverse electric (TE)-mode and transverse magnetic (TM)-mode photon production. Furthermore, with our approach we are able to calculate numerically the particle creation in a tuneable resonant superconducting cavity by the use of the generalized Robin boundary condition. We compare the numerical results with analytical predictions as well as a different numerical approach. Its extension to three dimensions is also straightforward.

Heat Transfer Analysis of Thermal Protection Structures for Hypersonic Vehicles

NASA Astrophysics Data System (ADS)

Zhou, Chen; Wang, Zhijin; Hou, Tianjiao

2017-11-01

This research aims to develop an analytical approach to study the heat transfer problem of thermal protection systems (TPS) for hypersonic vehicles. Laplace transform and integral method are used to describe the temperature distribution through the TPS subject to aerodynamic heating during flight. Time-dependent incident heat flux is also taken into account. Two different cases with heat flux and radiation boundary conditions are studied and discussed. The results are compared with those obtained by finite element analyses and show a good agreement. Although temperature profiles of such problems can be readily accessed via numerical simulations, analytical solutions give a greater insight into the physical essence of the heat transfer problem. Furthermore, with the analytical approach, rapid thermal analyses and even thermal optimization can be achieved during the preliminary TPS design.

Comparative study of soil erodibility and critical shear stress between loess and purple soils

NASA Astrophysics Data System (ADS)

Xing, Hang; Huang, Yu-han; Chen, Xiao-yan; Luo, Bang-lin; Mi, Hong-xing

2018-03-01

Loess and purple soils are two very important cultivated soils, with the former in the loess region and the latter in southern sub-tropical region of China, featured with high-risks of erosion, considerable differences of soil structures due to differences in mineral and nutrient compositions. Study on soil erodibility (Kr) and critical shear stress (τc) of these two soils is beneficial to predict soil erosion with such models as WEPP. In this study, rill erosion experimental data sets of the two soils are used for estimating their Kr and τc before they are compared to understand their differences of rill erosion behaviors. The maximum detachment rates of the loess and purple soils are calculated under different hydrodynamic conditions (flow rates: 2, 4, 8 L/min; slope gradients: 5°, 10°, 15°, 20°, 25°) through analytical and numerical methods respectively. Analytical method used the derivative of the function between sediment concentration and rill length to estimate potential detachment rates, at the rill beginning. Numerical method estimated potential detachment rates with the experimental data, at the rill beginning and 0.5 m location. The Kr and τc of these two soils are determined by the linear equation based on experimental data. Results show that the methods could well estimate the Kr and τc of these two soils as they remain basically unchanged under different hydrodynamic conditions. The Kr value of loess soil is about twice of the purple soil, whereas the τc is about half of that. The numerical results have good correlations with the analytical values. These results can be useful in modeling rill erosion processes of loess and purple soils.

Numerical Uncertainty Analysis for Computational Fluid Dynamics using Student T Distribution -- Application of CFD Uncertainty Analysis Compared to Exact Analytical Solution

NASA Technical Reports Server (NTRS)

Groves, Curtis E.; Ilie, marcel; Shallhorn, Paul A.

2014-01-01

Computational Fluid Dynamics (CFD) is the standard numerical tool used by Fluid Dynamists to estimate solutions to many problems in academia, government, and industry. CFD is known to have errors and uncertainties and there is no universally adopted method to estimate such quantities. This paper describes an approach to estimate CFD uncertainties strictly numerically using inputs and the Student-T distribution. The approach is compared to an exact analytical solution of fully developed, laminar flow between infinite, stationary plates. It is shown that treating all CFD input parameters as oscillatory uncertainty terms coupled with the Student-T distribution can encompass the exact solution.

Preliminary Upper Estimate of Peak Currents in Transcranial Magnetic Stimulation at Distant Locations from a TMS Coil

PubMed Central

Makarov, Sergey N.; Yanamadala, Janakinadh; Piazza, Matthew W.; Helderman, Alex M.; Thang, Niang S.; Burnham, Edward H.; Pascual-Leone, Alvaro

2016-01-01

Goals Transcranial magnetic stimulation (TMS) is increasingly used as a diagnostic and therapeutic tool for numerous neuropsychiatric disorders. The use of TMS might cause whole-body exposure to undesired induced currents in patients and TMS operators. The aim of the present study is to test and justify a simple analytical model known previously, which may be helpful as an upper estimate of eddy current density at a particular distant observation point for any body composition and any coil setup. Methods We compare the analytical solution with comprehensive adaptive mesh refinement-based FEM simulations of a detailed full-body human model, two coil types, five coil positions, about 100,000 observation points, and two distinct pulse rise times, thus providing a representative number of different data sets for comparison, while also using other numerical data. Results Our simulations reveal that, after a certain modification, the analytical model provides an upper estimate for the eddy current density at any location within the body. In particular, it overestimates the peak eddy currents at distant locations from a TMS coil by a factor of 10 on average. Conclusion The simple analytical model tested in the present study may be valuable as a rapid method to safely estimate levels of TMS currents at different locations within a human body. Significance At present, safe limits of general exposure to TMS electric and magnetic fields are an open subject, including fetal exposure for pregnant women. PMID:26685221

Analytic Closed-Form Solution of a Mixed Layer Model for Stratocumulus Clouds

NASA Astrophysics Data System (ADS)

Akyurek, Bengu Ozge

Stratocumulus clouds play an important role in climate cooling and are hard to predict using global climate and weather forecast models. Thus, previous studies in the literature use observations and numerical simulation tools, such as large-eddy simulation (LES), to solve the governing equations for the evolution of stratocumulus clouds. In contrast to the previous works, this work provides an analytic closed-form solution to the cloud thickness evolution of stratocumulus clouds in a mixed-layer model framework. With a focus on application over coastal lands, the diurnal cycle of cloud thickness and whether or not clouds dissipate are of particular interest. An analytic solution enables the sensitivity analysis of implicitly interdependent variables and extrema analysis of cloud variables that are hard to achieve using numerical solutions. In this work, the sensitivity of inversion height, cloud-base height, and cloud thickness with respect to initial and boundary conditions, such as Bowen ratio, subsidence, surface temperature, and initial inversion height, are studied. A critical initial cloud thickness value that can be dissipated pre- and post-sunrise is provided. Furthermore, an extrema analysis is provided to obtain the minima and maxima of the inversion height and cloud thickness within 24 h. The proposed solution is validated against LES results under the same initial and boundary conditions. Then, the proposed analytic framework is extended to incorporate multiple vertical columns that are coupled by advection through wind flow. This enables a bridge between the micro-scale and the mesoscale relations. The effect of advection on cloud evolution is studied and a sensitivity analysis is provided.

Boundary enhanced effects on the existence of quadratic solitons

NASA Astrophysics Data System (ADS)

Chen, Manna; Zhang, Ting; Li, Wenjie; Lu, Daquan; Guo, Qi; Hu, Wei

2018-05-01

We investigate, both analytically and numerically, the boundary enhanced effects exerted on the quadratic solitons consisting of fundamental waves and oscillatory second harmonics in the presence of boundary conditions. The nonlocal analogy predicts that the soliton for fundamental wave is supported by the balance between equivalent nonlinear confinement and diffraction (or dispersion). Under Snyder and Mitchell's strongly nonlocal approximation, we obtain the analytical soliton solutions both with and without the boundary conditions to show the impact of boundary conditions. We can distinguish explicitly the nonlinear confinement between the second harmonic mutual interaction and the enhanced effects caused by remote boundaries. Those boundary enhanced effects on the existence of solitons can be positive or negative, which depend on both sample size and nonlocal parameter. The piecewise existence regime of solitons can be explained analytically. The analytical soliton solutions are verified by the numerical ones and the discrepancy between them is also discussed.

Finite analytic numerical solution of heat transfer and flow past a square channel cavity

NASA Technical Reports Server (NTRS)

Chen, C.-J.; Obasih, K.

1982-01-01

A numerical solution of flow and heat transfer characteristics is obtained by the finite analytic method for a two dimensional laminar channel flow over a two-dimensional square cavity. The finite analytic method utilizes the local analytic solution in a small element of the problem region to form the algebraic equation relating an interior nodal value with its surrounding nodal values. Stable and rapidly converged solutions were obtained for Reynolds numbers ranging to 1000 and Prandtl number to 10. Streamfunction, vorticity and temperature profiles are solved. Local and mean Nusselt number are given. It is found that the separation streamlines between the cavity and channel flow are concave into the cavity at low Reynolds number and convex at high Reynolds number (Re greater than 100) and for square cavity the mean Nusselt number may be approximately correlated with Peclet number as Nu(m) = 0.365 Pe exp 0.2.

Nonlinear analysis for dual-frequency concurrent energy harvesting

NASA Astrophysics Data System (ADS)

Yan, Zhimiao; Lei, Hong; Tan, Ting; Sun, Weipeng; Huang, Wenhu

2018-05-01

The dual-frequency responses of the hybrid energy harvester undergoing the base excitation and galloping were analyzed numerically. In this work, an approximate dual-frequency analytical method is proposed for the nonlinear analysis of such a system. To obtain the approximate analytical solutions of the full coupled distributed-parameter model, the forcing interactions is first neglected. Then, the electromechanical decoupled governing equation is developed using the equivalent structure method. The hybrid mechanical response is finally separated to be the self-excited and forced responses for deriving the analytical solutions, which are confirmed by the numerical simulations of the full coupled model. The forced response has great impacts on the self-excited response. The boundary of Hopf bifurcation is analytically determined by the onset wind speed to galloping, which is linearly increased by the electrical damping. Quenching phenomenon appears when the increasing base excitation suppresses the galloping. The theoretical quenching boundary depends on the forced mode velocity. The quenching region increases with the base acceleration and electrical damping, but decreases with the wind speed. Superior to the base-excitation-alone case, the existence of the aerodynamic force protects the hybrid energy harvester at resonance from damages caused by the excessive large displacement. From the view of the harvested power, the hybrid system surpasses the base-excitation-alone system or the galloping-alone system. This study advances our knowledge on intrinsic nonlinear dynamics of the dual-frequency energy harvesting system by taking advantage of the analytical solutions.

Robustness of high-fidelity Rydberg gates with single-site addressability

NASA Astrophysics Data System (ADS)

Goerz, Michael H.; Halperin, Eli J.; Aytac, Jon M.; Koch, Christiane P.; Whaley, K. Birgitta

2014-09-01

Controlled-phase (cphase) gates can be realized with trapped neutral atoms by making use of the Rydberg blockade. Achieving the ultrahigh fidelities required for quantum computation with such Rydberg gates, however, is compromised by experimental inaccuracies in pulse amplitudes and timings, as well as by stray fields that cause fluctuations of the Rydberg levels. We report here a comparative study of analytic and numerical pulse sequences for the Rydberg cphase gate that specifically examines the robustness of the gate fidelity with respect to such experimental perturbations. Analytical pulse sequences of both simultaneous and stimulated Raman adiabatic passage (STIRAP) are found to be at best moderately robust under these perturbations. In contrast, optimal control theory is seen to allow generation of numerical pulses that are inherently robust within a predefined tolerance window. The resulting numerical pulse shapes display simple modulation patterns and can be rationalized in terms of an interference between distinct two-photon Rydberg excitation pathways. Pulses of such low complexity should be experimentally feasible, allowing gate fidelities of order 99.90-99.99% to be achievable under realistic experimental conditions.

Neoclassical tearing mode seeding by coupling with infernal modes in low-shear tokamaks

NASA Astrophysics Data System (ADS)

Kleiner, A.; Graves, J. P.; Brunetti, D.; Cooper, W. A.; Halpern, F. D.; Luciani, J.-F.; Lütjens, H.

2016-09-01

A numerical and an analytical study of the triggering of resistive MHD modes in tokamak plasmas with low magnetic shear core is presented. Flat q profiles give rise to fast growing pressure driven MHD modes, such as infernal modes. It has been shown that infernal modes drive fast growing islands on neighbouring rational surfaces. Numerical simulations of such instabilities in a MAST-like configuration are performed with the initial value stability code XTOR-2F in the resistive frame. The evolution of magnetic islands are computed from XTOR-2F simulations and an analytical model is developed based on Rutherford’s theory in combination with a model of resistive infernal modes. The parameter {{Δ }\\prime} is extended from the linear phase to the non-linear phase. Additionally, the destabilising contribution due to a helically perturbed bootstrap current is considered. Comparing the numerical XTOR-2F simulations to the model, we find that coupling has a strong destabilising effect on (neoclassical) tearing modes and is able to seed 2/1 magnetic islands in situations when the standard NTM theory predicts stability.

Analytical solutions by squeezing to the anisotropic Rabi model in the nonperturbative deep-strong-coupling regime

NASA Astrophysics Data System (ADS)

Zhang, Yu-Yu; Chen, Xiang-You

2017-12-01

An unexplored nonperturbative deep strong coupling (npDSC) achieved in superconducting circuits has been studied in the anisotropic Rabi model by the generalized squeezing rotating-wave approximation. Energy levels are evaluated analytically from the reformulated Hamiltonian and agree well with numerical ones in a wide range of coupling strength. Such improvement ascribes to deformation effects in the displaced-squeezed state presented by the squeezed momentum variance, which are omitted in previous displaced states. The atom population dynamics confirms the validity of our approach for the npDSC strength. Our approach offers the possibility to explore interesting phenomena analytically in the npDSC regime in qubit-oscillator experiments.

Analytic regularization of uniform cubic B-spline deformation fields.

PubMed

Shackleford, James A; Yang, Qi; Lourenço, Ana M; Shusharina, Nadya; Kandasamy, Nagarajan; Sharp, Gregory C

2012-01-01

Image registration is inherently ill-posed, and lacks a unique solution. In the context of medical applications, it is desirable to avoid solutions that describe physically unsound deformations within the patient anatomy. Among the accepted methods of regularizing non-rigid image registration to provide solutions applicable to medical practice is the penalty of thin-plate bending energy. In this paper, we develop an exact, analytic method for computing the bending energy of a three-dimensional B-spline deformation field as a quadratic matrix operation on the spline coefficient values. Results presented on ten thoracic case studies indicate the analytic solution is between 61-1371x faster than a numerical central differencing solution.

Analytical estimation show low depth-independent water loss due to vapor flux from deep aquifers

NASA Astrophysics Data System (ADS)

Selker, John S.

2017-06-01

Recent articles have provided estimates of evaporative flux from water tables in deserts that span 5 orders of magnitude. In this paper, we present an analytical calculation that indicates aquifer vapor flux to be limited to 0.01 mm/yr for sites where there is negligible recharge and the water table is well over 20 m below the surface. This value arises from the geothermal gradient, and therefore, is nearly independent of the actual depth of the aquifer. The value is in agreement with several numerical studies, but is 500 times lower than recently reported experimental values, and 100 times larger than an earlier analytical estimate.

Age-of-Air, Tape Recorder, and Vertical Transport Schemes

NASA Technical Reports Server (NTRS)

Lin, S.-J.; Einaudi, Franco (Technical Monitor)

2000-01-01

A numerical-analytic investigation of the impacts of vertical transport schemes on the model simulated age-of-air and the so-called 'tape recorder' will be presented using an idealized 1-D column transport model as well as a more realistic 3-D dynamical model. By comparing to the 'exact' solutions of 'age-of-air' and the 'tape recorder' obtainable in the 1-D setting, useful insight is gained on the impacts of numerical diffusion and dispersion of numerical schemes used in global models. Advantages and disadvantages of Eulerian, semi-Lagrangian, and Lagrangian transport schemes will be discussed. Vertical resolution requirement for numerical schemes as well as observing systems for capturing the fine details of the 'tape recorder' or any upward propagating wave-like structures can potentially be derived from the 1-D analytic model.

An analytical-numerical approach for parameter determination of a five-parameter single-diode model of photovoltaic cells and modules

NASA Astrophysics Data System (ADS)

Hejri, Mohammad; Mokhtari, Hossein; Azizian, Mohammad Reza; Söder, Lennart

2016-04-01

Parameter extraction of the five-parameter single-diode model of solar cells and modules from experimental data is a challenging problem. These parameters are evaluated from a set of nonlinear equations that cannot be solved analytically. On the other hand, a numerical solution of such equations needs a suitable initial guess to converge to a solution. This paper presents a new set of approximate analytical solutions for the parameters of a five-parameter single-diode model of photovoltaic (PV) cells and modules. The proposed solutions provide a good initial point which guarantees numerical analysis convergence. The proposed technique needs only a few data from the PV current-voltage characteristics, i.e. open circuit voltage Voc, short circuit current Isc and maximum power point current and voltage Im; Vm making it a fast and low cost parameter determination technique. The accuracy of the presented theoretical I-V curves is verified by experimental data.

Analytical, Numerical, and Experimental Investigation on a Non-Contact Method for the Measurements of Creep Properties of Ultra-High-Temperature Materials

NASA Technical Reports Server (NTRS)

Lee, Jonghyun; Hyers, Robert W.; Rogers, Jan R.; Rathz, Thomas J.; Choo, Hahn; Liaw, Peter

2006-01-01

Responsive access to space requires re-use of components such as rocket nozzles that operate at extremely high temperatures. For such applications, new ultra-hightemperature materials that can operate over 2,000 C are required. At the temperatures higher than the fifty percent of the melting temperature, the characterization of creep properties is indispensable. Since conventional methods for the measurement of creep is limited below 1,700 C, a new technique that can be applied at higher temperatures is strongly demanded. This research develops a non-contact method for the measurement of creep at the temperatures over 2,300 C. Using the electrostatic levitator in NASA MSFC, a spherical sample was rotated to cause creep deformation by centrifugal acceleration. The deforming sample was captured with a digital camera and analyzed to measure creep deformation. Numerical and analytical analyses have also been conducted to compare the experimental results. Analytical, numerical, and experimental results showed a good agreement with one another.

Onset of detachment in adhesive contact of an elastic half-space and flat-ended punches with non-circular shape: analytic estimates and comparison with numeric analysis

NASA Astrophysics Data System (ADS)

Li, Qiang; Argatov, Ivan; Popov, Valentin L.

2018-04-01

A recent paper by Popov, Pohrt and Li (PPL) in Friction investigated adhesive contacts of flat indenters in unusual shapes using numerical, analytical and experimental methods. Based on that paper, we analyze some special cases for which analytical solutions are known. As in the PPL paper, we consider adhesive contact in the Johnson-Kendall-Roberts approximation. Depending on the energy balance, different upper and lower estimates are obtained in terms of certain integral characteristics of the contact area. The special cases of an elliptical punch as well as a system of two circular punches are considered. Theoretical estimations for the first critical force (force at which the detachment process begins) are confirmed by numerical simulations using the adhesive boundary element method. It is shown that simpler approximations for the pull-off force, based both on the Holm radius of contact and the contact area, substantially overestimate the maximum adhesive force.

Numerical and analytical approaches to an advection-diffusion problem at small Reynolds number and large Péclet number

NASA Astrophysics Data System (ADS)

Fuller, Nathaniel J.; Licata, Nicholas A.

2018-05-01

Obtaining a detailed understanding of the physical interactions between a cell and its environment often requires information about the flow of fluid surrounding the cell. Cells must be able to effectively absorb and discard material in order to survive. Strategies for nutrient acquisition and toxin disposal, which have been evolutionarily selected for their efficacy, should reflect knowledge of the physics underlying this mass transport problem. Motivated by these considerations, in this paper we discuss the results from an undergraduate research project on the advection-diffusion equation at small Reynolds number and large Péclet number. In particular, we consider the problem of mass transport for a Stokesian spherical swimmer. We approach the problem numerically and analytically through a rescaling of the concentration boundary layer. A biophysically motivated first-passage problem for the absorption of material by the swimming cell demonstrates quantitative agreement between the numerical and analytical approaches. We conclude by discussing the connections between our results and the design of smart toxin disposal systems.

Visual Analytics of integrated Data Systems for Space Weather Purposes

NASA Astrophysics Data System (ADS)

Rosa, Reinaldo; Veronese, Thalita; Giovani, Paulo

Analysis of information from multiple data sources obtained through high resolution instrumental measurements has become a fundamental task in all scientific areas. The development of expert methods able to treat such multi-source data systems, with both large variability and measurement extension, is a key for studying complex scientific phenomena, especially those related to systemic analysis in space and environmental sciences. In this talk, we present a time series generalization introducing the concept of generalized numerical lattice, which represents a discrete sequence of temporal measures for a given variable. In this novel representation approach each generalized numerical lattice brings post-analytical data information. We define a generalized numerical lattice as a set of three parameters representing the following data properties: dimensionality, size and post-analytical measure (e.g., the autocorrelation, Hurst exponent, etc)[1]. From this representation generalization, any multi-source database can be reduced to a closed set of classified time series in spatiotemporal generalized dimensions. As a case study, we show a preliminary application in space science data, highlighting the possibility of a real time analysis expert system. In this particular application, we have selected and analyzed, using detrended fluctuation analysis (DFA), several decimetric solar bursts associated to X flare-classes. The association with geomagnetic activity is also reported. DFA method is performed in the framework of a radio burst automatic monitoring system. Our results may characterize the variability pattern evolution, computing the DFA scaling exponent, scanning the time series by a short windowing before the extreme event [2]. For the first time, the application of systematic fluctuation analysis for space weather purposes is presented. The prototype for visual analytics is implemented in a Compute Unified Device Architecture (CUDA) by using the K20 Nvidia graphics processing units (GPUs) to reduce the integrated analysis runtime. [1] Veronese et al. doi: 10.6062/jcis.2009.01.02.0021, 2010. [2] Veronese et al. doi:http://dx.doi.org/10.1016/j.jastp.2010.09.030, 2011.

Understanding the dynamics of superparamagnetic particles under the influence of high field gradient arrays

NASA Astrophysics Data System (ADS)

Barnsley, Lester C.; Carugo, Dario; Aron, Miles; Stride, Eleanor

2017-03-01

The aim of this study was to characterize the behaviour of superparamagnetic particles in magnetic drug targeting (MDT) schemes. A 3-dimensional mathematical model was developed, based on the analytical derivation of the trajectory of a magnetized particle suspended inside a fluid channel carrying laminar flow and in the vicinity of an external source of magnetic force. Semi-analytical expressions to quantify the proportion of captured particles, and their relative accumulation (concentration) as a function of distance along the wall of the channel were also derived. These were expressed in terms of a non-dimensional ratio of the relevant physical and physiological parameters corresponding to a given MDT protocol. The ability of the analytical model to assess magnetic targeting schemes was tested against numerical simulations of particle trajectories. The semi-analytical expressions were found to provide good first-order approximations for the performance of MDT systems in which the magnetic force is relatively constant over a large spatial range. The numerical model was then used to test the suitability of a range of different designs of permanent magnet assemblies for MDT. The results indicated that magnetic arrays that emit a strong magnetic force that varies rapidly over a confined spatial range are the most suitable for concentrating magnetic particles in a localized region. By comparison, commonly used magnet geometries such as button magnets and linear Halbach arrays result in distributions of accumulated particles that are less efficient for delivery. The trajectories predicted by the numerical model were verified experimentally by acoustically focusing magnetic microbeads flowing in a glass capillary channel, and optically tracking their path past a high field gradient Halbach array.

Numerical Nuclear Second Derivatives on a Computing Grid: Enabling and Accelerating Frequency Calculations on Complex Molecular Systems.

PubMed

Yang, Tzuhsiung; Berry, John F

2018-06-04

The computation of nuclear second derivatives of energy, or the nuclear Hessian, is an essential routine in quantum chemical investigations of ground and transition states, thermodynamic calculations, and molecular vibrations. Analytic nuclear Hessian computations require the resolution of costly coupled-perturbed self-consistent field (CP-SCF) equations, while numerical differentiation of analytic first derivatives has an unfavorable 6 N ( N = number of atoms) prefactor. Herein, we present a new method in which grid computing is used to accelerate and/or enable the evaluation of the nuclear Hessian via numerical differentiation: NUMFREQ@Grid. Nuclear Hessians were successfully evaluated by NUMFREQ@Grid at the DFT level as well as using RIJCOSX-ZORA-MP2 or RIJCOSX-ZORA-B2PLYP for a set of linear polyacenes with systematically increasing size. For the larger members of this group, NUMFREQ@Grid was found to outperform the wall clock time of analytic Hessian evaluation; at the MP2 or B2LYP levels, these Hessians cannot even be evaluated analytically. We also evaluated a 156-atom catalytically relevant open-shell transition metal complex and found that NUMFREQ@Grid is faster (7.7 times shorter wall clock time) and less demanding (4.4 times less memory requirement) than an analytic Hessian. Capitalizing on the capabilities of parallel grid computing, NUMFREQ@Grid can outperform analytic methods in terms of wall time, memory requirements, and treatable system size. The NUMFREQ@Grid method presented herein demonstrates how grid computing can be used to facilitate embarrassingly parallel computational procedures and is a pioneer for future implementations.

Modal method for Second Harmonic Generation in nanostructures

NASA Astrophysics Data System (ADS)

Héron, S.; Pardo, F.; Bouchon, P.; Pelouard, J.-L.; Haïdar, R.

2015-05-01

Nanophotonic devices show interesting features for nonlinear response enhancement but numerical tools are mandatory to fully determine their behaviour. To address this need, we present a numerical modal method dedicated to nonlinear optics calculations under the undepleted pump approximation. It is brie y explained in the frame of Second Harmonic Generation for both plane waves and focused beams. The nonlinear behaviour of selected nanostructures is then investigated to show comparison with existing analytical results and study the convergence of the code.

Analysis of a class of boundary value problems depending on left and right Caputo fractional derivatives

NASA Astrophysics Data System (ADS)

Antunes, Pedro R. S.; Ferreira, Rui A. C.

2017-07-01

In this work we study boundary value problems associated to a nonlinear fractional ordinary differential equation involving left and right Caputo derivatives. We discuss the regularity of the solutions of such problems and, in particular, give precise necessary conditions so that the solutions are C1([0, 1]). Taking into account our analytical results, we address the numerical solution of those problems by the augmented -RBF method. Several examples illustrate the good performance of the numerical method.

Modulation instability, Akhmediev breathers, and rogue waves in nonlinear fiber optics

NASA Astrophysics Data System (ADS)

Dudley, John M.; Genty, Go"ry; Dias, Frederic; Kibler, Bertrand; Akhmediev, Nail

2010-02-01

The development of the supercontinuum spectrum in the quasi-CW regime is studied analytically, numerically and experimentally. An interpretation in terms of localized periodic structures known as "Akhmediev Breathers" is proposed. Theory, numerical simulation and experiment are in excellent agreement. We also briefly consider the role of breather collisions in the presence of higher order dispersion and show that they lead to the formation of very large amplitude localized structures that may be analogous to the infamous oceanic rogue waves.

Gravitational Instabilities, Chondrule Formation, and the FU Orionis Phenomenon

NASA Astrophysics Data System (ADS)

Boley, Aaron C.; Durisen, Richard H.

2008-10-01

Using analytic arguments and numerical simulations, we examine whether chondrule formation and the FU Orionis phenomenon can be caused by the burstlike onset of gravitational instabilities (GIs) in dead zones. At least two scenarios for bursting dead zones can work, in principle. If the disk is on the verge of fragmentation, GI activation near r ~ 4-5 AU can produce chondrule-forming shocks, at least under extreme conditions. Mass fluxes are also high enough during the onset of GIs to suggest that the outburst is related to an FU Orionis phenomenon. This situation is demonstrated by numerical simulations. In contrast, as supported by analytic arguments, if the burst takes place close to r ~ 1 AU, then even low pitch angle spiral waves can create chondrule-producing shocks and outbursts. We also study the stability of the massive disks in our simulations against fragmentation and find that although disk evolution is sensitive to changes in opacity, the disks we study do not fragment, even at high resolution and even for extreme assumptions.

Analytical model of cracking due to rebar corrosion expansion in concrete considering the structure internal force

NASA Astrophysics Data System (ADS)

Lin, Xiangyue; Peng, Minli; Lei, Fengming; Tan, Jiangxian; Shi, Huacheng

2017-12-01

Based on the assumptions of uniform corrosion and linear elastic expansion, an analytical model of cracking due to rebar corrosion expansion in concrete was established, which is able to consider the structure internal force. And then, by means of the complex variable function theory and series expansion technology established by Muskhelishvili, the corresponding stress component functions of concrete around the reinforcement were obtained. Also, a comparative analysis was conducted between the numerical simulation model and present model in this paper. The results show that the calculation results of both methods were consistent with each other, and the numerical deviation was less than 10%, proving that the analytical model established in this paper is reliable.

Similarity solution of the Boussinesq equation

NASA Astrophysics Data System (ADS)

Lockington, D. A.; Parlange, J.-Y.; Parlange, M. B.; Selker, J.

Similarity transforms of the Boussinesq equation in a semi-infinite medium are available when the boundary conditions are a power of time. The Boussinesq equation is reduced from a partial differential equation to a boundary-value problem. Chen et al. [Trans Porous Media 1995;18:15-36] use a hodograph method to derive an integral equation formulation of the new differential equation which they solve by numerical iteration. In the present paper, the convergence of their scheme is improved such that numerical iteration can be avoided for all practical purposes. However, a simpler analytical approach is also presented which is based on Shampine's transformation of the boundary value problem to an initial value problem. This analytical approximation is remarkably simple and yet more accurate than the analytical hodograph approximations.

Analytical modeling of large amplitude free vibration of non-uniform beams carrying a both transversely and axially eccentric tip mass

NASA Astrophysics Data System (ADS)

Malaeke, Hasan; Moeenfard, Hamid

2016-03-01

The objective of this paper is to study large amplitude flexural-extensional free vibration of non-uniform cantilever beams carrying a both transversely and axially eccentric tip mass. The effects of variable axial force is also taken into account. Hamilton's principle is utilized to obtain the partial differential equations governing the nonlinear vibration of the system as well as the corresponding boundary conditions. A numerical finite difference scheme is proposed to find the natural frequencies and mode shapes of the system which is validated specifically for a beam with linearly varying cross section. Using a single mode approximation in conjunction with the Lagrange method, the governing equations are reduced to a set of two nonlinear ordinary differential equations in terms of end displacement components of the beam which are coupled due to the presence of the transverse eccentricity. These temporal coupled equations are then solved analytically using the multiple time scales perturbation technique. The obtained analytical results are compared with the numerical ones and excellent agreement is observed. The qualitative and quantitative knowledge resulting from this research is expected to enable the study of the effects of eccentric tip mass and non-uniformity on the large amplitude flexural-extensional vibration of beams for improved dynamic performance.

Reflection and diffraction corrections for nonlinear materials characterization by quasi-static pulse measurement

NASA Astrophysics Data System (ADS)

Nagy, Peter B.; Qu, Jianmin; Jacobs, Laurence J.

2014-02-01

A harmonic acoustic tone burst propagating through an elastic solid with quadratic nonlinearity produces not only a parallel burst of second harmonic but also an often neglected quasi-static pulse associated with the acoustic radiation-induced eigenstrain. Although initial analytical and experimental studies by Yost and Cantrell suggested that the pulse might have a right-angled triangular shape with the peak displacement at the leading edge being proportional to the length of the tone burst, more recent theoretical, analytical, numerical, and experimental studies proved that the pulse has a flat-top shape and the peak displacement is proportional to the propagation length. In this paper, analytical and numerical simulation results are presented to illustrate two types of finite-size effects. First, the finite axial dimension of the specimen cannot be simply accounted for by a linear reflection coefficient that neglects the nonlinear interaction between the combined incident and reflected fields. Second, the quasistatic pulse generated by a transducer of finite aperture suffers more severe divergence than both the fundamental and second harmonic pulses generated by the same transducer. These finite-size effects can make the top of the quasi-static pulse sloped rather than flat and therefore must be taken into consideration in the interpretation of experimental data.

Reflection and diffraction corrections for nonlinear materials characterization by quasi-static pulse measurement

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

Nagy, Peter B.; Qu, Jianmin; Jacobs, Laurence J.

A harmonic acoustic tone burst propagating through an elastic solid with quadratic nonlinearity produces not only a parallel burst of second harmonic but also an often neglected quasi-static pulse associated with the acoustic radiation-induced eigenstrain. Although initial analytical and experimental studies by Yost and Cantrell suggested that the pulse might have a right-angled triangular shape with the peak displacement at the leading edge being proportional to the length of the tone burst, more recent theoretical, analytical, numerical, and experimental studies proved that the pulse has a flat-top shape and the peak displacement is proportional to the propagation length. In thismore » paper, analytical and numerical simulation results are presented to illustrate two types of finite-size effects. First, the finite axial dimension of the specimen cannot be simply accounted for by a linear reflection coefficient that neglects the nonlinear interaction between the combined incident and reflected fields. Second, the quasistatic pulse generated by a transducer of finite aperture suffers more severe divergence than both the fundamental and second harmonic pulses generated by the same transducer. These finite-size effects can make the top of the quasi-static pulse sloped rather than flat and therefore must be taken into consideration in the interpretation of experimental data.« less

Magnetization-induced dynamics of a Josephson junction coupled to a nanomagnet

NASA Astrophysics Data System (ADS)

Ghosh, Roopayan; Maiti, Moitri; Shukrinov, Yury M.; Sengupta, K.

2017-11-01

We study the superconducting current of a Josephson junction (JJ) coupled to an external nanomagnet driven by a time-dependent magnetic field both without and in the presence of an external ac drive. We provide an analytic, albeit perturbative, solution for the Landau-Lifshitz (LL) equations governing the coupled JJ-nanomagnet system in the presence of a magnetic field with arbitrary time dependence oriented along the easy axis of the nanomagnet's magnetization and in the limit of weak dimensionless coupling ɛ0 between the JJ and the nanomagnet. We show the existence of Shapiro-type steps in the I -V characteristics of the JJ subjected to a voltage bias for a constant or periodically varying magnetic field and explore the effect of rotation of the magnetic field and the presence of an external ac drive on these steps. We support our analytic results with exact numerical solution of the LL equations. We also extend our results to dissipative nanomagnets by providing a perturbative solution to the Landau-Lifshitz-Gilbert (LLG) equations for weak dissipation. We study the fate of magnetization-induced Shapiro steps in the presence of dissipation both from our analytical results and via numerical solution of the coupled LLG equations. We discuss experiments which can test our theory.

Analytical, numerical, and experimental investigations on effective mechanical properties and performances of carbon nanotubes and nanotube based nanocomposites with novel three dimensional nanostructures

NASA Astrophysics Data System (ADS)

Askari, Davood

The theoretical objectives and accomplishment of this work are the analytical and numerical investigation of material properties and mechanical behavior of carbon nanotubes (CNTs) and nanotube nanocomposites when they are subjected to various loading conditions. First, the finite element method is employed to investigate numerically the effective Young's modulus and Poisson's ratio of a single-walled CNT. Next, the effects of chirality on the effective Young's modulus and Poisson's ratio are investigated and then variations of their effective coefficient of thermal expansions and effective thermal conductivities are studied for CNTs with different structural configurations. To study the influence of small vacancy defects on mechanical properties of CNTs, finite element analyses are performed and the behavior of CNTs with various structural configurations having different types of vacancy defects is studied. It is frequently reported that nano-materials are excellent candidates as reinforcements in nanocomposites to change or enhance material properties of polymers and their nanocomposites. Second, the inclusion of nano-materials can considerably improve electrical, thermal, and mechanical properties of the bonding agent, i.e., resin. Note that, materials atomic and molecular level do not usually show isotropic behaviour, rather they have orthotropic properties. Therefore, two-phase and three-phase cylindrically orthotropic composite models consisting of different constituents with orthotropic properties are developed and introduced in this work to analytically predict the effective mechanical properties and mechanical behavior of such structures when they are subjected to various external loading conditions. To verify the analytically obtained exact solutions, finite element analyses of identical cylindrical structures are also performed and then results are compared with those obtained analytically, and excellent agreement is achieved. The third part of this dissertation investigates the growth of vertically aligned, long, and high density arrays of CNTs and novel 3-D carbon nanotube nano-forests. A Chemical vapor deposition technique is used to grow radially aligned CNTs on various types of fibrous materials such as silicon carbide, carbon, Kevlar, and glass fibers and clothes that can be used for the fabrication of multifunctional high performing laminated nanocomposite structures. Using the CNTs nano-forest clothes, nanocomposite samples are prepared and tested giving promising results for the improvement of mechanical properties and performance of composites structures.

Transient well flow in vertically heterogeneous aquifers

NASA Astrophysics Data System (ADS)

Hemker, C. J.

1999-11-01

A solution for the general problem of computing well flow in vertically heterogeneous aquifers is found by an integration of both analytical and numerical techniques. The radial component of flow is treated analytically; the drawdown is a continuous function of the distance to the well. The finite-difference technique is used for the vertical flow component only. The aquifer is discretized in the vertical dimension and the heterogeneous aquifer is considered to be a layered (stratified) formation with a finite number of homogeneous sublayers, where each sublayer may have different properties. The transient part of the differential equation is solved with Stehfest's algorithm, a numerical inversion technique of the Laplace transform. The well is of constant discharge and penetrates one or more of the sublayers. The effect of wellbore storage on early drawdown data is taken into account. In this way drawdowns are found for a finite number of sublayers as a continuous function of radial distance to the well and of time since the pumping started. The model is verified by comparing results with published analytical and numerical solutions for well flow in homogeneous and heterogeneous, confined and unconfined aquifers. Instantaneous and delayed drainage of water from above the water table are considered, combined with the effects of partially penetrating and finite-diameter wells. The model is applied to demonstrate that the transient effects of wellbore storage in unconfined aquifers are less pronounced than previous numerical experiments suggest. Other applications of the presented solution technique are given for partially penetrating wells in heterogeneous formations, including a demonstration of the effect of decreasing specific storage values with depth in an otherwise homogeneous aquifer. The presented solution can be a powerful tool for the analysis of drawdown from pumping tests, because hydraulic properties of layered heterogeneous aquifer systems with partially penetrating wells may be estimated without the need to construct transient numerical models. A computer program based on the hybrid analytical-numerical technique is available from the author.

Concepts in Electromagnetic Scattering for Particulate-Systems Characterization

DTIC Science & Technology

2013-04-29

scientists attended and 19 presentations were given. Specific topics included remote sensing, polarimetry , analytic and numeric electromagnetic...presentations were given. Specific topics included remote sensing, polarimetry , ana- lytic and numeric electromagnetic theory, camouflage in nature

Lax-Wendroff and TVD finite volume methods for unidimensional thermomechanical numerical simulations of impacts on elastic-plastic solids

NASA Astrophysics Data System (ADS)

Heuzé, Thomas

2017-10-01

We present in this work two finite volume methods for the simulation of unidimensional impact problems, both for bars and plane waves, on elastic-plastic solid media within the small strain framework. First, an extension of Lax-Wendroff to elastic-plastic constitutive models with linear and nonlinear hardenings is presented. Second, a high order TVD method based on flux-difference splitting [1] and Superbee flux limiter [2] is coupled with an approximate elastic-plastic Riemann solver for nonlinear hardenings, and follows that of Fogarty [3] for linear ones. Thermomechanical coupling is accounted for through dissipation heating and thermal softening, and adiabatic conditions are assumed. This paper essentially focuses on one-dimensional problems since analytical solutions exist or can easily be developed. Accordingly, these two numerical methods are compared to analytical solutions and to the explicit finite element method on test cases involving discontinuous and continuous solutions. This allows to study in more details their respective performance during the loading, unloading and reloading stages. Particular emphasis is also paid to the accuracy of the computed plastic strains, some differences being found according to the numerical method used. Lax-Wendoff two-dimensional discretization of a one-dimensional problem is also appended at the end to demonstrate the extensibility of such numerical scheme to multidimensional problems.

Analytical model for the density distribution in the Io plasma torus

NASA Technical Reports Server (NTRS)

Mei, YI; Thorne, Richard M.; Bagenal, Fran

1995-01-01

An analytical model is developed for the diffusive equilibrium plasma density distribution in the Io plasma torus. The model has been employed successfully to follow the ray path of plasma waves in the multi-ion Jovian magnetosphere; it would also be valuable for other studies of the Io torus that require a smooth and continuous description of the plasma density and its gradients. Validity of the analytical treatment requires that the temperature of thermal electrons be much lower than the ion temperature and that superthermal electrons be much less abundant than the thermal electrons; these two conditions are satisfied in the warm outer region of the Io torus from L = 6 to L = 10. The analytical solutions agree well with exact numerical calculations for the most dense portion of the Io torus within 30 deg of the equator.

Long term evolution of planetary systems with a terrestrial planet and a giant planet.

NASA Astrophysics Data System (ADS)

Georgakarakos, Nikolaos; Dobbs-Dixon, Ian; Way, Michael J.

2017-06-01

We study the long term orbital evolution of a terrestrial planet under the gravitational perturbations of a giant planet. In particular, we are interested in situations where the two planets are in the same plane and are relatively close. We examine both possible configurations: the giant planet orbit being either outside or inside the orbit of the smaller planet. The perturbing potential is expanded to high orders and an analytical solution of the terrestrial planetary orbit is derived. The analytical estimates are then compared against results from the numerical integration of the full equations of motion and we find that the analytical solution works reasonably well. An interesting finding is that the new analytical estimates improve greatly the predictions for the timescales of the orbital evolution of the terrestrial planet compared to an octupole order expansion.

Untangling Slab Dynamics Using 3-D Numerical and Analytical Models

NASA Astrophysics Data System (ADS)

Holt, A. F.; Royden, L.; Becker, T. W.

2016-12-01

Increasingly sophisticated numerical models have enabled us to make significant strides in identifying the key controls on how subducting slabs deform. For example, 3-D models have demonstrated that subducting plate width, and the related strength of toroidal flow around the plate edge, exerts a strong control on both the curvature and the rate of migration of the trench. However, the results of numerical subduction models can be difficult to interpret, and many first order dynamics issues remain at least partially unresolved. Such issues include the dominant controls on trench migration, the interdependence of asthenospheric pressure and slab dynamics, and how nearby slabs influence each other's dynamics. We augment 3-D, dynamically evolving finite element models with simple, analytical force-balance models to distill the physics associated with subduction into more manageable parts. We demonstrate that for single, isolated subducting slabs much of the complexity of our fully numerical models can be encapsulated by simple analytical expressions. Rates of subduction and slab dip correlate strongly with the asthenospheric pressure difference across the subducting slab. For double subduction, an additional slab gives rise to more complex mantle pressure and flow fields, and significantly extends the range of plate kinematics (e.g., convergence rate, trench migration rate) beyond those present in single slab models. Despite these additional complexities, we show that much of the dynamics of such multi-slab systems can be understood using the physics illuminated by our single slab study, and that a force-balance method can be used to relate intra-plate stress to viscous pressure in the asthenosphere and coupling forces at plate boundaries. This method has promise for rapid modeling of large systems of subduction zones on a global scale.

A time-efficient implementation of Extended Kalman Filter for sequential orbit determination and a case study for onboard application

NASA Astrophysics Data System (ADS)

Tang, Jingshi; Wang, Haihong; Chen, Qiuli; Chen, Zhonggui; Zheng, Jinjun; Cheng, Haowen; Liu, Lin

2018-07-01

Onboard orbit determination (OD) is often used in space missions, with which mission support can be partially accomplished autonomously, with less dependency on ground stations. In major Global Navigation Satellite Systems (GNSS), inter-satellite link is also an essential upgrade in the future generations. To serve for autonomous operation, sequential OD method is crucial to provide real-time or near real-time solutions. The Extended Kalman Filter (EKF) is an effective and convenient sequential estimator that is widely used in onboard application. The filter requires the solutions of state transition matrix (STM) and the process noise transition matrix, which are always obtained by numerical integration. However, numerically integrating the differential equations is a CPU intensive process and consumes a large portion of the time in EKF procedures. In this paper, we present an implementation that uses the analytical solutions of these transition matrices to replace the numerical calculations. This analytical implementation is demonstrated and verified using a fictitious constellation based on selected medium Earth orbit (MEO) and inclined Geosynchronous orbit (IGSO) satellites. We show that this implementation performs effectively and converges quickly, steadily and accurately in the presence of considerable errors in the initial values, measurements and force models. The filter is able to converge within 2-4 h of flight time in our simulation. The observation residual is consistent with simulated measurement error, which is about a few centimeters in our scenarios. Compared to results implemented with numerically integrated STM, the analytical implementation shows results with consistent accuracy, while it takes only about half the CPU time to filter a 10-day measurement series. The future possible extensions are also discussed to fit in various missions.

Modeling of Compressible Flow with Friction and Heat Transfer Using the Generalized Fluid System Simulation Program (GFSSP)

NASA Technical Reports Server (NTRS)

Bandyopadhyay, Alak; Majumdar, Alok

2007-01-01

The present paper describes the verification and validation of a quasi one-dimensional pressure based finite volume algorithm, implemented in Generalized Fluid System Simulation Program (GFSSP), for predicting compressible flow with friction, heat transfer and area change. The numerical predictions were compared with two classical solutions of compressible flow, i.e. Fanno and Rayleigh flow. Fanno flow provides an analytical solution of compressible flow in a long slender pipe where incoming subsonic flow can be choked due to friction. On the other hand, Raleigh flow provides analytical solution of frictionless compressible flow with heat transfer where incoming subsonic flow can be choked at the outlet boundary with heat addition to the control volume. Nonuniform grid distribution improves the accuracy of numerical prediction. A benchmark numerical solution of compressible flow in a converging-diverging nozzle with friction and heat transfer has been developed to verify GFSSP's numerical predictions. The numerical predictions compare favorably in all cases.

Mathematical Simulation of Perturbations of Attack Angle of Asymmetric Nanosatellite Passing through Resonance

NASA Astrophysics Data System (ADS)

Lyubimov, V. V.; Kurkina, E. V.

2018-05-01

The authors consider the problem of a dynamic system passing through a low-order resonance, describing an uncontrolled atmospheric descent of an asymmetric nanosatellite in the Earth's atmosphere. The authors perform mathematical and numerical modeling of the motion of the nanosatellite with a small mass-aerodynamic asymmetry relative to the center of mass. The aim of the study is to obtain new reliable approximate analytical estimates of perturbations of the angle of attack of a nanosatellite passing through resonance at angles of attack of not more than 0.5π. By using the stationary phase method, the authors were able to investigate a discontinuous perturbation in the angle of attack of a nanosatellite passing through a resonance with two different nanosatellite designs. Comparison of the results of the numerical modeling and new approximate analytical estimates of the perturbation of the angle of attack confirms the reliability of the said estimates.

Implication of Two-Coupled Differential Van der Pol Duffing Oscillator in Weak Signal Detection

NASA Astrophysics Data System (ADS)

Peng, Hang-hang; Xu, Xue-mei; Yang, Bing-chu; Yin, Lin-zi

2016-04-01

The principle of the Van der Pol Duffing oscillator for state transition and for determining critical value is described, which has been studied to indicate that the application of the Van der Pol Duffing oscillator in weak signal detection is feasible. On the basis of this principle, an improved two-coupled differential Van der Pol Duffing oscillator is proposed which can identify signals under any frequency and ameliorate signal-to-noise ratio (SNR). The analytical methods of the proposed model and the construction of the proposed oscillator are introduced in detail. Numerical experiments on the properties of the proposed oscillator compared with those of the Van der Pol Duffing oscillator are carried out. Our numerical simulations have confirmed the analytical treatment. The results demonstrate that this novel oscillator has better detection performance than the Van der Pol Duffing oscillator.

Development of a computational testbed for numerical simulation of combustion instability

NASA Technical Reports Server (NTRS)

Grenda, Jeffrey; Venkateswaran, Sankaran; Merkle, Charles L.

1993-01-01

A synergistic hierarchy of analytical and computational fluid dynamic techniques is used to analyze three-dimensional combustion instabilities in liquid rocket engines. A mixed finite difference/spectral procedure is employed to study the effects of a distributed vaporization zone on standing and spinning instability modes within the chamber. Droplet atomization and vaporization are treated by a variety of classical models found in the literature. A multi-zone, linearized analytical solution is used to validate the accuracy of the numerical simulations at small amplitudes for a distributed vaporization region. This comparison indicates excellent amplitude and phase agreement under both stable and unstable operating conditions when amplitudes are small and proper grid resolution is used. As amplitudes get larger, expected nonlinearities are observed. The effect of liquid droplet temperature fluctuations was found to be of critical importance in driving the instabilities of the combustion chamber.

An elastic-plastic contact model for line contact structures

NASA Astrophysics Data System (ADS)

Zhu, Haibin; Zhao, Yingtao; He, Zhifeng; Zhang, Ruinan; Ma, Shaopeng

2018-06-01

Although numerical simulation tools are now very powerful, the development of analytical models is very important for the prediction of the mechanical behaviour of line contact structures for deeply understanding contact problems and engineering applications. For the line contact structures widely used in the engineering field, few analytical models are available for predicting the mechanical behaviour when the structures deform plastically, as the classic Hertz's theory would be invalid. Thus, the present study proposed an elastic-plastic model for line contact structures based on the understanding of the yield mechanism. A mathematical expression describing the global relationship between load history and contact width evolution of line contact structures was obtained. The proposed model was verified through an actual line contact test and a corresponding numerical simulation. The results confirmed that this model can be used to accurately predict the elastic-plastic mechanical behaviour of a line contact structure.

Predicting Upscaled Behavior of Aqueous Reactants in Heterogeneous Porous Media

NASA Astrophysics Data System (ADS)

Wright, E. E.; Hansen, S. K.; Bolster, D.; Richter, D. H.; Vesselinov, V. V.

2017-12-01

When modeling reactive transport, reaction rates are often overestimated due to the improper assumption of perfect mixing at the support scale of the transport model. In reality, fronts tend to form between participants in thermodynamically favorable reactions, leading to segregation of reactants into islands or fingers. When such a configuration arises, reactions are limited to the interface between the reactive solutes. Closure methods for estimating control-volume-effective reaction rates in terms of quantities defined at the control volume scale do not presently exist, but their development is crucial for effective field-scale modeling. We attack this problem through a combination of analytical and numerical means. Specifically, we numerically study reactive transport through an ensemble of realizations of two-dimensional heterogeneous porous media. We then employ regression analysis to calibrate an analytically-derived relationship between reaction rate and various dimensionless quantities representing conductivity-field heterogeneity and the respective strengths of diffusion, reaction and advection.

An analytical and numerical study of the Martian planetary boundary layer over slopes.

NASA Technical Reports Server (NTRS)

Blumsack, S. L.; Gierasch, P. J.; Wessel, W. R.

1973-01-01

A one-dimensional model of the Martian planetary boundary layer over sloping terrain is analyzed under a variety of conditions. Analytical results for the steady and diurnal components of the temperature and wind fields are found when a Boussinesq model with a Newtonian cooling law is considered. These results form a basis for understanding the numerical results which include more realistic representations for the heating and parametrizations for the eddy transfer of momentum and heat. The diurnal boundary layer thickness is determined primarily by radiative processes, and the amplitudes of the wind and temperature oscillations are found to depend in an important way on the latitude and slope magnitude. Typically, oscillations in the temperature of plus or minus 15 K and in the upslope wind of plus or minus 25 m/sec are found 1 km above a Martian slope of 0.005.

Collective effects and dynamics of non-adiabatic flame balls

NASA Astrophysics Data System (ADS)

D'Angelo, Yves; Joulin, Guy

2001-03-01

The dynamics of a homogeneous, polydisperse collection of non-adiabatic flame balls (FBs) is investigated by analytical/numerical means. A strongly temperature-dependent Arrhenius reaction rate is assumed, along with a light enough reactant characterized by a markedly less than unity Lewis number (Le). Combining activation-energy asymptotics with a mean-field type of treatment, the analysis yields a nonlinear integro-differential evolution equation (EE) for the FB population. The EE accounts for heat losses inside each FB and unsteadiness around it, as well as for its interactions with the entire FB population, namely mutual heating and faster (Le<1) consumption of the reactant pool. The initial FB number density and size distribution enter the EE explicitly. The latter is studied analytically at early times, then for small total FB number densities; it is subsequently solved numerically, yielding the whole population evolution and its lifetime. Generalizations and open questions relating to `spotty' turbulent combustion are finally evoked.

Simple estimation of linear 1+1 D tsunami run-up

NASA Astrophysics Data System (ADS)

Fuentes, M.; Campos, J. A.; Riquelme, S.

2016-12-01

An analytical expression is derived concerning the linear run-up for any given initial wave generated over a sloping bathymetry. Due to the simplicity of the linear formulation, complex transformations are unnecessay, because the shoreline motion is directly obtained in terms of the initial wave. This analytical result not only supports maximum run-up invariance between linear and non-linear theories, but also the time evolution of shoreline motion and velocity. The results exhibit good agreement with the non-linear theory. The present formulation also allows computing the shoreline motion numerically from a customised initial waveform, including non-smooth functions. This is useful for numerical tests, laboratory experiments or realistic cases in which the initial disturbance might be retrieved from seismic data rather than using a theoretical model. It is also shown that the real case studied is consistent with the field observations.

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

NASA Astrophysics Data System (ADS)

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

1990-04-01

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

Analytical studies on holographic superconductor in the probe limit

NASA Astrophysics Data System (ADS)

Peng, Yan; Liu, Guohua

2017-09-01

We investigate the holographic superconductor model constructed in the (2+1)-dimensional AdS soliton background in the probe limit. With analytical methods, we obtain the formula of critical phase transition points with respect to the scalar mass. We also generalize this formula to higher-dimensional space-time. We mention that these formulas are precise compared to numerical results. In addition, we find a correspondence between the value of the charged scalar field at the tip and the scalar operator at infinity around the phase transition points.

Stability of a Light Sail Riding on a Laser Beam

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

Manchester, Zachary; Loeb, Abraham, E-mail: zmanchester@seas.harvard.edu

2017-03-10

The stability of a light sail riding on a laser beam is analyzed both analytically and numerically. Conical sails on Gaussian beams, which have been studied in the past, are shown to be unstable without active control or additional mechanical modifications. A new architecture for a passively stable sail-and-beam configuration is proposed. The novel spherical shell design for the sail is capable of “beam riding” without the need for active feedback control. Full three-dimensional ray-tracing simulations are performed to verify our analytical results.

Satellite Orbit Under Influence of a Drag - Analytical Approach

NASA Astrophysics Data System (ADS)

Martinović, M. M.; Šegan, S. D.

2017-12-01

The report studies some changes in orbital elements of the artificial satellites of Earth under influence of atmospheric drag. In order to develop possibilities of applying the results in many future cases, an analytical interpretation of the orbital element perturbations is given via useful, but very long expressions. The development is based on the TD88 air density model, recently upgraded with some additional terms. Some expressions and formulae were developed by the computer algebra system Mathematica and tested in some hypothetical cases. The results have good agreement with iterative (numerical) approach.

An analytical approach to top predator interference on the dynamics of a food chain model

NASA Astrophysics Data System (ADS)

Senthamarai, R.; Vijayalakshmi, T.

2018-04-01

In this paper, a nonlinear mathematical model is proposed and analyzed to study of top predator interference on the dynamics of a food chain model. The mathematical model is formulated using the system of non-linear ordinary differential equations. In this model, there are three state dimensionless variables, viz, size of prey population x, size of intermediate predator y and size of top predator population z. The analytical results are compared with the numerical simulation using MATLAB software and satisfactory results are noticed.

General Analytical Procedure for Determination of Acidity Parameters of Weak Acids and Bases

PubMed Central

Pilarski, Bogusław; Kaliszan, Roman; Wyrzykowski, Dariusz; Młodzianowski, Janusz; Balińska, Agata

2015-01-01

The paper presents a new convenient, inexpensive, and reagent-saving general methodology for the determination of pK a values for components of the mixture of diverse chemical classes weak organic acids and bases in water solution, without the need to separate individual analytes. The data obtained from simple pH-metric microtitrations are numerically processed into reliable pK a values for each component of the mixture. Excellent agreement has been obtained between the determined pK a values and the reference literature data for compounds studied. PMID:25692072

General analytical procedure for determination of acidity parameters of weak acids and bases.

PubMed

Pilarski, Bogusław; Kaliszan, Roman; Wyrzykowski, Dariusz; Młodzianowski, Janusz; Balińska, Agata

2015-01-01

The paper presents a new convenient, inexpensive, and reagent-saving general methodology for the determination of pK a values for components of the mixture of diverse chemical classes weak organic acids and bases in water solution, without the need to separate individual analytes. The data obtained from simple pH-metric microtitrations are numerically processed into reliable pK a values for each component of the mixture. Excellent agreement has been obtained between the determined pK a values and the reference literature data for compounds studied.

High Rayleigh number convection in rectangular enclosures with differentially heated vertical walls and aspect ratios between zero and unity

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

Kassemi, S.A.

1988-04-01

High Rayleigh number convection in a rectangular cavity with insulated horizontal surfaces and differentially heated vertical walls was analyzed for an arbitrary aspect ratio smaller than or equal to unity. Unlike previous analytical studies, a systematic method of solution based on linearization technique and analytical iteration procedure was developed to obtain approximate closed-form solutions for a wide range of aspect ratios. The predicted velocity and temperature fields are shown to be in excellent agreement with available experimental and numerical data.

High Rayleigh number convection in rectangular enclosures with differentially heated vertical walls and aspect ratios between zero and unity

NASA Technical Reports Server (NTRS)

Kassemi, Siavash A.

1988-01-01

High Rayleigh number convection in a rectangular cavity with insulated horizontal surfaces and differentially heated vertical walls was analyzed for an arbitrary aspect ratio smaller than or equal to unity. Unlike previous analytical studies, a systematic method of solution based on linearization technique and analytical iteration procedure was developed to obtain approximate closed-form solutions for a wide range of aspect ratios. The predicted velocity and temperature fields are shown to be in excellent agreement with available experimental and numerical data.

Fractional Fourier transform of truncated elliptical Gaussian beams.

PubMed

Du, Xinyue; Zhao, Daomu

2006-12-20

Based on the fact that a hard-edged elliptical aperture can be expanded approximately as a finite sum of complex Gaussian functions in tensor form, an analytical expression for an elliptical Gaussian beam (EGB) truncated by an elliptical aperture and passing through a fractional Fourier transform system is derived by use of vector integration. The approximate analytical results provide more convenience for studying the propagation and transformation of truncated EGBs than the usual way by using the integral formula directly, and the efficiency of numerical calculation is significantly improved.

Studies of the effects of curvature on dilution jet mixing

NASA Technical Reports Server (NTRS)

Holdeman, James D.; Srinivasan, Ram; Reynolds, Robert S.; White, Craig D.

1992-01-01

An analytical program was conducted using both three-dimensional numerical and empirical models to investigate the effects of transition liner curvature on the mixing of jets injected into a confined crossflow. The numerical code is of the TEACH type with hybrid numerics; it uses the power-law and SIMPLER algorithms, an orthogonal curvilinear coordinate system, and an algebraic Reynolds stress turbulence model. From the results of the numerical calculations, an existing empirical model for the temperature field downstream of single and multiple rows of jets injected into a straight rectangular duct was extended to model the effects of curvature. Temperature distributions, calculated with both the numerical and empirical models, are presented to show the effects of radius of curvature and inner and outer wall injection for single and opposed rows of cool dilution jets injected into a hot mainstream flow.

Cylindrical optical resonators: fundamental properties and bio-sensing characteristics

NASA Astrophysics Data System (ADS)

Khozeymeh, Foroogh; Razaghi, Mohammad

2018-04-01

In this paper, detailed theoretical analysis of cylindrical resonators is demonstrated. As illustrated, these kinds of resonators can be used as optical bio-sensing devices. The proposed structure is analyzed using an analytical method based on Lam's approximation. This method is systematic and has simplified the tedious process of whispering-gallery mode (WGM) wavelength analysis in optical cylindrical biosensors. By this method, analysis of higher radial orders of high angular momentum WGMs has been possible. Using closed-form analytical equations, resonance wavelengths of higher radial and angular order WGMs of TE and TM polarization waves are calculated. It is shown that high angular momentum WGMs are more appropriate for bio-sensing applications. Some of the calculations are done using a numerical non-linear Newton method. A perfect match of 99.84% between the analytical and the numerical methods has been achieved. In order to verify the validity of the calculations, Meep simulations based on the finite difference time domain (FDTD) method are performed. In this case, a match of 96.70% between the analytical and FDTD results has been obtained. The analytical predictions are in good agreement with other experimental work (99.99% match). These results validate the proposed analytical modelling for the fast design of optical cylindrical biosensors. It is shown that by extending the proposed two-layer resonator structure analyzing scheme, it is possible to study a three-layer cylindrical resonator structure as well. Moreover, by this method, fast sensitivity optimization in cylindrical resonator-based biosensors has been possible. Sensitivity of the WGM resonances is analyzed as a function of the structural parameters of the cylindrical resonators. Based on the results, fourth radial order WGMs, with a resonator radius of 50 μm, display the most bulk refractive index sensitivity of 41.50 (nm/RIU).

Dispersive models describing mosquitoes’ population dynamics

NASA Astrophysics Data System (ADS)

Yamashita, W. M. S.; Takahashi, L. T.; Chapiro, G.

2016-08-01

The global incidences of dengue and, more recently, zica virus have increased the interest in studying and understanding the mosquito population dynamics. Understanding this dynamics is important for public health in countries where climatic and environmental conditions are favorable for the propagation of these diseases. This work is based on the study of nonlinear mathematical models dealing with the life cycle of the dengue mosquito using partial differential equations. We investigate the existence of traveling wave solutions using semi-analytical method combining dynamical systems techniques and numerical integration. Obtained solutions are validated through numerical simulations using finite difference schemes.

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

NASA Astrophysics Data System (ADS)

Kang, Ja-Young

2003-12-01

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

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

Stershic, Andrew J.; Dolbow, John E.; Moës, Nicolas

The Thick Level-Set (TLS) model is implemented to simulate brittle media undergoing dynamic fragmentation. This non-local model is discretized by the finite element method with damage represented as a continuous field over the domain. A level-set function defines the extent and severity of damage, and a length scale is introduced to limit the damage gradient. Numerical studies in one dimension demonstrate that the proposed method reproduces the rate-dependent energy dissipation and fragment length observations from analytical, numerical, and experimental approaches. In conclusion, additional studies emphasize the importance of appropriate bulk constitutive models and sufficient spatial resolution of the length scale.

Development of an analytical-numerical model to predict radiant emission or absorption

NASA Technical Reports Server (NTRS)

Wallace, Tim L.

1994-01-01

The development of an analytical-numerical model to predict radiant emission or absorption is discussed. A voigt profile is assumed to predict the spectral qualities of a singlet atomic transition line for atomic species of interest to the OPAD program. The present state of this model is described in each progress report required under contract. Model and code development is guided by experimental data where available. When completed, the model will be used to provide estimates of specie erosion rates from spectral data collected from rocket exhaust plumes or other sources.

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

NASA Technical Reports Server (NTRS)

Brown, R. L.

1979-01-01

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

Satellite recovery - Attitude dynamics of the targets

NASA Technical Reports Server (NTRS)

Cochran, J. E., Jr.; Lahr, B. S.

1986-01-01

The problems of categorizing and modeling the attitude dynamics of uncontrolled artificial earth satellites which may be targets in recovery attempts are addressed. Methods of classification presented are based on satellite rotational kinetic energy, rotational angular momentum and orbit and on the type of control present prior to the benign failure of the control system. The use of approximate analytical solutions and 'exact' numerical solutions to the equations governing satellite attitude motions to predict uncontrolled attitude motion is considered. Analytical and numerical results are presented for the evolution of satellite attitude motions after active control termination.

Spherical Pendulum Small Oscillations for Slewing Crane Motion

PubMed Central

Perig, Alexander V.; Stadnik, Alexander N.; Deriglazov, Alexander I.

2014-01-01

The present paper focuses on the Lagrange mechanics-based description of small oscillations of a spherical pendulum with a uniformly rotating suspension center. The analytical solution of the natural frequencies' problem has been derived for the case of uniform rotation of a crane boom. The payload paths have been found in the inertial reference frame fixed on earth and in the noninertial reference frame, which is connected with the rotating crane boom. The numerical amplitude-frequency characteristics of the relative payload motion have been found. The mechanical interpretation of the terms in Lagrange equations has been outlined. The analytical expression and numerical estimation for cable tension force have been proposed. The numerical computational results, which correlate very accurately with the experimental observations, have been shown. PMID:24526891

Numerical investigation of finite-volume effects for the HVP

NASA Astrophysics Data System (ADS)

Boyle, Peter; Gülpers, Vera; Harrison, James; Jüttner, Andreas; Portelli, Antonin; Sachrajda, Christopher

2018-03-01

It is important to correct for finite-volume (FV) effects in the presence of QED, since these effects are typically large due to the long range of the electromagnetic interaction. We recently made the first lattice calculation of electromagnetic corrections to the hadronic vacuum polarisation (HVP). For the HVP, an analytical derivation of FV corrections involves a two-loop calculation which has not yet been carried out. We instead calculate the universal FV corrections numerically, using lattice scalar QED as an effective theory. We show that this method gives agreement with known analytical results for scalar mass FV effects, before applying it to calculate FV corrections for the HVP. This method for numerical calculation of FV effects is also widely applicable to quantities beyond the HVP.

The liquid fuel jet in subsonic crossflow

NASA Technical Reports Server (NTRS)

Nguyen, T. T.; Karagozian, A. R.

1990-01-01

An analytical/numerical model is described which predicts the behavior of nonreacting and reacting liquid jets injected transversely into subsonic cross flow. The compressible flowfield about the elliptical jet cross section is solved at various locations along the jet trajectory by analytical means for free-stream local Mach number perpendicular to jet cross section smaller than 0.3 and by numerical means for free-stream local Mach number perpendicular to jet cross section in the range 0.3-1.0. External and internal boundary layers along the jet cross section are solved by integral and numerical methods, and the mass losses due to boundary layer shedding, evaporation, and combustion are calculated and incorporated into the trajectory calculation. Comparison of predicted trajectories is made with limited experimental observations.

Analytical solutions for coagulation and condensation kinetics of composite particles

NASA Astrophysics Data System (ADS)

Piskunov, Vladimir N.

2013-04-01

The processes of composite particles formation consisting of a mixture of different materials are essential for many practical problems: for analysis of the consequences of accidental releases in atmosphere; for simulation of precipitation formation in clouds; for description of multi-phase processes in chemical reactors and industrial facilities. Computer codes developed for numerical simulation of these processes require optimization of computational methods and verification of numerical programs. Kinetic equations of composite particle formation are given in this work in a concise form (impurity integrated). Coagulation, condensation and external sources associated with nucleation are taken into account. Analytical solutions were obtained in a number of model cases. The general laws for fraction redistribution of impurities were defined. The results can be applied to develop numerical algorithms considerably reducing the simulation effort, as well as to verify the numerical programs for calculation of the formation kinetics of composite particles in the problems of practical importance.

Stability of the fluid interface in a Hele-Shaw cell subjected to horizontal vibrations

NASA Astrophysics Data System (ADS)

Lyubimova, T. P.; Lyubimov, D. V.; Sadilov, E. S.; Popov, D. M.

2017-07-01

The stability of the horizontal interface of two immiscible viscous fluids in a Hele-Shaw cell subject to gravity and horizontal vibrations is studied. The problem is reduced to the generalized Hill equation, which is solved analytically by the multiple scale method and numerically. The long-wave instability, the resonance (parametric resonance) excitation of waves at finite frequencies of vibrations (for the first three resonances), and the limit of high-frequency vibrations are studied analytically under the assumption of small amplitudes of vibrations and small viscosity. For finite amplitudes of vibrations, finite wave numbers, and finite viscosity, the study is performed numerically. The influence of the specific natural control parameters and physical parameters of the system on its instability threshold is discussed. The results provide extension to other results [J. Bouchgl, S. Aniss, and M. Souhar, Phys. Rev. E 88, 023027 (2013), 10.1103/PhysRevE.88.023027], where the authors considered a similar problem but took into account viscosity in the basic state and did not consider it in the equations for perturbations.

Semi-analytical and Numerical Studies on the Flattened Brazilian Splitting Test Used for Measuring the Indirect Tensile Strength of Rocks

NASA Astrophysics Data System (ADS)

Huang, Y. G.; Wang, L. G.; Lu, Y. L.; Chen, J. R.; Zhang, J. H.

2015-09-01

Based on the two-dimensional elasticity theory, this study established a mechanical model under chordally opposing distributed compressive loads, in order to perfect the theoretical foundation of the flattened Brazilian splitting test used for measuring the indirect tensile strength of rocks. The stress superposition method was used to obtain the approximate analytic solutions of stress components inside the flattened Brazilian disk. These analytic solutions were then verified through a comparison with the numerical results of the finite element method (FEM). Based on the theoretical derivation, this research carried out a contrastive study on the effect of the flattened loading angles on the stress value and stress concentration degree inside the disk. The results showed that the stress concentration degree near the loading point and the ratio of compressive/tensile stress inside the disk dramatically decreased as the flattened loading angle increased, avoiding the crushing failure near-loading point of Brazilian disk specimens. However, only the tensile stress value and the tensile region were slightly reduced with the increase of the flattened loading angle. Furthermore, this study found that the optimal flattened loading angle was 20°-30°; flattened load angles that were too large or too small made it difficult to guarantee the central tensile splitting failure principle of the Brazilian splitting test. According to the Griffith strength failure criterion, the calculative formula of the indirect tensile strength of rocks was derived theoretically. This study obtained a theoretical indirect tensile strength that closely coincided with existing and experimental results. Finally, this paper simulated the fracture evolution process of rocks under different loading angles through the use of the finite element numerical software ANSYS. The modeling results showed that the Flattened Brazilian Splitting Test using the optimal loading angle could guarantee the tensile splitting failure initiated by a central crack.

Dynamic behaviour of a planar micro-beam loaded by a fluid-gap: Analytical and numerical approach in a high frequency range, benchmark solutions

NASA Astrophysics Data System (ADS)

Novak, A.; Honzik, P.; Bruneau, M.

2017-08-01

Miniaturized vibrating MEMS devices, active (receivers or emitters) or passive devices, and their use for either new applications (hearing, meta-materials, consumer devices,…) or metrological purposes under non-standard conditions, are involved today in several acoustic domains. More in-depth characterisation than the classical ones available until now are needed. In this context, the paper presents analytical and numerical approaches for describing the behaviour of three kinds of planar micro-beams of rectangular shape (suspended rigid or clamped elastic planar beam) loaded by a backing cavity or a fluid-gap, surrounded by very thin slits, and excited by an incident acoustic field. The analytical approach accounts for the coupling between the vibrating structure and the acoustic field in the backing cavity, the thermal and viscous diffusion processes in the boundary layers in the slits and the cavity, the modal behaviour for the vibrating structure, and the non-uniformity of the acoustic field in the backing cavity which is modelled in using an integral formulation with a suitable Green's function. Benchmark solutions are proposed in terms of beam motion (from which the sensitivity, input impedance, and pressure transfer function can be calculated). A numerical implementation (FEM) is handled against which the analytical results are tested.

Estimating statistical isotropy violation in CMB due to non-circular beam and complex scan in minutes

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

Pant, Nidhi; Das, Santanu; Mitra, Sanjit

Mild, unavoidable deviations from circular-symmetry of instrumental beams along with scan strategy can give rise to measurable Statistical Isotropy (SI) violation in Cosmic Microwave Background (CMB) experiments. If not accounted properly, this spurious signal can complicate the extraction of other SI violation signals (if any) in the data. However, estimation of this effect through exact numerical simulation is computationally intensive and time consuming. A generalized analytical formalism not only provides a quick way of estimating this signal, but also gives a detailed understanding connecting the leading beam anisotropy components to a measurable BipoSH characterisation of SI violation. In this paper,more » we provide an approximate generic analytical method for estimating the SI violation generated due to a non-circular (NC) beam and arbitrary scan strategy, in terms of the Bipolar Spherical Harmonic (BipoSH) spectra. Our analytical method can predict almost all the features introduced by a NC beam in a complex scan and thus reduces the need for extensive numerical simulation worth tens of thousands of CPU hours into minutes long calculations. As an illustrative example, we use WMAP beams and scanning strategy to demonstrate the easability, usability and efficiency of our method. We test all our analytical results against that from exact numerical simulations.« less

Effective description of domain wall strings

NASA Astrophysics Data System (ADS)

Rodrigues, Davi R.; Abanov, Ar.; Sinova, J.; Everschor-Sitte, K.

2018-04-01

The analysis of domain wall dynamics is often simplified to one-dimensional physics. For domain walls in thin films, more realistic approaches require the description as two-dimensional objects. This includes the study of vortices and curvatures along the domain walls as well as the influence of boundary effects. Here we provide a theory in terms of soft modes that allows us to analytically study the physics of extended domain walls and their stability. By considering irregularly shaped skyrmions as closed domain walls, we analyze their plasticity and compare their dynamics with those of circular skyrmions. Our theory directly provides an analytical description of the excitation modes of magnetic skyrmions, previously accessible only through sophisticated micromagnetic numerical calculations and spectral analysis. These analytical expressions provide the scaling behavior of the different physics on parameters that experiments can test.

Moisture Content and Migration Dynamics in Unsaturated Porous Media

NASA Technical Reports Server (NTRS)

Homsy, G. M.

1993-01-01

Fundamental studies of fluid mechanics and transport in partially saturated soils are presented. Solution of transient diffusion problems in support of the development of probes for the in-situ measurement of moisture content is given. Numerical and analytical methods are used to study the fundamental problem of meniscus and saturation front propagation in geometric models of porous media.

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

Vincenti, H.; Vay, J. -L.

Due to discretization effects and truncation to finite domains, many electromagnetic simulations present non-physical modifications of Maxwell's equations in space that may generate spurious signals affecting the overall accuracy of the result. Such modifications for instance occur when Perfectly Matched Layers (PMLs) are used at simulation domain boundaries to simulate open media. Another example is the use of arbitrary order Maxwell solver with domain decomposition technique that may under some condition involve stencil truncations at subdomain boundaries, resulting in small spurious errors that do eventually build up. In each case, a careful evaluation of the characteristics and magnitude of themore » errors resulting from these approximations, and their impact at any frequency and angle, requires detailed analytical and numerical studies. To this end, we present a general analytical approach that enables the evaluation of numerical discretization errors of fully three-dimensional arbitrary order finite-difference Maxwell solver, with arbitrary modification of the local stencil in the simulation domain. The analytical model is validated against simulations of domain decomposition technique and PMLs, when these are used with very high-order Maxwell solver, as well as in the infinite order limit of pseudo-spectral solvers. Results confirm that the new analytical approach enables exact predictions in each case. It also confirms that the domain decomposition technique can be used with very high-order Maxwell solver and a reasonably low number of guard cells with negligible effects on the whole accuracy of the simulation.« less

The serpentine optical waveguide: engineering the dispersion relations and the stopped light points.

PubMed

Scheuer, Jacob; Weiss, Ori

2011-06-06

We present a study a new type of optical slow-light structure comprising a serpentine shaped waveguide were the loops are coupled. The dispersion relation, group velocity and GVD are studied analytically using a transfer matrix method and numerically using finite difference time domain simulations. The structure exhibits zero group velocity points at the ends of the Brillouin zone, but also within the zone. The position of mid-zone zero group velocity point can be tuned by modifying the coupling coefficient between adjacent loops. Closed-form analytic expressions for the dispersion relations, group velocity and the mid-zone zero v(g) points are found and presented.

Entanglement entropy in the long-range Kitaev chain

NASA Astrophysics Data System (ADS)

Ares, Filiberto; Esteve, José G.; Falceto, Fernando; de Queiroz, Amilcar R.

2018-06-01

In this paper we complete the study on the asymptotic behavior of the entanglement entropy for Kitaev chains with long-range pairing. We discover that when the couplings decay with the distance with a critical exponent new properties for the asymptotic growth of the entropy appear. The coefficient of the leading term is not universal any more and the connection with conformal field theories is lost. We perform a numerical and analytical approach to the problem showing a perfect agreement. In order to carry out the analytical study, a technique for computing the asymptotic behavior of block Toeplitz determinants with discontinuous symbols has been developed.

Stabilization of the Rayleigh-Taylor instability in quantum magnetized plasmas

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

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

2012-07-15

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

Approximate analytic method for high-apogee twelve-hour orbits of artificial Earth's satellites

NASA Astrophysics Data System (ADS)

Vashkovyaka, M. A.; Zaslavskii, G. S.

2016-09-01

We propose an approach to the study of the evolution of high-apogee twelve-hour orbits of artificial Earth's satellites. We describe parameters of the motion model used for the artificial Earth's satellite such that the principal gravitational perturbations of the Moon and Sun, nonsphericity of the Earth, and perturbations from the light pressure force are approximately taken into account. To solve the system of averaged equations describing the evolution of the orbit parameters of an artificial satellite, we use both numeric and analytic methods. To select initial parameters of the twelve-hour orbit, we assume that the path of the satellite along the surface of the Earth is stable. Results obtained by the analytic method and by the numerical integration of the evolving system are compared. For intervals of several years, we obtain estimates of oscillation periods and amplitudes for orbital elements. To verify the results and estimate the precision of the method, we use the numerical integration of rigorous (not averaged) equations of motion of the artificial satellite: they take into account forces acting on the satellite substantially more completely and precisely. The described method can be applied not only to the investigation of orbit evolutions of artificial satellites of the Earth; it can be applied to the investigation of the orbit evolution for other planets of the Solar system provided that the corresponding research problem will arise in the future and the considered special class of resonance orbits of satellites will be used for that purpose.

Soliton polarization rotation in fiber lasers

NASA Astrophysics Data System (ADS)

Afanasjev, V. V.

1995-02-01

I have found the approximate analytical solution in explicit form for a vector soliton with an arbitrary component ratio. My solution describes the dependence of soliton intensity on polarization angle and also nonlinear polarization rotation. The analytical results agree well with the numerical simulations.

A free energy-based surface tension force model for simulation of multiphase flows by level-set method

NASA Astrophysics Data System (ADS)

Yuan, H. Z.; Chen, Z.; Shu, C.; Wang, Y.; Niu, X. D.; Shu, S.

2017-09-01

In this paper, a free energy-based surface tension force (FESF) model is presented for accurately resolving the surface tension force in numerical simulation of multiphase flows by the level set method. By using the analytical form of order parameter along the normal direction to the interface in the phase-field method and the free energy principle, FESF model offers an explicit and analytical formulation for the surface tension force. The only variable in this formulation is the normal distance to the interface, which can be substituted by the distance function solved by the level set method. On one hand, as compared to conventional continuum surface force (CSF) model in the level set method, FESF model introduces no regularized delta function, due to which it suffers less from numerical diffusions and performs better in mass conservation. On the other hand, as compared to the phase field surface tension force (PFSF) model, the evaluation of surface tension force in FESF model is based on an analytical approach rather than numerical approximations of spatial derivatives. Therefore, better numerical stability and higher accuracy can be expected. Various numerical examples are tested to validate the robustness of the proposed FESF model. It turns out that FESF model performs better than CSF model and PFSF model in terms of accuracy, stability, convergence speed and mass conservation. It is also shown in numerical tests that FESF model can effectively simulate problems with high density/viscosity ratio, high Reynolds number and severe topological interfacial changes.

Theoretical and Numerical Investigations on Shallow Tunnelling in Unsaturated Soils

NASA Astrophysics Data System (ADS)

Soranzo, Enrico; Wu, Wei

2013-04-01

Excavation of shallow tunnels with the New Austrian Tunnelling Method (NATM) requires proper assessing of the tunnel face stability, to enable an open-face excavation, and the estimation of the correspondent surface settlements. Soils in a partially saturated condition exhibit a higher cohesion than in a fully saturated state, which can be taken into account when assessing the stability of the tunnel face. For the assessment of the face support pressure, different methods are used in engineering practice, varying from simple empirical and analytical formulations to advanced finite element analysis. Such procedures can be modified to account for the unsaturated state of soils. In this study a method is presented to incorporate the effect of partial saturation in the numerical analysis. The results are then compared with a simple analytical formulation derived from parametric studies. As to the numerical analysis, the variation of cohesion and of Young's modulus with saturation can be considered when the water table lies below the tunnel in a soil exhibiting a certain capillary rise, so that the tunnel is driven in a partially saturated layer. The linear elastic model with Mohr-Coulomb failure criterion can be extended to partially saturated states and calibrated with triaxial tests on unsaturated. In order to model both positive and negative pore water pressure (suction), Bishop's effective stress is incorporated into Mohr-Coulomb's failure criterion. The effective stress parameter in Bishop's formulation is related to the degree of saturation as suggested by Fredlund. If a linear suction distribution is assumed, the degree of saturation can be calculated from the Soil Water Characteristic Curve (SWCC). Expressions exist that relate the Young's modulus of unsaturated soils to the net mean stress and the matric suction. The results of the numerical computation can be compared to Vermeer & Ruse's closed-form formula that expresses the limit support pressure of the tunnel face. The expression is derived from parametric studies and predicts stability of the tunnel face when negative values are returned, suggesting that open-face tunnelling can be performed. The formula can be modified to account for the variation of cohesion along the tunnel face. The results obtained from both the numerical analysis and the analytical formulation are well in agreement and show that the stability of the tunnel face can greatly benefit from the enhanced cohesion of partially saturated soils.

Analytical Solution of the Radiative Transfer Equation in a Thin Dusty Circumstellar Shell

NASA Astrophysics Data System (ADS)

Cruzalèbes, P.; Sacuto, S.

The radiative transfer equation can be solved analytically for optically thin shells. The solution leads to a semi-analytical expression of the visibility function, which can be compared to the numerical solution given by the DUSTY code. Best-fit model parameters are given using real measurements of ISO fluxes, ISI and VLTI-MIDI visibilities for 3 late-type stars.

Combined Numerical/Analytical Perturbation Solutions of the Navier-Stokes Equations for Aerodynamic Ejector/Mixer Nozzle Flows

NASA Technical Reports Server (NTRS)

DeChant, Lawrence Justin

1998-01-01

In spite of rapid advances in both scalar and parallel computational tools, the large number of variables involved in both design and inverse problems make the use of sophisticated fluid flow models impractical, With this restriction, it is concluded that an important family of methods for mathematical/computational development are reduced or approximate fluid flow models. In this study a combined perturbation/numerical modeling methodology is developed which provides a rigorously derived family of solutions. The mathematical model is computationally more efficient than classical boundary layer but provides important two-dimensional information not available using quasi-1-d approaches. An additional strength of the current methodology is its ability to locally predict static pressure fields in a manner analogous to more sophisticated parabolized Navier Stokes (PNS) formulations. To resolve singular behavior, the model utilizes classical analytical solution techniques. Hence, analytical methods have been combined with efficient numerical methods to yield an efficient hybrid fluid flow model. In particular, the main objective of this research has been to develop a system of analytical and numerical ejector/mixer nozzle models, which require minimal empirical input. A computer code, DREA Differential Reduced Ejector/mixer Analysis has been developed with the ability to run sufficiently fast so that it may be used either as a subroutine or called by an design optimization routine. Models are of direct use to the High Speed Civil Transport Program (a joint government/industry project seeking to develop an economically.viable U.S. commercial supersonic transport vehicle) and are currently being adopted by both NASA and industry. Experimental validation of these models is provided by comparison to results obtained from open literature and Limited Exclusive Right Distribution (LERD) sources, as well as dedicated experiments performed at Texas A&M. These experiments have been performed using a hydraulic/gas flow analog. Results of comparisons of DREA computations with experimental data, which include entrainment, thrust, and local profile information, are overall good. Computational time studies indicate that DREA provides considerably more information at a lower computational cost than contemporary ejector nozzle design models. Finally. physical limitations of the method, deviations from experimental data, potential improvements and alternative formulations are described. This report represents closure to the NASA Graduate Researchers Program. Versions of the DREA code and a user's guide may be obtained from the NASA Lewis Research Center.

Petermann I and II spot size: Accurate semi analytical description involving Nelder-Mead method of nonlinear unconstrained optimization and three parameter fundamental modal field

NASA Astrophysics Data System (ADS)

Roy Choudhury, Raja; Roy Choudhury, Arundhati; Kanti Ghose, Mrinal

2013-01-01

A semi-analytical model with three optimizing parameters and a novel non-Gaussian function as the fundamental modal field solution has been proposed to arrive at an accurate solution to predict various propagation parameters of graded-index fibers with less computational burden than numerical methods. In our semi analytical formulation the optimization of core parameter U which is usually uncertain, noisy or even discontinuous, is being calculated by Nelder-Mead method of nonlinear unconstrained minimizations as it is an efficient and compact direct search method and does not need any derivative information. Three optimizing parameters are included in the formulation of fundamental modal field of an optical fiber to make it more flexible and accurate than other available approximations. Employing variational technique, Petermann I and II spot sizes have been evaluated for triangular and trapezoidal-index fibers with the proposed fundamental modal field. It has been demonstrated that, the results of the proposed solution identically match with the numerical results over a wide range of normalized frequencies. This approximation can also be used in the study of doped and nonlinear fiber amplifier.

Modeling and control of flexible space platforms with articulated payloads

NASA Technical Reports Server (NTRS)

Graves, Philip C.; Joshi, Suresh M.

1989-01-01

The first steps in developing a methodology for spacecraft control-structure interaction (CSI) optimization are identification and classification of anticipated missions, and the development of tractable mathematical models in each mission class. A mathematical model of a generic large flexible space platform (LFSP) with multiple independently pointed rigid payloads is considered. The objective is not to develop a general purpose numerical simulation, but rather to develop an analytically tractable mathematical model of such composite systems. The equations of motion for a single payload case are derived, and are linearized about zero steady-state. The resulting model is then extended to include multiple rigid payloads, yielding the desired analytical form. The mathematical models developed clearly show the internal inertial/elastic couplings, and are therefore suitable for analytical and numerical studies. A simple decentralized control law is proposed for fine pointing the payloads and LFSP attitude control, and simulation results are presented for an example problem. The decentralized controller is shown to be adequate for the example problem chosen, but does not, in general, guarantee stability. A centralized dissipative controller is then proposed, requiring a symmetric form of the composite system equations. Such a controller guarantees robust closed loop stability despite unmodeled elastic dynamics and parameter uncertainties.

Computational Simulation of Acoustic Modes in Rocket Combustors

NASA Technical Reports Server (NTRS)

Harper, Brent (Technical Monitor); Merkle, C. L.; Sankaran, V.; Ellis, M.

2004-01-01

A combination of computational fluid dynamic analysis and analytical solutions is being used to characterize the dominant modes in liquid rocket engines in conjunction with laboratory experiments. The analytical solutions are based on simplified geometries and flow conditions and are used for careful validation of the numerical formulation. The validated computational model is then extended to realistic geometries and flow conditions to test the effects of various parameters on chamber modes, to guide and interpret companion laboratory experiments in simplified combustors, and to scale the measurements to engine operating conditions. In turn, the experiments are used to validate and improve the model. The present paper gives an overview of the numerical and analytical techniques along with comparisons illustrating the accuracy of the computations as a function of grid resolution. A representative parametric study of the effect of combustor mean flow Mach number and combustor aspect ratio on the chamber modes is then presented for both transverse and longitudinal modes. The results show that higher mean flow Mach numbers drive the modes to lower frequencies. Estimates of transverse wave mechanics in a high aspect ratio combustor are then contrasted with longitudinal modes in a long and narrow combustor to provide understanding of potential experimental simulations.

Exergy optimization in a steady moving bed heat exchanger.

PubMed

Soria-Verdugo, A; Almendros-Ibáñez, J A; Ruiz-Rivas, U; Santana, D

2009-04-01

This work provides an energy and exergy optimization analysis of a moving bed heat exchanger (MBHE). The exchanger is studied as a cross-flow heat exchanger where one of the phases is a moving granular medium. The optimal MBHE dimensions and the optimal particle diameter are obtained for a range of incoming fluid flow rates. The analyses are carried out over operation data of the exchanger obtained in two ways: a numerical simulation of the steady-state problem and an analytical solution of the simplified equations, neglecting the conduction terms. The numerical simulation considers, for the solid, the convection heat transfer to the fluid and the diffusion term in both directions, and for the fluid only the convection heat transfer to the solid. The results are compared with a well-known analytical solution (neglecting conduction effects) for the temperature distribution in the exchanger. Next, the analytical solution is used to derive an expression for the exergy destruction. The optimal length of the MBHE depends mainly on the flow rate and does not depend on particle diameter unless they become very small (thus increasing sharply the pressure drop). The exergy optimal length is always smaller than the thermal one, although the difference is itself small.

Numerical modeling of simultaneous tracer release and piscicide treatment for invasive species control in the Chicago Sanitary and Ship Canal, Chicago, Illinois

USGS Publications Warehouse

Zhu, Zhenduo; Motta, Davide; Jackson, P. Ryan; Garcia, Marcelo H.

2017-01-01

In December 2009, during a piscicide treatment targeting the invasive Asian carp in the Chicago Sanitary and Ship Canal, Rhodamine WT dye was released to track and document the transport and dispersion of the piscicide. In this study, two modeling approaches are presented to reproduce the advection and dispersion of the dye tracer (and piscicide), a one-dimensional analytical solution and a three-dimensional numerical model. The two approaches were compared with field measurements of concentration and their applicability is discussed. Acoustic Doppler current profiler measurements were used to estimate the longitudinal dispersion coefficients at ten cross sections, which were taken as reference for calibrating the longitudinal dispersion coefficient in the one-dimensional analytical solution. While the analytical solution is fast, relatively simple, and can fairly accurately predict the core of the observed concentration time series at points downstream, it does not capture the tail of the breakthrough curves. These tails are well reproduced by the three-dimensional model, because it accounts for the effects of dead zones and a power plant which withdraws nearly 80 % of the water from the canal for cooling purposes before returning it back to the canal.

The H I-to-H2 Transition in a Turbulent Medium

NASA Astrophysics Data System (ADS)

Bialy, Shmuel; Burkhart, Blakesley; Sternberg, Amiel

2017-07-01

We study the effect of density fluctuations induced by turbulence on the H I/H2 structure in photodissociation regions (PDRs) both analytically and numerically. We perform magnetohydrodynamic numerical simulations for both subsonic and supersonic turbulent gas and chemical H I/H2 balance calculations. We derive atomic-to-molecular density profiles and the H I column density probability density function (PDF) assuming chemical equilibrium. We find that, while the H I/H2 density profiles are strongly perturbed in turbulent gas, the mean H I column density is well approximated by the uniform-density analytic formula of Sternberg et al. The PDF width depends on (a) the radiation intensity-to-mean density ratio, (b) the sonic Mach number, and (c) the turbulence decorrelation scale, or driving scale. We derive an analytic model for the H I PDF and demonstrate how our model, combined with 21 cm observations, can be used to constrain the Mach number and driving scale of turbulent gas. As an example, we apply our model to observations of H I in the Perseus molecular cloud. We show that a narrow observed H I PDF may imply small-scale decorrelation, pointing to the potential importance of subcloud-scale turbulence driving.

Benchmarking fully analytic DFT force fields for vibrational spectroscopy: A study on halogenated compounds

NASA Astrophysics Data System (ADS)

Pietropolli Charmet, Andrea; Cornaton, Yann

2018-05-01

This work presents an investigation of the theoretical predictions yielded by anharmonic force fields having the cubic and quartic force constants are computed analytically by means of density functional theory (DFT) using the recursive scheme developed by M. Ringholm et al. (J. Comput. Chem. 35 (2014) 622). Different functionals (namely B3LYP, PBE, PBE0 and PW86x) and basis sets were used for calculating the anharmonic vibrational spectra of two halomethanes. The benchmark analysis carried out demonstrates the reliability and overall good performances offered by hybrid approaches, where the harmonic data obtained at the coupled cluster with single and double excitations level of theory augmented by a perturbational estimate of the effects of connected triple excitations, CCSD(T), are combined with the fully analytic higher order force constants yielded by DFT functionals. These methods lead to reliable and computationally affordable calculations of anharmonic vibrational spectra with an accuracy comparable to that yielded by hybrid force fields having the anharmonic force fields computed at second order Møller-Plesset perturbation theory (MP2) level of theory using numerical differentiation but without the corresponding potential issues related to computational costs and numerical errors.

An analytical model of prominence dynamics

NASA Astrophysics Data System (ADS)

Routh, Swati; Saha, Snehanshu; Bhat, Atul; Sundar, M. N.

2018-01-01

Solar prominences are magnetic structures incarcerating cool and dense gas in an otherwise hot solar corona. Prominences can be categorized as quiescent and active. Their origin and the presence of cool gas (∼104 K) within the hot (∼106K) solar corona remains poorly understood. The structure and dynamics of solar prominences was investigated in a large number of observational and theoretical (both analytical and numerical) studies. In this paper, an analytic model of quiescent solar prominence is developed and used to demonstrate that the prominence velocity increases exponentially, which means that some gas falls downward towards the solar surface, and that Alfvén waves are naturally present in the solar prominences. These theoretical predictions are consistent with the current observational data of solar quiescent prominences.

Analytical solution of groundwater flow in a sloping aquifer with stream-aquifer interaction.

NASA Astrophysics Data System (ADS)

Liu, X.; Zhan, H.

2017-12-01

This poster presents a new analytical solution to study water exchange, hydraulic head distribution and water flow in a stream-unconfined aquifer interaction system with a sloping bed and stream of varying heads in presence of two thin vertical sedimentary layers. The formation of a clogging bed of fine-grained sediments allows the interfaces among a sloping aquifer and two rivers as the third kind and Cauchy boundary conditions. The numerical solution of the corresponding nonlinear Boussinesq equation is also developed to compare the performance of the analytical solution. The effects of precipitation recharge, bed slope and stage variation rate of two rivers for water flow in the sloping aquifer are discussed in the results.

Analysis of silicon stress/strain relationships

NASA Technical Reports Server (NTRS)

Dillon, O.

1985-01-01

In the study of stress-strain relationships in silicon ribbon, numerous solutions were calculated for stresses, strain rates, and dislocation densities through the use of the Sumino model. It was concluded that many cases of failure of computer solutions to converge are analytical manifestations of shear bands (Luder's band) observed in experiments.

A Review of Meta-Analysis Packages in R

ERIC Educational Resources Information Center

Polanin, Joshua R.; Hennessy, Emily A.; Tanner-Smith, Emily E.

2017-01-01

Meta-analysis is a statistical technique that allows an analyst to synthesize effect sizes from multiple primary studies. To estimate meta-analysis models, the open-source statistical environment R is quickly becoming a popular choice. The meta-analytic community has contributed to this growth by developing numerous packages specific to…

Rotor-to-stator Partial Rubbing and Its Effects on Rotor Dynamic Response

NASA Technical Reports Server (NTRS)

Muszynska, Agnes; Franklin, Wesley D.; Hayashida, Robert D.

1991-01-01

Results from experimental and analytical studies on rotor to stationary element partial rubbings at several locations and their effects on rotor dynamic responses are presented. The mathematical model of a rubbing rotor is given. The computer program provides numerical results which agree with experimentally obtained rotor responses.

Development of an integrated BEM for hot fluid-structure interaction

NASA Technical Reports Server (NTRS)

Banerjee, P. K.; Dargush, G. F.

1989-01-01

The Boundary Element Method (BEM) is chosen as a basic analysis tool principally because the definition of quantities like fluxes, temperature, displacements, and velocities is very precise on a boundary base discretization scheme. One fundamental difficulty is, of course, that the entire analysis requires a very considerable amount of analytical work which is not present in other numerical methods. During the last 18 months all of this analytical work was completed and a two-dimensional, general purpose code was written. Some of the early results are described. It is anticipated that within the next two to three months almost all two-dimensional idealizations will be examined. It should be noted that the analytical work for the three-dimensional case has also been done and numerical implementation will begin next year.

LARC-1: a Los Alamos release calculation program for fission product transport in HTGRs during the LOFC accident

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

Carruthers, L.M.; Lee, C.E.

1976-10-01

The theoretical and numerical data base development of the LARC-1 code is described. Four analytical models of fission product release from an HTGR core during the loss of forced circulation accident are developed. Effects of diffusion, adsorption and evaporation of the metallics and precursors are neglected in this first LARC model. Comparison of the analytic models indicates that the constant release-renormalized model is adequate to describe the processes involved. The numerical data base for release constants, temperature modeling, fission product release rates, coated fuel particle failure fraction and aged coated fuel particle failure fractions is discussed. Analytic fits and graphicmore » displays for these data are given for the Ft. St. Vrain and GASSAR models.« less

Standards and Guidelines for Numerical Models for Tsunami Hazard Mitigation

NASA Astrophysics Data System (ADS)

Titov, V.; Gonzalez, F.; Kanoglu, U.; Yalciner, A.; Synolakis, C. E.

2006-12-01

An increased number of nations around the workd need to develop tsunami mitigation plans which invariably involve inundation maps for warning guidance and evacuation planning. There is the risk that inundation maps may be produced with older or untested methodology, as there are currently no standards for modeling tools. In the aftermath of the 2004 megatsunami, some models were used to model inundation for Cascadia events with results much larger than sediment records and existing state-of-the-art studies suggest leading to confusion among emergency management. Incorrectly assessing tsunami impact is hazardous, as recent events in 2006 in Tonga, Kythira, Greece and Central Java have suggested (Synolakis and Bernard, 2006). To calculate tsunami currents, forces and runup on coastal structures, and inundation of coastlines one must calculate the evolution of the tsunami wave from the deep ocean to its target site, numerically. No matter what the numerical model, validation (the process of ensuring that the model solves the parent equations of motion accurately) and verification (the process of ensuring that the model used represents geophysical reality appropriately) both are an essential. Validation ensures that the model performs well in a wide range of circumstances and is accomplished through comparison with analytical solutions. Verification ensures that the computational code performs well over a range of geophysical problems. A few analytic solutions have been validated themselves with laboratory data. Even fewer existing numerical models have been both validated with the analytical solutions and verified with both laboratory measurements and field measurements, thus establishing a gold standard for numerical codes for inundation mapping. While there is in principle no absolute certainty that a numerical code that has performed well in all the benchmark tests will also produce correct inundation predictions with any given source motions, validated codes reduce the level of uncertainty in their results to the uncertainty in the geophysical initial conditions. Further, when coupled with real--time free--field tsunami measurements from tsunameters, validated codes are the only choice for realistic forecasting of inundation; the consequences of failure are too ghastly to take chances with numerical procedures that have not been validated. We discuss a ten step process of benchmark tests for models used for inundation mapping. The associated methodology and algorithmes have to first be validated with analytical solutions, then verified with laboratory measurements and field data. The models need to be published in the scientific literature in peer-review journals indexed by ISI. While this process may appear onerous, it reflects our state of knowledge, and is the only defensible methodology when human lives are at stake. Synolakis, C.E., and Bernard, E.N, Tsunami science before and beyond Boxing Day 2004, Phil. Trans. R. Soc. A 364 1845, 2231--2263, 2005.

4D volcano gravimetry

USGS Publications Warehouse

Battaglia, Maurizio; Gottsmann, J.; Carbone, D.; Fernandez, J.

2008-01-01

Time-dependent gravimetric measurements can detect subsurface processes long before magma flow leads to earthquakes or other eruption precursors. The ability of gravity measurements to detect subsurface mass flow is greatly enhanced if gravity measurements are analyzed and modeled with ground-deformation data. Obtaining the maximum information from microgravity studies requires careful evaluation of the layout of network benchmarks, the gravity environmental signal, and the coupling between gravity changes and crustal deformation. When changes in the system under study are fast (hours to weeks), as in hydrothermal systems and restless volcanoes, continuous gravity observations at selected sites can help to capture many details of the dynamics of the intrusive sources. Despite the instrumental effects, mainly caused by atmospheric temperature, results from monitoring at Mt. Etna volcano show that continuous measurements are a powerful tool for monitoring and studying volcanoes.Several analytical and numerical mathematical models can beused to fit gravity and deformation data. Analytical models offer a closed-form description of the volcanic source. In principle, this allows one to readily infer the relative importance of the source parameters. In active volcanic sites such as Long Valley caldera (California, U.S.A.) and Campi Flegrei (Italy), careful use of analytical models and high-quality data sets has produced good results. However, the simplifications that make analytical models tractable might result in misleading volcanological inter-pretations, particularly when the real crust surrounding the source is far from the homogeneous/ isotropic assumption. Using numerical models allows consideration of more realistic descriptions of the sources and of the crust where they are located (e.g., vertical and lateral mechanical discontinuities, complex source geometries, and topography). Applications at Teide volcano (Tenerife) and Campi Flegrei demonstrate the importance of this more realistic description in gravity calculations. ?? 2008 Society of Exploration Geophysicists. All rights reserved.

Stomatal Opening: The Role of Cell-Wall Mechanical Anisotropy and Its Analytical Relations to the Bio-composite Characteristics

PubMed Central

Marom, Ziv; Shtein, Ilana; Bar-On, Benny

2017-01-01

Stomata are pores on the leaf surface, which are formed by a pair of curved, tubular guard cells; an increase in turgor pressure deforms the guard cells, resulting in the opening of the stomata. Recent studies employed numerical simulations, based on experimental data, to analyze the effects of various structural, chemical, and mechanical features of the guard cells on the stomatal opening characteristics; these studies all support the well-known qualitative observation that the mechanical anisotropy of the guard cells plays a critical role in stomatal opening. Here, we propose a computationally based analytical model that quantitatively establishes the relations between the degree of anisotropy of the guard cell, the bio-composite constituents of the cell wall, and the aperture and area of stomatal opening. The model introduces two non-dimensional key parameters that dominate the guard cell deformations—the inflation driving force and the anisotropy ratio—and it serves as a generic framework that is not limited to specific plant species. The modeling predictions are in line with a wide range of previous experimental studies, and its analytical formulation sheds new light on the relations between the structure, mechanics, and function of stomata. Moreover, the model provides an analytical tool to back-calculate the elastic characteristics of the matrix that composes the guard cell walls, which, to the best of our knowledge, cannot be probed by direct nano-mechanical experiments; indeed, the estimations of our model are in good agreement with recently published results of independent numerical optimization schemes. The emerging insights from the stomatal structure-mechanics “design guidelines” may promote the development of miniature, yet complex, multiscale composite actuation mechanisms for future engineering platforms. PMID:29312365

FAST MAGNETIC FIELD AMPLIFICATION IN THE EARLY UNIVERSE: GROWTH OF COLLISIONLESS PLASMA INSTABILITIES IN TURBULENT MEDIA

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

Falceta-Gonçalves, D.; Kowal, G.

2015-07-20

In this work we report on a numerical study of the cosmic magnetic field amplification due to collisionless plasma instabilities. The collisionless magnetohydrodynamic equations derived account for the pressure anisotropy that leads, in specific conditions, to the firehose and mirror instabilities. We study the time evolution of seed fields in turbulence under the influence of such instabilities. An approximate analytical time evolution of the magnetic field is provided. The numerical simulations and the analytical predictions are compared. We found that (i) amplification of the magnetic field was efficient in firehose-unstable turbulent regimes, but not in the mirror-unstable models; (ii) the growthmore » rate of the magnetic energy density is much faster than the turbulent dynamo; and (iii) the efficient amplification occurs at small scales. The analytical prediction for the correlation between the growth timescales and pressure anisotropy is confirmed by the numerical simulations. These results reinforce the idea that pressure anisotropies—driven naturally in a turbulent collisionless medium, e.g., the intergalactic medium, could efficiently amplify the magnetic field in the early universe (post-recombination era), previous to the collapse of the first large-scale gravitational structures. This mechanism, though fast for the small-scale fields (∼kpc scales), is unable to provide relatively strong magnetic fields at large scales. Other mechanisms that were not accounted for here (e.g., collisional turbulence once instabilities are quenched, velocity shear, or gravitationally induced inflows of gas into galaxies and clusters) could operate afterward to build up large-scale coherent field structures in the long time evolution.« less

Orbital theory in terms of KS elements with luni-solar perturbations

NASA Astrophysics Data System (ADS)

Sellamuthu, Harishkumar; Sharma, Ram

2016-07-01

Precise orbit computation of Earth orbiting satellites is essential for efficient mission planning of planetary exploration, navigation and satellite geodesy. The third-body perturbations of the Sun and the Moon predominantly affect the satellite motion in the high altitude and elliptical orbits, where the effect of atmospheric drag is negligible. The physics of the luni-solar gravity effect on Earth satellites have been studied extensively over the years. The combined luni-solar gravitational attraction will induce a cumulative effect on the dynamics of satellite orbits, which mainly oscillates the perigee altitude. Though accurate orbital parameters are computed by numerical integration with respect to complex force models, analytical theories are highly valued for the manifold of solutions restricted to relatively simple force models. During close approach, the classical equations of motion in celestial mechanics are almost singular and they are unstable for long-term orbit propagation. A new singularity-free analytical theory in terms of KS (Kustaanheimo and Stiefel) regular elements with respect to luni-solar perturbation is developed. These equations are regular everywhere and eccentric anomaly is the independent variable. Plataforma Solar de Almería (PSA) algorithm and a Fourier series algorithm are used to compute the accurate positions of the Sun and the Moon, respectively. Numerical studies are carried out for wide range of initial parameters and the analytical solutions are found to be satisfactory when compared with numerically integrated values. The symmetrical nature of the equations allows only two of the nine equations to be solved for computing the state vectors and the time. Only a change in the initial conditions is required to solve the other equations. This theory will find multiple applications including on-board software packages and for mission analysis purposes.

Computing the optimal path in stochastic dynamical systems

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

Bauver, Martha; Forgoston, Eric, E-mail: eric.forgoston@montclair.edu; Billings, Lora

2016-08-15

In stochastic systems, one is often interested in finding the optimal path that maximizes the probability of escape from a metastable state or of switching between metastable states. Even for simple systems, it may be impossible to find an analytic form of the optimal path, and in high-dimensional systems, this is almost always the case. In this article, we formulate a constructive methodology that is used to compute the optimal path numerically. The method utilizes finite-time Lyapunov exponents, statistical selection criteria, and a Newton-based iterative minimizing scheme. The method is applied to four examples. The first example is a two-dimensionalmore » system that describes a single population with internal noise. This model has an analytical solution for the optimal path. The numerical solution found using our computational method agrees well with the analytical result. The second example is a more complicated four-dimensional system where our numerical method must be used to find the optimal path. The third example, although a seemingly simple two-dimensional system, demonstrates the success of our method in finding the optimal path where other numerical methods are known to fail. In the fourth example, the optimal path lies in six-dimensional space and demonstrates the power of our method in computing paths in higher-dimensional spaces.« less

Exhaled breath condensate – from an analytical point of view

PubMed Central

Dodig, Slavica; Čepelak, Ivana

2013-01-01

Over the past three decades, the goal of many researchers is analysis of exhaled breath condensate (EBC) as noninvasively obtained sample. A total quality in laboratory diagnostic processes in EBC analysis was investigated: pre-analytical (formation, collection, storage of EBC), analytical (sensitivity of applied methods, standardization) and post-analytical (interpretation of results) phases. EBC analysis is still used as a research tool. Limitations referred to pre-analytical, analytical, and post-analytical phases of EBC analysis are numerous, e.g. low concentrations of EBC constituents, single-analyte methods lack in sensitivity, and multi-analyte has not been fully explored, and reference values are not established. When all, pre-analytical, analytical and post-analytical requirements are met, EBC biomarkers as well as biomarker patterns can be selected and EBC analysis can hopefully be used in clinical practice, in both, the diagnosis and in the longitudinal follow-up of patients, resulting in better outcome of disease. PMID:24266297

A semi-analytical study of positive corona discharge in wire-plane electrode configuration

NASA Astrophysics Data System (ADS)

Yanallah, K.; Pontiga, F.; Chen, J. H.

2013-08-01

Wire-to-plane positive corona discharge in air has been studied using an analytical model of two species (electrons and positive ions). The spatial distributions of electric field and charged species are obtained by integrating Gauss's law and the continuity equations of species along the Laplacian field lines. The experimental values of corona current intensity and applied voltage, together with Warburg's law, have been used to formulate the boundary condition for the electron density on the corona wire. To test the accuracy of the model, the approximate electric field distribution has been compared with the exact numerical solution obtained from a finite element analysis. A parametrical study of wire-to-plane corona discharge has then been undertaken using the approximate semi-analytical solutions. Thus, the spatial distributions of electric field and charged particles have been computed for different values of the gas pressure, wire radius and electrode separation. Also, the two dimensional distribution of ozone density has been obtained using a simplified plasma chemistry model. The approximate semi-analytical solutions can be evaluated in a negligible computational time, yet provide precise estimates of corona discharge variables.

The elastic and inelastic behavior of woven graphite fabric reinforced polyimide composites

NASA Astrophysics Data System (ADS)

Searles, Kevin H.

In many aerospace and conventional engineering applications, load-bearing composite structures are designed with the intent of being subjected to uniaxial stresses that are predominantly tensile or compressive. However, it is likely that biaxial and possibly triaxial states of stress will exist throughout the in-service life of the structure or component. The existing paradigm suggests that unidirectional tape materials are superior under uniaxial conditions since the vast majority of fibers lie in-plane and can be aligned to the loading axis. This may be true, but not without detriment to impact performance, interlaminar strength, strain to failure and complexity of part geometry. In circumstances where a sufficient balance of these properties is required, composites based on woven fabric reinforcements become attractive choices. In this thesis, the micro- and mesoscale elastic behavior of composites based on 8HS woven graphite fabric architectures and polyimide matrices is studied analytically and numerically. An analytical model is proposed to predict the composite elastic constants and is verified using numerical strain energy methods of equivalence. The model shows good agreement with the experiments and numerical strain energy equivalence. Lamina stresses generated numerically from in-plane shear loading show substantial shear and transverse normal stress concentrations in the transverse undulated tow which potentially leads to intralaminar damage. The macroscale inelastic behavior of the same composites is also studied experimentally and numerically. On an experimental basis, the biaxial and modified biaxial Iosipescu test methods are employed to study the weaker-mode shear and biaxial failure properties at room and elevated temperatures. On a numerical basis, the macroscale inelastic shear behavior of the composites is studied. Structural nonlinearities and material nonlinearities are identified and resolved. In terms of specimen-to-fixture interactions, load eccentricities, geometric (large strains and rotations) nonlinearities and boundary contact (friction) nonlinearities are explored. In terms of material nonlinearities, anisotropic plasticity and progressive damage are explored. A progressive damage criterion is proposed which accounts for the elastic strain energy densities in three directions. Of the types of nonlinearities studied, the nonlinear shear stress-strain behavior of the composites is principally from progressive intralaminar damage. Structural nonlinearities and elastoplastic deformation appear to be inconsequential.

Numerical Simulation and Experimental Study of a Dental Handpiece Air Turbine

NASA Astrophysics Data System (ADS)

Hsu, Chih-Neng; Chiang, Hsiao-Wei D.; Chang, Ya-Yi

2011-06-01

Dental air turbine handpieces have been widely used in clinical dentistry for over 30 years, however, little work has been reported on their performance. In dental air turbine handpieces, the types of flow channel and turbine blade shape can have very different designs. These different designs can have major influence on the torque, rotating speed, and power performance. This research is focused on the turbine blade and the flow channel designs. Using numerical simulation and experiments, the key design parameters which influence the performance of dental hand pieces can be studied. Three types of dental air turbine designs with different turbine blades, nozzle angles, nozzle flow channels, and shroud clearances were tested and analyzed. Very good agreement was demonstrated between the numerical simulation analyses and the experiments. Using the analytical model, parametric studies were performed to identify key design parameters.

Cocontraction of pairs of antagonistic muscles: analytical solution for planar static nonlinear optimization approaches.

PubMed

Herzog, W; Binding, P

1993-11-01

It has been stated in the literature that static, nonlinear optimization approaches cannot predict coactivation of pairs of antagonistic muscles; however, numerical solutions of such approaches have predicted coactivation of pairs of one-joint and multijoint antagonists. Analytical support for either finding is not available in the literature for systems containing more than one degree of freedom. The purpose of this study was to investigate analytically the possibility of cocontraction of pairs of antagonistic muscles using a static nonlinear optimization approach for a multidegree-of-freedom, two-dimensional system. Analytical solutions were found using the Karush-Kuhn-Tucker conditions, which were necessary and sufficient for optimality in this problem. The results show that cocontraction of pairs of one-joint antagonistic muscles is not possible, whereas cocontraction of pairs of multijoint antagonists is. These findings suggest that cocontraction of pairs of antagonistic muscles may be an "efficient" way to accomplish many movement tasks.

An efficient method for the calculation of mean extinction. I - The analyticity of the complex extinction efficiency of homogeneous spheres

NASA Astrophysics Data System (ADS)

Xing, Zhang-Fan; Greenberg, J. M.

1992-11-01

Results of an investigation of the analyticity of the complex extinction efficiency Q-tilde(ext) in different parameter domains are presented. In the size parameter domain, x = omega(a/c), numerical Hilbert transforms are used to study the analyticity properties of Q-tilde(ext) for homogeneous spheres. Q-tilde(ext) is found to be analytic in the entire lower complex x-tilde-plane when the refractive index, m, is fixed as a real constant (pure scattering) or infinity (perfect conductor); poles, however, appear in the left side of the lower complex x-tilde-plane as m becomes complex. The computation of the mean extinction produced by an extended size distribution of particles may be conveniently and accurately approximated using only a few values of the complex extinction evaluated in the complex plane.

Cross-stream diffusion under pressure-driven flow in microchannels with arbitrary aspect ratios: a phase diagram study using a three-dimensional analytical model

PubMed Central

Song, Hongjun; Wang, Yi; Pant, Kapil

2011-01-01

This article presents a three-dimensional analytical model to investigate cross-stream diffusion transport in rectangular microchannels with arbitrary aspect ratios under pressure-driven flow. The Fourier series solution to the three-dimensional convection–diffusion equation is obtained using a double integral transformation method and associated eigensystem calculation. A phase diagram derived from the dimensional analysis is presented to thoroughly interrogate the characteristics in various transport regimes and examine the validity of the model. The analytical model is verified against both experimental and numerical models in terms of the concentration profile, diffusion scaling law, and mixing efficiency with excellent agreement (with <0.5% relative error). Quantitative comparison against other prior analytical models in extensive parameter space is also performed, which demonstrates that the present model accommodates much broader transport regimes with significantly enhanced applicability. PMID:22247719

Cross-stream diffusion under pressure-driven flow in microchannels with arbitrary aspect ratios: a phase diagram study using a three-dimensional analytical model.

PubMed

Song, Hongjun; Wang, Yi; Pant, Kapil

2012-01-01

This article presents a three-dimensional analytical model to investigate cross-stream diffusion transport in rectangular microchannels with arbitrary aspect ratios under pressure-driven flow. The Fourier series solution to the three-dimensional convection-diffusion equation is obtained using a double integral transformation method and associated eigensystem calculation. A phase diagram derived from the dimensional analysis is presented to thoroughly interrogate the characteristics in various transport regimes and examine the validity of the model. The analytical model is verified against both experimental and numerical models in terms of the concentration profile, diffusion scaling law, and mixing efficiency with excellent agreement (with <0.5% relative error). Quantitative comparison against other prior analytical models in extensive parameter space is also performed, which demonstrates that the present model accommodates much broader transport regimes with significantly enhanced applicability.

Analytical drafting curves provide exact equations for plotted data

NASA Technical Reports Server (NTRS)

Stewart, R. B.

1967-01-01

Analytical drafting curves provide explicit mathematical expressions for any numerical data that appears in the form of graphical plots. The curves each have a reference coordinate axis system indicated on the curve as well as the mathematical equation from which the curve was generated.

Orthopositronium Lifetime: Analytic Results in O({alpha}) and O({alpha}{sup 3}ln{alpha})

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

Kniehl, Bernd A.; Kotikov, Anatoly V.; Veretin, Oleg L.

2008-11-07

We present the O({alpha}) and O({alpha}{sup 3}ln{alpha}) corrections to the total decay width of orthopositronium in closed analytic form, in terms of basic irrational numbers, which can be evaluated numerically to arbitrary precision.

Effect of coulomb spline on rotor dynamic response

NASA Technical Reports Server (NTRS)

Nataraj, C.; Nelson, H. D.; Arakere, N.

1985-01-01

A rigid rotor system coupled by a coulomb spline is modelled and analyzed by approximate analytical and numerical analytical methods. Expressions are derived for the variables of the resulting limit cycle and are shown to be quite accurate for a small departure from isotropy.

Analytical solutions for benchmarking cold regions subsurface water flow and energy transport models: one-dimensional soil thaw with conduction and advection

USGS Publications Warehouse

Kurylyk, Barret L.; McKenzie, Jeffrey M; MacQuarrie, Kerry T. B.; Voss, Clifford I.

2014-01-01

Numerous cold regions water flow and energy transport models have emerged in recent years. Dissimilarities often exist in their mathematical formulations and/or numerical solution techniques, but few analytical solutions exist for benchmarking flow and energy transport models that include pore water phase change. This paper presents a detailed derivation of the Lunardini solution, an approximate analytical solution for predicting soil thawing subject to conduction, advection, and phase change. Fifteen thawing scenarios are examined by considering differences in porosity, surface temperature, Darcy velocity, and initial temperature. The accuracy of the Lunardini solution is shown to be proportional to the Stefan number. The analytical solution results obtained for soil thawing scenarios with water flow and advection are compared to those obtained from the finite element model SUTRA. Three problems, two involving the Lunardini solution and one involving the classic Neumann solution, are recommended as standard benchmarks for future model development and testing.

Predicting playing frequencies for clarinets: A comparison between numerical simulations and simplified analytical formulas.

PubMed

Coyle, Whitney L; Guillemain, Philippe; Kergomard, Jean; Dalmont, Jean-Pierre

2015-11-01

When designing a wind instrument such as a clarinet, it can be useful to be able to predict the playing frequencies. This paper presents an analytical method to deduce these playing frequencies using the input impedance curve. Specifically there are two control parameters that have a significant influence on the playing frequency, the blowing pressure and reed opening. Four effects are known to alter the playing frequency and are examined separately: the flow rate due to the reed motion, the reed dynamics, the inharmonicity of the resonator, and the temperature gradient within the clarinet. The resulting playing frequencies for the first register of a particular professional level clarinet are found using the analytical formulas presented in this paper. The analytical predictions are then compared to numerically simulated results to validate the prediction accuracy. The main conclusion is that in general the playing frequency decreases above the oscillation threshold because of inharmonicity, then increases above the beating reed regime threshold because of the decrease of the flow rate effect.

How to Find a Bug in Ten Thousand Lines Transport Solver? Outline of Experiences from AN Advection-Diffusion Code Verification

NASA Astrophysics Data System (ADS)

Zamani, K.; Bombardelli, F.

2011-12-01

Almost all natural phenomena on Earth are highly nonlinear. Even simplifications to the equations describing nature usually end up being nonlinear partial differential equations. Transport (ADR) equation is a pivotal equation in atmospheric sciences and water quality. This nonlinear equation needs to be solved numerically for practical purposes so academicians and engineers thoroughly rely on the assistance of numerical codes. Thus, numerical codes require verification before they are utilized for multiple applications in science and engineering. Model verification is a mathematical procedure whereby a numerical code is checked to assure the governing equation is properly solved as it is described in the design document. CFD verification is not a straightforward and well-defined course. Only a complete test suite can uncover all the limitations and bugs. Results are needed to be assessed to make a distinction between bug-induced-defect and innate limitation of a numerical scheme. As Roache (2009) said, numerical verification is a state-of-the-art procedure. Sometimes novel tricks work out. This study conveys the synopsis of the experiences we gained during a comprehensive verification process which was done for a transport solver. A test suite was designed including unit tests and algorithmic tests. Tests were layered in complexity in several dimensions from simple to complex. Acceptance criteria defined for the desirable capabilities of the transport code such as order of accuracy, mass conservation, handling stiff source term, spurious oscillation, and initial shape preservation. At the begining, mesh convergence study which is the main craft of the verification is performed. To that end, analytical solution of ADR equation gathered. Also a new solution was derived. In the more general cases, lack of analytical solution could be overcome through Richardson Extrapolation and Manufactured Solution. Then, two bugs which were concealed during the mesh convergence study uncovered with the method of false injection and visualization of the results. Symmetry had dual functionality: there was a bug, which was hidden due to the symmetric nature of a test (it was detected afterward utilizing artificial false injection), on the other hand self-symmetry was used to design a new test, and in a case the analytical solution of the ADR equation was unknown. Assisting subroutines designed to check and post-process conservation of mass and oscillatory behavior. Finally, capability of the solver also checked for stiff reaction source term. The above test suite not only was a decent tool of error detection but also it provided a thorough feedback on the ADR solvers limitations. Such information is the crux of any rigorous numerical modeling for a modeler who deals with surface/subsurface pollution transport.

Exact Dispersion Study of an Asymmetric Thin Planar Slab Dielectric Waveguide without Computing {d^2}β/{d{k^2}} Numerically

NASA Astrophysics Data System (ADS)

Raghuwanshi, Sanjeev Kumar; Palodiya, Vikram

2017-08-01

Waveguide dispersion can be tailored but not the material dispersion. Hence, the total dispersion can be shifted at any desired band by adjusting the waveguide dispersion. Waveguide dispersion is proportional to {d^2}β/d{k^2} and need to be computed numerically. In this paper, we have tried to compute analytical expression for {d^2}β/d{k^2} in terms of {d^2}β/d{k^2} accurately with numerical technique, ≈ 10^{-5} decimal point. This constraint sometimes generates the error in calculation of waveguide dispersion. To formulate the problem we will use the graphical method. Our study reveals that we can compute the waveguide dispersion enough accurately for various modes by knowing - β only.

A deterministic particle method for one-dimensional reaction-diffusion equations

NASA Technical Reports Server (NTRS)

Mascagni, Michael

1995-01-01

We derive a deterministic particle method for the solution of nonlinear reaction-diffusion equations in one spatial dimension. This deterministic method is an analog of a Monte Carlo method for the solution of these problems that has been previously investigated by the author. The deterministic method leads to the consideration of a system of ordinary differential equations for the positions of suitably defined particles. We then consider the time explicit and implicit methods for this system of ordinary differential equations and we study a Picard and Newton iteration for the solution of the implicit system. Next we solve numerically this system and study the discretization error both analytically and numerically. Numerical computation shows that this deterministic method is automatically adaptive to large gradients in the solution.

Numerical prediction of the energy efficiency of the three-dimensional fish school using the discretized Adomian decomposition method

NASA Astrophysics Data System (ADS)

Lin, Yinwei

2018-06-01

A three-dimensional modeling of fish school performed by a modified Adomian decomposition method (ADM) discretized by the finite difference method is proposed. To our knowledge, few studies of the fish school are documented due to expensive cost of numerical computing and tedious three-dimensional data analysis. Here, we propose a simple model replied on the Adomian decomposition method to estimate the efficiency of energy saving of the flow motion of the fish school. First, the analytic solutions of Navier-Stokes equations are used for numerical validation. The influences of the distance between the side-by-side two fishes are studied on the energy efficiency of the fish school. In addition, the complete error analysis for this method is presented.

A strategy to determine operating parameters in tissue engineering hollow fiber bioreactors

PubMed Central

Shipley, RJ; Davidson, AJ; Chan, K; Chaudhuri, JB; Waters, SL; Ellis, MJ

2011-01-01

The development of tissue engineering hollow fiber bioreactors (HFB) requires the optimal design of the geometry and operation parameters of the system. This article provides a strategy for specifying operating conditions for the system based on mathematical models of oxygen delivery to the cell population. Analytical and numerical solutions of these models are developed based on Michaelis–Menten kinetics. Depending on the minimum oxygen concentration required to culture a functional cell population, together with the oxygen uptake kinetics, the strategy dictates the model needed to describe mass transport so that the operating conditions can be defined. If cmin ≫ Km we capture oxygen uptake using zero-order kinetics and proceed analytically. This enables operating equations to be developed that allow the user to choose the medium flow rate, lumen length, and ECS depth to provide a prescribed value of cmin. When , we use numerical techniques to solve full Michaelis–Menten kinetics and present operating data for the bioreactor. The strategy presented utilizes both analytical and numerical approaches and can be applied to any cell type with known oxygen transport properties and uptake kinetics. PMID:21370228

Modified symplectic schemes with nearly-analytic discrete operators for acoustic wave simulations

NASA Astrophysics Data System (ADS)

Liu, Shaolin; Yang, Dinghui; Lang, Chao; Wang, Wenshuai; Pan, Zhide

2017-04-01

Using a structure-preserving algorithm significantly increases the computational efficiency of solving wave equations. However, only a few explicit symplectic schemes are available in the literature, and the capabilities of these symplectic schemes have not been sufficiently exploited. Here, we propose a modified strategy to construct explicit symplectic schemes for time advance. The acoustic wave equation is transformed into a Hamiltonian system. The classical symplectic partitioned Runge-Kutta (PRK) method is used for the temporal discretization. Additional spatial differential terms are added to the PRK schemes to form the modified symplectic methods and then two modified time-advancing symplectic methods with all of positive symplectic coefficients are then constructed. The spatial differential operators are approximated by nearly-analytic discrete (NAD) operators, and we call the fully discretized scheme modified symplectic nearly analytic discrete (MSNAD) method. Theoretical analyses show that the MSNAD methods exhibit less numerical dispersion and higher stability limits than conventional methods. Three numerical experiments are conducted to verify the advantages of the MSNAD methods, such as their numerical accuracy, computational cost, stability, and long-term calculation capability.

An investigation of the information propagation and entropy transport aspects of Stirling machine numerical simulation

NASA Technical Reports Server (NTRS)

Goldberg, Louis F.

1992-01-01

Aspects of the information propagation modeling behavior of integral machine computer simulation programs are investigated in terms of a transmission line. In particular, the effects of pressure-linking and temporal integration algorithms on the amplitude ratio and phase angle predictions are compared against experimental and closed-form analytic data. It is concluded that the discretized, first order conservation balances may not be adequate for modeling information propagation effects at characteristic numbers less than about 24. An entropy transport equation suitable for generalized use in Stirling machine simulation is developed. The equation is evaluated by including it in a simulation of an incompressible oscillating flow apparatus designed to demonstrate the effect of flow oscillations on the enhancement of thermal diffusion. Numerical false diffusion is found to be a major factor inhibiting validation of the simulation predictions with experimental and closed-form analytic data. A generalized false diffusion correction algorithm is developed which allows the numerical results to match their analytic counterparts. Under these conditions, the simulation yields entropy predictions which satisfy Clausius' inequality.

Sound Emission of Rotor Induced Deformations of Generator Casings

NASA Technical Reports Server (NTRS)

Polifke, W.; Mueller, B.; Yee, H. C.; Mansour, Nagi (Technical Monitor)

2001-01-01

The casing of large electrical generators can be deformed slightly by the rotor's magnetic field. The sound emission produced by these periodic deformations, which could possibly exceed guaranteed noise emission limits, is analysed analytically and numerically. From the deformation of the casing, the normal velocity of the generator's surface is computed. Taking into account the corresponding symmetry, an analytical solution for the acoustic pressure outside the generator is round in terms of the Hankel function of second order. The normal velocity or the generator surface provides the required boundary condition for the acoustic pressure and determines the magnitude of pressure oscillations. For the numerical simulation, the nonlinear 2D Euler equations are formulated In a perturbation form for low Mach number Computational Aeroacoustics (CAA). The spatial derivatives are discretized by the classical sixth-order central interior scheme and a third-order boundary scheme. Spurious high frequency oscillations are damped by a characteristic-based artificial compression method (ACM) filter. The time derivatives are approximated by the classical 4th-order Runge-Kutta method. The numerical results are In excellent agreement with the analytical solution.

Complex Dynamics of Delay-Coupled Neural Networks

NASA Astrophysics Data System (ADS)

Mao, Xiaochen

2016-09-01

This paper reveals the complicated dynamics of a delay-coupled system that consists of a pair of sub-networks and multiple bidirectional couplings. Time delays are introduced into the internal connections and network-couplings, respectively. The stability and instability of the coupled network are discussed. The sufficient conditions for the existence of oscillations are given. Case studies of numerical simulations are given to validate the analytical results. Interesting and complicated neuronal activities are observed numerically, such as rest states, periodic oscillations, multiple switches of rest states and oscillations, and the coexistence of different types of oscillations.

Study of transient behavior of finned coil heat exchangers

NASA Technical Reports Server (NTRS)

Rooke, S. P.; Elissa, M. G.

1993-01-01

The status of research on the transient behavior of finned coil cross-flow heat exchangers using single phase fluids is reviewed. Applications with available analytical or numerical solutions are discussed. Investigation of water-to-air type cross-flow finned tube heat exchangers is examined through the use of simplified governing equations and an up-wind finite difference scheme. The degenerate case of zero air-side capacitance rate is compared with available exact solution. Generalization of the numerical model is discussed for application to multi-row multi-circuit heat exchangers.

Solitary wave solutions and their interactions for fully nonlinear water waves with surface tension in the generalized Serre equations

NASA Astrophysics Data System (ADS)

Dutykh, Denys; Hoefer, Mark; Mitsotakis, Dimitrios

2018-04-01

Some effects of surface tension on fully nonlinear, long, surface water waves are studied by numerical means. The differences between various solitary waves and their interactions in subcritical and supercritical surface tension regimes are presented. Analytical expressions for new peaked traveling wave solutions are presented in the dispersionless case of critical surface tension. Numerical experiments are performed using a high-accurate finite element method based on smooth cubic splines and the four-stage, classical, explicit Runge-Kutta method of order 4.

Stress-Strain State of a Combinational Soil Half-Space During Reconstruction

NASA Astrophysics Data System (ADS)

Prusov, D. E.

2014-03-01

A method for studying the stress-strain state of soil-retaining structures is proposed. It is based on the nonlinear theory of elasticity and plasticity of soils and allows for geometrical and physical nonlinearities. Numerical and analytical results on the stability of a retaining wall are compared. The influence of an inhomogeneous soil half-space on the stress-strain state of a deep-ditch wall is analyzed numerically. A scientific rationale for the redevelopment of densely built-up residential areas under adverse geological engineering conditions is recommended.

Realization of the Najafi-Golestanian microswimmer

NASA Astrophysics Data System (ADS)

Grosjean, Galien; Hubert, Maxime; Lagubeau, Guillaume; Vandewalle, Nicolas

2016-08-01

A paradigmatic microswimmer is the three-linked-spheres model, which follows a minimalist approach for propulsion by shape shifting. As such, it has been the subject of numerous analytical and numerical studies. In this Rapid Communication, an experimental three-linked-spheres swimmer is created by self-assembling ferromagnetic particles at an air-water interface. It is powered by a uniform oscillating magnetic field. A model, using two harmonic oscillators, reproduces the experimental findings. Because the model remains general, the same approach could be used to design a variety of efficient microswimmers.

Towards a wave-extraction method for numerical relativity. III. Analytical examples for the Beetle-Burko radiation scalar

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

Burko, Lior M.; Department of Physics, University of Alabama in Huntsville, Huntsville, Alabama 35899; Baumgarte, Thomas W.

2006-01-15

Beetle and Burko recently introduced a background-independent scalar curvature invariant for general relativity that carries information about the gravitational radiation in generic spacetimes, in cases where such radiation is incontrovertibly defined. In this paper we adopt a formalism that only uses spatial data as they are used in numerical relativity and compute the Beetle-Burko radiation scalar for a number of analytical examples, specifically linearized Einstein-Rosen cylindrical waves, linearized quadrupole waves, the Kerr spacetime, Bowen-York initial data, and the Kasner spacetime. These examples illustrate how the Beetle-Burko radiation scalar can be used to examine the gravitational wave content of numerically generatedmore » spacetimes, and how it may provide a useful diagnostic for initial data sets.« less

Information Based Numerical Practice.

DTIC Science & Technology

1987-02-01

characterization by comparative computational studies of various benchmark problems. See e.g. [MacNeal, Harder (1985)], [Robinson, Blackham (1981)] any...FOR NONADAPTIVE METHODS 2.1. THE QUADRATURE FORMULA The simplest example studied in detail in the literature is the problem of the optimal quadrature...formulae and the functional analytic prerequisites for the study of optimal formulae, we refer to the large monography (808 p) of [Sobolev (1974)]. Let us

A new approximation for pore pressure accumulation in marine sediment due to water waves

NASA Astrophysics Data System (ADS)

Jeng, D.-S.; Seymour, B. R.; Li, J.

2007-01-01

The residual mechanism of wave-induced pore water pressure accumulation in marine sediments is re-examined. An analytical approximation is derived using a linear relation for pore pressure generation in cyclic loading, and mistakes in previous solutions (Int. J. Numer. Anal. Methods Geomech. 2001; 25:885-907; J. Offshore Mech. Arctic Eng. (ASME) 1989; 111(1):1-11) are corrected. A numerical scheme is then employed to solve the case with a non-linear relation for pore pressure generation. Both analytical and numerical solutions are verified with experimental data (Laboratory and field investigation of wave-sediment interaction. Joseph H. Defrees Hydraulics Laboratory, School of Civil and Environmental Engineering, Cornell University, Ithaca, NY, 1983), and provide a better prediction of pore pressure accumulation than the previous solution (J. Offshore Mech. Arctic Eng. (ASME) 1989; 111(1):1-11). The parametric study concludes that the pore pressure accumulation and use of full non-linear relation of pore pressure become more important under the following conditions: (1) large wave amplitude, (2) longer wave period, (3) shallow water, (4) shallow soil and (5) softer soils with a low consolidation coefficient. Copyright

Approximate Solutions for Ideal Dam-Break Sediment-Laden Flows on Uniform Slopes

NASA Astrophysics Data System (ADS)

Ni, Yufang; Cao, Zhixian; Borthwick, Alistair; Liu, Qingquan

2018-04-01

Shallow water hydro-sediment-morphodynamic (SHSM) models have been applied increasingly widely in hydraulic engineering and geomorphological studies over the past few decades. Analytical and approximate solutions are usually sought to verify such models and therefore confirm their credibility. Dam-break flows are often evoked because such flows normally feature shock waves and contact discontinuities that warrant refined numerical schemes to solve. While analytical and approximate solutions to clear-water dam-break flows have been available for some time, such solutions are rare for sediment transport in dam-break flows. Here we aim to derive approximate solutions for ideal dam-break sediment-laden flows resulting from the sudden release of a finite volume of frictionless, incompressible water-sediment mixture on a uniform slope. The approximate solutions are presented for three typical sediment transport scenarios, i.e., pure advection, pure sedimentation, and concurrent entrainment and deposition. Although the cases considered in this paper are not real, the approximate solutions derived facilitate suitable benchmark tests for evaluating SHSM models, especially presently when shock waves can be numerically resolved accurately with a suite of finite volume methods, while the accuracy of the numerical solutions of contact discontinuities in sediment transport remains generally poorer.

Using nano- and micro-particles of silver in lignin analysis

Treesearch

Umesh P. Agarwal; Richard S. Reiner; Sally A. Ralph

2006-01-01

Although there are numerous techniques available to study lignin in its various states, they all have limitations and to extract most amount of information, a number of analytical techniques have to be jointly used. We have recently started applying a new approach to lignin analysis - namely using nano- and micro-particles of silver for study of native and residual...

A finite element model to study the effect of tissue anisotropy on ex vivo arterial shear wave elastography measurements

NASA Astrophysics Data System (ADS)

Shcherbakova, D. A.; Debusschere, N.; Caenen, A.; Iannaccone, F.; Pernot, M.; Swillens, A.; Segers, P.

2017-07-01

Shear wave elastography (SWE) is an ultrasound (US) diagnostic method for measuring the stiffness of soft tissues based on generated shear waves (SWs). SWE has been applied to bulk tissues, but in arteries it is still under investigation. Previously performed studies in arteries or arterial phantoms demonstrated the potential of SWE to measure arterial wall stiffness—a relevant marker in prediction of cardiovascular diseases. This study is focused on numerical modelling of SWs in ex vivo equine aortic tissue, yet based on experimental SWE measurements with the tissue dynamically loaded while rotating the US probe to investigate the sensitivity of SWE to the anisotropic structure. A good match with experimental shear wave group speed results was obtained. SWs were sensitive to the orthotropy and nonlinearity of the material. The model also allowed to study the nature of the SWs by performing 2D FFT-based and analytical phase analyses. A good match between numerical group velocities derived using the time-of-flight algorithm and derived from the dispersion curves was found in the cross-sectional and axial arterial views. The complexity of solving analytical equations for nonlinear orthotropic stressed plates was discussed.

Analytical, Numerical and Experimental Examination of Reinforced Composites Beams Covered with Carbon Fiber Reinforced Plastic

NASA Astrophysics Data System (ADS)

Kasimzade, A. A.; Tuhta, S.

2012-03-01

In the article, analytical, numerical (Finite Element Method) and experimental investigation results of beam that was strengthened with fiber reinforced plastic-FRP composite has been given as comparative, the effect of FRP wrapping number to the maximum load and moment capacity has been evaluated depending on this results. Carbon FRP qualitative dependences have been occurred between wrapping number and beam load and moment capacity for repair-strengthen the reinforced concrete beams with carbon fiber. Shown possibilities of application traditional known analysis programs, for the analysis of Carbon Fiber Reinforced Plastic (CFRP) strengthened structures.

Many-body delocalization in a strongly disordered system with long-range interactions: Finite-size scaling

NASA Astrophysics Data System (ADS)

Burin, Alexander L.

2015-03-01

Many-body localization in a disordered system of interacting spins coupled by the long-range interaction 1 /Rα is investigated combining analytical theory considering resonant interactions and a finite-size scaling of exact numerical solutions with number of spins N . The numerical results for a one-dimensional system are consistent with the general expectations of analytical theory for a d -dimensional system including the absence of localization in the infinite system at α <2 d and a universal scaling of a critical energy disordering Wc∝N2/d -α d .

Theory and Circuit Model for Lossy Coaxial Transmission Line

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

Genoni, T. C.; Anderson, C. N.; Clark, R. E.

2017-04-01

The theory of signal propagation in lossy coaxial transmission lines is revisited and new approximate analytic formulas for the line impedance and attenuation are derived. The accuracy of these formulas from DC to 100 GHz is demonstrated by comparison to numerical solutions of the exact field equations. Based on this analysis, a new circuit model is described which accurately reproduces the line response over the entire frequency range. Circuit model calculations are in excellent agreement with the numerical and analytic results, and with finite-difference-time-domain simulations which resolve the skindepths of the conducting walls.

Analytical formulation of selected activities of the remote manipulator system

NASA Technical Reports Server (NTRS)

Zimmerman, K. J.

1977-01-01

Existing analysis of Orbiter-RMS-Payload kinematics were surveyed, including equations dealing with the two body kinematics in the presence of a massless RMS and compares analytical explicit solutions with numerical solutions. For the following operational phases of the RMS numerical demonstration, problems are provided: (1) payload capture; (2) payload stowage and removal from cargo bay; and (3) payload deployment. The equation of motion provided accounted for RMS control forces and torque moments and could be extended to RMS flexibility and control loop simulation without increasing the degrees of freedom of the two body system.

Validation of a numerical method for interface-resolving simulation of multicomponent gas-liquid mass transfer and evaluation of multicomponent diffusion models

NASA Astrophysics Data System (ADS)

Woo, Mino; Wörner, Martin; Tischer, Steffen; Deutschmann, Olaf

2018-03-01

The multicomponent model and the effective diffusivity model are well established diffusion models for numerical simulation of single-phase flows consisting of several components but are seldom used for two-phase flows so far. In this paper, a specific numerical model for interfacial mass transfer by means of a continuous single-field concentration formulation is combined with the multicomponent model and effective diffusivity model and is validated for multicomponent mass transfer. For this purpose, several test cases for one-dimensional physical or reactive mass transfer of ternary mixtures are considered. The numerical results are compared with analytical or numerical solutions of the Maxell-Stefan equations and/or experimental data. The composition-dependent elements of the diffusivity matrix of the multicomponent and effective diffusivity model are found to substantially differ for non-dilute conditions. The species mole fraction or concentration profiles computed with both diffusion models are, however, for all test cases very similar and in good agreement with the analytical/numerical solutions or measurements. For practical computations, the effective diffusivity model is recommended due to its simplicity and lower computational costs.

Parametric resonance energy exchange and induction phenomenon in a one-dimensional nonlinear oscillator chain

NASA Astrophysics Data System (ADS)

Yoshimura, K.

2000-11-01

We study analytically the induction phenomenon in the Fermi-Pasta-Ulam β oscillator chain under initial conditions consisting of single mode excitation. Our study is based on the analytical computation of the largest characteristic exponent of an approximate version of the variational equation. The main results can be summarized as follows: (1) the energy density ɛ scaling of the induction time T is given by T~ɛ-1, and T becomes smaller for higher-frequency mode excitation; (2) there is a threshold energy density ɛc such that the induction time diverges when ɛ<ɛc=π2/6βN2, where N is the system size and β the nonlinearity parameter, and this expression for ɛc is correct in the limit N-->∞ (3) the threshold ɛc vanishes as ɛc~N-2 in the limit N-->∞ (4) the threshold ɛc does not depend on the mode number k that is excited in the initial condition; (5) the two modes k+/-m have the largest exponential growth rate, and m increases with increasing ɛ as m/N=3βɛ/π. The above analytical results are thoroughly verified in numerical experiments. Moreover, we discuss the energy exchange process after the induction period in some energy density regimes, based on the numerical results.

Effects of a finite aperture on the Inverse Born Approximation

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

Kogan, V.G.; Rose, J.H.

1983-01-01

One of the most important effects of complex part geometry is that the available entrance and exit angles for ultrasound are limited. We will present a study of the Inverse Born approximation in which we have data for incident (and exit) directions confined to a conical aperture. Modeling the direct problem by the Born Approximation, we obtained analytical results for (1) a weak spherical inclusion, and (2) a penny shaped crack (modeled by an oblate spheroid). General results are: (a) the value of the characteristic function ..gamma.. is constant in the interior of the flaw, but reduced in value; (b)more » the discontinuity at the boundary of the flaw occurs over the lighted portion of the flaw; (c) this discontinuity is contrasted by a region where ..gamma.. is negative; and (d) new non-physical discontinuities and non-analyticities appear in the reconstructed characteristic function. These general features also appear in numerical calculations which use as input strong scattering data from a spherical void and a flat penny shaped crack in Titanium. The numerical results can be straightforwardly interpreted in terms of the analytical calculation mentioned above, indicating that they will be useful in the study of realistic flaws. We conclude by discussing the stabilization of the aperture limited inversion problem and the removal of non-physical features in the reconstruction.« less

Axisymmetric capillary-gravity waves at the interface of two viscous, immiscible fluids - Initial value problem

NASA Astrophysics Data System (ADS)

Farsoiya, Palas Kumar; Dasgupta, Ratul

2017-11-01

When the interface between two radially unbounded, viscous fluids lying vertically in a stable configuration (denser fluid below) at rest, is perturbed, radially propagating capillary-gravity waves are formed which damp out with time. We study this process analytically using a recently developed linearised theory. For small amplitude initial perturbations, the analytical solution to the initial value problem, represented as a linear superposition of Bessel modes at time t = 0 , is found to agree very well with results obtained from direct numerical simulations of the Navier-Stokes equations, for a range of initial conditions. Our study extends the earlier work by John W. Miles who studied this initial value problem analytically, taking into account, a single viscous fluid only. Implications of this study for the mechanistic understanding of droplet impact into a deep pool, will be discussed. Some preliminary, qualitative comparison with experiments will also be presented. We thank SERB Dept. Science & Technology, Govt. of India, Grant No. EMR/2016/000830 for financial support.

Local box-counting dimensions of discrete quantum eigenvalue spectra: Analytical connection to quantum spectral statistics

NASA Astrophysics Data System (ADS)

Sakhr, Jamal; Nieminen, John M.

2018-03-01

Two decades ago, Wang and Ong, [Phys. Rev. A 55, 1522 (1997)], 10.1103/PhysRevA.55.1522 hypothesized that the local box-counting dimension of a discrete quantum spectrum should depend exclusively on the nearest-neighbor spacing distribution (NNSD) of the spectrum. In this Rapid Communication, we validate their hypothesis by deriving an explicit formula for the local box-counting dimension of a countably-infinite discrete quantum spectrum. This formula expresses the local box-counting dimension of a spectrum in terms of single and double integrals of the NNSD of the spectrum. As applications, we derive an analytical formula for Poisson spectra and closed-form approximations to the local box-counting dimension for spectra having Gaussian orthogonal ensemble (GOE), Gaussian unitary ensemble (GUE), and Gaussian symplectic ensemble (GSE) spacing statistics. In the Poisson and GOE cases, we compare our theoretical formulas with the published numerical data of Wang and Ong and observe excellent agreement between their data and our theory. We also study numerically the local box-counting dimensions of the Riemann zeta function zeros and the alternate levels of GOE spectra, which are often used as numerical models of spectra possessing GUE and GSE spacing statistics, respectively. In each case, the corresponding theoretical formula is found to accurately describe the numerically computed local box-counting dimension.

Stochastic coalescence in finite systems: an algorithm for the numerical solution of the multivariate master equation.

NASA Astrophysics Data System (ADS)

Alfonso, Lester; Zamora, Jose; Cruz, Pedro

2015-04-01

The stochastic approach to coagulation considers the coalescence process going in a system of a finite number of particles enclosed in a finite volume. Within this approach, the full description of the system can be obtained from the solution of the multivariate master equation, which models the evolution of the probability distribution of the state vector for the number of particles of a given mass. Unfortunately, due to its complexity, only limited results were obtained for certain type of kernels and monodisperse initial conditions. In this work, a novel numerical algorithm for the solution of the multivariate master equation for stochastic coalescence that works for any type of kernels and initial conditions is introduced. The performance of the method was checked by comparing the numerically calculated particle mass spectrum with analytical solutions obtained for the constant and sum kernels, with an excellent correspondence between the analytical and numerical solutions. In order to increase the speedup of the algorithm, software parallelization techniques with OpenMP standard were used, along with an implementation in order to take advantage of new accelerator technologies. Simulations results show an important speedup of the parallelized algorithms. This study was funded by a grant from Consejo Nacional de Ciencia y Tecnologia de Mexico SEP-CONACYT CB-131879. The authors also thanks LUFAC® Computacion SA de CV for CPU time and all the support provided.

Novel asymmetric representation method for solving the higher-order Ginzburg-Landau equation

PubMed Central

Wong, Pring; Pang, Lihui; Wu, Ye; Lei, Ming; Liu, Wenjun

2016-01-01

In ultrafast optics, optical pulses are generated to be of shorter pulse duration, which has enormous significance to industrial applications and scientific research. The ultrashort pulse evolution in fiber lasers can be described by the higher-order Ginzburg-Landau (GL) equation. However, analytic soliton solutions for this equation have not been obtained by use of existing methods. In this paper, a novel method is proposed to deal with this equation. The analytic soliton solution is obtained for the first time, and is proved to be stable against amplitude perturbations. Through the split-step Fourier method, the bright soliton solution is studied numerically. The analytic results here may extend the integrable methods, and could be used to study soliton dynamics for some equations in other disciplines. It may also provide the other way to obtain two-soliton solutions for higher-order GL equations. PMID:27086841

Analytical investigation of the faster-is-slower effect with a simplified phenomenological model

NASA Astrophysics Data System (ADS)

Suzuno, K.; Tomoeda, A.; Ueyama, D.

2013-11-01

We investigate the mechanism of the phenomenon called the “faster-is-slower”effect in pedestrian flow studies analytically with a simplified phenomenological model. It is well known that the flow rate is maximized at a certain strength of the driving force in simulations using the social force model when we consider the discharge of self-driven particles through a bottleneck. In this study, we propose a phenomenological and analytical model based on a mechanics-based modeling to reveal the mechanism of the phenomenon. We show that our reduced system, with only a few degrees of freedom, still has similar properties to the original many-particle system and that the effect comes from the competition between the driving force and the nonlinear friction from the model. Moreover, we predict the parameter dependences on the effect from our model qualitatively, and they are confirmed numerically by using the social force model.

Spatiotemporal Airy Ince-Gaussian wave packets in strongly nonlocal nonlinear media.

PubMed

Peng, Xi; Zhuang, Jingli; Peng, Yulian; Li, DongDong; Zhang, Liping; Chen, Xingyu; Zhao, Fang; Deng, Dongmei

2018-03-08

The self-accelerating Airy Ince-Gaussian (AiIG) and Airy helical Ince-Gaussian (AihIG) wave packets in strongly nonlocal nonlinear media (SNNM) are obtained by solving the strongly nonlocal nonlinear Schrödinger equation. For the first time, the propagation properties of three dimensional localized AiIG and AihIG breathers and solitons in the SNNM are demonstrated, these spatiotemporal wave packets maintain the self-accelerating and approximately non-dispersion properties in temporal dimension, periodically oscillating (breather state) or steady (soliton state) in spatial dimension. In particular, their numerical experiments of spatial intensity distribution, numerical simulations of spatiotemporal distribution, as well as the transverse energy flow and the angular momentum in SNNM are presented. Typical examples of the obtained solutions are based on the ratio between the input power and the critical power, the ellipticity and the strong nonlocality parameter. The comparisons of analytical solutions with numerical simulations and numerical experiments of the AiIG and AihIG optical solitons show that the numerical results agree well with the analytical solutions in the case of strong nonlocality.

Q-controlled amplitude modulation atomic force microscopy in liquids: An analysis

NASA Astrophysics Data System (ADS)

Hölscher, H.; Schwarz, U. D.

2006-08-01

An analysis of amplitude modulation atomic force microscopy in liquids is presented with respect to the application of the Q-Control technique. The equation of motion is solved by numerical and analytic methods with and without Q-Control in the presence of a simple model interaction force adequate for many liquid environments. In addition, the authors give an explicit analytical formula for the tip-sample indentation showing that higher Q factors reduce the tip-sample force. It is found that Q-Control suppresses unwanted deformations of the sample surface, leading to the enhanced image quality reported in several experimental studies.

Optimal design of piezoelectric transformers: a rational approach based on an analytical model and a deterministic global optimization.

PubMed

Pigache, Francois; Messine, Frédéric; Nogarede, Bertrand

2007-07-01

This paper deals with a deterministic and rational way to design piezoelectric transformers in radial mode. The proposed approach is based on the study of the inverse problem of design and on its reformulation as a mixed constrained global optimization problem. The methodology relies on the association of the analytical models for describing the corresponding optimization problem and on an exact global optimization software, named IBBA and developed by the second author to solve it. Numerical experiments are presented and compared in order to validate the proposed approach.

FIBER AND INTEGRATED OPTICS. FIBER WAVEGUIDE DEVICES: Stability of solitons in a two-mode fiber waveguide with a group velocity mismatch

NASA Astrophysics Data System (ADS)

Kivshar', Yu S.

1990-12-01

A study is reported of the stability of soliton pulses propagating in a two-mode fiber waveguide under conditions of a mismatch between the group velocities of the optical modes. An analytic explanation is proposed of the dependence of the threshold amplitude of an initial pulse, responsible for intermode locking of the pulses, on the mismatch between the group velocities. An analytically derived dependence is shown to be in good agreement with earlier numerical experiments. Decay of coupled intermode states of solitons due to dissipative losses is predicted.

Analytical studies of the Space Shuttle orbiter nose-gear tire

NASA Technical Reports Server (NTRS)

Noor, Ahmed K.; Tanner, John A.; Peters, Jeanne M.; Robinson, Martha P.

1991-01-01

A computational procedure is presented for evaluating the analytic sensitivity derivatives of the tire response with respect to material and geometrical properties of the tire. The tire is modeled by using a two-dimensional laminated anisotropic shell theory with the effects of variation in material and geometric parameters included. The computational procedure is applied to the case of the Space Shuttle orbiter nose-gear tire subjected to uniform inflation pressure. Numerical results are presented which show the sensitivity of the different tire response quantities to variations in the material characteristics of both the cord and rubber.

Theoretical Assessment of Compressibility Factor of Gases by Using Second Virial Coefficient

NASA Astrophysics Data System (ADS)

Mamedov, Bahtiyar A.; Somuncu, Elif; Askerov, Iskender M.

2018-01-01

We present a new analytical approximation for determining the compressibility factor of real gases at various temperature values. This algorithm is suitable for the accurate evaluation of the compressibility factor using the second virial coefficient with a Lennard-Jones (12-6) potential. Numerical examples are presented for the gases H2, N2, He, CO2, CH4 and air, and the results are compared with other studies in the literature. Our results showed good agreement with the data in the literature. The consistency of the results demonstrates the effectiveness of our analytical approximation for real gases.

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

Zúñiga-Segundo, Arturo; Juárez-Amaro, Raúl; Aguilar-Loreto, Omar

We study the atom–field interaction when the field is in a mixture of coherent states. We show that in this case it is possible to calculate analytically the field entropy for times of the order of twice the collapse time. Such analytical results are done with the help of numerical analysis. We also give an expression in terms of Chebyshev polynomials for power of density matrices. - Highlights: • We calculate the field entropy for times of the order of twice the collapse time. • We give a relation between powers of the density matrices of the subsystems. • Entropymore » operators for both subsystems are obtained.« less

Theoretical study on phase-locking of a radial array CO2 laser

NASA Astrophysics Data System (ADS)

Xu, Yonggen

2014-11-01

The phase-locking of the radial array CO2 laser (RAL) is introduced based on the injection-locking principle. The characteristic parameters of laser beams used in the phase-locking are described, and the coupling coefficient c00 between the injected mode and the eigenmode of RAL is calculated. The laser modes from RAL are the low-order Hermite Gaussian modes due to the diffraction loss. The analytical formula for the output beam through an ABCD optical system is derived according Collins formula. The numerical examples are given to illustrate our analytical results.

Development and application of a three dimensional numerical model for predicting pollutant and sediment transport using an Eulerian-Lagrangian marker particle technique

NASA Technical Reports Server (NTRS)

Pavish, D. L.; Spaulding, M. L.

1977-01-01

A computer coded Lagrangian marker particle in Eulerian finite difference cell solution to the three dimensional incompressible mass transport equation, Water Advective Particle in Cell Technique, WAPIC, was developed, verified against analytic solutions, and subsequently applied in the prediction of long term transport of a suspended sediment cloud resulting from an instantaneous dredge spoil release. Numerical results from WAPIC were verified against analytic solutions to the three dimensional incompressible mass transport equation for turbulent diffusion and advection of Gaussian dye releases in unbounded uniform and uniformly sheared uni-directional flow, and for steady-uniform plug channel flow. WAPIC was utilized to simulate an analytic solution for non-equilibrium sediment dropout from an initially vertically uniform particle distribution in one dimensional turbulent channel flow.

Revisiting analytical solutions for steady interface flow in subsea aquifers: Aquitard salinity effects

NASA Astrophysics Data System (ADS)

Werner, Adrian D.; Robinson, Neville I.

2018-06-01

Existing analytical solutions for the distribution of fresh groundwater in subsea aquifers presume that the overlying offshore aquitard, represented implicitly, contains seawater. Here, we consider the case where offshore fresh groundwater is the result of freshwater discharge from onshore aquifers, and neglect paleo-freshwater sources. A recent numerical modeling investigation, involving explicit simulation of the offshore aquitard, demonstrates that offshore aquitards more likely contain freshwater in areas of upward freshwater leakage to the sea. We integrate this finding into the existing analytical solutions by providing an alternative formulation for steady interface flow in subsea aquifers, whereby the salinity in the offshore aquitard can be chosen. The new solution, taking the aquitard salinity as that of freshwater, provides a closer match to numerical modeling results in which the aquitard is represented explicitly.

Distribution of Steps with Finite-Range Interactions: Analytic Approximations and Numerical Results

NASA Astrophysics Data System (ADS)

GonzáLez, Diego Luis; Jaramillo, Diego Felipe; TéLlez, Gabriel; Einstein, T. L.

2013-03-01

While most Monte Carlo simulations assume only nearest-neighbor steps interact elastically, most analytic frameworks (especially the generalized Wigner distribution) posit that each step elastically repels all others. In addition to the elastic repulsions, we allow for possible surface-state-mediated interactions. We investigate analytically and numerically how next-nearest neighbor (NNN) interactions and, more generally, interactions out to q'th nearest neighbor alter the form of the terrace-width distribution and of pair correlation functions (i.e. the sum over n'th neighbor distribution functions, which we investigated recently.[2] For physically plausible interactions, we find modest changes when NNN interactions are included and generally negligible changes when more distant interactions are allowed. We discuss methods for extracting from simulated experimental data the characteristic scale-setting terms in assumed potential forms.

Steering particles by breaking symmetries

NASA Astrophysics Data System (ADS)

Bet, Bram; Samin, Sela; Georgiev, Rumen; Burak Eral, Huseyin; van Roij, René

2018-06-01

We derive general equations of motions for highly-confined particles that perform quasi-two-dimensional motion in Hele-Shaw channels, which we solve analytically, aiming to derive design principles for self-steering particles. Based on symmetry properties of a particle, its equations of motion can be simplified, where we retrieve an earlier-known equation of motion for the orientation of dimer particles consisting of disks (Uspal et al 2013 Nat. Commun. 4), but now in full generality. Subsequently, these solutions are compared with particle trajectories that are obtained numerically. For mirror-symmetric particles, excellent agreement between the analytical and numerical solutions is found. For particles lacking mirror symmetry, the analytic solutions provide means to classify the motion based on particle geometry, while we find that taking the side-wall interactions into account is important to accurately describe the trajectories.

Nonlinear effects in a plain journal bearing. I - Analytical study. II - Results

NASA Technical Reports Server (NTRS)

Choy, F. K.; Braun, M. J.; Hu, Y.

1991-01-01

In the first part of this work, a numerical model is presented which couples the variable-property Reynolds equation with a rotor-dynamics model for the calculation of a plain journal bearing's nonlinear characteristics when working with a cryogenic fluid, LOX. The effects of load on the linear/nonlinear plain journal bearing characteristics are analyzed and presented in a parametric form. The second part of this work presents numerical results obtained for specific parametric-study input variables (lubricant inlet temperature, external load, angular rotational speed, and axial misalignment). Attention is given to the interrelations between pressure profiles and bearing linear and nonlinear characteristics.

A fully coupled method for massively parallel simulation of hydraulically driven fractures in 3-dimensions: FULLY COUPLED PARALLEL SIMULATION OF HYDRAULIC FRACTURES IN 3-D

DOE PAGES

Settgast, Randolph R.; Fu, Pengcheng; Walsh, Stuart D. C.; ...

2016-09-18

This study describes a fully coupled finite element/finite volume approach for simulating field-scale hydraulically driven fractures in three dimensions, using massively parallel computing platforms. The proposed method is capable of capturing realistic representations of local heterogeneities, layering and natural fracture networks in a reservoir. A detailed description of the numerical implementation is provided, along with numerical studies comparing the model with both analytical solutions and experimental results. The results demonstrate the effectiveness of the proposed method for modeling large-scale problems involving hydraulically driven fractures in three dimensions.

The Thick Level-Set model for dynamic fragmentation

DOE PAGES

Stershic, Andrew J.; Dolbow, John E.; Moës, Nicolas

2017-01-04

The Thick Level-Set (TLS) model is implemented to simulate brittle media undergoing dynamic fragmentation. This non-local model is discretized by the finite element method with damage represented as a continuous field over the domain. A level-set function defines the extent and severity of damage, and a length scale is introduced to limit the damage gradient. Numerical studies in one dimension demonstrate that the proposed method reproduces the rate-dependent energy dissipation and fragment length observations from analytical, numerical, and experimental approaches. In conclusion, additional studies emphasize the importance of appropriate bulk constitutive models and sufficient spatial resolution of the length scale.

A fully coupled method for massively parallel simulation of hydraulically driven fractures in 3-dimensions: FULLY COUPLED PARALLEL SIMULATION OF HYDRAULIC FRACTURES IN 3-D

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

Settgast, Randolph R.; Fu, Pengcheng; Walsh, Stuart D. C.

This study describes a fully coupled finite element/finite volume approach for simulating field-scale hydraulically driven fractures in three dimensions, using massively parallel computing platforms. The proposed method is capable of capturing realistic representations of local heterogeneities, layering and natural fracture networks in a reservoir. A detailed description of the numerical implementation is provided, along with numerical studies comparing the model with both analytical solutions and experimental results. The results demonstrate the effectiveness of the proposed method for modeling large-scale problems involving hydraulically driven fractures in three dimensions.

Effects of Energy Chirp on Echo-Enabled Harmonic Generation Free-Electron Lasers

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

Huang, Z.; Ratner, D.; Stupakov, G.

2009-02-23

We study effects of energy chirp on echo-enabled harmonic generation (EEHG). Analytical expressions are compared with numerical simulations for both harmonic and bunching factors. We also discuss the EEHG free-electron laser bandwidth increase due to an energy-modulated beam and its pulse length dependence on the electron energy chirp.

Classroom Learning and Achievement: How the Complexity of Classroom Interaction Impacts Students' Learning

ERIC Educational Resources Information Center

Podschuweit, Sören; Bernholt, Sascha; Brückmann, Maja

2016-01-01

Background: Complexity models have provided a suitable framework in various domains to assess students' educational achievement. Complexity is often used as the analytical focus when regarding learning outcomes, i.e. when analyzing written tests or problem-centered interviews. Numerous studies reveal negative correlations between the complexity of…

Mean motion resonances. [of asteroid belt structure

NASA Technical Reports Server (NTRS)

Froeschle, CL.; Greenberg, R.

1989-01-01

Recent research on the resonant structure of the asteroid belt is reviewed. The resonant mechanism is discussed, and analytical models for the study of mean motion resonances are examined. Numerical averaging methods and mapping methods are considered. It is shown how fresh insight can be obtained by means of a new semianalytical approach.

A Multi-State Factor-Analytic and Psychometric Meta-Analysis of Agricultural Mechanics Laboratory Management Competencies

ERIC Educational Resources Information Center

McKim, Billy R.; Saucier, P. Ryan

2012-01-01

For more than 20 years, the 50 agricultural mechanics laboratory management competencies identified by Johnson and Schumacher in 1989 have served as the basis for numerous needs assessments of secondary agriculture teachers. This study reevaluated Johnson and Schumacher's instrument, as modified by Saucier, Schumacher, Funkenbusch, Terry, and…

Using Meta-Regression to Explore Moderating Effects in Surveys of International Achievement

ERIC Educational Resources Information Center

Benton, Tom

2014-01-01

This article demonstrates how meta-analytic techniques, that have typically been used to synthesize findings across numerous studies, can also be applied to examine the reasons why relationships between background characteristics and outcomes may vary across different locations in a single multi-site survey. This application is particularly…

Feedback topology and XOR-dynamics in Boolean networks with varying input structure

NASA Astrophysics Data System (ADS)

Ciandrini, L.; Maffi, C.; Motta, A.; Bassetti, B.; Cosentino Lagomarsino, M.

2009-08-01

We analyze a model of fixed in-degree random Boolean networks in which the fraction of input-receiving nodes is controlled by the parameter γ . We investigate analytically and numerically the dynamics of graphs under a parallel XOR updating scheme. This scheme is interesting because it is accessible analytically and its phenomenology is at the same time under control and as rich as the one of general Boolean networks. We give analytical formulas for the dynamics on general graphs, showing that with a XOR-type evolution rule, dynamic features are direct consequences of the topological feedback structure, in analogy with the role of relevant components in Kauffman networks. Considering graphs with fixed in-degree, we characterize analytically and numerically the feedback regions using graph decimation algorithms (Leaf Removal). With varying γ , this graph ensemble shows a phase transition that separates a treelike graph region from one in which feedback components emerge. Networks near the transition point have feedback components made of disjoint loops, in which each node has exactly one incoming and one outgoing link. Using this fact, we provide analytical estimates of the maximum period starting from topological considerations.

A numerical test of the topographic bias

NASA Astrophysics Data System (ADS)

Sjöberg, L. E.; Joud, M. S. S.

2018-02-01

In 1962 A. Bjerhammar introduced the method of analytical continuation in physical geodesy, implying that surface gravity anomalies are downward continued into the topographic masses down to an internal sphere (the Bjerhammar sphere). The method also includes analytical upward continuation of the potential to the surface of the Earth to obtain the quasigeoid. One can show that also the common remove-compute-restore technique for geoid determination includes an analytical continuation as long as the complete density distribution of the topography is not known. The analytical continuation implies that the downward continued gravity anomaly and/or potential are/is in error by the so-called topographic bias, which was postulated by a simple formula of L E Sjöberg in 2007. Here we will numerically test the postulated formula by comparing it with the bias obtained by analytical downward continuation of the external potential of a homogeneous ellipsoid to an inner sphere. The result shows that the postulated formula holds: At the equator of the ellipsoid, where the external potential is downward continued 21 km, the computed and postulated topographic biases agree to less than a millimetre (when the potential is scaled to the unit of metre).

Feedback topology and XOR-dynamics in Boolean networks with varying input structure.

PubMed

Ciandrini, L; Maffi, C; Motta, A; Bassetti, B; Cosentino Lagomarsino, M

2009-08-01

We analyze a model of fixed in-degree random Boolean networks in which the fraction of input-receiving nodes is controlled by the parameter gamma. We investigate analytically and numerically the dynamics of graphs under a parallel XOR updating scheme. This scheme is interesting because it is accessible analytically and its phenomenology is at the same time under control and as rich as the one of general Boolean networks. We give analytical formulas for the dynamics on general graphs, showing that with a XOR-type evolution rule, dynamic features are direct consequences of the topological feedback structure, in analogy with the role of relevant components in Kauffman networks. Considering graphs with fixed in-degree, we characterize analytically and numerically the feedback regions using graph decimation algorithms (Leaf Removal). With varying gamma , this graph ensemble shows a phase transition that separates a treelike graph region from one in which feedback components emerge. Networks near the transition point have feedback components made of disjoint loops, in which each node has exactly one incoming and one outgoing link. Using this fact, we provide analytical estimates of the maximum period starting from topological considerations.

The PAC-MAN model: Benchmark case for linear acoustics in computational physics

NASA Astrophysics Data System (ADS)

Ziegelwanger, Harald; Reiter, Paul

2017-10-01

Benchmark cases in the field of computational physics, on the one hand, have to contain a certain complexity to test numerical edge cases and, on the other hand, require the existence of an analytical solution, because an analytical solution allows the exact quantification of the accuracy of a numerical simulation method. This dilemma causes a need for analytical sound field formulations of complex acoustic problems. A well known example for such a benchmark case for harmonic linear acoustics is the ;Cat's Eye model;, which describes the three-dimensional sound field radiated from a sphere with a missing octant analytically. In this paper, a benchmark case for two-dimensional (2D) harmonic linear acoustic problems, viz., the ;PAC-MAN model;, is proposed. The PAC-MAN model describes the radiated and scattered sound field around an infinitely long cylinder with a cut out sector of variable angular width. While the analytical calculation of the 2D sound field allows different angular cut-out widths and arbitrarily positioned line sources, the computational cost associated with the solution of this problem is similar to a 1D problem because of a modal formulation of the sound field in the PAC-MAN model.

Shock compression modeling of metallic single crystals: comparison of finite difference, steady wave, and analytical solutions

DOE PAGES

Lloyd, Jeffrey T.; Clayton, John D.; Austin, Ryan A.; ...

2015-07-10

Background: The shock response of metallic single crystals can be captured using a micro-mechanical description of the thermoelastic-viscoplastic material response; however, using a such a description within the context of traditional numerical methods may introduce a physical artifacts. Advantages and disadvantages of complex material descriptions, in particular the viscoplastic response, must be framed within approximations introduced by numerical methods. Methods: Three methods of modeling the shock response of metallic single crystals are summarized: finite difference simulations, steady wave simulations, and algebraic solutions of the Rankine-Hugoniot jump conditions. For the former two numerical techniques, a dislocation density based framework describes themore » rate- and temperature-dependent shear strength on each slip system. For the latter analytical technique, a simple (two-parameter) rate- and temperature-independent linear hardening description is necessarily invoked to enable simultaneous solution of the governing equations. For all models, the same nonlinear thermoelastic energy potential incorporating elastic constants of up to order 3 is applied. Results: Solutions are compared for plate impact of highly symmetric orientations (all three methods) and low symmetry orientations (numerical methods only) of aluminum single crystals shocked to 5 GPa (weak shock regime) and 25 GPa (overdriven regime). Conclusions: For weak shocks, results of the two numerical methods are very similar, regardless of crystallographic orientation. For strong shocks, artificial viscosity affects the finite difference solution, and effects of transverse waves for the lower symmetry orientations not captured by the steady wave method become important. The analytical solution, which can only be applied to highly symmetric orientations, provides reasonable accuracy with regards to prediction of most variables in the final shocked state but, by construction, does not provide insight into the shock structure afforded by the numerical methods.« less

Quantitative analysis of time-resolved microwave conductivity data

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

Reid, Obadiah G.; Moore, David T.; Li, Zhen

Flash-photolysis time-resolved microwave conductivity (fp-TRMC) is a versatile, highly sensitive technique for studying the complex photoconductivity of solution, solid, and gas-phase samples. The purpose of this paper is to provide a standard reference work for experimentalists interested in using microwave conductivity methods to study functional electronic materials, describing how to conduct and calibrate these experiments in order to obtain quantitative results. The main focus of the paper is on calculating the calibration factor, K, which is used to connect the measured change in microwave power absorption to the conductance of the sample. We describe the standard analytical formulae that havemore » been used in the past, and compare them to numerical simulations. This comparison shows that the most widely used analytical analysis of fp-TRMC data systematically under-estimates the transient conductivity by ~60%. We suggest a more accurate semi-empirical way of calibrating these experiments. However, we emphasize that the full numerical calculation is necessary to quantify both transient and steady-state conductance for arbitrary sample properties and geometry.« less

Quantitative analysis of time-resolved microwave conductivity data

DOE PAGES

Reid, Obadiah G.; Moore, David T.; Li, Zhen; ...

2017-11-10

Flash-photolysis time-resolved microwave conductivity (fp-TRMC) is a versatile, highly sensitive technique for studying the complex photoconductivity of solution, solid, and gas-phase samples. The purpose of this paper is to provide a standard reference work for experimentalists interested in using microwave conductivity methods to study functional electronic materials, describing how to conduct and calibrate these experiments in order to obtain quantitative results. The main focus of the paper is on calculating the calibration factor, K, which is used to connect the measured change in microwave power absorption to the conductance of the sample. We describe the standard analytical formulae that havemore » been used in the past, and compare them to numerical simulations. This comparison shows that the most widely used analytical analysis of fp-TRMC data systematically under-estimates the transient conductivity by ~60%. We suggest a more accurate semi-empirical way of calibrating these experiments. However, we emphasize that the full numerical calculation is necessary to quantify both transient and steady-state conductance for arbitrary sample properties and geometry.« less

Alpha particles diffusion due to charge changes

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

Clauser, C. F., E-mail: cesar.clauser@ib.edu.ar; Farengo, R.

2015-12-15

Alpha particles diffusion due to charge changes in a magnetized plasma is studied. Analytical calculations and numerical simulations are employed to show that this process can be very important in the pedestal-edge-SOL regions. This is the first study that presents clear evidence of the importance of atomic processes on the diffusion of alpha particles. A simple 1D model that includes inelastic collisions with plasma species, “cold” neutrals, and partially ionized species was employed. The code, which follows the exact particle orbits and includes the effect of inelastic collisions via a Monte Carlo type random process, runs on a graphic processormore » unit (GPU). The analytical and numerical results show excellent agreement when a uniform background (plasma and cold species) is assumed. The simulations also show that the gradients in the density of the plasma and cold species, which are large and opposite in the edge region, produce an inward flux of alpha particles. Calculations of the alpha particles flux reaching the walls or divertor plates should include these processes.« less

Development, integrated investigation, laboratory and in-flight testing of Chibis-M microsatellite ADCS

NASA Astrophysics Data System (ADS)

Ovchinnikov, M. Yu.; Ivanov, D. S.; Ivlev, N. A.; Karpenko, S. O.; Roldugin, D. S.; Tkachev, S. S.

2014-01-01

Design, analytical investigation, laboratory and in-flight testing of the attitude determination and control system (ADCS) of a microsatellites are considered. The system consists of three pairs of reaction wheels, three magnetorquers, a set of Sun sensors, a three-axis magnetometer and a control unit. The ADCS is designed for a small 10-50 kg LEO satellite. System development is accomplished in several steps: satellite dynamics preliminary study using asymptotical and numerical techniques, hardware and software design, laboratory testing of each actuator and sensor and the whole ADCS. Laboratory verification is carried out on the specially designed test-bench. In-flight ADCS exploitation results onboard the Russian microsatellite "Chibis-M" are presented. The satellite was developed, designed and manufactured by the Institute of Space Research of RAS. "Chibis-M" was launched by the "Progress-13M" cargo vehicle on January 25, 2012 after undocking from the International Space Station (ISS). This paper assess both the satellite and the ADCS mock-up dynamics. Analytical, numerical and laboratory study results are in good correspondence with in-flight data.

Effect of particle inertia on turbulence in a suspension.

PubMed

L'vov, Victor S; Ooms, Gijs; Pomyalov, Anna

2003-04-01

We propose a one-fluid analytical model for a turbulently flowing dilute suspension, based on a modified Navier-Stokes equation with a k-dependent effective density of suspension rho(eff)(k) and an additional damping term proportional, variant gamma(p)(k), representing the fluid-particle friction (described by Stokes law). The statistical description of turbulence within the model is simplified by a modification of the usual closure procedure based on the Richardson-Kolmogorov picture of turbulence with a differential approximation for the energy transfer term. The resulting ordinary differential equation for the energy budget is solved analytically for various important limiting cases and numerically in the general case. In the inertial interval of scales, we describe analytically two competing effects: the energy suppression due to the fluid-particle friction and the energy enhancement during the cascade process due to decrease of the effective density of the small-scale motions. An additional suppression or enhancement of the energy density may occur in the viscous subrange, caused by the variation of the extent of the inertial interval due to the combined effect of the fluid-particle friction and the decrease of the kinematic viscosity of the suspensions. The analytical description of the complicated interplay of these effects supported by numerical calculations is presented. Our findings allow one to rationalize the qualitative picture of the isotropic homogeneous turbulence of dilute suspensions as observed in direct numerical simulations.

Analytical solution and numerical simulation of the liquid nitrogen freezing-temperature field of a single pipe

NASA Astrophysics Data System (ADS)

Cai, Haibing; Xu, Liuxun; Yang, Yugui; Li, Longqi

2018-05-01

Artificial liquid nitrogen freezing technology is widely used in urban underground engineering due to its technical advantages, such as simple freezing system, high freezing speed, low freezing temperature, high strength of frozen soil, and absence of pollution. However, technical difficulties such as undefined range of liquid nitrogen freezing and thickness of frozen wall gradually emerge during the application process. Thus, the analytical solution of the freezing-temperature field of a single pipe is established considering the freezing temperature of soil and the constant temperature of freezing pipe wall. This solution is then applied in a liquid nitrogen freezing project. Calculation results show that the radius of freezing front of liquid nitrogen is proportional to the square root of freezing time. The radius of the freezing front also decreases with decreased the freezing temperature, and the temperature gradient of soil decreases with increased distance from the freezing pipe. The radius of cooling zone in the unfrozen area is approximately four times the radius of the freezing front. Meanwhile, the numerical simulation of the liquid nitrogen freezing-temperature field of a single pipe is conducted using the Abaqus finite-element program. Results show that the numerical simulation of soil temperature distribution law well agrees with the analytical solution, further verifies the reliability of the established analytical solution of the liquid nitrogen freezing-temperature field of a single pipe.

An explicit closed-form analytical solution for European options under the CGMY model

NASA Astrophysics Data System (ADS)

Chen, Wenting; Du, Meiyu; Xu, Xiang

2017-01-01

In this paper, we consider the analytical pricing of European path-independent options under the CGMY model, which is a particular type of pure jump Le´vy process, and agrees well with many observed properties of the real market data by allowing the diffusions and jumps to have both finite and infinite activity and variation. It is shown that, under this model, the option price is governed by a fractional partial differential equation (FPDE) with both the left-side and right-side spatial-fractional derivatives. In comparison to derivatives of integer order, fractional derivatives at a point not only involve properties of the function at that particular point, but also the information of the function in a certain subset of the entire domain of definition. This ;globalness; of the fractional derivatives has added an additional degree of difficulty when either analytical methods or numerical solutions are attempted. Albeit difficult, we still have managed to derive an explicit closed-form analytical solution for European options under the CGMY model. Based on our solution, the asymptotic behaviors of the option price and the put-call parity under the CGMY model are further discussed. Practically, a reliable numerical evaluation technique for the current formula is proposed. With the numerical results, some analyses of impacts of four key parameters of the CGMY model on European option prices are also provided.

Application of Semi-analytical Satellite Theory orbit propagator to orbit determination for space object catalog maintenance

NASA Astrophysics Data System (ADS)

Setty, Srinivas J.; Cefola, Paul J.; Montenbruck, Oliver; Fiedler, Hauke

2016-05-01

Catalog maintenance for Space Situational Awareness (SSA) demands an accurate and computationally lean orbit propagation and orbit determination technique to cope with the ever increasing number of observed space objects. As an alternative to established numerical and analytical methods, we investigate the accuracy and computational load of the Draper Semi-analytical Satellite Theory (DSST). The standalone version of the DSST was enhanced with additional perturbation models to improve its recovery of short periodic motion. The accuracy of DSST is, for the first time, compared to a numerical propagator with fidelity force models for a comprehensive grid of low, medium, and high altitude orbits with varying eccentricity and different inclinations. Furthermore, the run-time of both propagators is compared as a function of propagation arc, output step size and gravity field order to assess its performance for a full range of relevant use cases. For use in orbit determination, a robust performance of DSST is demonstrated even in the case of sparse observations, which is most sensitive to mismodeled short periodic perturbations. Overall, DSST is shown to exhibit adequate accuracy at favorable computational speed for the full set of orbits that need to be considered in space surveillance. Along with the inherent benefits of a semi-analytical orbit representation, DSST provides an attractive alternative to the more common numerical orbit propagation techniques.

Validation of the SINDA/FLUINT code using several analytical solutions

NASA Technical Reports Server (NTRS)

Keller, John R.

1995-01-01

The Systems Improved Numerical Differencing Analyzer and Fluid Integrator (SINDA/FLUINT) code has often been used to determine the transient and steady-state response of various thermal and fluid flow networks. While this code is an often used design and analysis tool, the validation of this program has been limited to a few simple studies. For the current study, the SINDA/FLUINT code was compared to four different analytical solutions. The thermal analyzer portion of the code (conduction and radiative heat transfer, SINDA portion) was first compared to two separate solutions. The first comparison examined a semi-infinite slab with a periodic surface temperature boundary condition. Next, a small, uniform temperature object (lumped capacitance) was allowed to radiate to a fixed temperature sink. The fluid portion of the code (FLUINT) was also compared to two different analytical solutions. The first study examined a tank filling process by an ideal gas in which there is both control volume work and heat transfer. The final comparison considered the flow in a pipe joining two infinite reservoirs of pressure. The results of all these studies showed that for the situations examined here, the SINDA/FLUINT code was able to match the results of the analytical solutions.

Coulomb double helical structure

NASA Astrophysics Data System (ADS)

Kamimura, Tetsuo; Ishihara, Osamu

2012-01-01

Structures of Coulomb clusters formed by dust particles in a plasma are studied by numerical simulation. Our study reveals the presence of various types of self-organized structures of a cluster confined in a prolate spheroidal electrostatic potential. The stable configurations depend on a prolateness parameter for the confining potential as well as on the number of dust particles in a cluster. One-dimensional string, two-dimensional zigzag structure and three-dimensional double helical structure are found as a result of the transition controlled by the prolateness parameter. The formation of stable double helical structures resulted from the transition associated with the instability of angular perturbations on double strings. Analytical perturbation study supports the findings of numerical simulations.

Analytical solutions of one-dimensional multispecies reactive transport in a permeable reactive barrier-aquifer system

NASA Astrophysics Data System (ADS)

Mieles, John; Zhan, Hongbin

2012-06-01

The permeable reactive barrier (PRB) remediation technology has proven to be more cost-effective than conventional pump-and-treat systems, and has demonstrated the ability to rapidly reduce the concentrations of specific chemicals of concern (COCs) by up to several orders of magnitude in some scenarios. This study derives new steady-state analytical solutions to multispecies reactive transport in a PRB-aquifer (dual domain) system. The advantage of the dual domain model is that it can account for the potential existence of natural degradation in the aquifer, when designing the required PRB thickness. The study focuses primarily on the steady-state analytical solutions of the tetrachloroethene (PCE) serial degradation pathway and secondly on the analytical solutions of the parallel degradation pathway. The solutions in this study can also be applied to other types of dual domain systems with distinct flow and transport properties. The steady-state analytical solutions are shown to be accurate and the numerical program RT3D is selected for comparison. The results of this study are novel in that the solutions provide improved modeling flexibility including: 1) every species can have unique first-order reaction rates and unique retardation factors, and 2) daughter species can be modeled with their individual input concentrations or solely as byproducts of the parent species. The steady-state analytical solutions exhibit a limitation that occurs when interspecies reaction rate factors equal each other, which result in undefined solutions. Excel spreadsheet programs were created to facilitate prompt application of the steady-state analytical solutions, for both the serial and parallel degradation pathways.