Dissipative nonlinear waves in a gravitating quantum fluid

NASA Astrophysics Data System (ADS)

Sahu, Biswajit; Sinha, Anjana; Roychoudhury, Rajkumar

2018-02-01

Nonlinear wave propagation is studied in a dissipative, self-gravitating Bose-Einstein condensate, starting from the Gross-Pitaevskii equation. In the absence of an exact analytical result, approximate methods like the linear analysis and perturbative approach are applied. The linear dispersion relation puts a restriction on the permissible range of the dissipation parameter. The waves get damped due to dissipation. The small amplitude analysis using reductive perturbation technique is found to yield a modified form of KdV equation, which is solved both analytically as well as numerically. Interestingly, the analytical and numerical plots match excellently with each other, in the realm of weak dissipation.

Rayleigh approximation to ground state of the Bose and Coulomb glasses

DOE PAGES

Ryan, S. D.; Mityushev, V.; Vinokur, V. M.; ...

2015-01-16

Glasses are rigid systems in which competing interactions prevent simultaneous minimization of local energies. This leads to frustration and highly degenerate ground states the nature and properties of which are still far from being thoroughly understood. We report an analytical approach based on the method of functional equations that allows us to construct the Rayleigh approximation to the ground state of a two-dimensional (2D) random Coulomb system with logarithmic interactions. We realize a model for 2D Coulomb glass as a cylindrical type II superconductor containing randomly located columnar defects (CD) which trap superconducting vortices induced by applied magnetic field. Ourmore » findings break ground for analytical studies of glassy systems, marking an important step towards understanding their properties.« less

Review and assessment of the database and numerical modeling for turbine heat transfer

NASA Technical Reports Server (NTRS)

Gladden, H. J.; Simoneau, R. J.

1988-01-01

The objectives of the HOST Turbine Heat Transfer subproject were to obtain a better understanding of the physics of the aerothermodynamic phenomena and to assess and improve the analytical methods used to predict the flow and heat transfer in high-temperature gas turbines. At the time the HOST project was initiated, an across-the-board improvement in turbine design technology was needed. A building-block approach was utilized and the research ranged from the study of fundamental phenomena and modeling to experiments in simulated real engine environments. Experimental research accounted for approximately 75 percent of the funding while the analytical efforts were approximately 25 percent. A healthy government/industry/university partnership, with industry providing almost half of the research, was created to advance the turbine heat transfer design technology base.

Evaluation of selected methods for determining streamflow during periods of ice effect

USGS Publications Warehouse

Melcher, Norwood B.; Walker, J.F.

1992-01-01

Seventeen methods for estimating ice-affected streamflow are evaluated for potential use with the U.S. Geological Survey streamflow-gaging station network. The methods evaluated were identified by written responses from U.S. Geological Survey field offices and by a comprehensive literature search. The methods selected and techniques used for applying the methods are described in this report. The methods are evaluated by comparing estimated results with data collected at three streamflow-gaging stations in Iowa during the winter of 1987-88. Discharge measurements were obtained at 1- to 5-day intervals during the ice-affected periods at the three stations to define an accurate baseline record. Discharge records were compiled for each method based on data available, assuming a 6-week field schedule. The methods are classified into two general categories-subjective and analytical--depending on whether individual judgment is necessary for method application. On the basis of results of the evaluation for the three Iowa stations, two of the subjective methods (discharge ratio and hydrographic-and-climatic comparison) were more accurate than the other subjective methods and approximately as accurate as the best analytical method. Three of the analytical methods (index velocity, adjusted rating curve, and uniform flow) could potentially be used at streamflow-gaging stations, where the need for accurate ice-affected discharge estimates justifies the expense of collecting additional field data. One analytical method (ice-adjustment factor) may be appropriate for use at stations with extremely stable stage-discharge ratings and measuring sections. Further research is needed to refine the analytical methods. The discharge-ratio and multiple-regression methods produce estimates of streamflow for varying ice conditions using information obtained from the existing U.S. Geological Survey streamflow-gaging network.

Reproduction of exact solutions of Lipkin model by nonlinear higher random-phase approximation

NASA Astrophysics Data System (ADS)

Terasaki, J.; Smetana, A.; Šimkovic, F.; Krivoruchenko, M. I.

2017-10-01

It is shown that the random-phase approximation (RPA) method with its nonlinear higher generalization, which was previously considered as approximation except for a very limited case, reproduces the exact solutions of the Lipkin model. The nonlinear higher RPA is based on an equation nonlinear on eigenvectors and includes many-particle-many-hole components in the creation operator of the excited states. We demonstrate the exact character of solutions analytically for the particle number N = 2 and numerically for N = 8. This finding indicates that the nonlinear higher RPA is equivalent to the exact Schrödinger equation.

Atmospheric guidance law for planar skip trajectories

NASA Technical Reports Server (NTRS)

Mease, K. D.; Mccreary, F. A.

1985-01-01

The applicability of an approximate, closed-form, analytical solution to the equations of motion, as a basis for a deterministic guidance law for controlling the in-plane motion during a skip trajectory, is investigated. The derivation of the solution by the method of matched asymptotic expansions is discussed. Specific issues that arise in the application of the solution to skip trajectories are addressed. Based on the solution, an explicit formula for the approximate energy loss due to an atmospheric pass is derived. A guidance strategy is proposed that illustrates the use of the approximate solution. A numerical example shows encouraging performance.

Photogrammetrically Measured Distortions of Composite Structure Microwave Reflectors at Approximately 90 K

NASA Technical Reports Server (NTRS)

Mule, Peter; Hill, Michael D.; Sampler, Henry P.

2000-01-01

The Microwave Anisotropy Probe (MAP) Observatory, scheduled for a fall 2000 launch, is designed to measure temperature fluctuations (anisotropy) and produce a high sensitivity and high spatial resolution (better than 0.3 deg.) map of the cosmic microwave background (CMB) radiation over the entire sky between 22 and 90 GHz. MAP utilizes back-to-back composite Gregorian telescopes supported on a composite truss structure to focus the microwave signals into 10 differential microwave receivers. Proper position and shape of the telescope reflectors at the operating temperature of approximately 90 K is a critical element to ensuring mission success. We describe the methods and analysis used to validate the in-flight position and shape predictions for the reflectors based on photogrammetric (PG) metrology data taken under vacuum with the reflectors at approximately 90 K. Contour maps showing reflector distortion analytical extrapolations were generated. The resulting reflector distortion data are shown to be crucial to the analytical assessment of the MAP instrument's microwave system in-flight performance.

Atomic density functional and diagram of structures in the phase field crystal model

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

Ankudinov, V. E., E-mail: vladimir@ankudinov.org; Galenko, P. K.; Kropotin, N. V.

2016-02-15

The phase field crystal model provides a continual description of the atomic density over the diffusion time of reactions. We consider a homogeneous structure (liquid) and a perfect periodic crystal, which are constructed from the one-mode approximation of the phase field crystal model. A diagram of 2D structures is constructed from the analytic solutions of the model using atomic density functionals. The diagram predicts equilibrium atomic configurations for transitions from the metastable state and includes the domains of existence of homogeneous, triangular, and striped structures corresponding to a liquid, a body-centered cubic crystal, and a longitudinal cross section of cylindricalmore » tubes. The method developed here is employed for constructing the diagram for the homogeneous liquid phase and the body-centered iron lattice. The expression for the free energy is derived analytically from density functional theory. The specific features of approximating the phase field crystal model are compared with the approximations and conclusions of the weak crystallization and 2D melting theories.« less

Development of an electrothermal vaporization ICP-MS method and assessment of its applicability to studies of the homogeneity of reference materials.

PubMed

Friese, K C; Grobecker, K H; Wätjen, U

2001-07-01

A method has been developed for measurement of the homogeneity of analyte distribution in powdered materials by use of electrothermal vaporization with inductively coupled plasma mass spectrometric (ETV-ICP-MS) detection. The method enabled the simultaneous determination of As, Cd, Cu, Fe, Mn, Pb, and Zn in milligram amounts of samples of biological origin. The optimized conditions comprised a high plasma power of 1,500 W, reduced aerosol transport flow, and heating ramps below 300 degrees C s(-1). A temperature ramp to 550 degrees C ensured effective pyrolysis of approximately 70% of the organic compounds without losses of analyte. An additional hold stage at 700 degrees C led to separation of most of the analyte signals from the evaporation of carbonaceous matrix compounds. The effect of time resolution of signal acquisition on the precision of the ETV measurements was investigated. An increase in the number of masses monitored up to 20 is possible with not more than 1% additional relative standard deviation of results caused by limited temporal resolution of the transient signals. Recording of signals from the nebulization of aqueous standards in each sample run enabled correction for drift of the sensitivity of the ETV-ICP-MS instrument. The applicability of the developed method to homogeneity studies was assessed by use of four certified reference materials. According to the best repeatability observed in these sample runs, the maximum contribution of the method to the standard deviation is approximately 5% to 6% for all the elements investigated.

Numerical and analytic models of spontaneous frequency sweeping for energetic particle-driven Alfven eigenmodes

NASA Astrophysics Data System (ADS)

Wang, Ge; Berk, H. L.

2011-10-01

The frequency chirping signal arising from spontaneous a toroidial Alfven eigenmode (TAE) excited by energetic particles is studied for both numerical and analytic models. The time-dependent numerical model is based on the 1D Vlasov equation. We use a sophisticated tracking method to lock onto the resonant structure to enable the chirping frequency to be nearly constant in the calculation frame. The accuracy of the adiabatic approximation is tested during the simulation which justifies the appropriateness of our analytic model. The analytic model uses the adiabatic approximation which allows us to solve the wave evolution equation in frequency space. Then, the resonant interactions between energetic particles and TAE yield predictions for the chirping rate, wave frequency and amplitudes vs. time. Here, an adiabatic invariant J is defined on the separatrix of a chirping mode to determine the region of confinement of the wave trapped distribution function. We examine the asymptotic behavior of the chirping signal for its long time evolution and find agreement in essential features with the results of the simulation. Work supported by Department of Energy contract DE-FC02-08ER54988.

Method development and validation for simultaneous quantification of 15 drugs of abuse and prescription drugs and 7 of their metabolites in whole blood relevant in the context of driving under the influence of drugs--usefulness of multi-analyte calibration.

PubMed

Steuer, Andrea E; Forss, Anna-Maria; Dally, Annika M; Kraemer, Thomas

2014-11-01

In the context of driving under the influence of drugs (DUID), not only common drugs of abuse may have an influence, but also medications with similar mechanisms of action. Simultaneous quantification of a variety of drugs and medications relevant in this context allows faster and more effective analyses. Therefore, multi-analyte approaches have gained more and more popularity in recent years. Usually, calibration curves for such procedures contain a mixture of all analytes, which might lead to mutual interferences. In this study we investigated whether the use of such mixtures leads to reliable results for authentic samples containing only one or two analytes. Five hundred microliters of whole blood were extracted by routine solid-phase extraction (SPE, HCX). Analysis was performed on an ABSciex 3200 QTrap instrument with ESI+ in scheduled MRM mode. The method was fully validated according to international guidelines including selectivity, recovery, matrix effects, accuracy and precision, stabilities, and limit of quantification. The selected SPE provided recoveries >60% for all analytes except 6-monoacetylmorphine (MAM) with coefficients of variation (CV) below 15% or 20% for quality controls (QC) LOW and HIGH, respectively. Ion suppression >30% was found for benzoylecgonine, hydrocodone, hydromorphone, MDA, oxycodone, and oxymorphone at QC LOW, however CVs were always below 10% (n=6 different whole blood samples). Accuracy and precision criteria were fulfilled for all analytes except for MAM. Systematic investigation of accuracy determined for QC MED in a multi-analyte mixture compared to samples containing only single analytes revealed no relevant differences for any analyte, indicating that a multi-analyte calibration is suitable for the presented method. Comparison of approximately 60 samples to a former GC-MS method showed good correlation. The newly validated method was successfully applied to more than 1600 routine samples and 3 proficiency tests. Copyright © 2014 Elsevier Ireland Ltd. All rights reserved.

Electron scattering intensities and Patterson functions of Skyrmions

NASA Astrophysics Data System (ADS)

Karliner, M.; King, C.; Manton, N. S.

2016-06-01

The scattering of electrons off nuclei is one of the best methods of probing nuclear structure. In this paper we focus on electron scattering off nuclei with spin and isospin zero within the Skyrme model. We consider two distinct methods and simplify our calculations by use of the Born approximation. The first method is to calculate the form factor of the spherically averaged Skyrmion charge density; the second uses the Patterson function to calculate the scattering intensity off randomly oriented Skyrmions, and spherically averages at the end. We compare our findings with experimental scattering data. We also find approximate analytical formulae for the first zero and first stationary point of a form factor.

BLUES function method in computational physics

NASA Astrophysics Data System (ADS)

Indekeu, Joseph O.; Müller-Nedebock, Kristian K.

2018-04-01

We introduce a computational method in physics that goes ‘beyond linear use of equation superposition’ (BLUES). A BLUES function is defined as a solution of a nonlinear differential equation (DE) with a delta source that is at the same time a Green’s function for a related linear DE. For an arbitrary source, the BLUES function can be used to construct an exact solution to the nonlinear DE with a different, but related source. Alternatively, the BLUES function can be used to construct an approximate piecewise analytical solution to the nonlinear DE with an arbitrary source. For this alternative use the related linear DE need not be known. The method is illustrated in a few examples using analytical calculations and numerical computations. Areas for further applications are suggested.

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.

Determination of synthetic ferric chelates used as fertilizers by liquid chromatography-electrospray/mass spectrometry in agricultural matrices.

PubMed

Alvarez-Fernández, Ana; Orera, Irene; Abadía, Javier; Abadía, Anunciación

2007-01-01

A high-performance liquid chromatography-electrospray ionization/mass spectrometry (time of flight) method has been developed for the simultaneous determination of synthetic Fe(III)-chelates used as fertilizers. Analytes included the seven major Fe(III)-chelates used in agriculture, Fe(III)-EDTA, Fe(III)-DTPA, Fe(III)-HEDTA, Fe(III)-CDTA, Fe(III)-o,oEDDHA, Fe(III)-o,pEDDHA, and Fe(III)-EDDHMA, and the method was validated using isotope labeled (57)Fe(III)-chelates as internal standards. Calibration curves had R values in the range 0.9962-0.9997. Limits of detection and quantification were in the ranges 3-164 and 14-945 pmol, respectively. Analyte concentrations could be determined between the limits of quantification and 25 muM (racemic and meso Fe(III)-o,oEDDHA and Fe(III)-EDDHMA) or 50 muM (Fe(III)-EDTA, Fe(III)-HEDTA, Fe(III)-DTPA, Fe(III)-CDTA and Fe(III)-o,pEDDHA). The average intraday repeatability values were approximately 0.5 and 5% for retention time and peak area, respectively, whereas the interday repeatability values were approximately 0.7 and 8% for retention time and peak area, respectively. The method was validated using four different agricultural matrices, including nutrient solution, irrigation water, soil solution, and plant xylem exudates, spiked with Fe(III)-chelate standards and their stable isotope-labeled corresponding chelates. Analyte recoveries found were in the ranges 92-101% (nutrient solution), 89-102% (irrigation water), 82-100% (soil solution), and 70-111% (plant xylem exudates). Recoveries depended on the analyte, with Fe(III)-EDTA and Fe(III)-DTPA showing the lowest recoveries (average values of 87 and 88%, respectively, for all agricultural matrices used), whereas for other analytes recoveries were between 91 and 101%. The method was also used to determine the real concentrations of Fe(III)-chelates in commercial fertilizers. Furthermore, the method is also capable of resolving two more synthetic Fe(III)-chelates, Fe(III)-EDDHSA and Fe(III)-EDDCHA, whose exact quantification is not currently possible because of lack of commercial standards.

Determination of Slope Safety Factor with Analytical Solution and Searching Critical Slip Surface with Genetic-Traversal Random Method

PubMed Central

2014-01-01

In the current practice, to determine the safety factor of a slope with two-dimensional circular potential failure surface, one of the searching methods for the critical slip surface is Genetic Algorithm (GA), while the method to calculate the slope safety factor is Fellenius' slices method. However GA needs to be validated with more numeric tests, while Fellenius' slices method is just an approximate method like finite element method. This paper proposed a new method to determine the minimum slope safety factor which is the determination of slope safety factor with analytical solution and searching critical slip surface with Genetic-Traversal Random Method. The analytical solution is more accurate than Fellenius' slices method. The Genetic-Traversal Random Method uses random pick to utilize mutation. A computer automatic search program is developed for the Genetic-Traversal Random Method. After comparison with other methods like slope/w software, results indicate that the Genetic-Traversal Random Search Method can give very low safety factor which is about half of the other methods. However the obtained minimum safety factor with Genetic-Traversal Random Search Method is very close to the lower bound solutions of slope safety factor given by the Ansys software. PMID:24782679

Evaluation of purity with its uncertainty value in high purity lead stick by conventional and electro-gravimetric methods

PubMed Central

2013-01-01

Background A conventional gravimetry and electro-gravimetry study has been carried out for the precise and accurate purity determination of lead (Pb) in high purity lead stick and for preparation of reference standard. Reference materials are standards containing a known amount of an analyte and provide a reference value to determine unknown concentrations or to calibrate analytical instruments. A stock solution of approximate 2 kg has been prepared after dissolving approximate 2 g of Pb stick in 5% ultra pure nitric acid. From the stock solution five replicates of approximate 50 g have been taken for determination of purity by each method. The Pb has been determined as PbSO4 by conventional gravimetry, as PbO2 by electro gravimetry. The percentage purity of the metallic Pb was calculated accordingly from PbSO4 and PbO2. Results On the basis of experimental observations it has been concluded that by conventional gravimetry and electro-gravimetry the purity of Pb was found to be 99.98 ± 0.24 and 99.97 ± 0.27 g/100 g and on the basis of Pb purity the concentration of reference standard solutions were found to be 1000.88 ± 2.44 and 1000.81 ± 2.68 mg kg-1 respectively with 95% confidence level (k = 2). The uncertainty evaluation has also been carried out in Pb determination following EURACHEM/GUM guidelines. The final analytical results quantifying uncertainty fulfills this requirement and gives a measure of the confidence level of the concerned laboratory. Conclusions Gravimetry is the most reliable technique in comparison to titremetry and instrumental method and the results of gravimetry are directly traceable to SI unit. Gravimetric analysis, if methods are followed carefully, provides for exceedingly precise analysis. In classical gravimetry the major uncertainties are due to repeatability but in electro-gravimetry several other factors also affect the final results. PMID:23800080

Evaluation of purity with its uncertainty value in high purity lead stick by conventional and electro-gravimetric methods.

PubMed

Singh, Nahar; Singh, Niranjan; Tripathy, S Swarupa; Soni, Daya; Singh, Khem; Gupta, Prabhat K

2013-06-26

A conventional gravimetry and electro-gravimetry study has been carried out for the precise and accurate purity determination of lead (Pb) in high purity lead stick and for preparation of reference standard. Reference materials are standards containing a known amount of an analyte and provide a reference value to determine unknown concentrations or to calibrate analytical instruments. A stock solution of approximate 2 kg has been prepared after dissolving approximate 2 g of Pb stick in 5% ultra pure nitric acid. From the stock solution five replicates of approximate 50 g have been taken for determination of purity by each method. The Pb has been determined as PbSO4 by conventional gravimetry, as PbO2 by electro gravimetry. The percentage purity of the metallic Pb was calculated accordingly from PbSO4 and PbO2. On the basis of experimental observations it has been concluded that by conventional gravimetry and electro-gravimetry the purity of Pb was found to be 99.98 ± 0.24 and 99.97 ± 0.27 g/100 g and on the basis of Pb purity the concentration of reference standard solutions were found to be 1000.88 ± 2.44 and 1000.81 ± 2.68 mg kg-1 respectively with 95% confidence level (k = 2). The uncertainty evaluation has also been carried out in Pb determination following EURACHEM/GUM guidelines. The final analytical results quantifying uncertainty fulfills this requirement and gives a measure of the confidence level of the concerned laboratory. Gravimetry is the most reliable technique in comparison to titremetry and instrumental method and the results of gravimetry are directly traceable to SI unit. Gravimetric analysis, if methods are followed carefully, provides for exceedingly precise analysis. In classical gravimetry the major uncertainties are due to repeatability but in electro-gravimetry several other factors also affect the final results.

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

Approximate bound-state solutions of the Dirac equation for the generalized yukawa potential plus the generalized tensor interaction

NASA Astrophysics Data System (ADS)

Ikot, Akpan N.; Maghsoodi, Elham; Hassanabadi, Hassan; Obu, Joseph A.

2014-05-01

In this paper, we obtain the approximate analytical bound-state solutions of the Dirac particle with the generalized Yukawa potential within the framework of spin and pseudospin symmetries for the arbitrary к state with a generalized tensor interaction. The generalized parametric Nikiforov-Uvarov method is used to obtain the energy eigenvalues and the corresponding wave functions in closed form. We also report some numerical results and present figures to show the effect of the tensor interaction.

Nanophotonic particle simulation and inverse design using artificial neural networks.

PubMed

Peurifoy, John; Shen, Yichen; Jing, Li; Yang, Yi; Cano-Renteria, Fidel; DeLacy, Brendan G; Joannopoulos, John D; Tegmark, Max; Soljačić, Marin

2018-06-01

We propose a method to use artificial neural networks to approximate light scattering by multilayer nanoparticles. We find that the network needs to be trained on only a small sampling of the data to approximate the simulation to high precision. Once the neural network is trained, it can simulate such optical processes orders of magnitude faster than conventional simulations. Furthermore, the trained neural network can be used to solve nanophotonic inverse design problems by using back propagation, where the gradient is analytical, not numerical.

Approximate Green's function methods for HZE transport in multilayered materials

NASA Technical Reports Server (NTRS)

Wilson, John W.; Badavi, Francis F.; Shinn, Judy L.; Costen, Robert C.

1993-01-01

A nonperturbative analytic solution of the high charge and energy (HZE) Green's function is used to implement a computer code for laboratory ion beam transport in multilayered materials. The code is established to operate on the Langley nuclear fragmentation model used in engineering applications. Computational procedures are established to generate linear energy transfer (LET) distributions for a specified ion beam and target for comparison with experimental measurements. The code was found to be highly efficient and compared well with the perturbation approximation.

Triangular dislocation: an analytical, artefact-free solution

NASA Astrophysics Data System (ADS)

Nikkhoo, Mehdi; Walter, Thomas R.

2015-05-01

Displacements and stress-field changes associated with earthquakes, volcanoes, landslides and human activity are often simulated using numerical models in an attempt to understand the underlying processes and their governing physics. The application of elastic dislocation theory to these problems, however, may be biased because of numerical instabilities in the calculations. Here, we present a new method that is free of artefact singularities and numerical instabilities in analytical solutions for triangular dislocations (TDs) in both full-space and half-space. We apply the method to both the displacement and the stress fields. The entire 3-D Euclidean space {R}3 is divided into two complementary subspaces, in the sense that in each one, a particular analytical formulation fulfils the requirements for the ideal, artefact-free solution for a TD. The primary advantage of the presented method is that the development of our solutions involves neither numerical approximations nor series expansion methods. As a result, the final outputs are independent of the scale of the input parameters, including the size and position of the dislocation as well as its corresponding slip vector components. Our solutions are therefore well suited for application at various scales in geoscience, physics and engineering. We validate the solutions through comparison to other well-known analytical methods and provide the MATLAB codes.

Analytic approximation for random muffin-tin alloys

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

Mills, R.; Gray, L.J.; Kaplan, T.

1983-03-15

The methods introduced in a previous paper under the name of ''traveling-cluster approximation'' (TCA) are applied, in a multiple-scattering approach, to the case of a random muffin-tin substitutional alloy. This permits the iterative part of a self-consistent calculation to be carried out entirely in terms of on-the-energy-shell scattering amplitudes. Off-shell components of the mean resolvent, needed for the calculation of spectral functions, are obtained by standard methods involving single-site scattering wave functions. The single-site TCA is just the usual coherent-potential approximation, expressed in a form particularly suited for iteration. A fixed-point theorem is proved for the general t-matrix TCA, ensuringmore » convergence upon iteration to a unique self-consistent solution with the physically essential Herglotz properties.« less

Solution of magnetic field and eddy current problem induced by rotating magnetic poles (abstract)

NASA Astrophysics Data System (ADS)

Liu, Z. J.; Low, T. S.

1996-04-01

The magnetic field and eddy current problems induced by rotating permanent magnet poles occur in electromagnetic dampers, magnetic couplings, and many other devices. Whereas numerical techniques, for example, finite element methods can be exploited to study various features of these problems, such as heat generation and drag torque development, etc., the analytical solution is always of interest to the designers since it helps them to gain the insight into the interdependence of the parameters involved and provides an efficient tool for designing. Some of the previous work showed that the solution of the eddy current problem due to the linearly moving magnet poles can give satisfactory approximation for the eddy current problem due to rotating fields. However, in many practical cases, especially when the number of magnet poles is small, there is significant effect of flux focusing due to the geometry. The above approximation can therefore lead to marked errors in the theoretical predictions of the device performance. Bernot et al. recently described an analytical solution in a polar coordinate system where the radial field is excited by a time-varying source. A discussion of an analytical solution of the magnetic field and eddy current problems induced by moving magnet poles in radial field machines will be given in this article. The theoretical predictions obtained from this method is compared with the results obtained from finite element calculations. The validity of the method is also checked by the comparison of the theoretical predictions and the measurements from a test machine. It is shown that the introduced solution leads to a significant improvement in the air gap field prediction as compared with the results obtained from the analytical solution that models the eddy current problems induced by linearly moving magnet poles.

An approximate analytical solution for interlaminar stresses in angle-ply laminates

NASA Technical Reports Server (NTRS)

Rose, Cheryl A.; Herakovich, Carl T.

1991-01-01

An improved approximate analytical solution for interlaminar stresses in finite width, symmetric, angle-ply laminated coupons subjected to axial loading is presented. The solution is based upon statically admissible stress fields which take into consideration local property mismatch effects and global equilibrium requirements. Unknown constants in the admissible stress states are determined through minimization of the complementary energy. Typical results are presented for through-the-thickness and interlaminar stress distributions for angle-ply laminates. It is shown that the results represent an improved approximate analytical solution for interlaminar stresses.

Finite element modeling of light propagation in fruit under illumination of continuous-wave beam

USDA-ARS?s Scientific Manuscript database

Spatially-resolved spectroscopy provides a means for measuring the optical properties of biological tissues, based on analytical solutions to diffusion approximation for semi-infinite media under the normal illumination of infinitely small size light beam. The method is, however, prone to error in m...

Finding Optimal Gains In Linear-Quadratic Control Problems

NASA Technical Reports Server (NTRS)

Milman, Mark H.; Scheid, Robert E., Jr.

1990-01-01

Analytical method based on Volterra factorization leads to new approximations for optimal control gains in finite-time linear-quadratic control problem of system having infinite number of dimensions. Circumvents need to analyze and solve Riccati equations and provides more transparent connection between dynamics of system and optimal gain.

Karhunen Loève approximation of random fields by generalized fast multipole methods

NASA Astrophysics Data System (ADS)

Schwab, Christoph; Todor, Radu Alexandru

2006-09-01

KL approximation of a possibly instationary random field a( ω, x) ∈ L2( Ω, d P; L∞( D)) subject to prescribed meanfield Ea(x)=∫a(ω,x) dP(ω) and covariance Va(x,x')=∫(a(ω,x)-Ea(x))(a(ω,x')-Ea(x')) dP(ω) in a polyhedral domain D⊂Rd is analyzed. We show how for stationary covariances Va( x, x') = ga(| x - x'|) with ga( z) analytic outside of z = 0, an M-term approximate KL-expansion aM( ω, x) of a( ω, x) can be computed in log-linear complexity. The approach applies in arbitrary domains D and for nonseparable covariances Ca. It involves Galerkin approximation of the KL eigenvalue problem by discontinuous finite elements of degree p ⩾ 0 on a quasiuniform, possibly unstructured mesh of width h in D, plus a generalized fast multipole accelerated Krylov-Eigensolver. The approximate KL-expansion aM( x, ω) of a( x, ω) has accuracy O(exp(- bM1/ d)) if ga is analytic at z = 0 and accuracy O( M- k/ d) if ga is Ck at zero. It is obtained in O( MN(log N) b) operations where N = O( h- d).

Application of the optimal homotopy asymptotic method to nonlinear Bingham fluid dampers

NASA Astrophysics Data System (ADS)

Marinca, Vasile; Ene, Remus-Daniel; Bereteu, Liviu

2017-10-01

Dynamic response time is an important feature for determining the performance of magnetorheological (MR) dampers in practical civil engineering applications. The objective of this paper is to show how to use the Optimal Homotopy Asymptotic Method (OHAM) to give approximate analytical solutions of the nonlinear differential equation of a modified Bingham model with non-viscous exponential damping. Our procedure does not depend upon small parameters and provides us with a convenient way to optimally control the convergence of the approximate solutions. OHAM is very efficient in practice for ensuring very rapid convergence of the solution after only one iteration and with a small number of steps.

Constructing analytic solutions on the Tricomi equation

NASA Astrophysics Data System (ADS)

Ghiasi, Emran Khoshrouye; Saleh, Reza

2018-04-01

In this paper, homotopy analysis method (HAM) and variational iteration method (VIM) are utilized to derive the approximate solutions of the Tricomi equation. Afterwards, the HAM is optimized to accelerate the convergence of the series solution by minimizing its square residual error at any order of the approximation. It is found that effect of the optimal values of auxiliary parameter on the convergence of the series solution is not negligible. Furthermore, the present results are found to agree well with those obtained through a closed-form equation available in the literature. To conclude, it is seen that the two are effective to achieve the solution of the partial differential equations.

Lead Slowing-Down Spectrometry Time Spectral Analysis for Spent Fuel Assay: FY11 Status Report

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

Kulisek, Jonathan A.; Anderson, Kevin K.; Bowyer, Sonya M.

2011-09-30

Developing a method for the accurate, direct, and independent assay of the fissile isotopes in bulk materials (such as used fuel) from next-generation domestic nuclear fuel cycles is a goal of the Office of Nuclear Energy, Fuel Cycle R&D, Material Protection and Control Technology (MPACT) Campaign. To meet this goal, MPACT supports a multi-institutional collaboration, of which PNNL is a part, to study the feasibility of Lead Slowing Down Spectroscopy (LSDS). This technique is an active nondestructive assay method that has the potential to provide independent, direct measurement of Pu and U isotopic masses in used fuel with an uncertaintymore » considerably lower than the approximately 10% typical of today's confirmatory assay methods. This document is a progress report for FY2011 PNNL analysis and algorithm development. Progress made by PNNL in FY2011 continues to indicate the promise of LSDS analysis and algorithms applied to used fuel. PNNL developed an empirical model based on calibration of the LSDS to responses generated from well-characterized used fuel. The empirical model, which accounts for self-shielding effects using empirical basis vectors calculated from the singular value decomposition (SVD) of a matrix containing the true self-shielding functions of the used fuel assembly models. The potential for the direct and independent assay of the sum of the masses of 239Pu and 241Pu to within approximately 3% over a wide used fuel parameter space was demonstrated. Also, in FY2011, PNNL continued to develop an analytical model. Such efforts included the addition of six more non-fissile absorbers in the analytical shielding function and the non-uniformity of the neutron flux across the LSDS assay chamber. A hybrid analytical-empirical approach was developed to determine the mass of total Pu (sum of the masses of 239Pu, 240Pu, and 241Pu), which is an important quantity in safeguards. Results using this hybrid method were of approximately the same accuracy as the pure empirical approach. In addition, total Pu with much better accuracy with the hybrid approach than the pure analytical approach. In FY2012, PNNL will continue efforts to optimize its empirical model and minimize its reliance on calibration data. In addition, PNNL will continue to develop an analytical model, considering effects such as neutron-scattering in the fuel and cladding, as well as neutrons streaming through gaps between fuel pins in the fuel assembly.« less

High-Contrast Gratings based Spoof Surface Plasmons

NASA Astrophysics Data System (ADS)

Li, Zhuo; Liu, Liangliang; Xu, Bingzheng; Ning, Pingping; Chen, Chen; Xu, Jia; Chen, Xinlei; Gu, Changqing; Qing, Quan

2016-02-01

In this work, we explore the existence of spoof surface plasmons (SSPs) supported by deep-subwavelength high-contrast gratings (HCGs) on a perfect electric conductor plane. The dispersion relation of the HCGs-based SSPs is derived analyt- ically by combining multimode network theory with rigorous mode matching method, which has nearly the same form with and can be degenerated into that of the SSPs arising from deep-subwavelength metallic gratings (MGs). Numerical simula- tions validate the analytical dispersion relation and an effective medium approximation is also presented to obtain the same analytical dispersion formula. This work sets up a unified theoretical framework for SSPs and opens up new vistas in surface plasmon optics.

Review and assessment of the database and numerical modeling for turbine heat transfer

NASA Technical Reports Server (NTRS)

Gladden, H. J.; Simoneau, R. J.

1989-01-01

The objectives of the NASA Hot Section Technology (HOST) Turbine Heat Transfer subproject were to obtain a better understanding of the physics of the aerothermodynamic phenomena and to assess and improve the analytical methods used to predict the flow and heat transfer in high-temperature gas turbines. At the time the HOST project was initiated, an across-the-board improvement in turbine design technology was needed. A building-block approach was utilized and the research ranged from the study of fundamental phenomena and modeling to experiments in simulated real engine environments. Experimental research accounted for approximately 75 percent of the funding while the analytical efforts were approximately 25 percent. A healthy government/industry/university partnership, with industry providing almost half of the research, was created to advance the turbine heat transfer design technology base.

A new frequency domain analytical solution of a cascade of diffusive channels for flood routing

NASA Astrophysics Data System (ADS)

Cimorelli, Luigi; Cozzolino, Luca; Della Morte, Renata; Pianese, Domenico; Singh, Vijay P.

2015-04-01

Simplified flood propagation models are often employed in practical applications for hydraulic and hydrologic analyses. In this paper, we present a new numerical method for the solution of the Linear Parabolic Approximation (LPA) of the De Saint Venant equations (DSVEs), accounting for the space variation of model parameters and the imposition of appropriate downstream boundary conditions. The new model is based on the analytical solution of a cascade of linear diffusive channels in the Laplace Transform domain. The time domain solutions are obtained using a Fourier series approximation of the Laplace Inversion formula. The new Inverse Laplace Transform Diffusive Flood Routing model (ILTDFR) can be used as a building block for the construction of real-time flood forecasting models or in optimization models, because it is unconditionally stable and allows fast and fairly precise computation.

An experimental and analytical method for approximate determination of the tilt rotor research aircraft rotor/wing download

NASA Technical Reports Server (NTRS)

Jordon, D. E.; Patterson, W.; Sandlin, D. R.

1985-01-01

The XV-15 Tilt Rotor Research Aircraft download phenomenon was analyzed. This phenomenon is a direct result of the two rotor wakes impinging on the wing upper surface when the aircraft is in the hover configuration. For this study the analysis proceeded along tow lines. First was a method whereby results from actual hover tests of the XV-15 aircraft were combined with drag coefficient results from wind tunnel tests of a wing that was representative of the aircraft wing. Second, an analytical method was used that modeled that airflow caused gy the two rotors. Formulas were developed in such a way that acomputer program could be used to calculate the axial velocities were then used in conjunction with the aforementioned wind tunnel drag coefficinet results to produce download values. An attempt was made to validate the analytical results by modeling a model rotor system for which direct download values were determinrd..

Generation of dark hollow beams by using a fractional radial Hilbert transform system

NASA Astrophysics Data System (ADS)

Xie, Qiansen; Zhao, Daomu

2007-07-01

The radial Hilbert transform has been extend to the fractional field, which could be called the fractional radial Hilbert transform (FRHT). Using edge-enhancement characteristics of this transform, we convert a Gaussian light beam into a variety of dark hollow beams (DHBs). Based on the fact that a hard-edged aperture can be expanded approximately as a finite sum of complex Gaussian functions, the analytical expression of a Gaussian beam passing through a FRHT system has been derived. As a numerical example, the properties of the DHBs with different fractional orders are illustrated graphically. The calculation results obtained by use of the analytical method and the integral method are also compared.

Analytical prediction of the heat transfer from a blood vessel near the skin surface when cooled by a symmetrical cooling strip

NASA Technical Reports Server (NTRS)

Chato, J. C.; Shitzer, A.

1971-01-01

An analytical method was developed to estimate the amount of heat extracted from an artery running close to the skin surface which is cooled in a symmetrical fashion by a cooling strip. The results indicate that the optimum width of a cooling strip is approximately three times the depth to the centerline of the artery. The heat extracted from an artery with such a strip is about 0.9 w/m-C which is too small to affect significantly the temperature of the blood flow through a main blood vessel, such as the carotid artery. The method is applicable to veins as well.

Universal analytical scattering form factor for shell-, core-shell, or homogeneous particles with continuously variable density profile shape.

PubMed

Foster, Tobias

2011-09-01

A novel analytical and continuous density distribution function with a widely variable shape is reported and used to derive an analytical scattering form factor that allows us to universally describe the scattering from particles with the radial density profile of homogeneous spheres, shells, or core-shell particles. Composed by the sum of two Fermi-Dirac distribution functions, the shape of the density profile can be altered continuously from step-like via Gaussian-like or parabolic to asymptotically hyperbolic by varying a single "shape parameter", d. Using this density profile, the scattering form factor can be calculated numerically. An analytical form factor can be derived using an approximate expression for the original Fermi-Dirac distribution function. This approximation is accurate for sufficiently small rescaled shape parameters, d/R (R being the particle radius), up to values of d/R ≈ 0.1, and thus captures step-like, Gaussian-like, and parabolic as well as asymptotically hyperbolic profile shapes. It is expected that this form factor is particularly useful in a model-dependent analysis of small-angle scattering data since the applied continuous and analytical function for the particle density profile can be compared directly with the density profile extracted from the data by model-free approaches like the generalized inverse Fourier transform method. © 2011 American Chemical Society

Analytic derivation of an approximate SU(3) symmetry inside the symmetry triangle of the interacting boson approximation model

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

Bonatsos, Dennis; Karampagia, S.; Casten, R. F.

2011-05-15

Using a contraction of the SU(3) algebra to the algebra of the rigid rotator in the large-boson-number limit of the interacting boson approximation (IBA) model, a line is found inside the symmetry triangle of the IBA, along which the SU(3) symmetry is preserved. The line extends from the SU(3) vertex to near the critical line of the first-order shape/phase transition separating the spherical and prolate deformed phases, and it lies within the Alhassid-Whelan arc of regularity, the unique valley of regularity connecting the SU(3) and U(5) vertices in the midst of chaotic regions. In addition to providing an explanation formore » the existence of the arc of regularity, the present line represents an example of an analytically determined approximate symmetry in the interior of the symmetry triangle of the IBA. The method is applicable to algebraic models possessing subalgebras amenable to contraction. This condition is equivalent to algebras in which the equilibrium ground state and its rotational band become energetically isolated from intrinsic excitations, as typified by deformed solutions to the IBA for large numbers of valence nucleons.« less

Inversion of the anomalous diffraction approximation for variable complex index of refraction near unity. [numerical tests for water-haze aerosol model

NASA Technical Reports Server (NTRS)

Smith, C. B.

1982-01-01

The Fymat analytic inversion method for retrieving a particle-area distribution function from anomalous diffraction multispectral extinction data and total area is generalized to the case of a variable complex refractive index m(lambda) near unity depending on spectral wavelength lambda. Inversion tests are presented for a water-haze aerosol model. An upper-phase shift limit of 5 pi/2 retrieved an accurate peak area distribution profile. Analytical corrections using both the total number and area improved the inversion.

Benchmark solutions for the galactic heavy-ion transport equations with energy and spatial coupling

NASA Technical Reports Server (NTRS)

Ganapol, Barry D.; Townsend, Lawrence W.; Lamkin, Stanley L.; Wilson, John W.

1991-01-01

Nontrivial benchmark solutions are developed for the galactic heavy ion transport equations in the straightahead approximation with energy and spatial coupling. Analytical representations of the ion fluxes are obtained for a variety of sources with the assumption that the nuclear interaction parameters are energy independent. The method utilizes an analytical LaPlace transform inversion to yield a closed form representation that is computationally efficient. The flux profiles are then used to predict ion dose profiles, which are important for shield design studies.

Stress Analysis of Beams with Shear Deformation of the Flanges

NASA Technical Reports Server (NTRS)

Kuhn, Paul

1937-01-01

This report discusses the fundamental action of shear deformation of the flanges on the basis of simplifying assumptions. The theory is developed to the point of giving analytical solutions for simple cases of beams and of skin-stringer panels under axial load. Strain-gage tests on a tension panel and on a beam corresponding to these simple cases are described and the results are compared with analytical results. For wing beams, an approximate method of applying the theory is given. As an alternative, the construction of a mechanical analyzer is advocated.

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.

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.

Analytical approximations to seawater optical phase functions of scattering

NASA Astrophysics Data System (ADS)

Haltrin, Vladimir I.

2004-11-01

This paper proposes a number of analytical approximations to the classic and recently measured seawater light scattering phase functions. The three types of analytical phase functions are derived: individual representations for 15 Petzold, 41 Mankovsky, and 91 Gulf of Mexico phase functions; collective fits to Petzold phase functions; and analytical representations that take into account dependencies between inherent optical properties of seawater. The proposed phase functions may be used for problems of radiative transfer, remote sensing, visibility and image propagation in natural waters of various turbidity.

On finding the analytic dependencies of the external field potential on the control function when optimizing the beam dynamics

NASA Astrophysics Data System (ADS)

Ovsyannikov, A. D.; Kozynchenko, S. A.; Kozynchenko, V. A.

2017-12-01

When developing a particle accelerator for generating the high-precision beams, the injection system design is of importance, because it largely determines the output characteristics of the beam. At the present paper we consider the injection systems consisting of electrodes with given potentials. The design of such systems requires carrying out simulation of beam dynamics in the electrostatic fields. For external field simulation we use the new approach, proposed by A.D. Ovsyannikov, which is based on analytical approximations, or finite difference method, taking into account the real geometry of the injection system. The software designed for solving the problems of beam dynamics simulation and optimization in the injection system for non-relativistic beams has been developed. Both beam dynamics and electric field simulations in the injection system which use analytical approach and finite difference method have been made and the results presented in this paper.

Verification of spatial and temporal pressure distributions in segmented solid rocket motors

NASA Technical Reports Server (NTRS)

Salita, Mark

1989-01-01

A wide variety of analytical tools are in use today to predict the history and spatial distributions of pressure in the combustion chambers of solid rocket motors (SRMs). Experimental and analytical methods are presented here that allow the verification of many of these predictions. These methods are applied to the redesigned space shuttle booster (RSRM). Girth strain-gage data is compared to the predictions of various one-dimensional quasisteady analyses in order to verify the axial drop in motor static pressure during ignition transients as well as quasisteady motor operation. The results of previous modeling of radial flows in the bore, slots, and around grain overhangs are supported by approximate analytical and empirical techniques presented here. The predictions of circumferential flows induced by inhibitor asymmetries, nozzle vectoring, and propellant slump are compared to each other and to subscale cold air and water tunnel measurements to ascertain their validity.

Exact exchange-correlation potentials of singlet two-electron systems

NASA Astrophysics Data System (ADS)

Ryabinkin, Ilya G.; Ospadov, Egor; Staroverov, Viktor N.

2017-10-01

We suggest a non-iterative analytic method for constructing the exchange-correlation potential, v XC ( r ) , of any singlet ground-state two-electron system. The method is based on a convenient formula for v XC ( r ) in terms of quantities determined only by the system's electronic wave function, exact or approximate, and is essentially different from the Kohn-Sham inversion technique. When applied to Gaussian-basis-set wave functions, the method yields finite-basis-set approximations to the corresponding basis-set-limit v XC ( r ) , whereas the Kohn-Sham inversion produces physically inappropriate (oscillatory and divergent) potentials. The effectiveness of the procedure is demonstrated by computing accurate exchange-correlation potentials of several two-electron systems (helium isoelectronic series, H2, H3 + ) using common ab initio methods and Gaussian basis sets.

Approximate series solution of multi-dimensional, time fractional-order (heat-like) diffusion equations using FRDTM.

PubMed

Singh, Brajesh K; Srivastava, Vineet K

2015-04-01

The main goal of this paper is to present a new approximate series solution of the multi-dimensional (heat-like) diffusion equation with time-fractional derivative in Caputo form using a semi-analytical approach: fractional-order reduced differential transform method (FRDTM). The efficiency of FRDTM is confirmed by considering four test problems of the multi-dimensional time fractional-order diffusion equation. FRDTM is a very efficient, effective and powerful mathematical tool which provides exact or very close approximate solutions for a wide range of real-world problems arising in engineering and natural sciences, modelled in terms of differential equations.

Approximate series solution of multi-dimensional, time fractional-order (heat-like) diffusion equations using FRDTM

PubMed Central

Singh, Brajesh K.; Srivastava, Vineet K.

2015-01-01

The main goal of this paper is to present a new approximate series solution of the multi-dimensional (heat-like) diffusion equation with time-fractional derivative in Caputo form using a semi-analytical approach: fractional-order reduced differential transform method (FRDTM). The efficiency of FRDTM is confirmed by considering four test problems of the multi-dimensional time fractional-order diffusion equation. FRDTM is a very efficient, effective and powerful mathematical tool which provides exact or very close approximate solutions for a wide range of real-world problems arising in engineering and natural sciences, modelled in terms of differential equations. PMID:26064639

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

Chhiber, R; Usmanov, AV; Matthaeus, WH

Simple estimates of the number of Coulomb collisions experienced by the interplanetary plasma to the point of observation, i.e., the “collisional age”, can be usefully employed in the study of non-thermal features of the solar wind. Usually these estimates are based on local plasma properties at the point of observation. Here we improve the method of estimation of the collisional age by employing solutions obtained from global three-dimensional magnetohydrodynamics simulations. This enables evaluation of the complete analytical expression for the collisional age without using approximations. The improved estimation of the collisional timescale is compared with turbulence and expansion timescales tomore » assess the relative importance of collisions. The collisional age computed using the approximate formula employed in previous work is compared with the improved simulation-based calculations to examine the validity of the simplified formula. We also develop an analytical expression for the evaluation of the collisional age and we find good agreement between the numerical and analytical results. Finally, we briefly discuss the implications for an improved estimation of collisionality along spacecraft trajectories, including Solar Probe Plus.« less

Efficient Band-to-Trap Tunneling Model Including Heterojunction Band Offset

DOE PAGES

Gao, Xujiao; Huang, Andy; Kerr, Bert

2017-10-25

In this paper, we present an efficient band-to-trap tunneling model based on the Schenk approach, in which an analytic density-of-states (DOS) model is developed based on the open boundary scattering method. The new model explicitly includes the effect of heterojunction band offset, in addition to the well-known field effect. Its analytic form enables straightforward implementation into TCAD device simulators. It is applicable to all one-dimensional potentials, which can be approximated to a good degree such that the approximated potentials lead to piecewise analytic wave functions with open boundary conditions. The model allows for simulating both the electric-field-enhanced and band-offset-enhanced carriermore » recombination due to the band-to-trap tunneling near the heterojunction in a heterojunction bipolar transistor (HBT). Simulation results of an InGaP/GaAs/GaAs NPN HBT show that the proposed model predicts significantly increased base currents, due to the hole-to-trap tunneling enhanced by the emitter-base junction band offset. Finally, the results compare favorably with experimental observation.« less

Accurate expressions for solar cell fill factors including series and shunt resistances

NASA Astrophysics Data System (ADS)

Green, Martin A.

2016-02-01

Together with open-circuit voltage and short-circuit current, fill factor is a key solar cell parameter. In their classic paper on limiting efficiency, Shockley and Queisser first investigated this factor's analytical properties showing, for ideal cells, it could be expressed implicitly in terms of the maximum power point voltage. Subsequently, fill factors usually have been calculated iteratively from such implicit expressions or from analytical approximations. In the absence of detrimental series and shunt resistances, analytical fill factor expressions have recently been published in terms of the Lambert W function available in most mathematical computing software. Using a recently identified perturbative relationship, exact expressions in terms of this function are derived in technically interesting cases when both series and shunt resistances are present but have limited impact, allowing a better understanding of their effect individually and in combination. Approximate expressions for arbitrary shunt and series resistances are then deduced, which are significantly more accurate than any previously published. A method based on the insights developed is also reported for deducing one-diode fits to experimental data.

Efficient Band-to-Trap Tunneling Model Including Heterojunction Band Offset

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

Gao, Xujiao; Huang, Andy; Kerr, Bert

In this paper, we present an efficient band-to-trap tunneling model based on the Schenk approach, in which an analytic density-of-states (DOS) model is developed based on the open boundary scattering method. The new model explicitly includes the effect of heterojunction band offset, in addition to the well-known field effect. Its analytic form enables straightforward implementation into TCAD device simulators. It is applicable to all one-dimensional potentials, which can be approximated to a good degree such that the approximated potentials lead to piecewise analytic wave functions with open boundary conditions. The model allows for simulating both the electric-field-enhanced and band-offset-enhanced carriermore » recombination due to the band-to-trap tunneling near the heterojunction in a heterojunction bipolar transistor (HBT). Simulation results of an InGaP/GaAs/GaAs NPN HBT show that the proposed model predicts significantly increased base currents, due to the hole-to-trap tunneling enhanced by the emitter-base junction band offset. Finally, the results compare favorably with experimental observation.« less

User's manual for the one-dimensional hypersonic experimental aero-thermodynamic (1DHEAT) data reduction code

NASA Technical Reports Server (NTRS)

Hollis, Brian R.

1995-01-01

A FORTRAN computer code for the reduction and analysis of experimental heat transfer data has been developed. This code can be utilized to determine heat transfer rates from surface temperature measurements made using either thin-film resistance gages or coaxial surface thermocouples. Both an analytical and a numerical finite-volume heat transfer model are implemented in this code. The analytical solution is based on a one-dimensional, semi-infinite wall thickness model with the approximation of constant substrate thermal properties, which is empirically corrected for the effects of variable thermal properties. The finite-volume solution is based on a one-dimensional, implicit discretization. The finite-volume model directly incorporates the effects of variable substrate thermal properties and does not require the semi-finite wall thickness approximation used in the analytical model. This model also includes the option of a multiple-layer substrate. Fast, accurate results can be obtained using either method. This code has been used to reduce several sets of aerodynamic heating data, of which samples are included in this report.

A finite element analysis of viscoelastically damped sandwich plates

NASA Astrophysics Data System (ADS)

Ma, B.-A.; He, J.-F.

1992-01-01

A finite element analysis associated with an asymptotic solution method for the harmonic flexural vibration of viscoelastically damped unsymmetrical sandwich plates is given. The element formulation is based on generalization of the discrete Kirchhoff theory (DKT) element formulation. The results obtained with the first order approximation of the asymptotic solution presented here are the same as those obtained by means of the modal strain energy (MSE) method. By taking more terms of the asymptotic solution, with successive calculations and use of the Padé approximants method, accuracy can be improved. The finite element computation has been verified by comparison with an analytical exact solution for rectangular plates with simply supported edges. Results for the same plates with clamped edges are also presented.

Analytical Plug-In Method for Kernel Density Estimator Applied to Genetic Neutrality Study

NASA Astrophysics Data System (ADS)

Troudi, Molka; Alimi, Adel M.; Saoudi, Samir

2008-12-01

The plug-in method enables optimization of the bandwidth of the kernel density estimator in order to estimate probability density functions (pdfs). Here, a faster procedure than that of the common plug-in method is proposed. The mean integrated square error (MISE) depends directly upon [InlineEquation not available: see fulltext.] which is linked to the second-order derivative of the pdf. As we intend to introduce an analytical approximation of [InlineEquation not available: see fulltext.], the pdf is estimated only once, at the end of iterations. These two kinds of algorithm are tested on different random variables having distributions known for their difficult estimation. Finally, they are applied to genetic data in order to provide a better characterisation in the mean of neutrality of Tunisian Berber populations.

Prediction of thermal cycling induced matrix cracking

NASA Technical Reports Server (NTRS)

Mcmanus, Hugh L.

1992-01-01

Thermal fatigue has been observed to cause matrix cracking in laminated composite materials. A method is presented to predict transverse matrix cracks in composite laminates subjected to cyclic thermal load. Shear lag stress approximations and a simple energy-based fracture criteria are used to predict crack densities as a function of temperature. Prediction of crack densities as a function of thermal cycling is accomplished by assuming that fatigue degrades the material's inherent resistance to cracking. The method is implemented as a computer program. A simple experiment provides data on progressive cracking of a laminate with decreasing temperature. Existing data on thermal fatigue is also used. Correlations of the analytical predictions to the data are very good. A parametric study using the analytical method is presented which provides insight into material behavior under cyclical thermal loads.

Analytical methods for determination of mycotoxins: a review.

PubMed

Turner, Nicholas W; Subrahmanyam, Sreenath; Piletsky, Sergey A

2009-01-26

Mycotoxins are small (MW approximately 700), toxic chemical products formed as secondary metabolites by a few fungal species that readily colonise crops and contaminate them with toxins in the field or after harvest. Ochratoxins and Aflatoxins are mycotoxins of major significance and hence there has been significant research on broad range of analytical and detection techniques that could be useful and practical. Due to the variety of structures of these toxins, it is impossible to use one standard technique for analysis and/or detection. Practical requirements for high-sensitivity analysis and the need for a specialist laboratory setting create challenges for routine analysis. Several existing analytical techniques, which offer flexible and broad-based methods of analysis and in some cases detection, have been discussed in this manuscript. There are a number of methods used, of which many are lab-based, but to our knowledge there seems to be no single technique that stands out above the rest, although analytical liquid chromatography, commonly linked with mass spectroscopy is likely to be popular. This review manuscript discusses (a) sample pre-treatment methods such as liquid-liquid extraction (LLE), supercritical fluid extraction (SFE), solid phase extraction (SPE), (b) separation methods such as (TLC), high performance liquid chromatography (HPLC), gas chromatography (GC), and capillary electrophoresis (CE) and (c) others such as ELISA. Further currents trends, advantages and disadvantages and future prospects of these methods have been discussed.

Comparison of approximate solutions to the phonon Boltzmann transport equation with the relaxation time approximation: Spherical harmonics expansions and the discrete ordinates method

NASA Astrophysics Data System (ADS)

Christenson, J. G.; Austin, R. A.; Phillips, R. J.

2018-05-01

The phonon Boltzmann transport equation is used to analyze model problems in one and two spatial dimensions, under transient and steady-state conditions. New, explicit solutions are obtained by using the P1 and P3 approximations, based on expansions in spherical harmonics, and are compared with solutions from the discrete ordinates method. For steady-state energy transfer, it is shown that analytic expressions derived using the P1 and P3 approximations agree quantitatively with the discrete ordinates method, in some cases for large Knudsen numbers, and always for Knudsen numbers less than unity. However, for time-dependent energy transfer, the PN solutions differ qualitatively from converged solutions obtained by the discrete ordinates method. Although they correctly capture the wave-like behavior of energy transfer at short times, the P1 and P3 approximations rely on one or two wave velocities, respectively, yielding abrupt, step-changes in temperature profiles that are absent when the angular dependence of the phonon velocities is captured more completely. It is shown that, with the gray approximation, the P1 approximation is formally equivalent to the so-called "hyperbolic heat equation." Overall, these results support the use of the PN approximation to find solutions to the phonon Boltzmann transport equation for steady-state conditions. Such solutions can be useful in the design and analysis of devices that involve heat transfer at nanometer length scales, where continuum-scale approaches become inaccurate.

Unbiased Scanning Method and Data Banking Approach Using Ultra-High Performance Liquid Chromatography Coupled with High-Resolution Mass Spectrometry for Quantitative Comparison of Metabolite Exposure in Plasma across Species Analyzed at Different Dates.

PubMed

Gao, Hongying; Deng, Shibing; Obach, R Scott

2015-12-01

An unbiased scanning methodology using ultra high-performance liquid chromatography coupled with high-resolution mass spectrometry was used to bank data and plasma samples for comparing the data generated at different dates. This method was applied to bank the data generated earlier in animal samples and then to compare the exposure to metabolites in animal versus human for safety assessment. With neither authentic standards nor prior knowledge of the identities and structures of metabolites, full scans for precursor ions and all ion fragments (AIF) were employed with a generic gradient LC method to analyze plasma samples at positive and negative polarity, respectively. In a total of 22 tested drugs and metabolites, 21 analytes were detected using this unbiased scanning method except that naproxen was not detected due to low sensitivity at negative polarity and interference at positive polarity; and 4'- or 5-hydroxy diclofenac was not separated by a generic UPLC method. Statistical analysis of the peak area ratios of the analytes versus the internal standard in five repetitive analyses over approximately 1 year demonstrated that the analysis variation was significantly different from sample instability. The confidence limits for comparing the exposure using peak area ratio of metabolites in animal plasma versus human plasma measured over approximately 1 year apart were comparable to the analysis undertaken side by side on the same days. These statistical analysis results showed it was feasible to compare data generated at different dates with neither authentic standards nor prior knowledge of the analytes.

Calibration Method to Eliminate Zeroth Order Effect in Lateral Shearing Interferometry

NASA Astrophysics Data System (ADS)

Fang, Chao; Xiang, Yang; Qi, Keqi; Chen, Dawei

2018-04-01

In this paper, a calibration method is proposed which eliminates the zeroth order effect in lateral shearing interferometry. An analytical expression of the calibration error function is deduced, and the relationship between the phase-restoration error and calibration error is established. The analytical results show that the phase-restoration error introduced by the calibration error is proportional to the phase shifting error and zeroth order effect. The calibration method is verified using simulations and experiments. The simulation results show that the phase-restoration error is approximately proportional to the phase shift error and zeroth order effect, when the phase shifting error is less than 2° and the zeroth order effect is less than 0.2. The experimental result shows that compared with the conventional method with 9-frame interferograms, the calibration method with 5-frame interferograms achieves nearly the same restoration accuracy.

Statistical analysis of loopy belief propagation in random fields

NASA Astrophysics Data System (ADS)

Yasuda, Muneki; Kataoka, Shun; Tanaka, Kazuyuki

2015-10-01

Loopy belief propagation (LBP), which is equivalent to the Bethe approximation in statistical mechanics, is a message-passing-type inference method that is widely used to analyze systems based on Markov random fields (MRFs). In this paper, we propose a message-passing-type method to analytically evaluate the quenched average of LBP in random fields by using the replica cluster variation method. The proposed analytical method is applicable to general pairwise MRFs with random fields whose distributions differ from each other and can give the quenched averages of the Bethe free energies over random fields, which are consistent with numerical results. The order of its computational cost is equivalent to that of standard LBP. In the latter part of this paper, we describe the application of the proposed method to Bayesian image restoration, in which we observed that our theoretical results are in good agreement with the numerical results for natural images.

A Higher Order Iterative Method for Computing the Drazin Inverse

PubMed Central

Soleymani, F.; Stanimirović, Predrag S.

2013-01-01

A method with high convergence rate for finding approximate inverses of nonsingular matrices is suggested and established analytically. An extension of the introduced computational scheme to general square matrices is defined. The extended method could be used for finding the Drazin inverse. The application of the scheme on large sparse test matrices alongside the use in preconditioning of linear system of equations will be presented to clarify the contribution of the paper. PMID:24222747

Thermal ionization of Cs Rydberg states

NASA Astrophysics Data System (ADS)

Glukhov, I. L.; Ovsiannikov, V. D.

2009-01-01

Rates Pnl of photoionization from Rydberg ns-, np-, nd-states of a valence electron in Cs, induced by black-body radiation, were calculated on the basis of the modified Fues model potential method. The numerical data were approximated with a three-term expression which reproduces in a simple analytical form the dependence of Pnl on the ambient temperature T and on the principal quantum number n. The comparison between approximate and exactly calculated values of the thermal ionization rate demonstrates the applicability of the proposed approximation for highly excited states with n from 20 to 100 in a wide temperature range of T from 100 to 10,000 K. We present coefficients of this approximation for the s-, p- and d-series of Rydberg states.

General properties and analytical approximations of photorefractive solitons

NASA Astrophysics Data System (ADS)

Geisler, A.; Homann, F.; Schmidt, H.-J.

2004-08-01

We investigate general properties of spatial 1-dimensional bright photorefractive solitons and discuss various analytical approximations for the soliton profile and the half width, both depending on an intensity parameter r. The case of dark solitons is also briefly addressed.

Unsteady Flow Simulation: A Numerical Challenge

DTIC Science & Technology

2003-03-01

drive to convergence the numerical unsteady term. The time marching procedure is based on the approximate implicit Newton method for systems of non...computed through analytical derivatives of S. The linear system stemming from equation (3) is solved at each integration step by the same iterative method...significant reduction of memory usage, thanks to the reduced dimensions of the linear system matrix during the implicit marching of the solution. The

Fuzzy and neural control

NASA Technical Reports Server (NTRS)

Berenji, Hamid R.

1992-01-01

Fuzzy logic and neural networks provide new methods for designing control systems. Fuzzy logic controllers do not require a complete analytical model of a dynamic system and can provide knowledge-based heuristic controllers for ill-defined and complex systems. Neural networks can be used for learning control. In this chapter, we discuss hybrid methods using fuzzy logic and neural networks which can start with an approximate control knowledge base and refine it through reinforcement learning.

Comptonization in Ultra-Strong Magnetic Fields: Numerical Solution to the Radiative Transfer Problem

NASA Technical Reports Server (NTRS)

Ceccobello, C.; Farinelli, R.; Titarchuk, L.

2014-01-01

We consider the radiative transfer problem in a plane-parallel slab of thermal electrons in the presence of an ultra-strong magnetic field (B approximately greater than B(sub c) approx. = 4.4 x 10(exp 13) G). Under these conditions, the magnetic field behaves like a birefringent medium for the propagating photons, and the electromagnetic radiation is split into two polarization modes, ordinary and extraordinary, that have different cross-sections. When the optical depth of the slab is large, the ordinary-mode photons are strongly Comptonized and the photon field is dominated by an isotropic component. Aims. The radiative transfer problem in strong magnetic fields presents many mathematical issues and analytical or numerical solutions can be obtained only under some given approximations. We investigate this problem both from the analytical and numerical point of view, provide a test of the previous analytical estimates, and extend these results with numerical techniques. Methods. We consider here the case of low temperature black-body photons propagating in a sub-relativistic temperature plasma, which allows us to deal with a semi-Fokker-Planck approximation of the radiative transfer equation. The problem can then be treated with the variable separation method, and we use a numerical technique to find solutions to the eigenvalue problem in the case of a singular kernel of the space operator. The singularity of the space kernel is the result of the strong angular dependence of the electron cross-section in the presence of a strong magnetic field. Results. We provide the numerical solution obtained for eigenvalues and eigenfunctions of the space operator, and the emerging Comptonization spectrum of the ordinary-mode photons for any eigenvalue of the space equation and for energies significantly lesser than the cyclotron energy, which is on the order of MeV for the intensity of the magnetic field here considered. Conclusions. We derived the specific intensity of the ordinary photons, under the approximation of large angle and large optical depth. These assumptions allow the equation to be treated using a diffusion-like approximation.

Drift-Free Position Estimation of Periodic or Quasi-Periodic Motion Using Inertial Sensors

PubMed Central

Latt, Win Tun; Veluvolu, Kalyana Chakravarthy; Ang, Wei Tech

2011-01-01

Position sensing with inertial sensors such as accelerometers and gyroscopes usually requires other aided sensors or prior knowledge of motion characteristics to remove position drift resulting from integration of acceleration or velocity so as to obtain accurate position estimation. A method based on analytical integration has previously been developed to obtain accurate position estimate of periodic or quasi-periodic motion from inertial sensors using prior knowledge of the motion but without using aided sensors. In this paper, a new method is proposed which employs linear filtering stage coupled with adaptive filtering stage to remove drift and attenuation. The prior knowledge of the motion the proposed method requires is only approximate band of frequencies of the motion. Existing adaptive filtering methods based on Fourier series such as weighted-frequency Fourier linear combiner (WFLC), and band-limited multiple Fourier linear combiner (BMFLC) are modified to combine with the proposed method. To validate and compare the performance of the proposed method with the method based on analytical integration, simulation study is performed using periodic signals as well as real physiological tremor data, and real-time experiments are conducted using an ADXL-203 accelerometer. Results demonstrate that the performance of the proposed method outperforms the existing analytical integration method. PMID:22163935

Environmental isotope investigation of groundwater flow in the Honey Lake Basin, California and Nevada

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

Rose, T.P.; Davisson, M.L.; Hudson, G.B.

The hydrology of Honey Lake Basin was studied using environmental isotope measurements of approximately 130 water samples collected during 1995 and 1996. The principal analytical methods included hydrogen, oxygen and carbon stable isotope ratio measurements, radiocarbon and tritium dating, and measurements of dissolved noble gas abundances.

Finite element modeling of light propagation in turbid media under illumination of a continuous-wave beam

USDA-ARS?s Scientific Manuscript database

Spatially-resolved spectroscopy provides a means for measuring the optical properties of biological tissues, based on analytical solutions to diffusion approximation for semi-infinite media under the normal illumination of infinitely small size light beam. The method is, however, prone to error in m...

Fatal overdose from Bendectin.

PubMed

Bayley, M; Walsh, F M; Valaske, M J

1975-05-01

A three-year old male ingested approximately 100 tablets of Bendectin. He developed tonic-clonic seizures followed by cardiac arrest. Toxicologic analysis yielded high levels of doxylamine, dicyclomine, and pyridoxine in blood, peritoneal fluid, and tissue homogenates. The antihistamine, doxylamine succinate appears to be the toxic constituent. Analytical methods used to document the case are herein described.

Quasi-analytical treatment of spatially averaged radiation transfer in complex terrain

NASA Astrophysics Data System (ADS)

Löwe, H.; Helbig, N.

2012-04-01

We provide a new quasi-analytical method to compute the topographic influence on the effective albedo of complex topography as required for meteorological, land-surface or climate models. We investigate radiative transfer in complex terrain via the radiosity equation on isotropic Gaussian random fields. Under controlled approximations we derive expressions for domain averages of direct, diffuse and terrain radiation and the sky view factor. Domain averaged quantities are related to a type of level-crossing probability of the random field which is approximated by longstanding results developed for acoustic scattering at ocean boundaries. This allows us to express all non-local horizon effects in terms of a local terrain parameter, namely the mean squared slope. Emerging integrals are computed numerically and fit formulas are given for practical purposes. As an implication of our approach we provide an expression for the effective albedo of complex terrain in terms of the sun elevation angle, mean squared slope, the area averaged surface albedo, and the direct-to-diffuse ratio of solar radiation. As an application, we compute the effective albedo for the Swiss Alps and discuss possible generalizations of the method.

Mixing of two co-directional Rayleigh surface waves in a nonlinear elastic material.

PubMed

Morlock, Merlin B; Kim, Jin-Yeon; Jacobs, Laurence J; Qu, Jianmin

2015-01-01

The mixing of two co-directional, initially monochromatic Rayleigh surface waves in an isotropic, homogeneous, and nonlinear elastic solid is investigated using analytical, finite element method, and experimental approaches. The analytical investigations show that while the horizontal velocity component can form a shock wave, the vertical velocity component can form a pulse independent of the specific ratios of the fundamental frequencies and amplitudes that are mixed. This analytical model is then used to simulate the development of the fundamentals, second harmonics, and the sum and difference frequency components over the propagation distance. The analytical model is further extended to include diffraction effects in the parabolic approximation. Finally, the frequency and amplitude ratios of the fundamentals are identified which provide maximum amplitudes of the second harmonics as well as of the sum and difference frequency components, to help guide effective material characterization; this approach should make it possible to measure the acoustic nonlinearity of a solid not only with the second harmonics, but also with the sum and difference frequency components. Results of the analytical investigations are then confirmed using the finite element method and the experimental feasibility of the proposed technique is validated for an aluminum specimen.

Analytical approximations for effective relative permeability in the capillary limit

NASA Astrophysics Data System (ADS)

Rabinovich, Avinoam; Li, Boxiao; Durlofsky, Louis J.

2016-10-01

We present an analytical method for calculating two-phase effective relative permeability, krjeff, where j designates phase (here CO2 and water), under steady state and capillary-limit assumptions. These effective relative permeabilities may be applied in experimental settings and for upscaling in the context of numerical flow simulations, e.g., for CO2 storage. An exact solution for effective absolute permeability, keff, in two-dimensional log-normally distributed isotropic permeability (k) fields is the geometric mean. We show that this does not hold for krjeff since log normality is not maintained in the capillary-limit phase permeability field (Kj=k·krj) when capillary pressure, and thus the saturation field, is varied. Nevertheless, the geometric mean is still shown to be suitable for approximating krjeff when the variance of ln⁡k is low. For high-variance cases, we apply a correction to the geometric average gas effective relative permeability using a Winsorized mean, which neglects large and small Kj values symmetrically. The analytical method is extended to anisotropically correlated log-normal permeability fields using power law averaging. In these cases, the Winsorized mean treatment is applied to the gas curves for cases described by negative power law exponents (flow across incomplete layers). The accuracy of our analytical expressions for krjeff is demonstrated through extensive numerical tests, using low-variance and high-variance permeability realizations with a range of correlation structures. We also present integral expressions for geometric-mean and power law average krjeff for the systems considered, which enable derivation of closed-form series solutions for krjeff without generating permeability realizations.

Partially Coherent Scattering in Stellar Chromospheres. Part 4; Analytic Wing Approximations

NASA Technical Reports Server (NTRS)

Gayley, K. G.

1993-01-01

Simple analytic expressions are derived to understand resonance-line wings in stellar chromospheres and similar astrophysical plasmas. The results are approximate, but compare well with accurate numerical simulations. The redistribution is modeled using an extension of the partially coherent scattering approximation (PCS) which we term the comoving-frame partially coherent scattering approximation (CPCS). The distinction is made here because Doppler diffusion is included in the coherent/noncoherent decomposition, in a form slightly improved from the earlier papers in this series.

Approximation methods of European option pricing in multiscale stochastic volatility model

NASA Astrophysics Data System (ADS)

Ni, Ying; Canhanga, Betuel; Malyarenko, Anatoliy; Silvestrov, Sergei

2017-01-01

In the classical Black-Scholes model for financial option pricing, the asset price follows a geometric Brownian motion with constant volatility. Empirical findings such as volatility smile/skew, fat-tailed asset return distributions have suggested that the constant volatility assumption might not be realistic. A general stochastic volatility model, e.g. Heston model, GARCH model and SABR volatility model, in which the variance/volatility itself follows typically a mean-reverting stochastic process, has shown to be superior in terms of capturing the empirical facts. However in order to capture more features of the volatility smile a two-factor, of double Heston type, stochastic volatility model is more useful as shown in Christoffersen, Heston and Jacobs [12]. We consider one modified form of such two-factor volatility models in which the volatility has multiscale mean-reversion rates. Our model contains two mean-reverting volatility processes with a fast and a slow reverting rate respectively. We consider the European option pricing problem under one type of the multiscale stochastic volatility model where the two volatility processes act as independent factors in the asset price process. The novelty in this paper is an approximating analytical solution using asymptotic expansion method which extends the authors earlier research in Canhanga et al. [5, 6]. In addition we propose a numerical approximating solution using Monte-Carlo simulation. For completeness and for comparison we also implement the semi-analytical solution by Chiarella and Ziveyi [11] using method of characteristics, Fourier and bivariate Laplace transforms.

An Approximate Solution to the Equation of Motion for Large-Angle Oscillations of the Simple Pendulum with Initial Velocity

ERIC Educational Resources Information Center

Johannessen, Kim

2010-01-01

An analytic approximation of the solution to the differential equation describing the oscillations of a simple pendulum at large angles and with initial velocity is discussed. In the derivation, a sinusoidal approximation has been applied, and an analytic formula for the large-angle period of the simple pendulum is obtained, which also includes…

Accelerating Electrostatic Surface Potential Calculation with Multiscale Approximation on Graphics Processing Units

PubMed Central

Anandakrishnan, Ramu; Scogland, Tom R. W.; Fenley, Andrew T.; Gordon, John C.; Feng, Wu-chun; Onufriev, Alexey V.

2010-01-01

Tools that compute and visualize biomolecular electrostatic surface potential have been used extensively for studying biomolecular function. However, determining the surface potential for large biomolecules on a typical desktop computer can take days or longer using currently available tools and methods. Two commonly used techniques to speed up these types of electrostatic computations are approximations based on multi-scale coarse-graining and parallelization across multiple processors. This paper demonstrates that for the computation of electrostatic surface potential, these two techniques can be combined to deliver significantly greater speed-up than either one separately, something that is in general not always possible. Specifically, the electrostatic potential computation, using an analytical linearized Poisson Boltzmann (ALPB) method, is approximated using the hierarchical charge partitioning (HCP) multiscale method, and parallelized on an ATI Radeon 4870 graphical processing unit (GPU). The implementation delivers a combined 934-fold speed-up for a 476,040 atom viral capsid, compared to an equivalent non-parallel implementation on an Intel E6550 CPU without the approximation. This speed-up is significantly greater than the 42-fold speed-up for the HCP approximation alone or the 182-fold speed-up for the GPU alone. PMID:20452792

Ergodicity of the Stochastic Nosé-Hoover Heat Bath

NASA Astrophysics Data System (ADS)

Wei Chung Lo,; Baowen Li,

2010-07-01

We numerically study the ergodicity of the stochastic Nosé-Hoover heat bath whose formalism is based on the Markovian approximation for the Nosé-Hoover equation [J. Phys. Soc. Jpn. 77 (2008) 103001]. The approximation leads to a Langevin-like equation driven by a fluctuating dissipative force and multiplicative Gaussian white noise. The steady state solution of the associated Fokker-Planck equation is the canonical distribution. We investigate the dynamics of this method for the case of (i) free particle, (ii) nonlinear oscillators and (iii) lattice chains. We derive the Fokker-Planck equation for the free particle and present approximate analytical solution for the stationary distribution in the context of the Markovian approximation. Numerical simulation results for nonlinear oscillators show that this method results in a Gaussian distribution for the particles velocity. We also employ the method as heat baths to study nonequilibrium heat flow in one-dimensional Fermi-Pasta-Ulam (FPU-β) and Frenkel-Kontorova (FK) lattices. The establishment of well-defined temperature profiles are observed only when the lattice size is large. Our results provide numerical justification for such Markovian approximation for classical single- and many-body systems.

Approximating natural connectivity of scale-free networks based on largest eigenvalue

NASA Astrophysics Data System (ADS)

Tan, S.-Y.; Wu, J.; Li, M.-J.; Lu, X.

2016-06-01

It has been recently proposed that natural connectivity can be used to efficiently characterize the robustness of complex networks. The natural connectivity has an intuitive physical meaning and a simple mathematical formulation, which corresponds to an average eigenvalue calculated from the graph spectrum. However, as a network model close to the real-world system that widely exists, the scale-free network is found difficult to obtain its spectrum analytically. In this article, we investigate the approximation of natural connectivity based on the largest eigenvalue in both random and correlated scale-free networks. It is demonstrated that the natural connectivity of scale-free networks can be dominated by the largest eigenvalue, which can be expressed asymptotically and analytically to approximate natural connectivity with small errors. Then we show that the natural connectivity of random scale-free networks increases linearly with the average degree given the scaling exponent and decreases monotonically with the scaling exponent given the average degree. Moreover, it is found that, given the degree distribution, the more assortative a scale-free network is, the more robust it is. Experiments in real networks validate our methods and results.

Numerical evaluation of electromagnetic fields due to dipole antennas in the presence of stratified media

NASA Technical Reports Server (NTRS)

Tsang, L.; Brown, R.; Kong, J. A.; Simmons, G.

1974-01-01

Two numerical methods are used to evaluate the integrals that express the em fields due to dipole antennas radiating in the presence of a stratified medium. The first method is a direct integration by means of Simpson's rule. The second method is indirect and approximates the kernel of the integral by means of the fast Fourier transform. In contrast to previous analytical methods that applied only to two-layer cases the numerical methods can be used for any arbitrary number of layers with general properties.

Application of the variational-asymptotical method to composite plates

NASA Technical Reports Server (NTRS)

Hodges, Dewey H.; Lee, Bok W.; Atilgan, Ali R.

1992-01-01

A method is developed for the 3D analysis of laminated plate deformation which is an extension of a variational-asymptotical method by Atilgan and Hodges (1991). Both methods are based on the treatment of plate deformation by splitting the 3D analysis into linear through-the-thickness analysis and 2D plate analysis. Whereas the first technique tackles transverse shear deformation in the second asymptotical approximation, the present method simplifies its treatment and restricts it to the first approximation. Both analytical techniques are applied to the linear cylindrical bending problem, and the strain and stress distributions are derived and compared with those of the exact solution. The present theory provides more accurate results than those of the classical laminated-plate theory for the transverse displacement of 2-, 3-, and 4-layer cross-ply laminated plates. The method can give reliable estimates of the in-plane strain and displacement distributions.

Precision of Sensitivity in the Design Optimization of Indeterminate Structures

NASA Technical Reports Server (NTRS)

Patnaik, Surya N.; Pai, Shantaram S.; Hopkins, Dale A.

2006-01-01

Design sensitivity is central to most optimization methods. The analytical sensitivity expression for an indeterminate structural design optimization problem can be factored into a simple determinate term and a complicated indeterminate component. Sensitivity can be approximated by retaining only the determinate term and setting the indeterminate factor to zero. The optimum solution is reached with the approximate sensitivity. The central processing unit (CPU) time to solution is substantially reduced. The benefit that accrues from using the approximate sensitivity is quantified by solving a set of problems in a controlled environment. Each problem is solved twice: first using the closed-form sensitivity expression, then using the approximation. The problem solutions use the CometBoards testbed as the optimization tool with the integrated force method as the analyzer. The modification that may be required, to use the stiffener method as the analysis tool in optimization, is discussed. The design optimization problem of an indeterminate structure contains many dependent constraints because of the implicit relationship between stresses, as well as the relationship between the stresses and displacements. The design optimization process can become problematic because the implicit relationship reduces the rank of the sensitivity matrix. The proposed approximation restores the full rank and enhances the robustness of the design optimization method.

A 2D Gaussian-Beam-Based Method for Modeling the Dichroic Surfaces of Quasi-Optical Systems

NASA Astrophysics Data System (ADS)

Elis, Kevin; Chabory, Alexandre; Sokoloff, Jérôme; Bolioli, Sylvain

2016-08-01

In this article, we propose an approach in the spectral domain to treat the interaction of a field with a dichroic surface in two dimensions. For a Gaussian beam illumination of the surface, the reflected and transmitted fields are approximated by one reflected and one transmitted Gaussian beams. Their characteristics are determined by means of a matching in the spectral domain, which requires a second-order approximation of the dichroic surface response when excited by plane waves. This approximation is of the same order as the one used in Gaussian beam shooting algorithm to model curved interfaces associated with lenses, reflector, etc. The method uses general analytical formulations for the GBs that depend either on a paraxial or far-field approximation. Numerical experiments are led to test the efficiency of the method in terms of accuracy and computation time. They include a parametric study and a case for which the illumination is provided by a horn antenna. For the latter, the incident field is firstly expressed as a sum of Gaussian beams by means of Gabor frames.

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.

Thermodynamic properties and static structure factor for a Yukawa fluid in the mean spherical approximation.

PubMed

Montes-Perez, J; Cruz-Vera, A; Herrera, J N

2011-12-01

This work presents the full analytic expressions for the thermodynamic properties and the static structure factor for a hard sphere plus 1-Yukawa fluid within the mean spherical approximation. To obtain these properties of the fluid type Yukawa analytically it was necessary to solve an equation of fourth order for the scaling parameter on a large scale. The physical root of this equation was determined by imposing physical conditions. The results of this work are obtained from seminal papers of Blum and Høye. We show that is not necessary the use the series expansion to solve the equation for the scaling parameter. We applied our theoretical result to find the thermodynamic and the static structure factor for krypton. Our results are in good agreement with those obtained in an experimental form or by simulation using the Monte Carlo method.

Exact PDF equations and closure approximations for advective-reactive transport

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

Venturi, D.; Tartakovsky, Daniel M.; Tartakovsky, Alexandre M.

2013-06-01

Mathematical models of advection–reaction phenomena rely on advective flow velocity and (bio) chemical reaction rates that are notoriously random. By using functional integral methods, we derive exact evolution equations for the probability density function (PDF) of the state variables of the advection–reaction system in the presence of random transport velocity and random reaction rates with rather arbitrary distributions. These PDF equations are solved analytically for transport with deterministic flow velocity and a linear reaction rate represented mathematically by a heterog eneous and strongly-correlated random field. Our analytical solution is then used to investigate the accuracy and robustness of the recentlymore » proposed large-eddy diffusivity (LED) closure approximation [1]. We find that the solution to the LED-based PDF equation, which is exact for uncorrelated reaction rates, is accurate even in the presence of strong correlations and it provides an upper bound of predictive uncertainty.« less

A hybrid analytical model for open-circuit field calculation of multilayer interior permanent magnet machines

NASA Astrophysics Data System (ADS)

Zhang, Zhen; Xia, Changliang; Yan, Yan; Geng, Qiang; Shi, Tingna

2017-08-01

Due to the complicated rotor structure and nonlinear saturation of rotor bridges, it is difficult to build a fast and accurate analytical field calculation model for multilayer interior permanent magnet (IPM) machines. In this paper, a hybrid analytical model suitable for the open-circuit field calculation of multilayer IPM machines is proposed by coupling the magnetic equivalent circuit (MEC) method and the subdomain technique. In the proposed analytical model, the rotor magnetic field is calculated by the MEC method based on the Kirchhoff's law, while the field in the stator slot, slot opening and air-gap is calculated by subdomain technique based on the Maxwell's equation. To solve the whole field distribution of the multilayer IPM machines, the coupled boundary conditions on the rotor surface are deduced for the coupling of the rotor MEC and the analytical field distribution of the stator slot, slot opening and air-gap. The hybrid analytical model can be used to calculate the open-circuit air-gap field distribution, back electromotive force (EMF) and cogging torque of multilayer IPM machines. Compared with finite element analysis (FEA), it has the advantages of faster modeling, less computation source occupying and shorter time consuming, and meanwhile achieves the approximate accuracy. The analytical model is helpful and applicable for the open-circuit field calculation of multilayer IPM machines with any size and pole/slot number combination.

First-order analytic propagation of satellites in the exponential atmosphere of an oblate planet

NASA Astrophysics Data System (ADS)

Martinusi, Vladimir; Dell'Elce, Lamberto; Kerschen, Gaëtan

2017-04-01

The paper offers the fully analytic solution to the motion of a satellite orbiting under the influence of the two major perturbations, due to the oblateness and the atmospheric drag. The solution is presented in a time-explicit form, and takes into account an exponential distribution of the atmospheric density, an assumption that is reasonably close to reality. The approach involves two essential steps. The first one concerns a new approximate mathematical model that admits a closed-form solution with respect to a set of new variables. The second step is the determination of an infinitesimal contact transformation that allows to navigate between the new and the original variables. This contact transformation is obtained in exact form, and afterwards a Taylor series approximation is proposed in order to make all the computations explicit. The aforementioned transformation accommodates both perturbations, improving the accuracy of the orbit predictions by one order of magnitude with respect to the case when the atmospheric drag is absent from the transformation. Numerical simulations are performed for a low Earth orbit starting at an altitude of 350 km, and they show that the incorporation of drag terms into the contact transformation generates an error reduction by a factor of 7 in the position vector. The proposed method aims at improving the accuracy of analytic orbit propagation and transforming it into a viable alternative to the computationally intensive numerical methods.

On the arbitrary l-wave solutions of the deformed hyperbolic manning-rosen potential including an improved approximation to the orbital centrifugal term

NASA Astrophysics Data System (ADS)

Xu, Chun-Long; Zhang, Min-Cang

2017-01-01

The arbitrary l-wave solutions to the Schrödinger equation for the deformed hyperbolic Manning-Rosen potential is investigated analytically by using the Nikiforov-Uvarov method, the centrifugal term is treated with an improved Greene and Aldrich's approximation scheme. The wavefunctions depend on the deformation parameter q, which is expressed in terms of the Jocobi polynomial or the hypergeometric function. The bound state energy is obtained, and the discrete spectrum is shown to be independent of the deformation parameter q.

Nanophotonic particle simulation and inverse design using artificial neural networks

PubMed Central

Peurifoy, John; Shen, Yichen; Jing, Li; Cano-Renteria, Fidel; DeLacy, Brendan G.; Joannopoulos, John D.; Tegmark, Max

2018-01-01

We propose a method to use artificial neural networks to approximate light scattering by multilayer nanoparticles. We find that the network needs to be trained on only a small sampling of the data to approximate the simulation to high precision. Once the neural network is trained, it can simulate such optical processes orders of magnitude faster than conventional simulations. Furthermore, the trained neural network can be used to solve nanophotonic inverse design problems by using back propagation, where the gradient is analytical, not numerical. PMID:29868640

Solution of the mean spherical approximation for polydisperse multi-Yukawa hard-sphere fluid mixture using orthogonal polynomial expansions

NASA Astrophysics Data System (ADS)

Kalyuzhnyi, Yurij V.; Cummings, Peter T.

2006-03-01

The Blum-Høye [J. Stat. Phys. 19 317 (1978)] solution of the mean spherical approximation for a multicomponent multi-Yukawa hard-sphere fluid is extended to a polydisperse multi-Yukawa hard-sphere fluid. Our extension is based on the application of the orthogonal polynomial expansion method of Lado [Phys. Rev. E 54, 4411 (1996)]. Closed form analytical expressions for the structural and thermodynamic properties of the model are presented. They are given in terms of the parameters that follow directly from the solution. By way of illustration the method of solution is applied to describe the thermodynamic properties of the one- and two-Yukawa versions of the model.

Relativistic symmetries in the Rosen—Morse potential and tensor interaction using the Nikiforov—Uvarov method

NASA Astrophysics Data System (ADS)

Sameer, M. Ikhdair; Majid, Hamzavi

2013-04-01

Approximate analytical bound-state solutions of the Dirac particle in the fields of attractive and repulsive Rosen—Morse (RM) potentials including the Coulomb-like tensor (CLT) potential are obtained for arbitrary spin-orbit quantum number κ. The Pekeris approximation is used to deal with the spin-orbit coupling terms κ (κ± 1)r-2. In the presence of exact spin and pseudospin (p-spin) symmetries, the energy eigenvalues and the corresponding normalized two-component wave functions are found by using the parametric generalization of the Nikiforov—Uvarov (NU) method. The numerical results show that the CLT interaction removes degeneracies between the spin and p-spin state doublets.

Approximate analytical relationships for linear optimal aeroelastic flight control laws

NASA Astrophysics Data System (ADS)

Kassem, Ayman Hamdy

1998-09-01

This dissertation introduces new methods to uncover functional relationships between design parameters of a contemporary control design technique and the resulting closed-loop properties. Three new methods are developed for generating such relationships through analytical expressions: the Direct Eigen-Based Technique, the Order of Magnitude Technique, and the Cost Function Imbedding Technique. Efforts concentrated on the linear-quadratic state-feedback control-design technique applied to an aeroelastic flight control task. For this specific application, simple and accurate analytical expressions for the closed-loop eigenvalues and zeros in terms of basic parameters such as stability and control derivatives, structural vibration damping and natural frequency, and cost function weights are generated. These expressions explicitly indicate how the weights augment the short period and aeroelastic modes, as well as the closed-loop zeros, and by what physical mechanism. The analytical expressions are used to address topics such as damping, nonminimum phase behavior, stability, and performance with robustness considerations, and design modifications. This type of knowledge is invaluable to the flight control designer and would be more difficult to formulate when obtained from numerical-based sensitivity analysis.

Full analytical solution of the bloch equation when using a hyperbolic-secant driving function.

PubMed

Zhang, Jinjin; Garwood, Michael; Park, Jang-Yeon

2017-04-01

The frequency-swept pulse known as the hyperbolic-secant (HS) pulse is popular in NMR for achieving adiabatic spin inversion. The HS pulse has also shown utility for achieving excitation and refocusing in gradient-echo and spin-echo sequences, including new ultrashort echo-time imaging (e.g., Sweep Imaging with Fourier Transform, SWIFT) and B 1 mapping techniques. To facilitate the analysis of these techniques, the complete theoretical solution of the Bloch equation, as driven by the HS pulse, was derived for an arbitrary state of initial magnetization. The solution of the Bloch-Riccati equation for transverse and longitudinal magnetization for an arbitrary initial state was derived analytically in terms of HS pulse parameters. The analytical solution was compared with the solutions using both the Runge-Kutta method and the small-tip approximation. The analytical solution was demonstrated on different initial states at different frequency offsets with/without a combination of HS pulses. Evolution of the transverse magnetization was influenced significantly by the choice of HS pulse parameters. The deviation of the magnitude of the transverse magnetization, as obtained by comparing the small-tip approximation to the analytical solution, was < 5% for flip angles < 30 °, but > 10% for the flip angles > 40 °. The derived analytical solution provides insights into the influence of HS pulse parameters on the magnetization evolution. Magn Reson Med 77:1630-1638, 2017. © 2016 International Society for Magnetic Resonance in Medicine. © 2016 International Society for Magnetic Resonance in Medicine.

Analytic representations of mK , FK, mη, and Fη in two loop S U (3 ) chiral perturbation theory

NASA Astrophysics Data System (ADS)

Ananthanarayan, B.; Bijnens, Johan; Friot, Samuel; Ghosh, Shayan

2018-06-01

In this work, we consider expressions for the masses and decay constants of the pseudoscalar mesons in S U (3 ) chiral perturbation theory. These involve sunset diagrams and their derivatives evaluated at p2=mP2 (P =π , K , η ). Recalling that there are three mass scales in this theory, mπ, mK and mη, there are instances when the finite part of the sunset diagrams do not admit an expression in terms of elementary functions, and have therefore been evaluated numerically in the past. In a recent publication, an expansion in the external momentum was performed to obtain approximate analytic expressions for mπ and Fπ, the pion mass and decay constant. We provide fully analytic exact expressions for mK and mη, the kaon and eta masses, and FK and Fη, the kaon and eta decay constants. These expressions, calculated using Mellin-Barnes methods, are in the form of double series in terms of two mass ratios. A numerical analysis of the results to evaluate the relative size of contributions coming from loops, chiral logarithms as well as phenomenological low-energy constants is presented. We also present a set of approximate analytic expressions for mK, FK, mη and Fη that facilitate comparisons with lattice results. Finally, we show how exact analytic expressions for mπ and Fπ may be obtained, the latter having been used in conjunction with the results for FK to produce a recently published analytic representation of FK/Fπ.

Towards the blackbox computation of magnetic exchange coupling parameters in polynuclear transition-metal complexes: theory, implementation, and application.

PubMed

Phillips, Jordan J; Peralta, Juan E

2013-05-07

We present a method for calculating magnetic coupling parameters from a single spin-configuration via analytic derivatives of the electronic energy with respect to the local spin direction. This method does not introduce new approximations beyond those found in the Heisenberg-Dirac Hamiltonian and a standard Kohn-Sham Density Functional Theory calculation, and in the limit of an ideal Heisenberg system it reproduces the coupling as determined from spin-projected energy-differences. Our method employs a generalized perturbative approach to constrained density functional theory, where exact expressions for the energy to second order in the constraints are obtained by analytic derivatives from coupled-perturbed theory. When the relative angle between magnetization vectors of metal atoms enters as a constraint, this allows us to calculate all the magnetic exchange couplings of a system from derivatives with respect to local spin directions from the high-spin configuration. Because of the favorable computational scaling of our method with respect to the number of spin-centers, as compared to the broken-symmetry energy-differences approach, this opens the possibility for the blackbox exploration of magnetic properties in large polynuclear transition-metal complexes. In this work we outline the motivation, theory, and implementation of this method, and present results for several model systems and transition-metal complexes with a variety of density functional approximations and Hartree-Fock.

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

McNamara, B.

Tandem and stellarator equilibria at high ..beta.. have proved hard to compute and the relaxation methods of Bauer et al., Chodura and Schluter, Hirshman, Strauss, and Pearlstein et al. have been slow to converge. This paper reports an extension of the low-..beta.. analytic method of Pearlstein, Kaiser, and Newcomb to arbitrary ..beta.. for tandem mirrors which converges in 10 to 20 iterations. Extensions of the method to stellarator equilibria are proposed and are very close to the analytic method of Johnson and Greene - the stellarator expansion. Most of the results of all these calculations can be adequately described bymore » low-..beta.. approximations since the MHD stability limits occur at low ..beta... The tandem mirror, having weak curvature and a long central cell, allows finite Larmor radius effects to eliminate most ballooning modes and offers the possibility of really high average ..beta... This is the interest in developing such three-dimensional numerical algorithms.« less

Derivation of phase functions from multiply scattered sunlight transmitted through a hazy atmosphere

NASA Technical Reports Server (NTRS)

Weinman, J. A.; Twitty, J. T.; Browning, S. R.; Herman, B. M.

1975-01-01

The intensity of sunlight multiply scattered in model atmospheres is derived from the equation of radiative transfer by an analytical small-angle approximation. The approximate analytical solutions are compared to rigorous numerical solutions of the same problem. Results obtained from an aerosol-laden model atmosphere are presented. Agreement between the rigorous and the approximate solutions is found to be within a few per cent. The analytical solution to the problem which considers an aerosol-laden atmosphere is then inverted to yield a phase function which describes a single scattering event at small angles. The effect of noisy data on the derived phase function is discussed.

Optimal clinical trial design based on a dichotomous Markov-chain mixed-effect sleep model.

PubMed

Steven Ernest, C; Nyberg, Joakim; Karlsson, Mats O; Hooker, Andrew C

2014-12-01

D-optimal designs for discrete-type responses have been derived using generalized linear mixed models, simulation based methods and analytical approximations for computing the fisher information matrix (FIM) of non-linear mixed effect models with homogeneous probabilities over time. In this work, D-optimal designs using an analytical approximation of the FIM for a dichotomous, non-homogeneous, Markov-chain phase advanced sleep non-linear mixed effect model was investigated. The non-linear mixed effect model consisted of transition probabilities of dichotomous sleep data estimated as logistic functions using piecewise linear functions. Theoretical linear and nonlinear dose effects were added to the transition probabilities to modify the probability of being in either sleep stage. D-optimal designs were computed by determining an analytical approximation the FIM for each Markov component (one where the previous state was awake and another where the previous state was asleep). Each Markov component FIM was weighted either equally or by the average probability of response being awake or asleep over the night and summed to derive the total FIM (FIM(total)). The reference designs were placebo, 0.1, 1-, 6-, 10- and 20-mg dosing for a 2- to 6-way crossover study in six dosing groups. Optimized design variables were dose and number of subjects in each dose group. The designs were validated using stochastic simulation/re-estimation (SSE). Contrary to expectations, the predicted parameter uncertainty obtained via FIM(total) was larger than the uncertainty in parameter estimates computed by SSE. Nevertheless, the D-optimal designs decreased the uncertainty of parameter estimates relative to the reference designs. Additionally, the improvement for the D-optimal designs were more pronounced using SSE than predicted via FIM(total). Through the use of an approximate analytic solution and weighting schemes, the FIM(total) for a non-homogeneous, dichotomous Markov-chain phase advanced sleep model was computed and provided more efficient trial designs and increased nonlinear mixed-effects modeling parameter precision.

Representing Functions in n Dimensions to Arbitrary Accuracy

NASA Technical Reports Server (NTRS)

Scotti, Stephen J.

2007-01-01

A method of approximating a scalar function of n independent variables (where n is a positive integer) to arbitrary accuracy has been developed. This method is expected to be attractive for use in engineering computations in which it is necessary to link global models with local ones or in which it is necessary to interpolate noiseless tabular data that have been computed from analytic functions or numerical models in n-dimensional spaces of design parameters.

Reemission spectra and inelastic processes at interaction of attosecond and shorter duration electromagnetic pulses with atoms

NASA Astrophysics Data System (ADS)

Makarov, D. N.; Matveev, V. I.

2017-01-01

Inelastic processes and the reemission of attosecond and shorter electromagnetic pulses by atoms have been considered within the analytical solution of the Schrödinger equation in the sudden perturbation approximation. A method of calculations with the exact inclusion of spatial inhomogeneity of the field of an ultrashort pulse and the momenta of photons in the reemission processes has been developed. The probabilities of inelastic processes and spectra of reemission of ultrashort electromagnetic pulses by one- and many-electron atoms have been calculated. The results have been presented in the form of analytical formulas.

A simplified analytical solution for thermal response of a one-dimensional, steady state transpiration cooling system in radiative and convective environment

NASA Technical Reports Server (NTRS)

Kubota, H.

1976-01-01

A simplified analytical method for calculation of thermal response within a transpiration-cooled porous heat shield material in an intense radiative-convective heating environment is presented. The essential assumptions of the radiative and convective transfer processes in the heat shield matrix are the two-temperature approximation and the specified radiative-convective heatings of the front surface. Sample calculations for porous silica with CO2 injection are presented for some typical parameters of mass injection rate, porosity, and material thickness. The effect of these parameters on the cooling system is discussed.

Analytical solutions of the Dirac equation under Hellmann–Frost–Musulin potential

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

Onate, C.A., E-mail: oaclems14@physicist.net; Onyeaju, M.C.; Ikot, A.N.

2016-12-15

The approximate analytical solutions of the Dirac equation with Hellmann–Frost–Musulin potential have been studied by using the generalized parametric Nikiforov–Uvarov (NU) method for arbitrary spin–orbit quantum number k under the spin and pseudospin symmetries. The Hellmann–Frost–Musulin potential is a superposition potential that consists of Yukawa potential, Coulomb potential, and Frost–Musulin potential. As a particular case, we found the energy levels of the non-relativistic limit of the spin symmetry. The energy equation of Yukawa potential, Coulomb potential, Hellmann potential and Frost–Musulin potential are obtained. Energy values are generated for some diatomic molecules.

Analytical and numerical treatment of the heat conduction equation obtained via time-fractional distributed-order heat conduction law

NASA Astrophysics Data System (ADS)

Želi, Velibor; Zorica, Dušan

2018-02-01

Generalization of the heat conduction equation is obtained by considering the system of equations consisting of the energy balance equation and fractional-order constitutive heat conduction law, assumed in the form of the distributed-order Cattaneo type. The Cauchy problem for system of energy balance equation and constitutive heat conduction law is treated analytically through Fourier and Laplace integral transform methods, as well as numerically by the method of finite differences through Adams-Bashforth and Grünwald-Letnikov schemes for approximation derivatives in temporal domain and leap frog scheme for spatial derivatives. Numerical examples, showing time evolution of temperature and heat flux spatial profiles, demonstrate applicability and good agreement of both methods in cases of multi-term and power-type distributed-order heat conduction laws.

Numerical simulation of KdV equation by finite difference method

NASA Astrophysics Data System (ADS)

Yokus, A.; Bulut, H.

2018-05-01

In this study, the numerical solutions to the KdV equation with dual power nonlinearity by using the finite difference method are obtained. Discretize equation is presented in the form of finite difference operators. The numerical solutions are secured via the analytical solution to the KdV equation with dual power nonlinearity which is present in the literature. Through the Fourier-Von Neumann technique and linear stable, we have seen that the FDM is stable. Accuracy of the method is analyzed via the L2 and L_{∞} norm errors. The numerical, exact approximations and absolute error are presented in tables. We compare the numerical solutions with the exact solutions and this comparison is supported with the graphic plots. Under the choice of suitable values of parameters, the 2D and 3D surfaces for the used analytical solution are plotted.

Modeling Sound Propagation Through Non-Axisymmetric Jets

NASA Technical Reports Server (NTRS)

Leib, Stewart J.

2014-01-01

A method for computing the far-field adjoint Green's function of the generalized acoustic analogy equations under a locally parallel mean flow approximation is presented. The method is based on expanding the mean-flow-dependent coefficients in the governing equation and the scalar Green's function in truncated Fourier series in the azimuthal direction and a finite difference approximation in the radial direction in circular cylindrical coordinates. The combined spectral/finite difference method yields a highly banded system of algebraic equations that can be efficiently solved using a standard sparse system solver. The method is applied to test cases, with mean flow specified by analytical functions, corresponding to two noise reduction concepts of current interest: the offset jet and the fluid shield. Sample results for the Green's function are given for these two test cases and recommendations made as to the use of the method as part of a RANS-based jet noise prediction code.

Development of Finite Elements for Two-Dimensional Structural Analysis Using the Integrated Force Method

NASA Technical Reports Server (NTRS)

Kaljevic, Igor; Patnaik, Surya N.; Hopkins, Dale A.

1996-01-01

The Integrated Force Method has been developed in recent years for the analysis of structural mechanics problems. This method treats all independent internal forces as unknown variables that can be calculated by simultaneously imposing equations of equilibrium and compatibility conditions. In this paper a finite element library for analyzing two-dimensional problems by the Integrated Force Method is presented. Triangular- and quadrilateral-shaped elements capable of modeling arbitrary domain configurations are presented. The element equilibrium and flexibility matrices are derived by discretizing the expressions for potential and complementary energies, respectively. The displacement and stress fields within the finite elements are independently approximated. The displacement field is interpolated as it is in the standard displacement method, and the stress field is approximated by using complete polynomials of the correct order. A procedure that uses the definitions of stress components in terms of an Airy stress function is developed to derive the stress interpolation polynomials. Such derived stress fields identically satisfy the equations of equilibrium. Moreover, the resulting element matrices are insensitive to the orientation of local coordinate systems. A method is devised to calculate the number of rigid body modes, and the present elements are shown to be free of spurious zero-energy modes. A number of example problems are solved by using the present library, and the results are compared with corresponding analytical solutions and with results from the standard displacement finite element method. The Integrated Force Method not only gives results that agree well with analytical and displacement method results but also outperforms the displacement method in stress calculations.

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.

On the design of recursive digital filters

NASA Technical Reports Server (NTRS)

Shenoi, K.; Narasimha, M. J.; Peterson, A. M.

1976-01-01

A change of variables is described which transforms the problem of designing a recursive digital filter to that of approximation by a ratio of polynomials on a finite interval. Some analytic techniques for the design of low-pass filters are presented, illustrating the use of the transformation. Also considered are methods for the design of phase equalizers.

Flame characteristics for fires in southern fuels

Treesearch

Ralph M. Nelson

1980-01-01

A flame model and analytical method are used to derive forest fire flame characteristics. Approximate solutions are used to express flame lengths, angles, heights, and tip velocities of headfires and calm-air fires in terms of fire intensity. Equations are compared with data from low-intensity controlled burns in southern fuels and with data from the literature.

Approximate analysis of thermal convection in a crystal-growth cell for Spacelab 3

NASA Technical Reports Server (NTRS)

Dressler, R. F.

1982-01-01

The transient and steady thermal convection in microgravity is described. The approach is applicable to many three dimensional flows in containers of various shapes with various thermal gradients imposed. The method employs known analytical solutions to two dimensional thermal flows in simpler geometries, and does not require recourse to numerical calculations by computer.

APPROXIMATE AND ANALYTICAL SOLUTIONS FOR SOLUTE TRANSPORT FROM AN INJECTION WELL INTO A SINGLE FRACTURE

EPA Science Inventory

In dealing with problems related to land-based nuclear waste management, a number of analytical and approximate solutions were developed to quantify radionuclide transport through fractures contained in the porous formation. t has been reported that by treating the radioactive de...

An analytical solution of Richards' equation providing the physical basis of SCS curve number method and its proportionality relationship

NASA Astrophysics Data System (ADS)

Hooshyar, Milad; Wang, Dingbao

2016-08-01

The empirical proportionality relationship, which indicates that the ratio of cumulative surface runoff and infiltration to their corresponding potentials are equal, is the basis of the extensively used Soil Conservation Service Curve Number (SCS-CN) method. The objective of this paper is to provide the physical basis of the SCS-CN method and its proportionality hypothesis from the infiltration excess runoff generation perspective. To achieve this purpose, an analytical solution of Richards' equation is derived for ponded infiltration in shallow water table environment under the following boundary conditions: (1) the soil is saturated at the land surface; and (2) there is a no-flux boundary which moves downward. The solution is established based on the assumptions of negligible gravitational effect, constant soil water diffusivity, and hydrostatic soil moisture profile between the no-flux boundary and water table. Based on the derived analytical solution, the proportionality hypothesis is a reasonable approximation for rainfall partitioning at the early stage of ponded infiltration in areas with a shallow water table for coarse textured soils.

Simultaneous determination of fluoride, chloride, sulfate, phosphate, monofluorophosphate, glycerophosphate, sorbate, and saccharin in gargles by ion chromatography.

PubMed

Zhang, Yan-zhen; Zhou, Yan-chun; Liu, Li; Zhu, Yan

2007-07-01

Simple, reliable and sensitive analytical methods to determine anticariogenic agents, preservatives, and artificial sweeteners contained in commercial gargles are necessary for evaluating their effectiveness, safety, and quality. An ion chromatography (IC) method has been described to analyze simultaneously eight anions including fluoride, chloride, sulfate, phosphate, monofluorophosphate, glycerophosphate (anticariogenic agents), sorbate (a preservative), and saccharin (an artificial sweetener) in gargles. In this IC system, we applied a mobile phased gradient elution with KOH, separation by IonPac AS18 columns, and suppressed conductivity detection. Optimized analytical conditions were further evaluated for accuracy. The relative standard deviations (RSDs) of the inter-day's retention time and peak area of all species were less than 0.938% and 8.731%, respectively, while RSDs of 5-day retention time and peak area were less than 1.265% and 8.934%, respectively. The correlation coefficients for targeted analytes ranged from 0.999 7 to 1.000 0. The spiked recoveries for the anions were 90% approximately 102.5%. We concluded that the method can be applied for comprehensive evaluation of commercial gargles.

Analytical treatment of self-phase-modulation beyond the slowly varying envelope approximation

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

Syrchin, M.S.; Zheltikov, A.M.; International Laser Center, M.V. Lomonosov Moscow State University, 119899 Moscow

Analytical treatment of the self-phase-modulation of an ultrashort light pulse is extended beyond the slowly varying envelope approximation. The resulting wave equation is modified to include corrections to self-phase-modulation due to higher-order spatial and temporal derivatives. Analytical solutions are found in the limiting regimes of high nonlinearities and very short pulses. Our results reveal features that can significantly impact both pulse shape and the evolution of the phase.

An Incompressible, Depth-Averaged Lattice Boltzmann Method for Liquid Flow in Microfluidic Devices with Variable Aperture

DOE PAGES

Laleian, Artin; Valocchi, Albert J.; Werth, Charles J.

2015-11-24

Two-dimensional (2D) pore-scale models have successfully simulated microfluidic experiments of aqueous-phase flow with mixing-controlled reactions in devices with small aperture. A standard 2D model is not generally appropriate when the presence of mineral precipitate or biomass creates complex and irregular three-dimensional (3D) pore geometries. We modify the 2D lattice Boltzmann method (LBM) to incorporate viscous drag from the top and bottom microfluidic device (micromodel) surfaces, typically excluded in a 2D model. Viscous drag from these surfaces can be approximated by uniformly scaling a steady-state 2D velocity field at low Reynolds number. We demonstrate increased accuracy by approximating the viscous dragmore » with an analytically-derived body force which assumes a local parabolic velocity profile across the micromodel depth. Accuracy of the generated 2D velocity field and simulation permeability have not been evaluated in geometries with variable aperture. We obtain permeabilities within approximately 10% error and accurate streamlines from the proposed 2D method relative to results obtained from 3D simulations. Additionally, the proposed method requires a CPU run time approximately 40 times less than a standard 3D method, representing a significant computational benefit for permeability calculations.« less

Multiple Scattering Effects on Pulse Propagation in Optically Turbid Media.

NASA Astrophysics Data System (ADS)

Joelson, Bradley David

The effects of multiple scattering in a optically turbid media is examined for an impulse solution to the radiative transfer equation for a variety of geometries and phase functions. In regions where the complexities of the phase function proved too cumbersome for analytic methods Monte Carlo techniques were developed to describe the entire scalar radiance distribution. The determination of a general spread function is strongly dependent on geometry and particular regions where limits can be placed on the variables of the problem. Hence, the general spread function is first simplified by considering optical regions which reduce the complexity of the variable dependence. First, in the small-angle limit we calculate some contracted spread functions along with their moments and then use Monte Carlo techniques to establish the limitations imposed by the small-angle approximation in planar geometry. The point spread function (PSF) for a spherical geometry is calculated for the full angular spread in the forward direction of ocean waters using Monte Carlo methods in the optically thin and moderate depths and analytic methods in the diffusion domain. The angular dependence of the PSF for various ocean waters is examined for a range of optical parameters. The analytic method used in the diffusion calculation is justified by examining the angular dependence of the radiance of a impulse solution in a planar geometry for a prolongated Henyey-Greenstein phase function of asymmetry factor approximately equal to that of the ocean phase functions. The Legendre moments of the radiance are examined in order to examine the viability of the diffusion approximation which assumes a linearly anisotropic angular distribution for the radiance. A realistic lidar calculation is performed for a variety of ocean waters to determine the effects of multiple scattering on the determination of the speed of sound by using the range gated frequency spectrum of the lidar signal. It is shown that the optical properties of the ocean help to ensure single scatter form for the frequency spectra of the lidar signal. This spectra can then be used to compute the speed of sound and backscatter probability.

Analytical and numerical treatment of drift-tearing modes in plasma slab

NASA Astrophysics Data System (ADS)

Mirnov, V. V.; Hegna, C. C.; Sovinec, C. R.; Howell, E. C.

2016-10-01

Two-fluid corrections to linear tearing modes includes 1) diamagnetic drifts that reduce the growth rate and 2) electron and ion decoupling on short scales that can lead to fast reconnection. We have recently developed an analytical model that includes effects 1) and 2) and important contribution from finite electron parallel thermal conduction. Both the tendencies 1) and 2) are confirmed by an approximate analytic dispersion relation that is derived using a perturbative approach of small ion-sound gyroradius ρs. This approach is only valid at the beginning of the transition from the collisional to semi-collisional regimes. Further analytical and numerical work is performed to cover the full interval of ρs connecting these two limiting cases. Growth rates are computed from analytic theory with a shooting method. They match the resistive MHD regime with the dispersion relations known at asymptotically large ion-sound gyroradius. A comparison between this analytical treatment and linear numerical simulations using the NIMROD code with cold ions and hot electrons in plasma slab is reported. The material is based on work supported by the U.S. DOE and NSF.

Optical properties of electrohydrodynamic convection patterns: rigorous and approximate methods.

PubMed

Bohley, Christian; Heuer, Jana; Stannarius, Ralf

2005-12-01

We analyze the optical behavior of two-dimensionally periodic structures that occur in electrohydrodynamic convection (EHC) patterns in nematic sandwich cells. These structures are anisotropic, locally uniaxial, and periodic on the scale of micrometers. For the first time, the optics of these structures is investigated with a rigorous method. The method used for the description of the electromagnetic waves interacting with EHC director patterns is a numerical approach that discretizes directly the Maxwell equations. It works as a space-grid-time-domain method and computes electric and magnetic fields in time steps. This so-called finite-difference-time-domain (FDTD) method is able to generate the fields with arbitrary accuracy. We compare this rigorous method with earlier attempts based on ray-tracing and analytical approximations. Results of optical studies of EHC structures made earlier based on ray-tracing methods are confirmed for thin cells, when the spatial periods of the pattern are sufficiently large. For the treatment of small-scale convection structures, the FDTD method is without alternatives.

Modeling and Analysis of Structural Dynamics for a One-Tenth Scale Model NGST Sunshield

NASA Technical Reports Server (NTRS)

Johnston, John; Lienard, Sebastien; Brodeur, Steve (Technical Monitor)

2001-01-01

New modeling and analysis techniques have been developed for predicting the dynamic behavior of the Next Generation Space Telescope (NGST) sunshield. The sunshield consists of multiple layers of pretensioned, thin-film membranes supported by deployable booms. Modeling the structural dynamic behavior of the sunshield is a challenging aspect of the problem due to the effects of membrane wrinkling. A finite element model of the sunshield was developed using an approximate engineering approach, the cable network method, to account for membrane wrinkling effects. Ground testing of a one-tenth scale model of the NGST sunshield were carried out to provide data for validating the analytical model. A series of analyses were performed to predict the behavior of the sunshield under the ground test conditions. Modal analyses were performed to predict the frequencies and mode shapes of the test article and transient response analyses were completed to simulate impulse excitation tests. Comparison was made between analytical predictions and test measurements for the dynamic behavior of the sunshield. In general, the results show good agreement with the analytical model correctly predicting the approximate frequency and mode shapes for the significant structural modes.

A discussion on validity of the diffusion theory by Monte Carlo method

NASA Astrophysics Data System (ADS)

Peng, Dong-qing; Li, Hui; Xie, Shusen

2008-12-01

Diffusion theory was widely used as a basis of the experiments and methods in determining the optical properties of biological tissues. A simple analytical solution could be obtained easily from the diffusion equation after a series of approximations. Thus, a misinterpret of analytical solution would be made: while the effective attenuation coefficient of several semi-infinite bio-tissues were the same, the distribution of light fluence in the tissues would be the same. In order to assess the validity of knowledge above, depth resolved internal fluence of several semi-infinite biological tissues which have the same effective attenuation coefficient were simulated with wide collimated beam in the paper by using Monte Carlo method in different condition. Also, the influence of bio-tissue refractive index on the distribution of light fluence was discussed in detail. Our results showed that, when the refractive index of several bio-tissues which had the same effective attenuation coefficient were the same, the depth resolved internal fluence would be the same; otherwise, the depth resolved internal fluence would be not the same. The change of refractive index of tissue would have affection on the light depth distribution in tissue. Therefore, the refractive index is an important optical property of tissue, and should be taken in account while using the diffusion approximation theory.

Limited Rationality and Its Quantification Through the Interval Number Judgments With Permutations.

PubMed

Liu, Fang; Pedrycz, Witold; Zhang, Wei-Guo

2017-12-01

The relative importance of alternatives expressed in terms of interval numbers in the fuzzy analytic hierarchy process aims to capture the uncertainty experienced by decision makers (DMs) when making a series of comparisons. Under the assumption of full rationality, the judgements of DMs in the typical analytic hierarchy process could be consistent. However, since the uncertainty in articulating the opinions of DMs is unavoidable, the interval number judgements are associated with the limited rationality. In this paper, we investigate the concept of limited rationality by introducing interval multiplicative reciprocal comparison matrices. By analyzing the consistency of interval multiplicative reciprocal comparison matrices, it is observed that the interval number judgements are inconsistent. By considering the permutations of alternatives, the concepts of approximation-consistency and acceptable approximation-consistency of interval multiplicative reciprocal comparison matrices are proposed. The exchange method is designed to generate all the permutations. A novel method of determining the interval weight vector is proposed under the consideration of randomness in comparing alternatives, and a vector of interval weights is determined. A new algorithm of solving decision making problems with interval multiplicative reciprocal preference relations is provided. Two numerical examples are carried out to illustrate the proposed approach and offer a comparison with the methods available in the literature.

Bayesian inference based on stationary Fokker-Planck sampling.

PubMed

Berrones, Arturo

2010-06-01

A novel formalism for bayesian learning in the context of complex inference models is proposed. The method is based on the use of the stationary Fokker-Planck (SFP) approach to sample from the posterior density. Stationary Fokker-Planck sampling generalizes the Gibbs sampler algorithm for arbitrary and unknown conditional densities. By the SFP procedure, approximate analytical expressions for the conditionals and marginals of the posterior can be constructed. At each stage of SFP, the approximate conditionals are used to define a Gibbs sampling process, which is convergent to the full joint posterior. By the analytical marginals efficient learning methods in the context of artificial neural networks are outlined. Offline and incremental bayesian inference and maximum likelihood estimation from the posterior are performed in classification and regression examples. A comparison of SFP with other Monte Carlo strategies in the general problem of sampling from arbitrary densities is also presented. It is shown that SFP is able to jump large low-probability regions without the need of a careful tuning of any step-size parameter. In fact, the SFP method requires only a small set of meaningful parameters that can be selected following clear, problem-independent guidelines. The computation cost of SFP, measured in terms of loss function evaluations, grows linearly with the given model's dimension.

On the Accuracy of Double Scattering Approximation for Atmospheric Polarization Computations

NASA Technical Reports Server (NTRS)

Korkin, Sergey V.; Lyapustin, Alexei I.; Marshak, Alexander L.

2011-01-01

Interpretation of multi-angle spectro-polarimetric data in remote sensing of atmospheric aerosols require fast and accurate methods of solving the vector radiative transfer equation (VRTE). The single and double scattering approximations could provide an analytical framework for the inversion algorithms and are relatively fast, however accuracy assessments of these approximations for the aerosol atmospheres in the atmospheric window channels have been missing. This paper provides such analysis for a vertically homogeneous aerosol atmosphere with weak and strong asymmetry of scattering. In both cases, the double scattering approximation gives a high accuracy result (relative error approximately 0.2%) only for the low optical path - 10(sup -2) As the error rapidly grows with optical thickness, a full VRTE solution is required for the practical remote sensing analysis. It is shown that the scattering anisotropy is not important at low optical thicknesses neither for reflected nor for transmitted polarization components of radiation.

Phonon dispersion on Ag (100) surface: A modified analytic embedded atom method study

NASA Astrophysics Data System (ADS)

Xiao-Jun, Zhang; Chang-Le, Chen

2016-01-01

Within the harmonic approximation, the analytic expression of the dynamical matrix is derived based on the modified analytic embedded atom method (MAEAM) and the dynamics theory of surface lattice. The surface phonon dispersions along three major symmetry directions , and X¯M¯ are calculated for the clean Ag (100) surface by using our derived formulas. We then discuss the polarization and localization of surface modes at points X¯ and M¯ by plotting the squared polarization vectors as a function of the layer index. The phonon frequencies of the surface modes calculated by MAEAM are compared with the available experimental and other theoretical data. It is found that the present results are generally in agreement with the referenced experimental or theoretical results, with a maximum deviation of 10.4%. The agreement shows that the modified analytic embedded atom method is a reasonable many-body potential model to quickly describe the surface lattice vibration. It also lays a significant foundation for studying the surface lattice vibration in other metals. Project supported by the National Natural Science Foundation of China (Grant Nos. 61471301 and 61078057), the Scientific Research Program Funded by Shaanxi Provincial Education Department, China (Grant No. 14JK1301), and the Specialized Research Fund for the Doctoral Program of Higher Education, China (Grant No. 20126102110045).

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

Modal element method for potential flow in non-uniform ducts: Combining closed form analysis with CFD

NASA Technical Reports Server (NTRS)

Baumeister, Kenneth J.; Baumeister, Joseph F.

1994-01-01

An analytical procedure is presented, called the modal element method, that combines numerical grid based algorithms with eigenfunction expansions developed by separation of variables. A modal element method is presented for solving potential flow in a channel with two-dimensional cylindrical like obstacles. The infinite computational region is divided into three subdomains; the bounded finite element domain, which is characterized by the cylindrical obstacle and the surrounding unbounded uniform channel entrance and exit domains. The velocity potential is represented approximately in the grid based domain by a finite element solution and is represented analytically by an eigenfunction expansion in the uniform semi-infinite entrance and exit domains. The calculated flow fields are in excellent agreement with exact analytical solutions. By eliminating the grid surrounding the obstacle, the modal element method reduces the numerical grid size, employs a more precise far field boundary condition, as well as giving theoretical insight to the interaction of the obstacle with the mean flow. Although the analysis focuses on a specific geometry, the formulation is general and can be applied to a variety of problems as seen by a comparison to companion theories in aeroacoustics and electromagnetics.

Sample displacement chromatography as a method for purification of proteins and peptides from complex mixtures

PubMed Central

Gajdosik, Martina Srajer; Clifton, James; Josic, Djuro

2012-01-01

Sample displacement chromatography (SDC) in reversed-phase and ion-exchange modes was introduced approximately twenty years ago. This method takes advantage of relative binding affinities of components in a sample mixture. During loading, there is a competition among different sample components for the sorption on the surface of the stationary phase. SDC was first used for the preparative purification of proteins. Later, it was demonstrated that this kind of chromatography can also be performed in ion-exchange, affinity and hydrophobic-interaction mode. It has also been shown that SDC can be performed on monoliths and membrane-based supports in both analytical and preparative scale. Recently, SDC in ion-exchange and hydrophobic interaction mode was also employed successfully for the removal of trace proteins from monoclonal antibody preparations and for the enrichment of low abundance proteins from human plasma. In this review, the principals of SDC are introduced, and the potential for separation of proteins and peptides in micro-analytical, analytical and preparative scale is discussed. PMID:22520159

Analytical modeling of electron energy loss spectroscopy of graphene: Ab initio study versus extended hydrodynamic model.

PubMed

Djordjević, Tijana; Radović, Ivan; Despoja, Vito; Lyon, Keenan; Borka, Duško; Mišković, Zoran L

2018-01-01

We present an analytical modeling of the electron energy loss (EEL) spectroscopy data for free-standing graphene obtained by scanning transmission electron microscope. The probability density for energy loss of fast electrons traversing graphene under normal incidence is evaluated using an optical approximation based on the conductivity of graphene given in the local, i.e., frequency-dependent form derived by both a two-dimensional, two-fluid extended hydrodynamic (eHD) model and an ab initio method. We compare the results for the real and imaginary parts of the optical conductivity in graphene obtained by these two methods. The calculated probability density is directly compared with the EEL spectra from three independent experiments and we find very good agreement, especially in the case of the eHD model. Furthermore, we point out that the subtraction of the zero-loss peak from the experimental EEL spectra has a strong influence on the analytical model for the EEL spectroscopy data. Copyright © 2017 Elsevier B.V. All rights reserved.

The NIST Quantitative Infrared Database

PubMed Central

Chu, P. M.; Guenther, F. R.; Rhoderick, G. C.; Lafferty, W. J.

1999-01-01

With the recent developments in Fourier transform infrared (FTIR) spectrometers it is becoming more feasible to place these instruments in field environments. As a result, there has been enormous increase in the use of FTIR techniques for a variety of qualitative and quantitative chemical measurements. These methods offer the possibility of fully automated real-time quantitation of many analytes; therefore FTIR has great potential as an analytical tool. Recently, the U.S. Environmental Protection Agency (U.S.EPA) has developed protocol methods for emissions monitoring using both extractive and open-path FTIR measurements. Depending upon the analyte, the experimental conditions and the analyte matrix, approximately 100 of the hazardous air pollutants (HAPs) listed in the 1990 U.S.EPA Clean Air Act amendment (CAAA) can be measured. The National Institute of Standards and Technology (NIST) has initiated a program to provide quality-assured infrared absorption coefficient data based on NIST prepared primary gas standards. Currently, absorption coefficient data has been acquired for approximately 20 of the HAPs. For each compound, the absorption coefficient spectrum was calculated using nine transmittance spectra at 0.12 cm−1 resolution and the Beer’s law relationship. The uncertainties in the absorption coefficient data were estimated from the linear regressions of the transmittance data and considerations of other error sources such as the nonlinear detector response. For absorption coefficient values greater than 1 × 10−4 μmol/mol)−1 m−1 the average relative expanded uncertainty is 2.2 %. This quantitative infrared database is currently an ongoing project at NIST. Additional spectra will be added to the database as they are acquired. Our current plans include continued data acquisition of the compounds listed in the CAAA, as well as the compounds that contribute to global warming and ozone depletion.

An extension of the Derrida-Lebowitz-Speer-Spohn equation

NASA Astrophysics Data System (ADS)

Bordenave, Charles; Germain, Pierre; Trogdon, Thomas

2015-12-01

We show how the derivation of the Derrida-Lebowitz-Speer-Spohn equation can be prolonged to obtain a new equation, generalizing the models obtained in the paper by these authors. We then investigate its properties from both an analytical and numerical perspective. Specifically, a numerical method is presented to approximate solutions of the prolonged equation. Using this method, we investigate the relationship between the solutions of the prolonged equation and the Tracy-Widom GOE distribution.

Approximate analytic expression for the Skyrmions crystal

NASA Astrophysics Data System (ADS)

Grandi, Nicolás; Sturla, Mauricio

2018-01-01

We find approximate solutions for the two-dimensional nonlinear Σ-model with Dzyalioshinkii-Moriya term, representing magnetic Skyrmions. They are built in an analytic form, by pasting different approximate solutions found in different regions of space. We verify that our construction reproduces the phenomenology known from numerical solutions and Monte Carlo simulations, giving rise to a Skyrmion lattice at an intermediate range of magnetic field, flanked by spiral and spin-polarized phases for low and high magnetic fields, respectively.

On the singular perturbations for fractional differential equation.

PubMed

Atangana, Abdon

2014-01-01

The goal of this paper is to examine the possible extension of the singular perturbation differential equation to the concept of fractional order derivative. To achieve this, we presented a review of the concept of fractional calculus. We make use of the Laplace transform operator to derive exact solution of singular perturbation fractional linear differential equations. We make use of the methodology of three analytical methods to present exact and approximate solution of the singular perturbation fractional, nonlinear, nonhomogeneous differential equation. These methods are including the regular perturbation method, the new development of the variational iteration method, and the homotopy decomposition method.

Critical evaluation of connectivity-based point of care testing systems of glucose in a hospital environment.

PubMed

Floré, Katelijne M J; Fiers, Tom; Delanghe, Joris R

2008-01-01

In recent years a number of point of care testing (POCT) glucometers were introduced on the market. We investigated the analytical variability (lot-to-lot variation, calibration error, inter-instrument and inter-operator variability) of glucose POCT systems in a university hospital environment and compared these results with the analytical needs required for tight glucose monitoring. The reference hexokinase method was compared to different POCT systems based on glucose oxidase (blood gas instruments) or glucose dehydrogenase (handheld glucometers). Based upon daily internal quality control data, total errors were calculated for the various glucose methods and the analytical variability of the glucometers was estimated. The total error of the glucometers exceeded by far the desirable analytical specifications (based on a biological variability model). Lot-to-lot variation, inter-instrument variation and inter-operator variability contributed approximately equally to total variance. As in a hospital environment, distribution of hematocrit values is broad, converting blood glucose into plasma values using a fixed factor further increases variance. The percentage of outliers exceeded the ISO 15197 criteria in a broad glucose concentration range. Total analytical variation of handheld glucometers is larger than expected. Clinicians should be aware that the variability of glucose measurements obtained by blood gas instruments is lower than results obtained with handheld glucometers on capillary blood.

Application of adjusted data in calculating fission-product decay energies and spectra

NASA Astrophysics Data System (ADS)

George, D. C.; Labauve, R. J.; England, T. R.

1982-06-01

The code ADENA, which approximately calculates fussion-product beta and gamma decay energies and spectra in 19 or fewer energy groups from a mixture of U235 and Pu239 fuels, is described. The calculation uses aggregate, adjusted data derived from a combination of several experiments and summation results based on the ENDF/B-V fission product file. The method used to obtain these adjusted data and the method used by ADENA to calculate fission-product decay energy with an absorption correction are described, and an estimate of the uncertainty of the ADENA results is given. Comparisons of this approximate method are made to experimental measurements, to the ANSI/ANS 5.1-1979 standard, and to other calculational methods. A listing of the complete computer code (ADENA) is contained in an appendix. Included in the listing are data statements containing the adjusted data in the form of parameters to be used in simple analytic functions.

An approximate theoretical method for modeling the static thrust performance of non-axisymmetric two-dimensional convergent-divergent nozzles. M.S. Thesis - George Washington Univ.

NASA Technical Reports Server (NTRS)

Hunter, Craig A.

1995-01-01

An analytical/numerical method has been developed to predict the static thrust performance of non-axisymmetric, two-dimensional convergent-divergent exhaust nozzles. Thermodynamic nozzle performance effects due to over- and underexpansion are modeled using one-dimensional compressible flow theory. Boundary layer development and skin friction losses are calculated using an approximate integral momentum method based on the classic karman-Polhausen solution. Angularity effects are included with these two models in a computational Nozzle Performance Analysis Code, NPAC. In four different case studies, results from NPAC are compared to experimental data obtained from subscale nozzle testing to demonstrate the capabilities and limitations of the NPAC method. In several cases, the NPAC prediction matched experimental gross thrust efficiency data to within 0.1 percent at a design NPR, and to within 0.5 percent at off-design conditions.

Speed-up of the volumetric method of moments for the approximate RCS of large arbitrary-shaped dielectric targets

NASA Astrophysics Data System (ADS)

Moreno, Javier; Somolinos, Álvaro; Romero, Gustavo; González, Iván; Cátedra, Felipe

2017-08-01

A method for the rigorous computation of the electromagnetic scattering of large dielectric volumes is presented. One goal is to simplify the analysis of large dielectric targets with translational symmetries taken advantage of their Toeplitz symmetry. Then, the matrix-fill stage of the Method of Moments is efficiently obtained because the number of coupling terms to compute is reduced. The Multilevel Fast Multipole Method is applied to solve the problem. Structured meshes are obtained efficiently to approximate the dielectric volumes. The regular mesh grid is achieved by using parallelepipeds whose centres have been identified as internal to the target. The ray casting algorithm is used to classify the parallelepiped centres. It may become a bottleneck when too many points are evaluated in volumes defined by parametric surfaces, so a hierarchical algorithm is proposed to minimize the number of evaluations. Measurements and analytical results are included for validation purposes.

Comparison of the iterated equation of motion approach and the density matrix formalism for the quantum Rabi model

NASA Astrophysics Data System (ADS)

Kalthoff, Mona; Keim, Frederik; Krull, Holger; Uhrig, Götz S.

2017-05-01

The density matrix formalism and the equation of motion approach are two semi-analytical methods that can be used to compute the non-equilibrium dynamics of correlated systems. While for a bilinear Hamiltonian both formalisms yield the exact result, for any non-bilinear Hamiltonian a truncation is necessary. Due to the fact that the commonly used truncation schemes differ for these two methods, the accuracy of the obtained results depends significantly on the chosen approach. In this paper, both formalisms are applied to the quantum Rabi model. This allows us to compare the approximate results and the exact dynamics of the system and enables us to discuss the accuracy of the approximations as well as the advantages and the disadvantages of both methods. It is shown to which extent the results fulfill physical requirements for the observables and which properties of the methods lead to unphysical results.

Accelerating electrostatic surface potential calculation with multi-scale approximation on graphics processing units.

PubMed

Anandakrishnan, Ramu; Scogland, Tom R W; Fenley, Andrew T; Gordon, John C; Feng, Wu-chun; Onufriev, Alexey V

2010-06-01

Tools that compute and visualize biomolecular electrostatic surface potential have been used extensively for studying biomolecular function. However, determining the surface potential for large biomolecules on a typical desktop computer can take days or longer using currently available tools and methods. Two commonly used techniques to speed-up these types of electrostatic computations are approximations based on multi-scale coarse-graining and parallelization across multiple processors. This paper demonstrates that for the computation of electrostatic surface potential, these two techniques can be combined to deliver significantly greater speed-up than either one separately, something that is in general not always possible. Specifically, the electrostatic potential computation, using an analytical linearized Poisson-Boltzmann (ALPB) method, is approximated using the hierarchical charge partitioning (HCP) multi-scale method, and parallelized on an ATI Radeon 4870 graphical processing unit (GPU). The implementation delivers a combined 934-fold speed-up for a 476,040 atom viral capsid, compared to an equivalent non-parallel implementation on an Intel E6550 CPU without the approximation. This speed-up is significantly greater than the 42-fold speed-up for the HCP approximation alone or the 182-fold speed-up for the GPU alone. Copyright (c) 2010 Elsevier Inc. All rights reserved.

Decoupling the NLO-coupled QED⊗QCD, DGLAP evolution equations, using Laplace transform method

NASA Astrophysics Data System (ADS)

Mottaghizadeh, Marzieh; Eslami, Parvin; Taghavi-Shahri, Fatemeh

2017-05-01

We analytically solved the QED⊗QCD-coupled DGLAP evolution equations at leading order (LO) quantum electrodynamics (QED) and next-to-leading order (NLO) quantum chromodynamics (QCD) approximations, using the Laplace transform method and then computed the proton structure function in terms of the unpolarized parton distribution functions. Our analytical solutions for parton densities are in good agreement with those from CT14QED (1.2952 < Q2 < 1010) (Ref. 6) global parametrizations and APFEL (A PDF Evolution Library) (2 < Q2 < 108) (Ref. 4). We also compared the proton structure function, F2p(x,Q2), with the experimental data released by the ZEUS and H1 collaborations at HERA. There is a nice agreement between them in the range of low and high x and Q2.

On the line-shape analysis of Compton profiles and its application to neutron scattering

NASA Astrophysics Data System (ADS)

Romanelli, G.; Krzystyniak, M.

2016-05-01

Analytical properties of Compton profiles are used in order to simplify the analysis of neutron Compton scattering experiments. In particular, the possibility to fit the difference of Compton profiles is discussed as a way to greatly decrease the level of complexity of the data treatment, making the analysis easier, faster and more robust. In the context of the novel method proposed, two mathematical models describing the shapes of differenced Compton profiles are discussed: the simple Gaussian approximation for harmonic and isotropic local potential, and an analytical Gauss-Hermite expansion for an anharmonic or anisotropic potential. The method is applied to data collected by VESUVIO spectrometer at ISIS neutron and muon pulsed source (UK) on Copper and Aluminium samples at ambient and low temperatures.

Multi-level methods and approximating distribution functions

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

Wilson, D., E-mail: daniel.wilson@dtc.ox.ac.uk; Baker, R. E.

2016-07-15

Biochemical reaction networks are often modelled using discrete-state, continuous-time Markov chains. System statistics of these Markov chains usually cannot be calculated analytically and therefore estimates must be generated via simulation techniques. There is a well documented class of simulation techniques known as exact stochastic simulation algorithms, an example of which is Gillespie’s direct method. These algorithms often come with high computational costs, therefore approximate stochastic simulation algorithms such as the tau-leap method are used. However, in order to minimise the bias in the estimates generated using them, a relatively small value of tau is needed, rendering the computational costs comparablemore » to Gillespie’s direct method. The multi-level Monte Carlo method (Anderson and Higham, Multiscale Model. Simul. 10:146–179, 2012) provides a reduction in computational costs whilst minimising or even eliminating the bias in the estimates of system statistics. This is achieved by first crudely approximating required statistics with many sample paths of low accuracy. Then correction terms are added until a required level of accuracy is reached. Recent literature has primarily focussed on implementing the multi-level method efficiently to estimate a single system statistic. However, it is clearly also of interest to be able to approximate entire probability distributions of species counts. We present two novel methods that combine known techniques for distribution reconstruction with the multi-level method. We demonstrate the potential of our methods using a number of examples.« less

Eigen solutions and entropic system for Hellmann potential in the presence of the Schrödinger equation

NASA Astrophysics Data System (ADS)

Onate, C. A.; Onyeaju, M. C.; Ikot, A. N.; Ebomwonyi, O.

2017-11-01

By using the supersymmetric approach, we studied the approximate analytic solutions of the three-dimensional Schrödinger equation with the Hellmann potential by applying a suitable approximation scheme to the centrifugal term. The solutions of other useful potentials, such as Coulomb potential and Yukawa potential, are obtained by transformation of variables from the Hellmann potential. Finally, we calculated the Tsallis entropy and Rényi entropy both in position and momentum spaces under the Hellmann potential using integral method. The effects of these entropies on the angular momentum quantum number are investigated in detail.

The D-dimensional non-relativistic particle in the Scarf Trigonometry plus Non-Central Rosen-Morse Potentials

NASA Astrophysics Data System (ADS)

Deta, U. A.; Lestari, N. A.; Yantidewi, M.; Suparmi, A.; Cari, C.

2018-03-01

The D-Dimensional Non-Relativistic Particle Properties in the Scarf Trigonometry plus Non-Central Rosen-Morse Potentials was investigated using an analytical method. The bound state energy is given approximately in the closed form. The approximate wave function for arbitrary l-state in D-dimensions are expressed in the form of generalised Jacobi Polynomials. The energy spectra of the particle are increased when the dimensions are higher. The relationship between the orbital number in each dimension is recursive. The special case in 3 dimensions is given to the ground state.

Reducing microwave absorption with fast frequency modulation.

PubMed

Qin, Juehang; Hubler, A

2017-05-01

We study the response of a two-level quantum system to a chirp signal, using both numerical and analytical methods. The numerical method is based on numerical solutions of the Schrödinger solution of the two-level system, while the analytical method is based on an approximate solution of the same equations. We find that when two-level systems are perturbed by a chirp signal, the peak population of the initially unpopulated state exhibits a high sensitivity to frequency modulation rate. We also find that the aforementioned sensitivity depends on the strength of the forcing, and weaker forcings result in a higher sensitivity, where the frequency modulation rate required to produce the same reduction in peak population would be lower. We discuss potential applications of this result in the field of microwave power transmission, as it shows applying fast frequency modulation to transmitted microwaves used for power transmission could decrease unintended absorption of microwaves by organic tissue.

A Method for Calculating Viscosity and Thermal Conductivity of a Helium-Xenon Gas Mixture

NASA Technical Reports Server (NTRS)

Johnson, Paul K.

2006-01-01

A method for calculating viscosity and thermal conductivity of a helium-xenon (He-Xe) gas mixture was employed, and results were compared to AiResearch (part of Honeywell) analytical data. The method of choice was that presented by Hirschfelder with Singh's third-order correction factor applied to thermal conductivity. Values for viscosity and thermal conductivity were calculated over a temperature range of 400 to 1200 K for He-Xe gas mixture molecular weights of 20.183, 39.94, and 83.8 kg/kmol. First-order values for both transport properties were in good agreement with AiResearch analytical data. Third-order-corrected thermal conductivity values were all greater than AiResearch data, but were considered to be a better approximation of thermal conductivity because higher-order effects of mass and temperature were taken into consideration. Viscosity, conductivity, and Prandtl number were then compared to experimental data presented by Taylor.

Analytic energy gradients for the coupled-cluster singles and doubles method with the density-fitting approximation

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

Bozkaya, Uğur, E-mail: ugur.bozkaya@hacettepe.edu.tr; Department of Chemistry, Atatürk University, Erzurum 25240; Sherrill, C. David

2016-05-07

An efficient implementation is presented for analytic gradients of the coupled-cluster singles and doubles (CCSD) method with the density-fitting approximation, denoted DF-CCSD. Frozen core terms are also included. When applied to a set of alkanes, the DF-CCSD analytic gradients are significantly accelerated compared to conventional CCSD for larger molecules. The efficiency of our DF-CCSD algorithm arises from the acceleration of several different terms, which are designated as the “gradient terms”: computation of particle density matrices (PDMs), generalized Fock-matrix (GFM), solution of the Z-vector equation, formation of the relaxed PDMs and GFM, back-transformation of PDMs and GFM to the atomic orbitalmore » (AO) basis, and evaluation of gradients in the AO basis. For the largest member of the alkane set (C{sub 10}H{sub 22}), the computational times for the gradient terms (with the cc-pVTZ basis set) are 2582.6 (CCSD) and 310.7 (DF-CCSD) min, respectively, a speed up of more than 8-folds. For gradient related terms, the DF approach avoids the usage of four-index electron repulsion integrals. Based on our previous study [U. Bozkaya, J. Chem. Phys. 141, 124108 (2014)], our formalism completely avoids construction or storage of the 4-index two-particle density matrix (TPDM), using instead 2- and 3-index TPDMs. The DF approach introduces negligible errors for equilibrium bond lengths and harmonic vibrational frequencies.« less

Identification of multiple leaks in pipeline: Linearized model, maximum likelihood, and super-resolution localization

NASA Astrophysics Data System (ADS)

Wang, Xun; Ghidaoui, Mohamed S.

2018-07-01

This paper considers the problem of identifying multiple leaks in a water-filled pipeline based on inverse transient wave theory. The analytical solution to this problem involves nonlinear interaction terms between the various leaks. This paper shows analytically and numerically that these nonlinear terms are of the order of the leak sizes to the power two and; thus, negligible. As a result of this simplification, a maximum likelihood (ML) scheme that identifies leak locations and leak sizes separately is formulated and tested. It is found that the ML estimation scheme is highly efficient and robust with respect to noise. In addition, the ML method is a super-resolution leak localization scheme because its resolvable leak distance (approximately 0.15λmin , where λmin is the minimum wavelength) is below the Nyquist-Shannon sampling theorem limit (0.5λmin). Moreover, the Cramér-Rao lower bound (CRLB) is derived and used to show the efficiency of the ML scheme estimates. The variance of the ML estimator approximates the CRLB proving that the ML scheme belongs to class of best unbiased estimator of leak localization methods.

On the exactness of effective Floquet Hamiltonians employed in solid-state NMR spectroscopy

NASA Astrophysics Data System (ADS)

Garg, Rajat; Ramachandran, Ramesh

2017-05-01

Development of theoretical models based on analytic theory has remained an active pursuit in molecular spectroscopy for its utility both in the design of experiments as well as in the interpretation of spectroscopic data. In particular, the role of "Effective Hamiltonians" in the evolution of theoretical frameworks is well known across all forms of spectroscopy. Nevertheless, a constant revalidation of the approximations employed in the theoretical frameworks is necessitated by the constant improvements on the experimental front in addition to the complexity posed by the systems under study. Here in this article, we confine our discussion to the derivation of effective Floquet Hamiltonians based on the contact transformation procedure. While the importance of the effective Floquet Hamiltonians in the qualitative description of NMR experiments has been realized in simpler cases, its extension in quantifying spectral data deserves a cautious approach. With this objective, the validity of the approximations employed in the derivation of the effective Floquet Hamiltonians is re-examined through a comparison with exact numerical methods under differing experimental conditions. The limitations arising from the existing analytic methods are outlined along with remedial measures for improving the accuracy of the derived effective Floquet Hamiltonians.

Exact analytical formulae for linearly distributed vortex and source sheets in uence computation in 2D vortex methods

NASA Astrophysics Data System (ADS)

Kuzmina, K. S.; Marchevsky, I. K.; Ryatina, E. P.

2017-11-01

We consider the methodology of numerical schemes development for two-dimensional vortex method. We describe two different approaches to deriving integral equation for unknown vortex sheet intensity. We simulate the velocity of the surface line of an airfoil as the influence of attached vortex and source sheets. We consider a polygonal approximation of the airfoil and assume intensity distributions of free and attached vortex sheets and attached source sheet to be approximated with piecewise constant or piecewise linear (continuous or discontinuous) functions. We describe several specific numerical schemes that provide different accuracy and have a different computational cost. The study shows that a Galerkin-type approach to solving boundary integral equation requires computing several integrals and double integrals over the panels. We obtain exact analytical formulae for all the necessary integrals, which makes it possible to raise significantly the accuracy of vortex sheet intensity computation and improve the quality of velocity and vorticity field representation, especially in proximity to the surface line of the airfoil. All the formulae are written down in the invariant form and depend only on the geometric relationship between the positions of the beginnings and ends of the panels.

Analytical and numerical construction of vertical periodic orbits about triangular libration points based on polynomial expansion relations among directions

NASA Astrophysics Data System (ADS)

Qian, Ying-Jing; Yang, Xiao-Dong; Zhai, Guan-Qiao; Zhang, Wei

2017-08-01

Innovated by the nonlinear modes concept in the vibrational dynamics, the vertical periodic orbits around the triangular libration points are revisited for the Circular Restricted Three-body Problem. The ζ -component motion is treated as the dominant motion and the ξ and η -component motions are treated as the slave motions. The slave motions are in nature related to the dominant motion through the approximate nonlinear polynomial expansions with respect to the ζ -position and ζ -velocity during the one of the periodic orbital motions. By employing the relations among the three directions, the three-dimensional system can be transferred into one-dimensional problem. Then the approximate three-dimensional vertical periodic solution can be analytically obtained by solving the dominant motion only on ζ -direction. To demonstrate the effectiveness of the proposed method, an accuracy study was carried out to validate the polynomial expansion (PE) method. As one of the applications, the invariant nonlinear relations in polynomial expansion form are used as constraints to obtain numerical solutions by differential correction. The nonlinear relations among the directions provide an alternative point of view to explore the overall dynamics of periodic orbits around libration points with general rules.

Properties of atomic pairs produced in the collision of Bose-Einstein condensates

NASA Astrophysics Data System (ADS)

Ziń, Paweł; Wasak, Tomasz

2018-04-01

During a collision of Bose-Einstein condensates correlated pairs of atoms are emitted. The scattered massive particles, in analogy to photon pairs in quantum optics, might be used in the violation of Bell's inequalities, demonstration of Einstein-Podolsky-Rosen correlations, or sub-shot-noise atomic interferometry. Usually, a theoretical description of the collision relies either on stochastic numerical methods or on analytical treatments involving various approximations. Here, we investigate elastic scattering of atoms from colliding elongated Bose-Einstein condensates within the Bogoliubov method, carefully controlling performed approximations at every stage of the analysis. We derive expressions for the one- and two-particle correlation functions. The obtained formulas, which relate the correlation functions to the condensate wave function, are convenient for numerical calculations. We employ the variational approach for condensate wave functions to obtain analytical expressions for the correlation functions, whose properties we analyze in detail. We also present a useful semiclassical model of the process and compare its results with the quantum one. The results are relevant for recent experiments with excited helium atoms, as well as for planned experiments aimed at investigating the nonclassicality of the system.

An Accurate Analytic Approximation for Light Scattering by Non-absorbing Spherical Aerosol Particles

NASA Astrophysics Data System (ADS)

Lewis, E. R.

2017-12-01

The scattering of light by particles in the atmosphere is a ubiquitous and important phenomenon, with applications to numerous fields of science and technology. The problem of scattering of electromagnetic radiation by a uniform spherical particle can be solved by the method of Mie and Debye as a series of terms depending on the size parameter, x=2πr/λ, and the complex index of refraction, m. However, this solution does not provide insight into the dependence of the scattering on the radius of the particle, the wavelength, or the index of refraction, or how the scattering varies with relative humidity. Van de Hulst demonstrated that the scattering efficiency (the scattering cross section divided by the geometric cross section) of a non-absorbing sphere, over a wide range of particle sizes of atmospheric importance, depends not on x and m separately, but on the quantity 2x(m-1); this is the basis for the anomalous diffraction approximation. Here an analytic approximation for the scattering efficiency of a non-absorbing spherical particle is presented in terms of this new quantity that is accurate over a wide range of particle sizes of atmospheric importance and which readily displays the dependences of the scattering efficiency on particle radius, index of refraction, and wavelength. For an aerosol for which the particle size distribution is parameterized as a gamma function, this approximation also yields analytical results for the scattering coefficient and for the Ångström exponent, with the dependences of scattering properties on wavelength and index of refraction clearly displayed. This approximation provides insight into the dependence of light scattering properties on factors such as relative humidity, readily enables conversion of scattering from one index of refraction to another, and demonstrates the conditions under which the aerosol index (the product of the aerosol optical depth and the Ångström exponent) is a useful proxy for the number of cloud condensation nuclei.

Size-dependent resonance frequencies of cantilevered and bridged nanosensors

NASA Astrophysics Data System (ADS)

Shi, W.; Zou, J.; Lee, K. Y.; Li, X. F.

2018-03-01

This paper studies transverse vibration of nanoscale cantilevered and bridged sensors carrying a nanoparticle. The nanoscale sensors are modelled as Euler-Bernoulli beams with surface effect and nanoparticle as a concentrated mass. Frequency equations of cantilevered and bridged beam-mass system are derived and exact resonance frequencies are calculated. An alternative Fredholm integral equation method is used to obtain an approximate explicit expression for the fundamental frequency for both cases. A comparison between the approximate and analytical results is made and the approximation accuracy is satisfactory. The influences of the residual surface stress, surface elasticity, and attached mass on the resonance frequencies and mode shapes are discussed. These results are useful to illustrate the surface phenomena and are helpful to design micro-/nano-mechanical sensors.

On-line focusing of flavin derivatives using Dynamic pH junction-sweeping capillary electrophoresis with laser-induced fluorescence detection.

PubMed

Britz-McKibbin, Philip; Otsuka, Koji; Terabe, Shigeru

2002-08-01

Simple yet effective methods to enhance concentration sensitivity is needed for capillary electrophoresis (CE) to become a practical method to analyze trace levels of analytes in real samples. In this report, the development of a novel on-line preconcentration technique combining dynamic pH junction and sweeping modes of focusing is applied to the sensitive and selective analysis of three flavin derivatives: riboflavin, flavin mononucleotide (FMN) and flavin adenine dinucleotide (FAD). Picomolar (pM) detectability of flavins by CE with laser-induced fluorescence (LIF) detection is demonstrated through effective focusing of large sample volumes (up to 22% capillary length) using a dual pH junction-sweeping focusing mode. This results in greater than a 1,200-fold improvement in sensitivity relative to conventional injection methods, giving a limit of detection (S/N = 3) of approximately 4.0 pM for FAD and FMN. Flavin focusing is examined in terms of analyte mobility dependence on buffer pH, borate complexation and SDS interaction. Dynamic pH junction-sweeping extends on-line focusing to both neutral (hydrophobic) and weakly acidic (hydrophilic) species and is considered useful in cases when either conventional sweeping or dynamic pH junction techniques used alone are less effective for certain classes of analytes. Enhanced focusing performance by this hyphenated method was demonstrated by greater than a 4-fold reduction in flavin bandwidth, as compared to either sweeping or dynamic pH junction, reflected by analyte detector bandwidths <0.20 cm. Novel on-line focusing strategies are required to improve sensitivity in CE, which may be applied toward more effective biochemical analysis methods for diverse types of analytes.

A Machine-Learning-Driven Sky Model.

PubMed

Satylmys, Pynar; Bashford-Rogers, Thomas; Chalmers, Alan; Debattista, Kurt

2017-01-01

Sky illumination is responsible for much of the lighting in a virtual environment. A machine-learning-based approach can compactly represent sky illumination from both existing analytic sky models and from captured environment maps. The proposed approach can approximate the captured lighting at a significantly reduced memory cost and enable smooth transitions of sky lighting to be created from a small set of environment maps captured at discrete times of day. The author's results demonstrate accuracy close to the ground truth for both analytical and capture-based methods. The approach has a low runtime overhead, so it can be used as a generic approach for both offline and real-time applications.

Analytical expressions for the correlation function of a hard sphere dimer fluid

NASA Astrophysics Data System (ADS)

Kim, Soonho; Chang, Jaeeon; Kim, Hwayong

A closed form expression is given for the correlation function of a hard sphere dimer fluid. A set of integral equations is obtained from Wertheim's multidensity Ornstein-Zernike integral equation theory with Percus-Yevick approximation. Applying the Laplace transformation method to the integral equations and then solving the resulting equations algebraically, the Laplace transforms of the individual correlation functions are obtained. By the inverse Laplace transformation, the radial distribution function (RDF) is obtained in closed form out to 3D (D is the segment diameter). The analytical expression for the RDF of the hard dimer should be useful in developing the perturbation theory of dimer fluids.

Analytical expression for the correlation function of a hard sphere chain fluid

NASA Astrophysics Data System (ADS)

Chang, Jaeeon; Kim, Hwayong

A closed form expression is given for the correlation function of flexible hard sphere chain fluid. A set of integral equations obtained from Wertheim's multidensity Ornstein-Zernike integral equation theory with the polymer Percus-Yevick ideal chain approximation is considered. Applying the Laplace transformation method to the integral equations and then solving the resulting equations algebraically, the Laplace transforms of individual correlation functions are obtained. By inverse Laplace transformation the inter- and intramolecular radial distribution functions (RDFs) are obtained in closed forms up to 3D(D is segment diameter). These analytical expressions for the RDFs would be useful in developing the perturbation theory of chain fluids.

A Galerkin method for linear PDE systems in circular geometries with structural acoustic applications

NASA Technical Reports Server (NTRS)

Smith, Ralph C.

1994-01-01

A Galerkin method for systems of PDE's in circular geometries is presented with motivating problems being drawn from structural, acoustic, and structural acoustic applications. Depending upon the application under consideration, piecewise splines or Legendre polynomials are used when approximating the system dynamics with modifications included to incorporate the analytic solution decay near the coordinate singularity. This provides an efficient method which retains its accuracy throughout the circular domain without degradation at singularity. Because the problems under consideration are linear or weakly nonlinear with constant or piecewise constant coefficients, transform methods for the problems are not investigated. While the specific method is developed for the two dimensional wave equations on a circular domain and the equation of transverse motion for a thin circular plate, examples demonstrating the extension of the techniques to a fully coupled structural acoustic system are used to illustrate the flexibility of the method when approximating the dynamics of more complex systems.

Methods for describing the electromagnetic properties of silver and gold nanoparticles.

PubMed

Zhao, Jing; Pinchuk, Anatoliy O; McMahon, Jeffrey M; Li, Shuzhou; Ausman, Logan K; Atkinson, Ariel L; Schatz, George C

2008-12-01

This Account provides an overview of the methods that are currently being used to study the electromagnetics of silver and gold nanoparticles, with an emphasis on the determination of extinction and surface-enhanced Raman scattering (SERS) spectra. These methods have proven to be immensely useful in recent years for interpreting a wide range of nanoscience experiments and providing the capability to describe optical properties of particles up to several hundred nanometers in dimension, including arbitrary particle structures and complex dielectric environments (adsorbed layers of molecules, nearby metal films, and other particles). While some of the methods date back to Mie's celebrated work a century ago, others are still at the forefront of algorithm development in computational electromagnetics. This Account gives a qualitative description of the physical and mathematical basis behind the most commonly used methods, including both analytical and numerical methods, as well as representative results of applications that are relevant to current experiments. The analytical methods that we discuss are either derived from Mie theory for spheres or from the quasistatic (Gans) model as applied to spheres and spheroids. In this discussion, we describe the use of Mie theory to determine electromagnetic contributions to SERS enhancements that include for retarded dipole emission effects, and the use of the quasistatic approximation for spheroidal particles interacting with dye adsorbate layers. The numerical methods include the discrete dipole approximation (DDA), the finite difference time domain (FDTD) method, and the finite element method (FEM) based on Whitney forms. We discuss applications such as using DDA to describe the interaction of two gold disks to define electromagnetic hot spots, FDTD for light interacting with metal wires that go from particle-like plasmonic response to the film-like transmission as wire dimension is varied, and FEM studies of electromagnetic fields near cubic particles.

Precise analytic approximations for the Bessel function J1 (x)

NASA Astrophysics Data System (ADS)

Maass, Fernando; Martin, Pablo

2018-03-01

Precise and straightforward analytic approximations for the Bessel function J1 (x) have been found. Power series and asymptotic expansions have been used to determine the parameters of the approximation, which is as a bridge between both expansions, and it is a combination of rational and trigonometric functions multiplied with fractional powers of x. Here, several improvements with respect to the so called Multipoint Quasirational Approximation technique have been performed. Two procedures have been used to determine the parameters of the approximations. The maximum absolute errors are in both cases smaller than 0.01. The zeros of the approximation are also very precise with less than 0.04 per cent for the first one. A second approximation has been also determined using two more parameters, and in this way the accuracy has been increased to less than 0.001.

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.

A rapid method for the determination of some antihypertensive and antipyretic drugs by thermometric titrimetry.

PubMed

Abbasi, U M; Chand, F; Bhanger, M I; Memon, S A

1986-02-01

A simple and rapid method is described for the direct thermometric determination of milligram amounts of methyl dopa, propranolol hydrochloride, 1-phenyl-3-methylpyrazolone (MPP) and 2,3-dimethyl-1-phenylpyrazol-5-one (phenazone) in the presence of excipients. The compounds are reacted with N'-bromosuccinimide and the heat of reaction is used to determine the end-point of the titration. The time required is approximately 2 min, and the accuracy is analytically acceptable.

On Chaotic Behavior of Temperature Distribution in a Heat Exchanger

NASA Astrophysics Data System (ADS)

Bagyalakshmi, Morachan; Gangadharan, Saisundarakrishnan; Ganesh, Madhu

The objective of this paper is to introduce the notion of fractional derivatives in the energy equations and to study the chaotic nature of the temperature distribution in a heat exchanger with variation of temperature dependent transport properties. The governing fractional partial differential equations are transformed to a set of recurrence relations using fractional differential transform method and solved using inverse transform. The approximate analytical solution obtained by the proposed method has good agreement with the existing results.

Designing stellarator coils by a modified Newton method using FOCUS

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

Zhu, Caoxiang; Hudson, Stuart R.; Song, Yuntao

To find the optimal coils for stellarators, nonlinear optimization algorithms are applied in existing coil design codes. However, none of these codes have used the information from the second-order derivatives. In this paper, we present a modified Newton method in the recently developed code FOCUS. The Hessian matrix is calculated with analytically derived equations. Its inverse is approximated by a modified Cholesky factorization and applied in the iterative scheme of a classical Newton method. Using this method, FOCUS is able to recover the W7-X modular coils starting from a simple initial guess. Results demonstrate significant advantages.

Designing stellarator coils by a modified Newton method using FOCUS

NASA Astrophysics Data System (ADS)

Zhu, Caoxiang; Hudson, Stuart R.; Song, Yuntao; Wan, Yuanxi

2018-06-01

To find the optimal coils for stellarators, nonlinear optimization algorithms are applied in existing coil design codes. However, none of these codes have used the information from the second-order derivatives. In this paper, we present a modified Newton method in the recently developed code FOCUS. The Hessian matrix is calculated with analytically derived equations. Its inverse is approximated by a modified Cholesky factorization and applied in the iterative scheme of a classical Newton method. Using this method, FOCUS is able to recover the W7-X modular coils starting from a simple initial guess. Results demonstrate significant advantages.

Designing stellarator coils by a modified Newton method using FOCUS

DOE PAGES

Zhu, Caoxiang; Hudson, Stuart R.; Song, Yuntao; ...

2018-03-22

To find the optimal coils for stellarators, nonlinear optimization algorithms are applied in existing coil design codes. However, none of these codes have used the information from the second-order derivatives. In this paper, we present a modified Newton method in the recently developed code FOCUS. The Hessian matrix is calculated with analytically derived equations. Its inverse is approximated by a modified Cholesky factorization and applied in the iterative scheme of a classical Newton method. Using this method, FOCUS is able to recover the W7-X modular coils starting from a simple initial guess. Results demonstrate significant advantages.

Targeted Analyte Detection by Standard Addition Improves Detection Limits in MALDI Mass Spectrometry

PubMed Central

Eshghi, Shadi Toghi; Li, Xingde; Zhang, Hui

2014-01-01

Matrix-assisted laser desorption/ionization has proven an effective tool for fast and accurate determination of many molecules. However, the detector sensitivity and chemical noise compromise the detection of many invaluable low-abundance molecules from biological and clinical samples. To challenge this limitation, we developed a targeted analyte detection (TAD) technique. In TAD, the target analyte is selectively elevated by spiking a known amount of that analyte into the sample, thereby raising its concentration above the noise level, where we take advantage of the improved sensitivity to detect the presence of the endogenous analyte in the sample. We assessed TAD on three peptides in simple and complex background solutions with various exogenous analyte concentrations in two MALDI matrices. TAD successfully improved the limit of detection (LOD) of target analytes when the target peptides were added to the sample in a concentration close to optimum concentration. The optimum exogenous concentration was estimated through a quantitative method to be approximately equal to the original LOD for each target. Also, we showed that TAD could achieve LOD improvements on an average of 3-fold in a simple and 2-fold in a complex sample. TAD provides a straightforward assay to improve the LOD of generic target analytes without the need for costly hardware modifications. PMID:22877355

Targeted analyte detection by standard addition improves detection limits in matrix-assisted laser desorption/ionization mass spectrometry.

PubMed

Toghi Eshghi, Shadi; Li, Xingde; Zhang, Hui

2012-09-18

Matrix-assisted laser desorption/ionization (MALDI) has proven an effective tool for fast and accurate determination of many molecules. However, the detector sensitivity and chemical noise compromise the detection of many invaluable low-abundance molecules from biological and clinical samples. To challenge this limitation, we developed a targeted analyte detection (TAD) technique. In TAD, the target analyte is selectively elevated by spiking a known amount of that analyte into the sample, thereby raising its concentration above the noise level, where we take advantage of the improved sensitivity to detect the presence of the endogenous analyte in the sample. We assessed TAD on three peptides in simple and complex background solutions with various exogenous analyte concentrations in two MALDI matrices. TAD successfully improved the limit of detection (LOD) of target analytes when the target peptides were added to the sample in a concentration close to optimum concentration. The optimum exogenous concentration was estimated through a quantitative method to be approximately equal to the original LOD for each target. Also, we showed that TAD could achieve LOD improvements on an average of 3-fold in a simple and 2-fold in a complex sample. TAD provides a straightforward assay to improve the LOD of generic target analytes without the need for costly hardware modifications.

Approximate message passing for nonconvex sparse regularization with stability and asymptotic analysis

NASA Astrophysics Data System (ADS)

Sakata, Ayaka; Xu, Yingying

2018-03-01

We analyse a linear regression problem with nonconvex regularization called smoothly clipped absolute deviation (SCAD) under an overcomplete Gaussian basis for Gaussian random data. We propose an approximate message passing (AMP) algorithm considering nonconvex regularization, namely SCAD-AMP, and analytically show that the stability condition corresponds to the de Almeida-Thouless condition in spin glass literature. Through asymptotic analysis, we show the correspondence between the density evolution of SCAD-AMP and the replica symmetric (RS) solution. Numerical experiments confirm that for a sufficiently large system size, SCAD-AMP achieves the optimal performance predicted by the replica method. Through replica analysis, a phase transition between replica symmetric and replica symmetry breaking (RSB) region is found in the parameter space of SCAD. The appearance of the RS region for a nonconvex penalty is a significant advantage that indicates the region of smooth landscape of the optimization problem. Furthermore, we analytically show that the statistical representation performance of the SCAD penalty is better than that of \

Substructure Versus Property-Level Dispersed Modes Calculation

NASA Technical Reports Server (NTRS)

Stewart, Eric C.; Peck, Jeff A.; Bush, T. Jason; Fulcher, Clay W.

2016-01-01

This paper calculates the effect of perturbed finite element mass and stiffness values on the eigenvectors and eigenvalues of the finite element model. The structure is perturbed in two ways: at the "subelement" level and at the material property level. In the subelement eigenvalue uncertainty analysis the mass and stiffness of each subelement is perturbed by a factor before being assembled into the global matrices. In the property-level eigenvalue uncertainty analysis all material density and stiffness parameters of the structure are perturbed modified prior to the eigenvalue analysis. The eigenvalue and eigenvector dispersions of each analysis (subelement and property-level) are also calculated using an analytical sensitivity approximation. Two structural models are used to compare these methods: a cantilevered beam model, and a model of the Space Launch System. For each structural model it is shown how well the analytical sensitivity modes approximate the exact modes when the uncertainties are applied at the subelement level and at the property level.

Peculiarities of the momentum distribution functions of strongly correlated charged fermions

NASA Astrophysics Data System (ADS)

Larkin, A. S.; Filinov, V. S.; Fortov, V. E.

2018-01-01

New numerical version of the Wigner approach to quantum thermodynamics of strongly coupled systems of particles has been developed for extreme conditions, when analytical approximations based on different kinds of perturbation theories cannot be applied. An explicit analytical expression of the Wigner function has been obtained in linear and harmonic approximations. Fermi statistical effects are accounted for by effective pair pseudopotential depending on coordinates, momenta and degeneracy parameter of particles and taking into account Pauli blocking of fermions. A new quantum Monte-Carlo method for calculations of average values of arbitrary quantum operators has been developed. Calculations of the momentum distribution functions and the pair correlation functions of degenerate ideal Fermi gas have been carried out for testing the developed approach. Comparison of the obtained momentum distribution functions of strongly correlated Coulomb systems with the Maxwell-Boltzmann and the Fermi distributions shows the significant influence of interparticle interaction both at small momenta and in high energy quantum ‘tails’.

Protecting a quantum state from environmental noise by an incompatible finite-time measurement

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

Brasil, Carlos Alexandre; Castro, L. A. de; Napolitano, R. d. J.

We show that measurements of finite duration performed on an open two-state system can protect the initial state from a phase-noisy environment, provided the measured observable does not commute with the perturbing interaction. When the measured observable commutes with the environmental interaction, the finite-duration measurement accelerates the rate of decoherence induced by the phase noise. For the description of the measurement of an observable that is incompatible with the interaction between system and environment, we have found an approximate analytical expression, valid at zero temperature and weak coupling with the measuring device. We have tested the validity of the analyticalmore » predictions against an exact numerical approach, based on the superoperator-splitting method, that confirms the protection of the initial state of the system. When the coupling between the system and the measuring apparatus increases beyond the range of validity of the analytical approximation, the initial state is still protected by the finite-time measurement, according with the exact numerical calculations.« less

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

Dechant, Lawrence J.

Wave packet analysis provides a connection between linear small disturbance theory and subsequent nonlinear turbulent spot flow behavior. The traditional association between linear stability analysis and nonlinear wave form is developed via the method of stationary phase whereby asymptotic (simplified) mean flow solutions are used to estimate dispersion behavior and stationary phase approximation are used to invert the associated Fourier transform. The resulting process typically requires nonlinear algebraic equations inversions that can be best performed numerically, which partially mitigates the value of the approximation as compared to a more complete, e.g. DNS or linear/nonlinear adjoint methods. To obtain a simpler,more » closed-form analytical result, the complete packet solution is modeled via approximate amplitude (linear convected kinematic wave initial value problem) and local sinusoidal (wave equation) expressions. Significantly, the initial value for the kinematic wave transport expression follows from a separable variable coefficient approximation to the linearized pressure fluctuation Poisson expression. The resulting amplitude solution, while approximate in nature, nonetheless, appears to mimic many of the global features, e.g. transitional flow intermittency and pressure fluctuation magnitude behavior. A low wave number wave packet models also recover meaningful auto-correlation and low frequency spectral behaviors.« less

The singular values of the imbedding operators of some classes of analytic functions of several variables

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

Parfenov, O.G.

1994-12-25

We discuss three results. The first exhibits the order of decrease of the s-values as a function of the CR-dimension of a compact set on which we approximate the class of analytic functions being studied. The second is an asymptotic formula for the case when the domain of analyticity and the compact set are Reinhart domains. The third is the computation of the s-values of a special operator that is of interest for approximation theory on one-dimensional manifolds.

Anharmonic phonon-phonon scattering modeling of three-dimensional atomistic transport: An efficient quantum treatment

NASA Astrophysics Data System (ADS)

Lee, Y.; Bescond, M.; Logoteta, D.; Cavassilas, N.; Lannoo, M.; Luisier, M.

2018-05-01

We propose an efficient method to quantum mechanically treat anharmonic interactions in the atomistic nonequilibrium Green's function simulation of phonon transport. We demonstrate that the so-called lowest-order approximation, implemented through a rescaling technique and analytically continued by means of the Padé approximants, can be used to accurately model third-order anharmonic effects. Although the paper focuses on a specific self-energy, the method is applicable to a very wide class of physical interactions. We apply this approach to the simulation of anharmonic phonon transport in realistic Si and Ge nanowires with uniform or discontinuous cross sections. The effect of increasing the temperature above 300 K is also investigated. In all the considered cases, we are able to obtain a good agreement with the routinely adopted self-consistent Born approximation, at a remarkably lower computational cost. In the more complicated case of high temperatures (≫300 K), we find that the first-order Richardson extrapolation applied to the sequence of the Padé approximants N -1 /N results in a significant acceleration of the convergence.

Approximate analytical solutions to the double-stance dynamics of the lossy spring-loaded inverted pendulum.

PubMed

Shahbazi, Mohammad; Saranlı, Uluç; Babuška, Robert; Lopes, Gabriel A D

2016-12-05

This paper introduces approximate time-domain solutions to the otherwise non-integrable double-stance dynamics of the 'bipedal' spring-loaded inverted pendulum (B-SLIP) in the presence of non-negligible damping. We first introduce an auxiliary system whose behavior under certain conditions is approximately equivalent to the B-SLIP in double-stance. Then, we derive approximate solutions to the dynamics of the new system following two different methods: (i) updated-momentum approach that can deal with both the lossy and lossless B-SLIP models, and (ii) perturbation-based approach following which we only derive a solution to the lossless case. The prediction performance of each method is characterized via a comprehensive numerical analysis. The derived representations are computationally very efficient compared to numerical integrations, and, hence, are suitable for online planning, increasing the autonomy of walking robots. Two application examples of walking gait control are presented. The proposed solutions can serve as instrumental tools in various fields such as control in legged robotics and human motion understanding in biomechanics.

Theory and simulation of the time-dependent rate coefficients of diffusion-influenced reactions.

PubMed Central

Zhou, H X; Szabo, A

1996-01-01

A general formalism is developed for calculating the time-dependent rate coefficient k(t) of an irreversible diffusion-influenced reaction. This formalism allows one to treat most factors that affect k(t), including rotational Brownian motion and conformational gating of reactant molecules and orientation constraint for product formation. At long times k(t) is shown to have the asymptotic expansion k(infinity)[1 + k(infinity) (pie Dt)-1/2 /4 pie D + ...], where D is the relative translational diffusion constant. An approximate analytical method for calculating k(t) is presented. This is based on the approximation that the probability density of the reactant pair in the reactive region keeps the equilibrium distribution but with a decreasing amplitude. The rate coefficient then is determined by the Green function in the absence of chemical reaction. Within the framework of this approximation, two general relations are obtained. The first relation allows the rate coefficient for an arbitrary amplitude of the reactivity to be found if the rate coefficient for one amplitude of the reactivity is known. The second relation allows the rate coefficient in the presence of conformational gating to be found from that in the absence of conformational gating. The ratio k(t)/k(0) is shown to be the survival probability of the reactant pair at time t starting from an initial distribution that is localized in the reactive region. This relation forms the basis of the calculation of k(t) through Brownian dynamics simulations. Two simulation procedures involving the propagation of nonreactive trajectories initiated only from the reactive region are described and illustrated on a model system. Both analytical and simulation results demonstrate the accuracy of the equilibrium-distribution approximation method. PMID:8913584

Modeling systematic errors: polychromatic sources of Beer-Lambert deviations in HPLC/UV and nonchromatographic spectrophotometric assays.

PubMed

Galli, C

2001-07-01

It is well established that the use of polychromatic radiation in spectrophotometric assays leads to excursions from the Beer-Lambert limit. This Note models the resulting systematic error as a function of assay spectral width, slope of molecular extinction coefficient, and analyte concentration. The theoretical calculations are compared with recent experimental results; a parameter is introduced which can be used to estimate the magnitude of the systematic error in both chromatographic and nonchromatographic spectrophotometric assays. It is important to realize that the polychromatic radiation employed in common laboratory equipment can yield assay errors up to approximately 4%, even at absorption levels generally considered 'safe' (i.e. absorption <1). Thus careful consideration of instrumental spectral width, analyte concentration, and slope of molecular extinction coefficient is required to ensure robust analytical methods.

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.

Closed-form analytical solutions of high-temperature heat pipe startup and frozen startup limitation

NASA Technical Reports Server (NTRS)

Cao, Y.; Faghri, A.

1992-01-01

Previous numerical and experimental studies indicate that the high-temperature heat pipe startup process is characterized by a moving hot zone with relatively sharp fronts. Based on the above observation, a flat-front model for an approximate analytical solution is proposed. A closed-form solution related to the temperature distribution in the hot zone and the hot zone length as a function of time are obtained. The analytical results agree well with the corresponding experimental data, and provide a quick prediction method for the heat pipe startup performance. Finally, a heat pipe limitation related to the frozen startup process is identified, and an explicit criterion for the high-temperature heat pipe startup is derived. The frozen startup limit identified in this paper provides a fundamental guidance for high-temperature heat pipe design.

Analytic approximations to the modon dispersion relation. [in oceanography

NASA Technical Reports Server (NTRS)

Boyd, J. P.

1981-01-01

Three explicit analytic approximations are given to the modon dispersion relation developed by Flierl et al. (1980) to describe Gulf Stream rings and related phenomena in the oceans and atmosphere. The solutions are in the form of k(q), and are developed in the form of a power series in q for small q, an inverse power series in 1/q for large q, and a two-point Pade approximant. The low order Pade approximant is shown to yield a solution for the dispersion relation with a maximum relative error for the lowest branch of the function equal to one in 700 in the q interval zero to infinity.

Trace level detection of compounds related to the chemical weapons convention by 1H-detected 13C NMR spectroscopy executed with a sensitivity-enhanced, cryogenic probehead.

PubMed

Cullinan, David B; Hondrogiannis, George; Henderson, Terry J

2008-04-15

Two-dimensional 1H-13C HSQC (heteronuclear single quantum correlation) and fast-HMQC (heteronuclear multiple quantum correlation) pulse sequences were implemented using a sensitivity-enhanced, cryogenic probehead for detecting compounds relevant to the Chemical Weapons Convention present in complex mixtures. The resulting methods demonstrated exceptional sensitivity for detecting the analytes at trace level concentrations. 1H-13C correlations of target analytes at < or = 25 microg/mL were easily detected in a sample where the 1H solvent signal was approximately 58,000-fold more intense than the analyte 1H signals. The problem of overlapping signals typically observed in conventional 1H spectroscopy was essentially eliminated, while 1H and 13C chemical shift information could be derived quickly and simultaneously from the resulting spectra. The fast-HMQC pulse sequences generated magnitude mode spectra suitable for detailed analysis in approximately 4.5 h and can be used in experiments to efficiently screen a large number of samples. The HSQC pulse sequences, on the other hand, required roughly twice the data acquisition time to produce suitable spectra. These spectra, however, were phase-sensitive, contained considerably more resolution in both dimensions, and proved to be superior for detecting analyte 1H-13C correlations. Furthermore, a HSQC spectrum collected with a multiplicity-edited pulse sequence provided additional structural information valuable for identifying target analytes. The HSQC pulse sequences are ideal for collecting high-quality data sets with overnight acquisitions and logically follow the use of fast-HMQC pulse sequences to rapidly screen samples for potential target analytes. Use of the pulse sequences considerably improves the performance of NMR spectroscopy as a complimentary technique for the screening, identification, and validation of chemical warfare agents and other small-molecule analytes present in complex mixtures and environmental samples.

Quantum Bell inequalities from macroscopic locality

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

Yang, Tzyh Haur; Sheridan, Lana; Navascues, Miguel

2011-02-15

We propose a method to generate analytical quantum Bell inequalities based on the principle of macroscopic locality. By imposing locality over binary processings of virtual macroscopic intensities, we establish a correspondence between Bell inequalities and quantum Bell inequalities in bipartite scenarios with dichotomic observables. We discuss how to improve the latter approximation and how to extend our ideas to scenarios with more than two outcomes per setting.

Phonon scattering in nanoscale systems: lowest order expansion of the current and power expressions

NASA Astrophysics Data System (ADS)

Paulsson, Magnus; Frederiksen, Thomas; Brandbyge, Mads

2006-04-01

We use the non-equilibrium Green's function method to describe the effects of phonon scattering on the conductance of nano-scale devices. Useful and accurate approximations are developed that both provide (i) computationally simple formulas for large systems and (ii) simple analytical models. In addition, the simple models can be used to fit experimental data and provide physical parameters.

Transport of phase space densities through tetrahedral meshes using discrete flow mapping

NASA Astrophysics Data System (ADS)

Bajars, Janis; Chappell, David J.; Søndergaard, Niels; Tanner, Gregor

2017-01-01

Discrete flow mapping was recently introduced as an efficient ray based method determining wave energy distributions in complex built up structures. Wave energy densities are transported along ray trajectories through polygonal mesh elements using a finite dimensional approximation of a ray transfer operator. In this way the method can be viewed as a smoothed ray tracing method defined over meshed surfaces. Many applications require the resolution of wave energy distributions in three-dimensional domains, such as in room acoustics, underwater acoustics and for electromagnetic cavity problems. In this work we extend discrete flow mapping to three-dimensional domains by propagating wave energy densities through tetrahedral meshes. The geometric simplicity of the tetrahedral mesh elements is utilised to efficiently compute the ray transfer operator using a mixture of analytic and spectrally accurate numerical integration. The important issue of how to choose a suitable basis approximation in phase space whilst maintaining a reasonable computational cost is addressed via low order local approximations on tetrahedral faces in the position coordinate and high order orthogonal polynomial expansions in momentum space.

Methods in the study of discrete upper hybrid waves

NASA Astrophysics Data System (ADS)

Yoon, P. H.; Ye, S.; Labelle, J.; Weatherwax, A. T.; Menietti, J. D.

2007-11-01

Naturally occurring plasma waves characterized by fine frequency structure or discrete spectrum, detected by satellite, rocket-borne instruments, or ground-based receivers, can be interpreted as eigenmodes excited and trapped in field-aligned density structures. This paper overviews various theoretical methods to study such phenomena for a one-dimensional (1-D) density structure. Among the various methods are parabolic approximation, eikonal matching, eigenfunction matching, and full numerical solution based upon shooting method. Various approaches are compared against the full numerical solution. Among the analytic methods it is found that the eigenfunction matching technique best approximates the actual numerical solution. The analysis is further extended to 2-D geometry. A detailed comparative analysis between the eigenfunction matching and fully numerical methods is carried out for the 2-D case. Although in general the two methods compare favorably, significant differences are also found such that for application to actual observations it is prudent to employ the fully numerical method. Application of the methods developed in the present paper to actual geophysical problems will be given in a companion paper.

Ribozyme-mediated signal augmentation on a mass-sensitive biosensor.

PubMed

Knudsen, Scott M; Lee, Joonhyung; Ellington, Andrew D; Savran, Cagri A

2006-12-20

Mass-based detection methods such as the quartz crystal microbalance (QCM) offer an attractive option to label-based methods; however the sensitivity is generally lower by comparison. In particular, low-molecular-weight analytes can be difficult to detect based on mass addition alone. In this communication, we present the use of effector-dependent ribozymes (aptazymes) as reagents for augmenting small ligand detection on a mass-sensitive device. Two distinct aptazymes were chosen: an L1-ligase-based aptazyme (L1-Rev), which is activated by a small peptide (MW approximately 2.4 kDa) from the HIV-1 Rev protein, and a hammerhead cleavase-based aptazyme (HH-theo3) activated by theophylline (MW = 180 Da). Aptazyme activity was observed in real time, and low-molecular-weight analyte detection has been successfully demonstrated with both aptazymes.

Decoupling of the Leading Order DGLAP Evolution Equation with Spin Dependent Structure Functions

NASA Astrophysics Data System (ADS)

Azadbakht, F. Teimoury; Boroun, G. R.

2018-02-01

We propose an analytical solution for DGLAP evolution equations with polarized splitting functions at the Leading Order (LO) approximation based on the Laplace transform method. It is shown that the DGLAP evolution equations can be decoupled completely into two second order differential equations which then are solved analytically by using the initial conditions δ FS(x,Q2)=F[partial δ FS0(x), δ FS0(x)] and {δ G}(x,Q2)=G[partial δ G0(x), δ G0(x)]. We used this method to obtain the polarized structure function of the proton as well as the polarized gluon distribution function inside the proton and compared the numerical results with experimental data of COMPASS, HERMES, and AAC'08 Collaborations. It was found that there is a good agreement between our predictions and the experiments.

A finite element study of the EIDI system. [Electro-Impulse De-Icing System

NASA Technical Reports Server (NTRS)

Khatkhate, A. A.; Scavuzzo, R. J.; Chu, M. L.

1988-01-01

This paper presents a method for modeling the structural dynamics of an Electro-Impulse De-Icing System, using finite element analyses procedures. A guideline for building a representative finite element model is discussed. Modeling was done initially using four noded cubic elements, four noded isoparametric plate elements and eight noded isoparametric shell elements. Due to the size of the problem and due to the underestimation of shear stress results when compared to previous analytical work an approximate model was created to predict possible areas of shedding of ice. There appears to be good agreement with the test data provided by The Boeing Commercial Airplane Company. Thus these initial results of this method were found to be encouraging. Additional analytical work and comparison with experiment is needed in order to completely evaluate this approach.

Analytical method for the effects of the asteroid belt on planetary orbits

NASA Technical Reports Server (NTRS)

Mayo, A. P.

1979-01-01

Analytic expressions are derived for the perturbation of planetary orbits due to a thick constant-density asteroid belt. The derivations include extensions and adaptations of Plakhov's (1968) analytic expressions for the perturbations in five of the orbital elements for closed orbits around Saturn's rings. The equations of Plakhov are modified to include the effect of ring thickness, and additional equations are derived for the perturbations in the sixth orbital element, the mean anomaly. The gravitational potential and orbital perturbations are derived for the asteroid belt with and without thickness, and for a hoop approximation to the belt. The procedures are also applicable to Saturn's rings and the newly discovered rings of Uranus. The effects of the asteroid belt thickness on the gravitational potential coefficients and the orbital motions are demonstrated. Comparisons between the Mars orbital perturbations obtained by using the analytic expressions and those obtained by numerical integration are discussed. The effects of the asteroid belt on earth-based ranging to Mars are also demonstrated.

Review and assessment of the HOST turbine heat transfer program

NASA Technical Reports Server (NTRS)

Gladden, Herbert J.

1988-01-01

The objectives of the HOST Turbine Heat Transfer subproject were to obtain a better understanding of the physics of the aerothermodynamic phenomena occurring in high-performance gas turbine engines and to assess and improve the analytical methods used to predict the fluid dynamics and heat transfer phenomena. At the time the HOST project was initiated, an across-the-board improvement in turbine design technology was needed. Therefore, a building-block approach was utilized, with research ranging from the study of fundamental phenomena and analytical modeling to experiments in simulated real-engine environments. Experimental research accounted for 75 percent of the project, and analytical efforts accounted for approximately 25 percent. Extensive experimental datasets were created depicting the three-dimensional flow field, high free-stream turbulence, boundary-layer transition, blade tip region heat transfer, film cooling effects in a simulated engine environment, rough-wall cooling enhancement in a rotating passage, and rotor-stator interaction effects. In addition, analytical modeling of these phenomena was initiated using boundary-layer assumptions as well as Navier-Stokes solutions.

On the Singular Perturbations for Fractional Differential Equation

PubMed Central

Atangana, Abdon

2014-01-01

The goal of this paper is to examine the possible extension of the singular perturbation differential equation to the concept of fractional order derivative. To achieve this, we presented a review of the concept of fractional calculus. We make use of the Laplace transform operator to derive exact solution of singular perturbation fractional linear differential equations. We make use of the methodology of three analytical methods to present exact and approximate solution of the singular perturbation fractional, nonlinear, nonhomogeneous differential equation. These methods are including the regular perturbation method, the new development of the variational iteration method, and the homotopy decomposition method. PMID:24683357

Implementing a Matrix-free Analytical Jacobian to Handle Nonlinearities in Models of 3D Lithospheric Deformation

NASA Astrophysics Data System (ADS)

Kaus, B.; Popov, A.

2015-12-01

The analytical expression for the Jacobian is a key component to achieve fast and robust convergence of the nonlinear Newton-Raphson iterative solver. Accomplishing this task in practice often requires a significant algebraic effort. Therefore it is quite common to use a cheap alternative instead, for example by approximating the Jacobian with a finite difference estimation. Despite its simplicity it is a relatively fragile and unreliable technique that is sensitive to the scaling of the residual and unknowns, as well as to the perturbation parameter selection. Unfortunately no universal rule can be applied to provide both a robust scaling and a perturbation. The approach we use here is to derive the analytical Jacobian for the coupled set of momentum, mass, and energy conservation equations together with the elasto-visco-plastic rheology and a marker in cell/staggered finite difference method. The software project LaMEM (Lithosphere and Mantle Evolution Model) is primarily developed for the thermo-mechanically coupled modeling of the 3D lithospheric deformation. The code is based on a staggered grid finite difference discretization in space, and uses customized scalable solvers form PETSc library to efficiently run on the massively parallel machines (such as IBM Blue Gene/Q). Currently LaMEM relies on the Jacobian-Free Newton-Krylov (JFNK) nonlinear solver, which approximates the Jacobian-vector product using a simple finite difference formula. This approach never requires an assembled Jacobian matrix and uses only the residual computation routine. We use an approximate Jacobian (Picard) matrix to precondition the Krylov solver with the Galerkin geometric multigrid. Because of the inherent problems of the finite difference Jacobian estimation, this approach doesn't always result in stable convergence. In this work we present and discuss a matrix-free technique in which the Jacobian-vector product is replaced by analytically-derived expressions and compare results with those obtained with a finite difference approximation of the Jacobian. This project is funded by ERC Starting Grant 258830 and computer facilities were provided by Jülich supercomputer center (Germany).

PPM mixtures of formaldehyde in gas cylinders: Stability and analysis

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

Wong, K.C.; Miller, S.B.; Patterson, L.M.

1999-07-01

Scott Specialty Gases has been successful in producing stable calibration gases of formaldehyde at low concentration. Critical to this success has been the development of a treatment process for high pressure aluminum cylinders. Formaldehyde cylinders having concentrations of 20ppm and 4ppm were found to show only small decline in concentrations over a period of approximately 12 months. Since no NIST traceable formaldehyde standards (or Standard Reference Material) are available, all Scott's formaldehyde cylinders were originally certified by traditional impinger method. This method involves an extremely tedious purification procedure for 2,4-dinitrophenylhydrazine (2,4-DNPH). A modified version of the impinger method has beenmore » developed and does not require extensive reagent purification for formaldehyde analysis. Extremely low formaldehyde blanks have been obtained with the modified method. The HPLC conditions in the original method were used for chromatographic separations. The modified method results in a lower analytical uncertainty for the formaldehyde standard mixtures. Consequently, it is possible to discern small differences between analytical results that are important for stability study.« less

Analytical estimation on divergence and flutter vibrations of symmetrical three-phase induction stator via field-synchronous coordinates

NASA Astrophysics Data System (ADS)

Xia, Ying; Wang, Shiyu; Sun, Wenjia; Xiu, Jie

2017-01-01

The electromagnetically induced parametric vibration of the symmetrical three-phase induction stator is examined. While it can be analyzed by an approximate analytical or numerical method, more accurate and simple analytical method is desirable. This work proposes a new method based on the field-synchronous coordinates. A mechanical-electromagnetic coupling model is developed under this frame such that a time-invariant governing equation with gyroscopic term can be developed. With the general vibration theory, the eigenvalue is formulated; the transition curves between the stable and unstable regions, and response are all determined as closed-form expressions of basic mechanical-electromagnetic parameters. The dependence of these parameters on the instability behaviors is demonstrated. The results imply that the divergence and flutter instabilities can occur even for symmetrical motors with balanced, constant amplitude and sinusoidal voltage. To verify the analytical predictions, this work also builds up a time-variant model of the same system under the conventional inertial frame. The Floquét theory is employed to predict the parametric instability and the numerical integration is used to obtain the parametric response. The parametric instability and response are both well compared against those under the field-synchronous coordinates. The proposed field-synchronous coordinates allows a quick estimation on the electromagnetically induced vibration. The convenience offered by the body-fixed coordinates is discussed across various fields.

Sampling and sample processing in pesticide residue analysis.

PubMed

Lehotay, Steven J; Cook, Jo Marie

2015-05-13

Proper sampling and sample processing in pesticide residue analysis of food and soil have always been essential to obtain accurate results, but the subject is becoming a greater concern as approximately 100 mg test portions are being analyzed with automated high-throughput analytical methods by agrochemical industry and contract laboratories. As global food trade and the importance of monitoring increase, the food industry and regulatory laboratories are also considering miniaturized high-throughput methods. In conjunction with a summary of the symposium "Residues in Food and Feed - Going from Macro to Micro: The Future of Sample Processing in Residue Analytical Methods" held at the 13th IUPAC International Congress of Pesticide Chemistry, this is an opportune time to review sampling theory and sample processing for pesticide residue analysis. If collected samples and test portions do not adequately represent the actual lot from which they came and provide meaningful results, then all costs, time, and efforts involved in implementing programs using sophisticated analytical instruments and techniques are wasted and can actually yield misleading results. This paper is designed to briefly review the often-neglected but crucial topic of sample collection and processing and put the issue into perspective for the future of pesticide residue analysis. It also emphasizes that analysts should demonstrate the validity of their sample processing approaches for the analytes/matrices of interest and encourages further studies on sampling and sample mass reduction to produce a test portion.

Error analysis regarding the calculation of nonlinear force-free field

NASA Astrophysics Data System (ADS)

Liu, S.; Zhang, H. Q.; Su, J. T.

2012-02-01

Magnetic field extrapolation is an alternative method to study chromospheric and coronal magnetic fields. In this paper, two semi-analytical solutions of force-free fields (Low and Lou in Astrophys. J. 352:343, 1990) have been used to study the errors of nonlinear force-free (NLFF) fields based on force-free factor α. Three NLFF fields are extrapolated by approximate vertical integration (AVI) Song et al. (Astrophys. J. 649:1084, 2006), boundary integral equation (BIE) Yan and Sakurai (Sol. Phys. 195:89, 2000) and optimization (Opt.) Wiegelmann (Sol. Phys. 219:87, 2004) methods. Compared with the first semi-analytical field, it is found that the mean values of absolute relative standard deviations (RSD) of α along field lines are about 0.96-1.19, 0.63-1.07 and 0.43-0.72 for AVI, BIE and Opt. fields, respectively. While for the second semi-analytical field, they are about 0.80-1.02, 0.67-1.34 and 0.33-0.55 for AVI, BIE and Opt. fields, respectively. As for the analytical field, the calculation error of <| RSD|> is about 0.1˜0.2. It is also found that RSD does not apparently depend on the length of field line. These provide the basic estimation on the deviation of extrapolated field obtained by proposed methods from the real force-free field.

Nonlinear core deflection in injection molding

NASA Astrophysics Data System (ADS)

Poungthong, P.; Giacomin, A. J.; Saengow, C.; Kolitawong, C.; Liao, H.-C.; Tseng, S.-C.

2018-05-01

Injection molding of thin slender parts is often complicated by core deflection. This deflection is caused by molten plastics race tracking through the slit between the core and the rigid cavity wall. The pressure of this liquid exerts a lateral force of the slender core causing the core to bend, and this bending is governed by a nonlinear fifth order ordinary differential equation for the deflection that is not directly in the position along the core. Here we subject this differential equation to 6 sets of boundary conditions, corresponding to 6 commercial core constraints. For each such set of boundary conditions, we develop an explicit approximate analytical solution, including both a linear term and a nonlinear term. By comparison with finite difference solutions, we find our new analytical solutions to be accurate. We then use these solutions to derive explicit analytical approximations for maximum deflections and for the core position of these maximum deflections. Our experiments on the base-gated free-tip boundary condition agree closely with our new explicit approximate analytical solution.

Modeling of phonon scattering in n-type nanowire transistors using one-shot analytic continuation technique

NASA Astrophysics Data System (ADS)

Bescond, Marc; Li, Changsheng; Mera, Hector; Cavassilas, Nicolas; Lannoo, Michel

2013-10-01

We present a one-shot current-conserving approach to model the influence of electron-phonon scattering in nano-transistors using the non-equilibrium Green's function formalism. The approach is based on the lowest order approximation (LOA) to the current and its simplest analytic continuation (LOA+AC). By means of a scaling argument, we show how both LOA and LOA+AC can be easily obtained from the first iteration of the usual self-consistent Born approximation (SCBA) algorithm. Both LOA and LOA+AC are then applied to model n-type silicon nanowire field-effect-transistors and are compared to SCBA current characteristics. In this system, the LOA fails to describe electron-phonon scattering, mainly because of the interactions with acoustic phonons at the band edges. In contrast, the LOA+AC still well approximates the SCBA current characteristics, thus demonstrating the power of analytic continuation techniques. The limits of validity of LOA+AC are also discussed, and more sophisticated and general analytic continuation techniques are suggested for more demanding cases.

Indexing and retrieving point and region objects

NASA Astrophysics Data System (ADS)

Ibrahim, Azzam T.; Fotouhi, Farshad A.

1996-03-01

R-tree and its variants are examples of spatial data structures for paged-secondary memory. To process a query, these structures require multiple path traversals. In this paper, we present a new image access method, SB+-tree which requires a single path traversal to process a query. Also, SB+-tree will allow commercial databases an access method for spatial objects without a major change, since most commercial databases already support B+-tree as an access method for text data. The SB+-tree can be used for zero and non-zero size data objects. Non-zero size objects are approximated by their minimum bounding rectangles (MBRs). The number of SB+-trees generated is dependent upon the number of dimensions of the approximation of the object. The structure supports efficient spatial operations such as regions-overlap, distance and direction. In this paper, we experimentally and analytically demonstrate the superiority of SB+-tree over R-tree.

Preliminary design methods for fiber reinforced composite structures employing a personal computer

NASA Technical Reports Server (NTRS)

Eastlake, C. N.

1986-01-01

The objective of this project was to develop a user-friendly interactive computer program to be used as an analytical tool by structural designers. Its intent was to do preliminary, approximate stress analysis to help select or verify sizing choices for composite structural members. The approach to the project was to provide a subroutine which uses classical lamination theory to predict an effective elastic modulus for a laminate of arbitrary material and ply orientation. This effective elastic modulus can then be used in a family of other subroutines which employ the familiar basic structural analysis methods for isotropic materials. This method is simple and convenient to use but only approximate, as is appropriate for a preliminary design tool which will be subsequently verified by more sophisticated analysis. Additional subroutines have been provided to calculate laminate coefficient of thermal expansion and to calculate ply-by-ply strains within a laminate.

A third-order approximation method for three-dimensional wheel-rail contact

NASA Astrophysics Data System (ADS)

Negretti, Daniele

2012-03-01

Multibody train analysis is used increasingly by railway operators whenever a reliable and time-efficient method to evaluate the contact between wheel and rail is needed; particularly, the wheel-rail contact is one of the most important aspects that affects a reliable and time-efficient vehicle dynamics computation. The focus of the approach proposed here is to carry out such tasks by means of online wheel-rail elastic contact detection. In order to improve efficiency and save time, a main analytical approach is used for the definition of wheel and rail surfaces as well as for contact detection, then a final numerical evaluation is used to locate contact. The final numerical procedure consists in finding the zeros of a nonlinear function in a single variable. The overall method is based on the approximation of the wheel surface, which does not influence the contact location significantly, as shown in the paper.

A finite difference Davidson procedure to sidestep full ab initio hessian calculation: Application to characterization of stationary points and transition state searches

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

Sharada, Shaama Mallikarjun; Bell, Alexis T., E-mail: mhg@bastille.cchem.berkeley.edu, E-mail: bell@cchem.berkeley.edu; Head-Gordon, Martin, E-mail: mhg@bastille.cchem.berkeley.edu, E-mail: bell@cchem.berkeley.edu

2014-04-28

The cost of calculating nuclear hessians, either analytically or by finite difference methods, during the course of quantum chemical analyses can be prohibitive for systems containing hundreds of atoms. In many applications, though, only a few eigenvalues and eigenvectors, and not the full hessian, are required. For instance, the lowest one or two eigenvalues of the full hessian are sufficient to characterize a stationary point as a minimum or a transition state (TS), respectively. We describe here a method that can eliminate the need for hessian calculations for both the characterization of stationary points as well as searches for saddlemore » points. A finite differences implementation of the Davidson method that uses only first derivatives of the energy to calculate the lowest eigenvalues and eigenvectors of the hessian is discussed. This method can be implemented in conjunction with geometry optimization methods such as partitioned-rational function optimization (P-RFO) to characterize stationary points on the potential energy surface. With equal ease, it can be combined with interpolation methods that determine TS guess structures, such as the freezing string method, to generate approximate hessian matrices in lieu of full hessians as input to P-RFO for TS optimization. This approach is shown to achieve significant cost savings relative to exact hessian calculation when applied to both stationary point characterization as well as TS optimization. The basic reason is that the present approach scales one power of system size lower since the rate of convergence is approximately independent of the size of the system. Therefore, the finite-difference Davidson method is a viable alternative to full hessian calculation for stationary point characterization and TS search particularly when analytical hessians are not available or require substantial computational effort.« less

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

Sakaguchi, Hidetsugu; Malomed, Boris A.

We study ordinary solitons and gap solitons (GS's) in the framework of the one-dimensional Gross-Pitaevskii equation (GPE) with a combination of both linear and nonlinear lattice potentials. The main points of the analysis are the effects of (in)commensurability between the lattices, the development of analytical methods, viz., the variational approximation (VA) for narrow ordinary solitons and various forms of the averaging method for broad solitons of both types, and also the study of the mobility of the solitons. Under the direct commensurability (equal periods of the lattices, L{sub lin}=L{sub nonlin}), the family of ordinary solitons is similar to its counterpartmore » in the GPE without external potentials. In the case of the subharmonic commensurability with L{sub lin}=(1/2)L{sub nonlin}, or incommensurability, there is an existence threshold for the ordinary solitons and the scaling relation between their amplitude and width is different from that in the absence of the potentials. GS families demonstrate a bistability unless the direct commensurability takes place. Specific scaling relations are found for them as well. Ordinary solitons can be readily set in motion by kicking. GS's are also mobile and feature inelastic collisions. The analytical approximations are shown to be quite accurate, predicting correct scaling relations for the soliton families in different cases. The stability of the ordinary solitons is fully determined by the Vakhitov-Kolokolov (VK) criterion (i.e., a negative slope in the dependence between the solitons's chemical potential mu and norm N). The stability of GS families obeys an inverted ('anti-VK') criterion dmu/dN>0, which is explained by the approximation based on the averaging method. The present system provides for the unique possibility to check the anti-VK criterion, as mu(N) dependencies for GS's feature turning points except in the case of direct commensurability.« less

Characterization of Perovskite Oxide/Semiconductor Heterostructures

NASA Astrophysics Data System (ADS)

Walker, Phillip

The tools developed for the use of investigating dynamical systems have provided critical understanding to a wide range of physical phenomena. Here these tools are used to gain further insight into scalar transport, and how it is affected by mixing. The aim of this research is to investigate the efficiency of several different partitioning methods which demarcate flow fields into dynamically distinct regions, and the correlation of finite-time statistics from the advection-diffusion equation to these regions. For autonomous systems, invariant manifold theory can be used to separate the system into dynamically distinct regions. Despite there being no equivalent method for nonautonomous systems, a similar analysis can be done. Systems with general time dependencies must resort to using finite-time transport barriers for partitioning; these barriers are the edges of Lagrangian coherent structures (LCS), the analog to the stable and unstable manifolds of invariant manifold theory. Using the coherent structures of a flow to analyze the statistics of trapping, flight, and residence times, the signature of anomalous diffusion are obtained. This research also investigates the use of linear models for approximating the elements of the covariance matrix of nonlinear flows, and then applying the covariance matrix approximation over coherent regions. The first and second-order moments can be used to fully describe an ensemble evolution in linear systems, however there is no direct method for nonlinear systems. The problem is only compounded by the fact that the moments for nonlinear flows typically don't have analytic representations, therefore direct numerical simulations would be needed to obtain the moments throughout the domain. To circumvent these many computations, the nonlinear system is approximated as many linear systems for which analytic expressions for the moments exist. The parameters introduced in the linear models are obtained locally from the nonlinear deformation tensor.

An analytic, approximate method for modeling steady, three-dimensional flow to partially penetrating wells

NASA Astrophysics Data System (ADS)

Bakker, Mark

2001-05-01

An analytic, approximate solution is derived for the modeling of three-dimensional flow to partially penetrating wells. The solution is written in terms of a correction on the solution for a fully penetrating well and is obtained by dividing the aquifer up, locally, in a number of aquifer layers. The resulting system of differential equations is solved by application of the theory for multiaquifer flow. The presented approach has three major benefits. First, the solution may be applied to any groundwater model that can simulate flow to a fully penetrating well; the solution may be superimposed onto the solution for the fully penetrating well to simulate the local three-dimensional drawdown and flow field. Second, the approach is applicable to isotropic, anisotropic, and stratified aquifers and to both confined and unconfined flow. Third, the solution extends over a small area around the well only; outside this area the three-dimensional effect of the partially penetrating well is negligible, and no correction to the fully penetrating well is needed. A number of comparisons are made to existing three-dimensional, analytic solutions, including radial confined and unconfined flow and a well in a uniform flow field. It is shown that a subdivision in three layers is accurate for many practical cases; very accurate solutions are obtained with more layers.

An approach to estimate spatial distribution of analyte within cells using spectrally-resolved fluorescence microscopy.

PubMed

Sharma, Dharmendar Kumar; Irfanullah, Mir; Basu, Santanu Kumar; Madhu, Sheri; De, Suman; Jadhav, Sameer; Ravikanth, Mangalampalli; Chowdhury, Arindam

2017-01-18

While fluorescence microscopy has become an essential tool amongst chemists and biologists for the detection of various analyte within cellular environments, non-uniform spatial distribution of sensors within cells often restricts extraction of reliable information on relative abundance of analytes in different subcellular regions. As an alternative to existing sensing methodologies such as ratiometric or FRET imaging, where relative proportion of analyte with respect to the sensor can be obtained within cells, we propose a methodology using spectrally-resolved fluorescence microscopy, via which both the relative abundance of sensor as well as their relative proportion with respect to the analyte can be simultaneously extracted for local subcellular regions. This method is exemplified using a BODIPY sensor, capable of detecting mercury ions within cellular environments, characterized by spectral blue-shift and concurrent enhancement of emission intensity. Spectral emission envelopes collected from sub-microscopic regions allowed us to compare the shift in transition energies as well as integrated emission intensities within various intracellular regions. Construction of a 2D scatter plot using spectral shifts and emission intensities, which depend on the relative amount of analyte with respect to sensor and the approximate local amounts of the probe, respectively, enabled qualitative extraction of relative abundance of analyte in various local regions within a single cell as well as amongst different cells. Although the comparisons remain semi-quantitative, this approach involving analysis of multiple spectral parameters opens up an alternative way to extract spatial distribution of analyte in heterogeneous systems. The proposed method would be especially relevant for fluorescent probes that undergo relatively nominal shift in transition energies compared to their emission bandwidths, which often restricts their usage for quantitative ratiometric imaging in cellular media due to strong cross-talk between energetically separated detection channels.

An approach to estimate spatial distribution of analyte within cells using spectrally-resolved fluorescence microscopy

NASA Astrophysics Data System (ADS)

Sharma, Dharmendar Kumar; Irfanullah, Mir; Basu, Santanu Kumar; Madhu, Sheri; De, Suman; Jadhav, Sameer; Ravikanth, Mangalampalli; Chowdhury, Arindam

2017-03-01

While fluorescence microscopy has become an essential tool amongst chemists and biologists for the detection of various analyte within cellular environments, non-uniform spatial distribution of sensors within cells often restricts extraction of reliable information on relative abundance of analytes in different subcellular regions. As an alternative to existing sensing methodologies such as ratiometric or FRET imaging, where relative proportion of analyte with respect to the sensor can be obtained within cells, we propose a methodology using spectrally-resolved fluorescence microscopy, via which both the relative abundance of sensor as well as their relative proportion with respect to the analyte can be simultaneously extracted for local subcellular regions. This method is exemplified using a BODIPY sensor, capable of detecting mercury ions within cellular environments, characterized by spectral blue-shift and concurrent enhancement of emission intensity. Spectral emission envelopes collected from sub-microscopic regions allowed us to compare the shift in transition energies as well as integrated emission intensities within various intracellular regions. Construction of a 2D scatter plot using spectral shifts and emission intensities, which depend on the relative amount of analyte with respect to sensor and the approximate local amounts of the probe, respectively, enabled qualitative extraction of relative abundance of analyte in various local regions within a single cell as well as amongst different cells. Although the comparisons remain semi-quantitative, this approach involving analysis of multiple spectral parameters opens up an alternative way to extract spatial distribution of analyte in heterogeneous systems. The proposed method would be especially relevant for fluorescent probes that undergo relatively nominal shift in transition energies compared to their emission bandwidths, which often restricts their usage for quantitative ratiometric imaging in cellular media due to strong cross-talk between energetically separated detection channels. Dedicated to Professor Kankan Bhattacharyya.

Computational Aspects of Sensitivity Calculations in Linear Transient Structural Analysis. Ph.D. Thesis

NASA Technical Reports Server (NTRS)

Greene, William H.

1989-01-01

A study has been performed focusing on the calculation of sensitivities of displacements, velocities, accelerations, and stresses in linear, structural, transient response problems. One significant goal was to develop and evaluate sensitivity calculation techniques suitable for large-order finite element analyses. Accordingly, approximation vectors such as vibration mode shapes are used to reduce the dimensionality of the finite element model. Much of the research focused on the accuracy of both response quantities and sensitivities as a function of number of vectors used. Two types of sensitivity calculation techniques were developed and evaluated. The first type of technique is an overall finite difference method where the analysis is repeated for perturbed designs. The second type of technique is termed semianalytical because it involves direct, analytical differentiation of the equations of motion with finite difference approximation of the coefficient matrices. To be computationally practical in large-order problems, the overall finite difference methods must use the approximation vectors from the original design in the analyses of the perturbed models.

Gaussian process regression for sensor networks under localization uncertainty

USGS Publications Warehouse

Jadaliha, M.; Xu, Yunfei; Choi, Jongeun; Johnson, N.S.; Li, Weiming

2013-01-01

In this paper, we formulate Gaussian process regression with observations under the localization uncertainty due to the resource-constrained sensor networks. In our formulation, effects of observations, measurement noise, localization uncertainty, and prior distributions are all correctly incorporated in the posterior predictive statistics. The analytically intractable posterior predictive statistics are proposed to be approximated by two techniques, viz., Monte Carlo sampling and Laplace's method. Such approximation techniques have been carefully tailored to our problems and their approximation error and complexity are analyzed. Simulation study demonstrates that the proposed approaches perform much better than approaches without considering the localization uncertainty properly. Finally, we have applied the proposed approaches on the experimentally collected real data from a dye concentration field over a section of a river and a temperature field of an outdoor swimming pool to provide proof of concept tests and evaluate the proposed schemes in real situations. In both simulation and experimental results, the proposed methods outperform the quick-and-dirty solutions often used in practice.

Ballistic-diffusive approximation for the thermal dynamics of metallic nanoparticles in nanocomposite materials

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

Shirdel-Havar, A. H., E-mail: Amir.hushang.shirdel@gmail.com; Masoudian Saadabad, R.

2015-03-21

Based on ballistic-diffusive approximation, a method is presented to model heat transfer in nanocomposites containing metal nanoparticles. This method provides analytical expression for the temperature dynamics of metallic nanoparticles embedded in a dielectric medium. In this study, nanoparticles are considered as spherical shells, so that Boltzmann equation is solved using ballistic-diffusive approximation to calculate the electron and lattice thermal dynamics in gold nanoparticles, while thermal exchange between the particles is taken into account. The model was used to investigate the influence of particle size and metal concentration of the medium on the electron and lattice thermal dynamics. It is shownmore » that these two parameters are crucial in determining the nanocomposite thermal behavior. Our results showed that the heat transfer rate from nanoparticles to the matrix decreases as the nanoparticle size increases. On the other hand, increasing the metal concentration of the medium can also decrease the heat transfer rate.« less

Swiss-cheese models and the Dyer-Roeder approximation

NASA Astrophysics Data System (ADS)

Fleury, Pierre

2014-06-01

In view of interpreting the cosmological observations precisely, especially when they involve narrow light beams, it is crucial to understand how light propagates in our statistically homogeneous, clumpy, Universe. Among the various approaches to tackle this issue, Swiss-cheese models propose an inhomogeneous spacetime geometry which is an exact solution of Einstein's equation, while the Dyer-Roeder approximation deals with inhomogeneity in an effective way. In this article, we demonstrate that the distance-redshift relation of a certain class of Swiss-cheese models is the same as the one predicted by the Dyer-Roeder approach, at a well-controlled level of approximation. Both methods are therefore equivalent when applied to the interpretation of, e.g., supernova obervations. The proof relies on completely analytical arguments, and is illustrated by numerical results.

Charge noise in quantum dot qubits: beyond the Markovian approximation.

NASA Astrophysics Data System (ADS)

Yang, Yuan-Chi; Friesen, Mark; Coppersmith, S. N.

Charge noise is a limiting factor in the performance of semiconductor quantum dot qubits, including both spin and charge qubits. In this work, we develop an analytical formalism for treating semiclassical noise beyond the Markovian approximation, which allows us to investigate noise models relevant for quantum dots, such as 1 / f noise. We apply our methods to both charge qubits and quantum dot hybrid qubits, and study the effects of charge noise on single-qubit rotations in these systems. The formalism is also directly applicable to the case of strong microwave driving, for which the rotating wave approximation breaks down. This work was supported in part by ARO (W911NF-12-0607) and ONR (N00014-15-1-0029), and the University of Wisconsin-Madison.

Effectiveness of “Thin-Layer” and “Effective Medium” Approximations in Numerical Simulation of Dielectric Spectra of Biological Cell Suspensions

NASA Astrophysics Data System (ADS)

Asami, Koji

2010-12-01

There are a few concerns in dielectric modeling of biological cells by the finite-element method (FEM) to simulate their dielectric spectra. Cells possess thin plasma membranes and membrane-bound intracellular organelles, requiring extra fine meshes and considerable computational tasks in the simulation. To solve the problems, the “thin-layer” approximation (TLA) and the “effective medium” approximation (EMA) were adopted. TLA deals with the membrane as an interface of the specific membrane impedance, and therefore it is not necessary to divide the membrane region. EMA regards the composite cytoplasm as an effective homogeneous phase whose dielectric properties are calculated separately. It was proved that TLA and EMA were both useful for greatly reducing computational tasks while accurately coinciding with analytical solutions.

Thin airfoil theory based on approximate solution of the transonic flow equation

NASA Technical Reports Server (NTRS)

Spreiter, John R; Alksne, Alberta Y

1957-01-01

A method is presented for the approximate solution of the nonlinear equations transonic flow theory. Solutions are found for two-dimensional flows at a Mach number of 1 and for purely subsonic and purely supersonic flows. Results are obtained in closed analytic form for a large and significant class of nonlifting airfoils. At a Mach number of 1 general expressions are given for the pressure distribution on an airfoil of specified geometry and for the shape of an airfoil having a prescribed pressure distribution. Extensive comparisons are made with available data, particularly for a Mach number of 1, and with existing solutions.

CHClF2 /F-22/ in the earth's atmosphere

NASA Technical Reports Server (NTRS)

Rasmussen, R. A.; Khalil, M. A. K.; Penkett, S. A.; Prosser, N. J. D.

1980-01-01

Recent global measurements of CHClF2 (F-22) are reported. Originally, GC/MS techniques were used to obtain these data. Since then, significant advances using an O2-doped electron capture detector have been made in the analytical techniques, so that F-22 can be measured by EC/GC methods at ambient concentrations. The atmospheric burden of F-22 calculated from these measurements (average mixing ratio, mid-1979, approximately 45 pptv) is considerably greater than that expected from the estimates of direct industrial emissions (average mixing ratio, mid-1979, approximately 30 pptv). This difference is probably due to underestimates of F-22 emissions.

The Surface Density Distribution in the Solar Nebula

NASA Technical Reports Server (NTRS)

Davis, Sanford S.

2004-01-01

The commonly used minimum mass power law representation of the pre-solar nebula is reanalyzed using a new cumulative-mass-model. This model predicts a smoother surface density approximation compared with methods based on direct computation of surface density. The density is quantified using two independent analytical formulations. First, a best-fit transcendental function is applied directly to the basic planetary data. Next a solution to the time-dependent disk evolution equation is parametrically adapted to the solar nebula data. The latter model is shown to be a good approximation to the finite-size early Solar Nebula, and by extension to other extra solar protoplanetary disks.

Application of geometric approximation to the CPMG experiment: Two- and three-site exchange.

PubMed

Chao, Fa-An; Byrd, R Andrew

2017-04-01

The Carr-Purcell-Meiboom-Gill (CPMG) experiment is one of the most classical and well-known relaxation dispersion experiments in NMR spectroscopy, and it has been successfully applied to characterize biologically relevant conformational dynamics in many cases. Although the data analysis of the CPMG experiment for the 2-site exchange model can be facilitated by analytical solutions, the data analysis in a more complex exchange model generally requires computationally-intensive numerical analysis. Recently, a powerful computational strategy, geometric approximation, has been proposed to provide approximate numerical solutions for the adiabatic relaxation dispersion experiments where analytical solutions are neither available nor feasible. Here, we demonstrate the general potential of geometric approximation by providing a data analysis solution of the CPMG experiment for both the traditional 2-site model and a linear 3-site exchange model. The approximate numerical solution deviates less than 0.5% from the numerical solution on average, and the new approach is computationally 60,000-fold more efficient than the numerical approach. Moreover, we find that accurate dynamic parameters can be determined in most cases, and, for a range of experimental conditions, the relaxation can be assumed to follow mono-exponential decay. The method is general and applicable to any CPMG RD experiment (e.g. N, C', C α , H α , etc.) The approach forms a foundation of building solution surfaces to analyze the CPMG experiment for different models of 3-site exchange. Thus, the geometric approximation is a general strategy to analyze relaxation dispersion data in any system (biological or chemical) if the appropriate library can be built in a physically meaningful domain. Published by Elsevier Inc.

H2+, HeH and H2: Approximating potential curves, calculating rovibrational states

NASA Astrophysics Data System (ADS)

Olivares-Pilón, Horacio; Turbiner, Alexander V.

2018-06-01

Analytic consideration of the Bohr-Oppenheimer (BO) potential curves for diatomic molecules is proposed: accurate analytic interpolation for a potential curve consistent with its rovibrational spectra is found. It is shown that in the BO approximation for four lowest electronic states 1 sσg and 2 pσu, 2 pπu and 3 dπg of H2+, the ground state X2Σ+ of HeH and the two lowest states 1 Σg+ and 3 Σu+ of H2, the potential curves can be analytically interpolated in full range of internuclear distances R with not less than 4-5-6 s.d. Approximation based on matching the Laurant-type expansion at small R and a combination of the multipole expansion with one-instanton type contribution at large distances R is given by two-point Padé approximant. The position of minimum, when exists, is predicted within 1% or better. For the molecular ion H2+ in the Lagrange mesh method, the spectra of vibrational, rotational and rovibrational states (ν , L) associated with 1 sσg and 2 pσu, 2 pπu and 3 dπg potential curves are calculated. In general, it coincides with spectra found via numerical solution of the Schrödinger equation (when available) within six s.d. It is shown that 1 sσg curve contains 19 vibrational states (ν , 0) , while 2 pσu curve contains a single one (0 , 0) and 2 pπu state contains 12 vibrational states (ν , 0) . In general, 1 sσg electronic curve contains 420 rovibrational states, which increases up to 423 when we are beyond BO approximation. For the state 2 pσu the total number of rovibrational states (all with ν = 0) is equal to 3, within or beyond Bohr-Oppenheimer approximation. As for the state 2 pπu within the Bohr-Oppenheimer approximation the total number of the rovibrational bound states is equal to 284. The state 3 dπg is repulsive, no rovibrational state is found. It is confirmed in Lagrange mesh formalism the statement that the ground state potential curve of the heteronuclear molecule HeH does not support rovibrational states. Accurate analytical expression for the potential curves of the hydrogen molecule H2 for the states 1Σg+ and 3 Σu+ is presented. The ground state 1 Σg+ contains 15 vibrational states (ν , 0) , ν = 0- 14. In general, this state supports 301 rovibrational states. The potential curve of the state 3Σu+ has a shallow minimum: it does not support any rovibrational state, it is repulsive.

Analytical methods involving separation techniques for determination of low-molecular-weight biothiols in human plasma and blood.

PubMed

Isokawa, Muneki; Kanamori, Takahiro; Funatsu, Takashi; Tsunoda, Makoto

2014-08-01

Low-molecular-weight biothiols such as homocysteine, cysteine, and glutathione are metabolites of the sulfur cycle and play important roles in biological processes such as the antioxidant defense network, methionine cycle, and protein synthesis. Thiol concentrations in human plasma and blood are related to diseases such as cardiovascular disease, neurodegenerative disease, and cancer. The concentrations of homocysteine, cysteine, and glutathione in plasma samples from healthy human subjects are approximately in the range of 5-15, 200-300, and 1-5 μM, respectively. Glutathione concentration in the whole blood is in the millimolar range. Measurement of biothiol levels in plasma and blood is thought to be important for understanding the physiological roles and biomarkers for certain diseases. This review summarizes the relationship of biothiols with certain disease as well as pre-analytical treatment and analytical methods for determination of biothiols in human plasma and blood by using high-performance liquid chromatography and capillary electrophoresis coupled with ultraviolet, fluorescence, or chemiluminescence detection; or mass spectrometry. Copyright © 2014 Elsevier B.V. All rights reserved.

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.

Transmission and reflection of terahertz plasmons in two-dimensional plasmonic devices

DOE PAGES

Sydoruk, Oleksiy; Choonee, Kaushal; Dyer, Gregory Conrad

2015-03-10

We found that plasmons in two-dimensional semiconductor devices will be reflected by discontinuities, notably, junctions between gated and non-gated electron channels. The transmitted and reflected plasmons can form spatially- and frequency-varying signals, and their understanding is important for the design of terahertz detectors, oscillators, and plasmonic crystals. Using mode decomposition, we studied terahertz plasmons incident on a junction between a gated and a nongated channel. The plasmon reflection and transmission coefficients were found numerically and analytically and studied between 0.3 and 1 THz for a range of electron densities. At higher frequencies, we could describe the plasmons by a simplifiedmore » model of channels in homogeneous dielectrics, for which the analytical approximations were accurate. At low frequencies, however, the full geometry and mode spectrum had to be taken into account. Moreover, the results agreed with simulations by the finite-element method. As a result, mode decomposition thus proved to be a powerful method for plasmonic devices, combining the rigor of complete solutions of Maxwell's equations with the convenience of analytical expressions.« less

Estimation of the curvature of the solid liquid interface during Bridgman crystal growth

NASA Astrophysics Data System (ADS)

Barat, Catherine; Duffar, Thierry; Garandet, Jean-Paul

1998-11-01

An approximate solution for the solid/liquid interface curvature due to the crucible effect in crystal growth is derived from simple heat flux considerations. The numerical modelling of the problem carried out with the help of the finite element code FIDAP supports the predictions of our analytical expression and allows to identify its range of validity. Experimental interface curvatures, measured in gallium antimonide samples grown by the vertical Bridgman method, are seen to compare satisfactorily to analytical and numerical results. Other literature data are also in fair agreement with the predictions of our models in the case where the amount of heat carried by the crucible is small compared to the overall heat flux.

Analytical solution for a class of network dynamics with mechanical and financial applications

NASA Astrophysics Data System (ADS)

Krejčí, P.; Lamba, H.; Melnik, S.; Rachinskii, D.

2014-09-01

We show that for a certain class of dynamics at the nodes the response of a network of any topology to arbitrary inputs is defined in a simple way by its response to a monotone input. The nodes may have either a discrete or continuous set of states and there is no limit on the complexity of the network. The results provide both an efficient numerical method and the potential for accurate analytic approximation of the dynamics on such networks. As illustrative applications, we introduce a quasistatic mechanical model with objects interacting via frictional forces and a financial market model with avalanches and critical behavior that are generated by momentum trading strategies.

Boundary Electron and Beta Dosimetry-Quantification of the Effects of Dissimilar Media on Absorbed Dose

NASA Astrophysics Data System (ADS)

Nunes, Josane C.

1991-02-01

This work quantifies the changes effected in electron absorbed dose to a soft-tissue equivalent medium when part of this medium is replaced by a material that is not soft -tissue equivalent. That is, heterogeneous dosimetry is addressed. Radionuclides which emit beta particles are the electron sources of primary interest. They are used in brachytherapy and in nuclear medicine: for example, beta -ray applicators made with strontium-90 are employed in certain ophthalmic treatments and iodine-131 is used to test thyroid function. More recent medical procedures under development and which involve beta radionuclides include radioimmunotherapy and radiation synovectomy; the first is a cancer modality and the second deals with the treatment of rheumatoid arthritis. In addition, the possibility of skin surface contamination exists whenever there is handling of radioactive material. Determination of absorbed doses in the examples of the preceding paragraph requires considering boundaries of interfaces. Whilst the Monte Carlo method can be applied to boundary calculations, for routine work such as in clinical situations, or in other circumstances where doses need to be determined quickly, analytical dosimetry would be invaluable. Unfortunately, few analytical methods for boundary beta dosimetry exist. Furthermore, the accuracy of results from both Monte Carlo and analytical methods has to be assessed. Although restricted to one radionuclide, phosphorus -32, the experimental data obtained in this work serve several purposes, one of which is to provide standards against which calculated results can be tested. The experimental data also contribute to the relatively sparse set of published boundary dosimetry data. At the same time, they may be useful in developing analytical boundary dosimetry methodology. The first application of the experimental data is demonstrated. Results from two Monte Carlo codes and two analytical methods, which were developed elsewhere, are compared with experimental data. Monte Carlo results compare satisfactory with experimental results for the boundaries considered. The agreement with experimental results for air interfaces is of particular interest because of discrepancies reported previously by another investigator who used data obtained from a different experimental technique. Results from one of the analytical methods differ significantly from the experimental data obtained here. The second analytical method provided data which approximate experimental results to within 30%. This is encouraging but it remains to be determined whether this method performs equally well for other source energies.

An analytical technique for approximating unsteady aerodynamics in the time domain

NASA Technical Reports Server (NTRS)

Dunn, H. J.

1980-01-01

An analytical technique is presented for approximating unsteady aerodynamic forces in the time domain. The order of elements of a matrix Pade approximation was postulated, and the resulting polynomial coefficients were determined through a combination of least squares estimates for the numerator coefficients and a constrained gradient search for the denominator coefficients which insures stable approximating functions. The number of differential equations required to represent the aerodynamic forces to a given accuracy tends to be smaller than that employed in certain existing techniques where the denominator coefficients are chosen a priori. Results are shown for an aeroelastic, cantilevered, semispan wing which indicate a good fit to the aerodynamic forces for oscillatory motion can be achieved with a matrix Pade approximation having fourth order numerator and second order denominator polynomials.

Analytical approximations to the Hotelling trace for digital x-ray detectors

NASA Astrophysics Data System (ADS)

Clarkson, Eric; Pineda, Angel R.; Barrett, Harrison H.

2001-06-01

The Hotelling trace is the signal-to-noise ratio for the ideal linear observer in a detection task. We provide an analytical approximation for this figure of merit when the signal is known exactly and the background is generated by a stationary random process, and the imaging system is an ideal digital x-ray detector. This approximation is based on assuming that the detector is infinite in extent. We test this approximation for finite-size detectors by comparing it to exact calculations using matrix inversion of the data covariance matrix. After verifying the validity of the approximation under a variety of circumstances, we use it to generate plots of the Hotelling trace as a function of pairs of parameters of the system, the signal and the background.

Application of vector-valued rational approximations to the matrix eigenvalue problem and connections with Krylov subspace methods

NASA Technical Reports Server (NTRS)

Sidi, Avram

1992-01-01

Let F(z) be a vectored-valued function F: C approaches C sup N, which is analytic at z=0 and meromorphic in a neighborhood of z=0, and let its Maclaurin series be given. We use vector-valued rational approximation procedures for F(z) that are based on its Maclaurin series in conjunction with power iterations to develop bona fide generalizations of the power method for an arbitrary N X N matrix that may be diagonalizable or not. These generalizations can be used to obtain simultaneously several of the largest distinct eigenvalues and the corresponding invariant subspaces, and present a detailed convergence theory for them. In addition, it is shown that the generalized power methods of this work are equivalent to some Krylov subspace methods, among them the methods of Arnoldi and Lanczos. Thus, the theory provides a set of completely new results and constructions for these Krylov subspace methods. This theory suggests at the same time a new mode of usage for these Krylov subspace methods that were observed to possess computational advantages over their common mode of usage.

Exact and approximate Fourier rebinning algorithms for the solution of the data truncation problem in 3-D PET.

PubMed

Bouallègue, Fayçal Ben; Crouzet, Jean-François; Comtat, Claude; Fourcade, Marjolaine; Mohammadi, Bijan; Mariano-Goulart, Denis

2007-07-01

This paper presents an extended 3-D exact rebinning formula in the Fourier space that leads to an iterative reprojection algorithm (iterative FOREPROJ), which enables the estimation of unmeasured oblique projection data on the basis of the whole set of measured data. In first approximation, this analytical formula also leads to an extended Fourier rebinning equation that is the basis for an approximate reprojection algorithm (extended FORE). These algorithms were evaluated on numerically simulated 3-D positron emission tomography (PET) data for the solution of the truncation problem, i.e., the estimation of the missing portions in the oblique projection data, before the application of algorithms that require complete projection data such as some rebinning methods (FOREX) or 3-D reconstruction algorithms (3DRP or direct Fourier methods). By taking advantage of all the 3-D data statistics, the iterative FOREPROJ reprojection provides a reliable alternative to the classical FOREPROJ method, which only exploits the low-statistics nonoblique data. It significantly improves the quality of the external reconstructed slices without loss of spatial resolution. As for the approximate extended FORE algorithm, it clearly exhibits limitations due to axial interpolations, but will require clinical studies with more realistic measured data in order to decide on its pertinence.

Virtual-source diffusion approximation for enhanced near-field modeling of photon-migration in low-albedo medium.

PubMed

Jia, Mengyu; Chen, Xueying; Zhao, Huijuan; Cui, Shanshan; Liu, Ming; Liu, Lingling; Gao, Feng

2015-01-26

Most analytical methods for describing light propagation in turbid medium exhibit low effectiveness in the near-field of a collimated source. Motivated by the Charge Simulation Method in electromagnetic theory as well as the established discrete source based modeling, we herein report on an improved explicit model for a semi-infinite geometry, referred to as "Virtual Source" (VS) diffuse approximation (DA), to fit for low-albedo medium and short source-detector separation. In this model, the collimated light in the standard DA is analogously approximated as multiple isotropic point sources (VS) distributed along the incident direction. For performance enhancement, a fitting procedure between the calculated and realistic reflectances is adopted in the near-field to optimize the VS parameters (intensities and locations). To be practically applicable, an explicit 2VS-DA model is established based on close-form derivations of the VS parameters for the typical ranges of the optical parameters. This parameterized scheme is proved to inherit the mathematical simplicity of the DA approximation while considerably extending its validity in modeling the near-field photon migration in low-albedo medium. The superiority of the proposed VS-DA method to the established ones is demonstrated in comparison with Monte-Carlo simulations over wide ranges of the source-detector separation and the medium optical properties.

Semi-preparative HPLC preparation and HPTLC quantification of tetrahydroamentoflavone as marker in Semecarpus anacardium and its polyherbal formulations.

PubMed

Aravind, S G; Arimboor, Ranjith; Rangan, Meena; Madhavan, Soumya N; Arumughan, C

2008-11-04

Application of modern scientific knowledge coupled with sensitive analytical technique is important for the quality evaluation and standardization of polyherbal formulations. Semecarpus anacardium, an important medicinal plant with wide medicinal properties, is frequently used in a large number of traditional herbal preparations. Tetrahydroamentoflavone (THA), a major bioactive biflavonoid was selected as a chemical marker of S. anacardium and RP-semi-preparative HPLC conditions were optimized for the isolation of tetrahydroamentoflavone. HPTLC analytical method was developed for the fingerprinting of S. anacardium flavonoids and quantification of tetrahydroamentoflavone. The method was validated in terms of their linearity, LOD, LOQ, precision and accuracy and compared with RP-HPLC-DAD method. The methods were demonstrated for the chemical fingerprinting of S. anacardium plant parts and some commercial polyherbal formulations and the amount of tetrahydroamentoflavone was quantified. HPTLC analysis showed that S. anacardium seed contained approximately 10 g kg(-1) of tetrahydroamentoflavone. The methods were able to identify and quantify tetrahydroamentoflavone from complex mixtures of phytochemicals and could be extended to the marker-based standardization of polyherbal formulations, containing S. anacardium.

Rapid fusion method for the determination of refractory thorium and uranium isotopes in soil samples

DOE PAGES

Maxwell, Sherrod L.; Hutchison, Jay B.; McAlister, Daniel R.

2015-02-14

Recently, approximately 80% of participating laboratories failed to accurately determine uranium isotopes in soil samples in the U.S Department of Energy Mixed Analyte Performance Evaluation Program (MAPEP) Session 30, due to incomplete dissolution of refractory particles in the samples. Failing laboratories employed acid dissolution methods, including hydrofluoric acid, to recover uranium from the soil matrix. The failures illustrate the importance of rugged soil dissolution methods for the accurate measurement of analytes in the sample matrix. A new rapid fusion method has been developed by the Savannah River National Laboratory (SRNL) to prepare 1-2 g soil sample aliquots very quickly, withmore » total dissolution of refractory particles. Soil samples are fused with sodium hydroxide at 600 ºC in zirconium crucibles to enable complete dissolution of the sample. Uranium and thorium are separated on stacked TEVA and TRU extraction chromatographic resin cartridges, prior to isotopic measurements by alpha spectrometry on cerium fluoride microprecipitation sources. Plutonium can also be separated and measured using this method. Batches of 12 samples can be prepared for measurement in <5 hours.« less

Rapid analysis of aminoglycoside antibiotics in bovine tissues using disposable pipette extraction and ultrahigh performance liquid chromatography-tandem mass spectrometry.

PubMed

Lehotay, Steven J; Mastovska, Katerina; Lightfield, Alan R; Nuñez, Alberto; Dutko, Terry; Ng, Chilton; Bluhm, Louis

2013-10-25

A high-throughput qualitative screening and identification method for 9 aminoglycosides of regulatory interest has been developed, validated, and implemented for bovine kidney, liver, and muscle tissues. The method involves extraction at previously validated conditions, cleanup using disposable pipette extraction, and analysis by a 3 min ultrahigh-performance liquid chromatography-tandem mass spectrometry (UHPLC-MS/MS) method. The drug analytes include neomycin, streptomycin, dihydrosptreptomycin, and spectinomycin, which have residue tolerances in bovine in the US, and kanamicin, gentamicin, apramycin, amikacin, and hygromycin, which do not have US tolerances established in bovine tissues. Tobramycin was used as an internal standard. An additional drug, paromomycin also was validated in the method, but it was dropped during implementation due to conversion of neomycin into paromomycin. Proposed fragmentation patterns for the monitored ions of each analyte were elucidated with the aid of high resolution MS using a quadrupole-time-of-flight instrument. Recoveries from spiking experiments at regulatory levels of concern showed that all analytes averaged 70-120% recoveries in all tissues, except hygromycin averaged 61% recovery. Lowest calibrated levels were as low as 0.005 μg/g in matrix extracts, which approximately corresponded to the limit of detection for screening purposes. Drug identifications at levels <0.05 μg/g were made in spiked and/or real samples for all analytes and tissues tested. Analyses of 60 samples from 20 slaughtered cattle previously screened positive for aminoglycosides showed that this method worked well in practice. The UHPLC-MS/MS method has several advantages compared to the previous microbial inhibition screening assay, especially for distinguishing individual drugs from a mixture and improving identification of gentamicin in tissue samples. Published by Elsevier B.V.

Utility perspective on USEPA analytical methods program redirection

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

Koch, B.; Davis, M.K.; Krasner, S.W.

1996-11-01

The Metropolitan Water District of Southern California (Metropolitan) is a public, municipal corporation, created by the State of California, which wholesales supplemental water trough 27 member agencies (cities and water districts). Metropolitan serves nearly 16 million people in an area along the coastal plain of Southern California that covers approximately 5200 square miles. Water deliveries have averaged up to 2.5 million acre-feet per year. Metropolitan`s Water Quality Laboratory (WQL) conducts compliance monitoring of its source and finished drinking waters for chemical and microbial constituents. The laboratory maintains certification of a large number and variety of analytical procedures. The WQL operatesmore » in a 17,000-square-foot facility. The equipment is state-of-the-art analytical instrumentation. The staff consists of 40 professional chemists and microbiologists whose experience and expertise are extensive and often highly specialized. The staff turnover is very low, and the laboratory is consistently, efficiently, and expertly run.« less

On the analytical modeling of the nonlinear vibrations of pretensioned space structures

NASA Technical Reports Server (NTRS)

Housner, J. M.; Belvin, W. K.

1983-01-01

Pretensioned structures are receiving considerable attention as candidate large space structures. A typical example is a hoop-column antenna. The large number of preloaded members requires efficient analytical methods for concept validation and design. Validation through analyses is especially important since ground testing may be limited due to gravity effects and structural size. The present investigation has the objective to present an examination of the analytical modeling of pretensioned members undergoing nonlinear vibrations. Two approximate nonlinear analysis are developed to model general structural arrangements which include beam-columns and pretensioned cables attached to a common nucleus, such as may occur at a joint of a pretensioned structure. Attention is given to structures undergoing nonlinear steady-state oscillations due to sinusoidal excitation forces. Three analyses, linear, quasi-linear, and nonlinear are conducted and applied to study the response of a relatively simple cable stiffened structure.

Analytic double product integrals for all-frequency relighting.

PubMed

Wang, Rui; Pan, Minghao; Chen, Weifeng; Ren, Zhong; Zhou, Kun; Hua, Wei; Bao, Hujun

2013-07-01

This paper presents a new technique for real-time relighting of static scenes with all-frequency shadows from complex lighting and highly specular reflections from spatially varying BRDFs. The key idea is to depict the boundaries of visible regions using piecewise linear functions, and convert the shading computation into double product integrals—the integral of the product of lighting and BRDF on visible regions. By representing lighting and BRDF with spherical Gaussians and approximating their product using Legendre polynomials locally in visible regions, we show that such double product integrals can be evaluated in an analytic form. Given the precomputed visibility, our technique computes the visibility boundaries on the fly at each shading point, and performs the analytic integral to evaluate the shading color. The result is a real-time all-frequency relighting technique for static scenes with dynamic, spatially varying BRDFs, which can generate more accurate shadows than the state-of-the-art real-time PRT methods.

Study of a vibrating plate: comparison between experimental (ESPI) and analytical results

NASA Astrophysics Data System (ADS)

Romero, G.; Alvarez, L.; Alanís, E.; Nallim, L.; Grossi, R.

2003-07-01

Real-time electronic speckle pattern interferometry (ESPI) was used for tuning and visualization of natural frequencies of a trapezoidal plate. The plate was excited to resonant vibration by a sinusoidal acoustical source, which provided a continuous range of audio frequencies. Fringe patterns produced during the time-average recording of the vibrating plate—corresponding to several resonant frequencies—were registered. From these interferograms, calculations of vibrational amplitudes by means of zero-order Bessel functions were performed in some particular cases. The system was also studied analytically. The analytical approach developed is based on the Rayleigh-Ritz method and on the use of non-orthogonal right triangular co-ordinates. The deflection of the plate is approximated by a set of beam characteristic orthogonal polynomials generated by using the Gram-Schmidt procedure. A high degree of correlation between computational analysis and experimental results was observed.

Linearized semiclassical initial value time correlation functions with maximum entropy analytic continuation.

PubMed

Liu, Jian; Miller, William H

2008-09-28

The maximum entropy analytic continuation (MEAC) method is used to extend the range of accuracy of the linearized semiclassical initial value representation (LSC-IVR)/classical Wigner approximation for real time correlation functions. LSC-IVR provides a very effective "prior" for the MEAC procedure since it is very good for short times, exact for all time and temperature for harmonic potentials (even for correlation functions of nonlinear operators), and becomes exact in the classical high temperature limit. This combined MEAC+LSC/IVR approach is applied here to two highly nonlinear dynamical systems, a pure quartic potential in one dimensional and liquid para-hydrogen at two thermal state points (25 and 14 K under nearly zero external pressure). The former example shows the MEAC procedure to be a very significant enhancement of the LSC-IVR for correlation functions of both linear and nonlinear operators, and especially at low temperature where semiclassical approximations are least accurate. For liquid para-hydrogen, the LSC-IVR is seen already to be excellent at T=25 K, but the MEAC procedure produces a significant correction at the lower temperature (T=14 K). Comparisons are also made as to how the MEAC procedure is able to provide corrections for other trajectory-based dynamical approximations when used as priors.

A fast method to compute Three-Dimensional Infrared Radiative Transfer in non scattering medium

NASA Astrophysics Data System (ADS)

Makke, Laurent; Musson-Genon, Luc; Carissimo, Bertrand

2014-05-01

The Atmospheric Radiation field has seen the development of more accurate and faster methods to take into account absoprtion in participating media. Radiative fog appears with clear sky condition due to a significant cooling during the night, so scattering is left out. Fog formation modelling requires accurate enough method to compute cooling rates. Thanks to High Performance Computing, multi-spectral approach of Radiative Transfer Equation resolution is most often used. Nevertheless, the coupling of three-dimensionnal radiative transfer with fluid dynamics is very detrimental to the computational cost. To reduce the time spent in radiation calculations, the following method uses analytical absorption functions fitted by Sasamori (1968) on Yamamoto's charts (Yamamoto,1956) to compute a local linear absorption coefficient. By averaging radiative properties, this method eliminates the spectral integration. For an isothermal atmosphere, analytical calculations lead to an explicit formula between emissivities functions and linear absorption coefficient. In the case of cooling to space approximation, this analytical expression gives very accurate results compared to correlated k-distribution. For non homogeneous paths, we propose a two steps algorithm. One-dimensional radiative quantities and linear absorption coefficient are computed by a two-flux method. Then, three-dimensional RTE under the grey medium assumption is solved with the DOM. Comparisons with measurements of radiative quantities during ParisFOG field (2006) shows the cability of this method to handle strong vertical variations of pressure/temperature and gases concentrations.

Numerically stable formulas for a particle-based explicit exponential integrator

NASA Astrophysics Data System (ADS)

Nadukandi, Prashanth

2015-05-01

Numerically stable formulas are presented for the closed-form analytical solution of the X-IVAS scheme in 3D. This scheme is a state-of-the-art particle-based explicit exponential integrator developed for the particle finite element method. Algebraically, this scheme involves two steps: (1) the solution of tangent curves for piecewise linear vector fields defined on simplicial meshes and (2) the solution of line integrals of piecewise linear vector-valued functions along these tangent curves. Hence, the stable formulas presented here have general applicability, e.g. exact integration of trajectories in particle-based (Lagrangian-type) methods, flow visualization and computer graphics. The Newton form of the polynomial interpolation definition is used to express exponential functions of matrices which appear in the analytical solution of the X-IVAS scheme. The divided difference coefficients in these expressions are defined in a piecewise manner, i.e. in a prescribed neighbourhood of removable singularities their series approximations are computed. An optimal series approximation of divided differences is presented which plays a critical role in this methodology. At least ten significant decimal digits in the formula computations are guaranteed to be exact using double-precision floating-point arithmetic. The worst case scenarios occur in the neighbourhood of removable singularities found in fourth-order divided differences of the exponential function.

Establishing the traceability of a uranyl nitrate solution to a standard reference material

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

Jackson, C.H.; Clark, J.P.

1978-01-01

A uranyl nitrate solution for use as a Working Calibration and Test Material (WCTM) was characterized, using a statistically designed procedure to document traceability to National Bureau of Standards Reference Material (SPM-960). A Reference Calibration and Test Material (PCTM) was prepared from SRM-960 uranium metal to approximate the acid and uranium concentration of the WCTM. This solution was used in the characterization procedure. Details of preparing, handling, and packaging these solutions are covered. Two outside laboratories, each having measurement expertise using a different analytical method, were selected to measure both solutions according to the procedure for characterizing the WCTM. Twomore » different methods were also used for the in-house characterization work. All analytical results were tested for statistical agreement before the WCTM concentration and limit of error values were calculated. A concentration value was determined with a relative limit of error (RLE) of approximately 0.03% which was better than the target RLE of 0.08%. The use of this working material eliminates the expense of using SRMs to fulfill traceability requirements for uranium measurements on this type material. Several years' supply of uranyl nitrate solution with NBS traceability was produced. The cost of this material was less than 10% of an equal quantity of SRM-960 uranium metal.« less

Versions of the collocation and least squares method for solving biharmonic equations in non-canonical domains

NASA Astrophysics Data System (ADS)

Belyaev, V. A.; Shapeev, V. P.

2017-10-01

New versions of the collocations and least squares method of high-order accuracy are proposed and implemented for the numerical solution of the boundary value problems for the biharmonic equation in non-canonical domains. The solution of the biharmonic equation is used for simulating the stress-strain state of an isotropic plate under the action of transverse load. The differential problem is projected into a space of fourth-degree polynomials by the CLS method. The boundary conditions for the approximate solution are put down exactly on the boundary of the computational domain. The versions of the CLS method are implemented on the grids which are constructed in two different ways. It is shown that the approximate solution of problems converges with high order. Thus it matches with high accuracy with the analytical solution of the test problems in the case of known solution in the numerical experiments on the convergence of the solution of various problems on a sequence of grids.

Using radiance predicted by the P3 approximation in a spherical geometry to predict tissue optical properties

NASA Astrophysics Data System (ADS)

Dickey, Dwayne J.; Moore, Ronald B.; Tulip, John

2001-01-01

For photodynamic therapy of solid tumors, such as prostatic carcinoma, to be achieved, an accurate model to predict tissue parameters and light dose must be found. Presently, most analytical light dosimetry models are fluence based and are not clinically viable for tissue characterization. Other methods of predicting optical properties, such as Monet Carlo, are accurate but far too time consuming for clinical application. However, radiance predicted by the P3-Approximation, an anaylitical solution to the transport equation, may be a viable and accurate alternative. The P3-Approximation accurately predicts optical parameters in intralipid/methylene blue based phantoms in a spherical geometry. The optical parameters furnished by the radiance, when introduced into fluence predicted by both P3- Approximation and Grosjean Theory, correlate well with experimental data. The P3-Approximation also predicts the optical properties of prostate tissue, agreeing with documented optical parameters. The P3-Approximation could be the clinical tool necessary to facilitate PDT of solid tumors because of the limited number of invasive measurements required and the speed in which accurate calculations can be performed.

Vibrations of a thin cylindrical shell stiffened by rings with various stiffness

NASA Astrophysics Data System (ADS)

Nesterchuk, G. A.

2018-05-01

The problem of vibrations of a thin-walled elastic cylindrical shell reinforced by frames of different rigidity is investigated. The solution for the case of the clamped shell edges was obtained by asymptotic methods and refined by the finite element method. Rings with zero eccentricity and stiffness varying along the generatrix of the shell cylinder are considered. Varying the optimal coefficients of the distribution functions of the rigidity of the frames and finding more precise parameters makes it possible to find correction factors for analytical formulas of approximate calculation.

A binomial stochastic kinetic approach to the Michaelis-Menten mechanism

NASA Astrophysics Data System (ADS)

Lente, Gábor

2013-05-01

This Letter presents a new method that gives an analytical approximation of the exact solution of the stochastic Michaelis-Menten mechanism without computationally demanding matrix operations. The method is based on solving the deterministic rate equations and then using the results as guiding variables of calculating probability values using binomial distributions. This principle can be generalized to a number of different kinetic schemes and is expected to be very useful in the evaluation of measurements focusing on the catalytic activity of one or a few individual enzyme molecules.

Geometry of Thin Nematic Elastomer Sheets

NASA Astrophysics Data System (ADS)

Aharoni, Hillel; Sharon, Eran; Kupferman, Raz

A thin sheet of nematic elastomer attains 3D configurations depending on the nematic director field upon heating. In this talk we describe the intrinsic geometry of such a sheet, and derive an expression for the metric induced by general smooth nematic director fields. Furthermore, we investigate the reverse problem of constructing a director field that induces a specified 2D geometry. We provide an explicit analytical recipe for constructing any surface of revolution using this method. We demonstrate how the design of an arbitrary 2D geometry is accessible using approximate numerical methods.

Partial Wave Dispersion Relations: Application to Electron-Atom Scattering

NASA Technical Reports Server (NTRS)

Temkin, A.; Drachman, Richard J.

1999-01-01

In this Letter we propose the use of partial wave dispersion relations (DR's) as the way of solving the long-standing problem of correctly incorporating exchange in a valid DR for electron-atom scattering. In particular a method is given for effectively calculating the contribution of the discontinuity and/or poles of the partial wave amplitude which occur in the negative E plane. The method is successfully tested in three cases: (i) the analytically solvable exponential potential, (ii) the Hartree potential, and (iii) the S-wave exchange approximation for electron-hydrogen scattering.

An analytic performance model of disk arrays and its application

NASA Technical Reports Server (NTRS)

Lee, Edward K.; Katz, Randy H.

1991-01-01

As disk arrays become widely used, tools for understanding and analyzing their performance become increasingly important. In particular, performance models can be invaluable in both configuring and designing disk arrays. Accurate analytic performance models are desirable over other types of models because they can be quickly evaluated, are applicable under a wide range of system and workload parameters, and can be manipulated by a range of mathematical techniques. Unfortunately, analytical performance models of disk arrays are difficult to formulate due to the presence of queuing and fork-join synchronization; a disk array request is broken up into independent disk requests which must all complete to satisfy the original request. We develop, validate, and apply an analytic performance model for disk arrays. We derive simple equations for approximating their utilization, response time, and throughput. We then validate the analytic model via simulation and investigate the accuracy of each approximation used in deriving the analytical model. Finally, we apply the analytical model to derive an equation for the optimal unit of data striping in disk arrays.

A study of the temporal robustness of the growing global container-shipping network

PubMed Central

Wang, Nuo; Wu, Nuan; Dong, Ling-ling; Yan, Hua-kun; Wu, Di

2016-01-01

Whether they thrive as they grow must be determined for all constantly expanding networks. However, few studies have focused on this important network feature or the development of quantitative analytical methods. Given the formation and growth of the global container-shipping network, we proposed the concept of network temporal robustness and quantitative method. As an example, we collected container liner companies’ data at two time points (2004 and 2014) and built a shipping network with ports as nodes and routes as links. We thus obtained a quantitative value of the temporal robustness. The temporal robustness is a significant network property because, for the first time, we can clearly recognize that the shipping network has become more vulnerable to damage over the last decade: When the node failure scale reached 50% of the entire network, the temporal robustness was approximately −0.51% for random errors and −12.63% for intentional attacks. The proposed concept and analytical method described in this paper are significant for other network studies. PMID:27713549

Subtracting infrared renormalons from Wilson coefficients: Uniqueness and power dependences on ΛQCD

NASA Astrophysics Data System (ADS)

Mishima, Go; Sumino, Yukinari; Takaura, Hiromasa

2017-06-01

In the context of operator product expansion (OPE) and using the large-β0 approximation, we propose a method to define Wilson coefficients free from uncertainties due to IR renormalons. We first introduce a general observable X (Q2) with an explicit IR cutoff, and then we extract a genuine UV contribution XUV as a cutoff-independent part. XUV includes power corrections ˜(ΛQCD2/Q2)n which are independent of renormalons. Using the integration-by-regions method, we observe that XUV coincides with the leading Wilson coefficient in OPE and also clarify that the power corrections originate from UV region. We examine scheme dependence of XUV and single out a specific scheme favorable in terms of analytical properties. Our method would be optimal with respect to systematicity, analyticity and stability. We test our formulation with the examples of the Adler function, QCD force between Q Q ¯, and R -ratio in e+e- collision.

Two-dimensional fracture analysis of piezoelectric material based on the scaled boundary node method

NASA Astrophysics Data System (ADS)

Shen-Shen, Chen; Juan, Wang; Qing-Hua, Li

2016-04-01

A scaled boundary node method (SBNM) is developed for two-dimensional fracture analysis of piezoelectric material, which allows the stress and electric displacement intensity factors to be calculated directly and accurately. As a boundary-type meshless method, the SBNM employs the moving Kriging (MK) interpolation technique to an approximate unknown field in the circumferential direction and therefore only a set of scattered nodes are required to discretize the boundary. As the shape functions satisfy Kronecker delta property, no special techniques are required to impose the essential boundary conditions. In the radial direction, the SBNM seeks analytical solutions by making use of analytical techniques available to solve ordinary differential equations. Numerical examples are investigated and satisfactory solutions are obtained, which validates the accuracy and simplicity of the proposed approach. Project supported by the National Natural Science Foundation of China (Grant Nos. 11462006 and 21466012), the Foundation of Jiangxi Provincial Educational Committee, China (Grant No. KJLD14041), and the Foundation of East China Jiaotong University, China (Grant No. 09130020).

Probabilistic dual heuristic programming-based adaptive critic

NASA Astrophysics Data System (ADS)

Herzallah, Randa

2010-02-01

Adaptive critic (AC) methods have common roots as generalisations of dynamic programming for neural reinforcement learning approaches. Since they approximate the dynamic programming solutions, they are potentially suitable for learning in noisy, non-linear and non-stationary environments. In this study, a novel probabilistic dual heuristic programming (DHP)-based AC controller is proposed. Distinct to current approaches, the proposed probabilistic (DHP) AC method takes uncertainties of forward model and inverse controller into consideration. Therefore, it is suitable for deterministic and stochastic control problems characterised by functional uncertainty. Theoretical development of the proposed method is validated by analytically evaluating the correct value of the cost function which satisfies the Bellman equation in a linear quadratic control problem. The target value of the probabilistic critic network is then calculated and shown to be equal to the analytically derived correct value. Full derivation of the Riccati solution for this non-standard stochastic linear quadratic control problem is also provided. Moreover, the performance of the proposed probabilistic controller is demonstrated on linear and non-linear control examples.

Stationary responses of a Rayleigh viscoelastic system with zero barrier impacts under external random excitation.

PubMed

Wang, Deli; Xu, Wei; Zhao, Xiangrong

2016-03-01

This paper aims to deal with the stationary responses of a Rayleigh viscoelastic system with zero barrier impacts under external random excitation. First, the original stochastic viscoelastic system is converted to an equivalent stochastic system without viscoelastic terms by approximately adding the equivalent stiffness and damping. Relying on the means of non-smooth transformation of state variables, the above system is replaced by a new system without an impact term. Then, the stationary probability density functions of the system are observed analytically through stochastic averaging method. By considering the effects of the biquadratic nonlinear damping coefficient and the noise intensity on the system responses, the effectiveness of the theoretical method is tested by comparing the analytical results with those generated from Monte Carlo simulations. Additionally, it does deserve attention that some system parameters can induce the occurrence of stochastic P-bifurcation.

Green's functions in equilibrium and nonequilibrium from real-time bold-line Monte Carlo

NASA Astrophysics Data System (ADS)

Cohen, Guy; Gull, Emanuel; Reichman, David R.; Millis, Andrew J.

2014-03-01

Green's functions for the Anderson impurity model are obtained within a numerically exact formalism. We investigate the limits of analytical continuation for equilibrium systems, and show that with real time methods even sharp high-energy features can be reliably resolved. Continuing to an Anderson impurity in a junction, we evaluate two-time correlation functions, spectral properties, and transport properties, showing how the correspondence between the spectral function and the differential conductance breaks down when nonequilibrium effects are taken into account. Finally, a long-standing dispute regarding this model has involved the voltage splitting of the Kondo peak, an effect which was predicted over a decade ago by approximate analytical methods but never successfully confirmed by numerics. We settle the issue by demonstrating in an unbiased manner that this splitting indeed occurs. Yad Hanadiv-Rothschild Foundation, TG-DMR120085, TG-DMR130036, NSF CHE-1213247, NSF DMR 1006282, DOE ER 46932.

Self-excited oscillation and monostable operation of a bistable light emitting diode (BILED)

NASA Astrophysics Data System (ADS)

Okumura, K.; Ogawa, Y.; Ito, H.; Inaba, H.

1983-07-01

A new simple opto-electronic bistable device has been obtained by combining a light emitting diode (LED) and a photodetector (PD) with electronic feedback using a broad bandpass filter. This has interesting dynamic characteristics which are expected to have such various applications as optical oscillators, optical pulse generators and optical pulsewidth modulators. The dynamic characteristics are represented by second-order nonlinear differential equations. In the analyses of these nonlinear systems, instead of numerical analyses with a computer, an approximate analytical method devised for this purpose has been used. This method has been used for investigating the characteristics of the proposed device quantitatively. These include the frequency of oscillations, pulsewidths and hysteresis. The results of the analyses agree approximately with experimentally observed values, thus the dynamic characteristics of the proposed device can be explained.

Efficient time-dependent density functional theory approximations for hybrid density functionals: analytical gradients and parallelization.

PubMed

Petrenko, Taras; Kossmann, Simone; Neese, Frank

2011-02-07

In this paper, we present the implementation of efficient approximations to time-dependent density functional theory (TDDFT) within the Tamm-Dancoff approximation (TDA) for hybrid density functionals. For the calculation of the TDDFT/TDA excitation energies and analytical gradients, we combine the resolution of identity (RI-J) algorithm for the computation of the Coulomb terms and the recently introduced "chain of spheres exchange" (COSX) algorithm for the calculation of the exchange terms. It is shown that for extended basis sets, the RIJCOSX approximation leads to speedups of up to 2 orders of magnitude compared to traditional methods, as demonstrated for hydrocarbon chains. The accuracy of the adiabatic transition energies, excited state structures, and vibrational frequencies is assessed on a set of 27 excited states for 25 molecules with the configuration interaction singles and hybrid TDDFT/TDA methods using various basis sets. Compared to the canonical values, the typical error in transition energies is of the order of 0.01 eV. Similar to the ground-state results, excited state equilibrium geometries differ by less than 0.3 pm in the bond distances and 0.5° in the bond angles from the canonical values. The typical error in the calculated excited state normal coordinate displacements is of the order of 0.01, and relative error in the calculated excited state vibrational frequencies is less than 1%. The errors introduced by the RIJCOSX approximation are, thus, insignificant compared to the errors related to the approximate nature of the TDDFT methods and basis set truncation. For TDDFT/TDA energy and gradient calculations on Ag-TB2-helicate (156 atoms, 2732 basis functions), it is demonstrated that the COSX algorithm parallelizes almost perfectly (speedup ~26-29 for 30 processors). The exchange-correlation terms also parallelize well (speedup ~27-29 for 30 processors). The solution of the Z-vector equations shows a speedup of ~24 on 30 processors. The parallelization efficiency for the Coulomb terms can be somewhat smaller (speedup ~15-25 for 30 processors), but their contribution to the total calculation time is small. Thus, the parallel program completes a Becke3-Lee-Yang-Parr energy and gradient calculation on the Ag-TB2-helicate in less than 4 h on 30 processors. We also present the necessary extension of the Lagrangian formalism, which enables the calculation of the TDDFT excited state properties in the frozen-core approximation. The algorithms described in this work are implemented into the ORCA electronic structure system.

Algorithms and analytical solutions for rapidly approximating long-term dispersion from line and area sources

NASA Astrophysics Data System (ADS)

Barrett, Steven R. H.; Britter, Rex E.

Predicting long-term mean pollutant concentrations in the vicinity of airports, roads and other industrial sources are frequently of concern in regulatory and public health contexts. Many emissions are represented geometrically as ground-level line or area sources. Well developed modelling tools such as AERMOD and ADMS are able to model dispersion from finite (i.e. non-point) sources with considerable accuracy, drawing upon an up-to-date understanding of boundary layer behaviour. Due to mathematical difficulties associated with line and area sources, computationally expensive numerical integration schemes have been developed. For example, some models decompose area sources into a large number of line sources orthogonal to the mean wind direction, for which an analytical (Gaussian) solution exists. Models also employ a time-series approach, which involves computing mean pollutant concentrations for every hour over one or more years of meteorological data. This can give rise to computer runtimes of several days for assessment of a site. While this may be acceptable for assessment of a single industrial complex, airport, etc., this level of computational cost precludes national or international policy assessments at the level of detail available with dispersion modelling. In this paper, we extend previous work [S.R.H. Barrett, R.E. Britter, 2008. Development of algorithms and approximations for rapid operational air quality modelling. Atmospheric Environment 42 (2008) 8105-8111] to line and area sources. We introduce approximations which allow for the development of new analytical solutions for long-term mean dispersion from line and area sources, based on hypergeometric functions. We describe how these solutions can be parameterized from a single point source run from an existing advanced dispersion model, thereby accounting for all processes modelled in the more costly algorithms. The parameterization method combined with the analytical solutions for long-term mean dispersion are shown to produce results several orders of magnitude more efficiently with a loss of accuracy small compared to the absolute accuracy of advanced dispersion models near sources. The method can be readily incorporated into existing dispersion models, and may allow for additional computation time to be expended on modelling dispersion processes more accurately in future, rather than on accounting for source geometry.

A non-grey analytical model for irradiated atmospheres. II. Analytical vs. numerical solutions

NASA Astrophysics Data System (ADS)

Parmentier, Vivien; Guillot, Tristan; Fortney, Jonathan J.; Marley, Mark S.

2015-02-01

Context. The recent discovery and characterization of the diversity of the atmospheres of exoplanets and brown dwarfs calls for the development of fast and accurate analytical models. Aims: We wish to assess the goodness of the different approximations used to solve the radiative transfer problem in irradiated atmospheres analytically, and we aim to provide a useful tool for a fast computation of analytical temperature profiles that remains correct over a wide range of atmospheric characteristics. Methods: We quantify the accuracy of the analytical solution derived in paper I for an irradiated, non-grey atmosphere by comparing it to a state-of-the-art radiative transfer model. Then, using a grid of numerical models, we calibrate the different coefficients of our analytical model for irradiated solar-composition atmospheres of giant exoplanets and brown dwarfs. Results: We show that the so-called Eddington approximation used to solve the angular dependency of the radiation field leads to relative errors of up to ~5% on the temperature profile. For grey or semi-grey atmospheres (i.e., when the visible and thermal opacities, respectively, can be considered independent of wavelength), we show that the presence of a convective zone has a limited effect on the radiative atmosphere above it and leads to modifications of the radiative temperature profile of approximately ~2%. However, for realistic non-grey planetary atmospheres, the presence of a convective zone that extends to optical depths smaller than unity can lead to changes in the radiative temperature profile on the order of 20% or more. When the convective zone is located at deeper levels (such as for strongly irradiated hot Jupiters), its effect on the radiative atmosphere is again on the same order (~2%) as in the semi-grey case. We show that the temperature inversion induced by a strong absorber in the optical, such as TiO or VO is mainly due to non-grey thermal effects reducing the ability of the upper atmosphere to cool down rather than an enhanced absorption of the stellar light as previously thought. Finally, we provide a functional form for the coefficients of our analytical model for solar-composition giant exoplanets and brown dwarfs. This leads to fully analytical pressure-temperature profiles for irradiated atmospheres with a relative accuracy better than 10% for gravities between 2.5 m s-2 and 250 m s-2 and effective temperatures between 100 K and 3000 K. This is a great improvement over the commonly used Eddington boundary condition. A FORTRAN implementation of the analytical model is only available at the CDS via anonymous ftp to http://cdsarc.u-strasbg.fr (ftp://130.79.128.5) or via http://cdsarc.u-strasbg.fr/viz-bin/qcat?J/A+A/574/A35 or at http://www.oca.eu/parmentier/nongrey.Appendix A is available in electronic form at http://www.aanda.org

On Bayesian Testing of Additive Conjoint Measurement Axioms Using Synthetic Likelihood.

PubMed

Karabatsos, George

2018-06-01

This article introduces a Bayesian method for testing the axioms of additive conjoint measurement. The method is based on an importance sampling algorithm that performs likelihood-free, approximate Bayesian inference using a synthetic likelihood to overcome the analytical intractability of this testing problem. This new method improves upon previous methods because it provides an omnibus test of the entire hierarchy of cancellation axioms, beyond double cancellation. It does so while accounting for the posterior uncertainty that is inherent in the empirical orderings that are implied by these axioms, together. The new method is illustrated through a test of the cancellation axioms on a classic survey data set, and through the analysis of simulated data.

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.

Approximation of optimal filter for Ornstein-Uhlenbeck process with quantised discrete-time observation

NASA Astrophysics Data System (ADS)

Bania, Piotr; Baranowski, Jerzy

2018-02-01

Quantisation of signals is a ubiquitous property of digital processing. In many cases, it introduces significant difficulties in state estimation and in consequence control. Popular approaches either do not address properly the problem of system disturbances or lead to biased estimates. Our intention was to find a method for state estimation for stochastic systems with quantised and discrete observation, that is free of the mentioned drawbacks. We have formulated a general form of the optimal filter derived by a solution of Fokker-Planck equation. We then propose the approximation method based on Galerkin projections. We illustrate the approach for the Ornstein-Uhlenbeck process, and derive analytic formulae for the approximated optimal filter, also extending the results for the variant with control. Operation is illustrated with numerical experiments and compared with classical discrete-continuous Kalman filter. Results of comparison are substantially in favour of our approach, with over 20 times lower mean squared error. The proposed filter is especially effective for signal amplitudes comparable to the quantisation thresholds. Additionally, it was observed that for high order of approximation, state estimate is very close to the true process value. The results open the possibilities of further analysis, especially for more complex processes.

A closed form, physical optics expression for the radar cross section of a perfectly conducting flat plate over a dielectric half-space

NASA Astrophysics Data System (ADS)

Anastassiu, Hristos T.

2003-04-01

The physical optics approximation is employed in the derivation of a closed form expression for the radar cross section (RCS) of a flat, perfectly conducting plate of various shapes, located over a dielectric, possibly lossy half-space. The half-space is assumed to lie in the far field region of the plate. The well-known "four-path model" is invoked in a first-order approximation of the half-space contribution to the scattering mechanisms. Numerical results are compared to a reference, Moment Method solution, and the agreement is investigated, to assess the accuracy of the approximations used. The analytical expressions derived can facilitate very fast RCS calculations for realistic scatterers, such as ships in a sea environment, or aircraft flying low over the ground.

Revisiting the effective medium approximation in all-dielectric subwavelength multilayers: Breakdown and rebuilding

NASA Astrophysics Data System (ADS)

Lei, Xinrui; Mao, Lei; Lu, Yonghua; Wang, Pei

2017-07-01

Here, we present a comprehensive analysis of the effective medium approximation (EMA) breakdown in all-dielectric deep-subwavelength multilayers made of alternating layers by means of the transfer matrix method. We demonstrated that the approximation is invalid at the vicinity of the effective medium's critical angle for total internal reflection and obtained an analytical criterion for the breakdown of the EMA, which depends on the layer thickness, the incident angle, and the permittivity difference between the alternate layers. We rebuilt the EMA by adding higher-order correction onto the traditional effective permittivity. Furthermore, we found that the EMA breakdown that arises from the boundary effect cannot be repaired in the traditional homogenization strategy with only one layer of effective medium. By adding an artificial matched layer after the conventional effective layer, the boundary effect breakdown was neatly removed.

Exact and approximate solutions for transient squeezing flow

NASA Astrophysics Data System (ADS)

Lang, Ji; Santhanam, Sridhar; Wu, Qianhong

2017-10-01

In this paper, we report two novel theoretical approaches to examine a fast-developing flow in a thin fluid gap, which is widely observed in industrial applications and biological systems. The problem is featured by a very small Reynolds number and Strouhal number, making the fluid convective acceleration negligible, while its local acceleration is not. We have developed an exact solution for this problem which shows that the flow starts with an inviscid limit when the viscous effect has no time to appear and is followed by a subsequent developing flow, in which the viscous effect continues to penetrate into the entire fluid gap. An approximate solution is also developed using a boundary layer integral method. This solution precisely captures the general behavior of the transient fluid flow process and agrees very well with the exact solution. We also performed numerical simulation using Ansys-CFX. Excellent agreement between the analytical and the numerical solutions is obtained, indicating the validity of the analytical approaches. The study presented herein fills the gap in the literature and will have a broad impact on industrial and biomedical applications.

Yarkovsky-O'Keefe-Radzievskii-Paddack effect on tumbling objects

NASA Astrophysics Data System (ADS)

Breiter, S.; Rożek, A.; Vokrouhlický, D.

2011-11-01

A semi-analytical model of the Yarkovsky-O'Keefe-Radzievskii-Paddack (YORP) effect on an asteroid spin in a non-principal axis rotation state is developed. The model describes the spin-state evolution in Deprit-Elipe variables, first-order averaged with respect to rotation and Keplerian orbital motion. Assuming zero conductivity, the YORP torque is represented by spherical harmonic series with vectorial coefficients, allowing us to use any degree and order of approximation. Within the quadrupole approximation of the illumination function we find the same first integrals involving rotational momentum, obliquity and dynamical inertia that were obtained by Cicaló & Scheeres. The integrals do not exist when higher degree terms of the illumination function are included, and then the asymptotic states known from Vokrouhlický et al. appear. This resolves an apparent contradiction between earlier results. Averaged equations of motion admit stable and unstable limit cycle solutions that were not previously detected. Non-averaged numerical integration by the Taylor series method for an exemplary shape of 3103 Eger is in good agreement with the semi-analytical theory.

Modeling Three-Dimensional Flow in Confined Aquifers by Superposition of Both Two- and Three-Dimensional Analytic Functions

NASA Astrophysics Data System (ADS)

Haitjema, Henk M.

1985-10-01

A technique is presented to incorporate three-dimensional flow in a Dupuit-Forchheimer model. The method is based on superposition of approximate analytic solutions to both two- and three-dimensional flow features in a confined aquifer of infinite extent. Three-dimensional solutions are used in the domain of interest, while farfield conditions are represented by two-dimensional solutions. Approximate three- dimensional solutions have been derived for a partially penetrating well and a shallow creek. Each of these solutions satisfies the condition that no flow occurs across the confining layers of the aquifer. Because of this condition, the flow at some distance of a three-dimensional feature becomes nearly horizontal. Consequently, remotely from a three-dimensional feature, its three-dimensional solution is replaced by a corresponding two-dimensional one. The latter solution is trivial as compared to its three-dimensional counterpart, and its use greatly enhances the computational efficiency of the model. As an example, the flow is modeled between a partially penetrating well and a shallow creek that occur in a regional aquifer system.

High-beta analytic equilibria in circular, elliptical, and D-shaped large aspect ratio axisymmetric configurations with poloidal and toroidal flows

NASA Astrophysics Data System (ADS)

López, O. E.; Guazzotto, L.

2017-03-01

The Grad-Shafranov-Bernoulli system of equations is a single fluid magnetohydrodynamical description of axisymmetric equilibria with mass flows. Using a variational perturbative approach [E. Hameiri, Phys. Plasmas 20, 024504 (2013)], analytic approximations for high-beta equilibria in circular, elliptical, and D-shaped cross sections in the high aspect ratio approximation are found, which include finite toroidal and poloidal flows. Assuming a polynomial dependence of the free functions on the poloidal flux, the equilibrium problem is reduced to an inhomogeneous Helmholtz partial differential equation (PDE) subject to homogeneous Dirichlet conditions. An application of the Green's function method leads to a closed form for the circular solution and to a series solution in terms of Mathieu functions for the elliptical case, which is valid for arbitrary elongations. To extend the elliptical solution to a D-shaped domain, a boundary perturbation in terms of the triangularity is used. A comparison with the code FLOW [L. Guazzotto et al., Phys. Plasmas 11(2), 604-614 (2004)] is presented for relevant scenarios.

Analytical mass formula and nuclear surface properties in the ETF approximation. Part II: asymmetric nuclei

NASA Astrophysics Data System (ADS)

Aymard, François; Gulminelli, Francesca; Margueron, Jérôme

2016-08-01

We have recently addressed the problem of the determination of the nuclear surface energy for symmetric nuclei in the framework of the extended Thomas-Fermi (ETF) approximation using Skyrme functionals. We presently extend this formalism to the case of asymmetric nuclei and the question of the surface symmetry energy. We propose an approximate expression for the diffuseness and the surface energy. These quantities are analytically related to the parameters of the energy functional. In particular, the influence of the different equation of state parameters can be explicitly quantified. Detailed analyses of the different energy components (local/non-local, isoscalar/isovector, surface/curvature and higher order) are also performed. Our analytical solution of the ETF integral improves previous models and leads to a precision of better than 200 keV per nucleon in the determination of the nuclear binding energy for dripline nuclei.

Chemical Sensor Array Response Modeling Using Quantitative Structure-Activity Relationships Technique

NASA Astrophysics Data System (ADS)

Shevade, Abhijit V.; Ryan, Margaret A.; Homer, Margie L.; Zhou, Hanying; Manfreda, Allison M.; Lara, Liana M.; Yen, Shiao-Pin S.; Jewell, April D.; Manatt, Kenneth S.; Kisor, Adam K.

We have developed a Quantitative Structure-Activity Relationships (QSAR) based approach to correlate the response of chemical sensors in an array with molecular descriptors. A novel molecular descriptor set has been developed; this set combines descriptors of sensing film-analyte interactions, representing sensor response, with a basic analyte descriptor set commonly used in QSAR studies. The descriptors are obtained using a combination of molecular modeling tools and empirical and semi-empirical Quantitative Structure-Property Relationships (QSPR) methods. The sensors under investigation are polymer-carbon sensing films which have been exposed to analyte vapors at parts-per-million (ppm) concentrations; response is measured as change in film resistance. Statistically validated QSAR models have been developed using Genetic Function Approximations (GFA) for a sensor array for a given training data set. The applicability of the sensor response models has been tested by using it to predict the sensor activities for test analytes not considered in the training set for the model development. The validated QSAR sensor response models show good predictive ability. The QSAR approach is a promising computational tool for sensing materials evaluation and selection. It can also be used to predict response of an existing sensing film to new target analytes.

Realistic Analytical Polyhedral MRI Phantoms

PubMed Central

Ngo, Tri M.; Fung, George S. K.; Han, Shuo; Chen, Min; Prince, Jerry L.; Tsui, Benjamin M. W.; McVeigh, Elliot R.; Herzka, Daniel A.

2015-01-01

Purpose Analytical phantoms have closed form Fourier transform expressions and are used to simulate MRI acquisitions. Existing 3D analytical phantoms are unable to accurately model shapes of biomedical interest. It is demonstrated that polyhedral analytical phantoms have closed form Fourier transform expressions and can accurately represent 3D biomedical shapes. Theory The derivations of the Fourier transform of a polygon and polyhedron are presented. Methods The Fourier transform of a polyhedron was implemented and its accuracy in representing faceted and smooth surfaces was characterized. Realistic anthropomorphic polyhedral brain and torso phantoms were constructed and their use in simulated 3D/2D MRI acquisitions was described. Results Using polyhedra, the Fourier transform of faceted shapes can be computed to within machine precision. Smooth surfaces can be approximated with increasing accuracy by increasing the number of facets in the polyhedron; the additional accumulated numerical imprecision of the Fourier transform of polyhedra with many faces remained small. Simulations of 3D/2D brain and 2D torso cine acquisitions produced realistic reconstructions free of high frequency edge aliasing as compared to equivalent voxelized/rasterized phantoms. Conclusion Analytical polyhedral phantoms are easy to construct and can accurately simulate shapes of biomedical interest. PMID:26479724

Rapid, specific determination of iodine and iodide by combined solid-phase extraction/diffuse reflectance spectroscopy

NASA Technical Reports Server (NTRS)

Arena, Matteo P.; Porter, Marc D.; Fritz, James S.

2002-01-01

A new, rapid methodology for trace analysis using solid-phase extraction is described. The two-step methodology is based on the concentration of an analyte onto a membrane disk and on the determination by diffuse reflectance spectroscopy of the amount of analyte extracted on the disk surface. This method, which is adaptable to a wide range of analytes, has been used for monitoring ppm levels of iodine and iodide in spacecraft water. Iodine is used as a biocide in spacecraft water. For these determinations, a water sample is passed through a membrane disk by means of a 10-mL syringe that is attached to a disk holder assembly. The disk, which is a polystyrene-divinylbenzene composite, is impregnated with poly(vinylpyrrolidone) (PVP), which exhaustively concentrates iodine as a yellow iodine-PVP complex. The amount of concentrated iodine is then determined in only 2 s by using a hand-held diffuse reflectance spectrometer by comparing the result with a calibration curve based on the Kubelka-Munk function. The same general procedure can be used to determine iodide levels after its facile and exhaustive oxidation to iodine by peroxymonosulfate (i.e., Oxone reagent). For samples containing both analytes, a two-step procedure can be used in which the iodide concentration is calculated from the difference in iodine levels before and after treatment of the sample with peroxymonosulfate. With this methodology, iodine and iodide levels in the 0.1-5.0 ppm range can be determined with a total workup time of approximately 60 s with a RSD of approximately 6%.

Stability analysis of magnetized neutron stars - a semi-analytic approach

NASA Astrophysics Data System (ADS)

Herbrik, Marlene; Kokkotas, Kostas D.

2017-04-01

We implement a semi-analytic approach for stability analysis, addressing the ongoing uncertainty about stability and structure of neutron star magnetic fields. Applying the energy variational principle, a model system is displaced from its equilibrium state. The related energy density variation is set up analytically, whereas its volume integration is carried out numerically. This facilitates the consideration of more realistic neutron star characteristics within the model compared to analytical treatments. At the same time, our method retains the possibility to yield general information about neutron star magnetic field and composition structures that are likely to be stable. In contrast to numerical studies, classes of parametrized systems can be studied at once, finally constraining realistic configurations for interior neutron star magnetic fields. We apply the stability analysis scheme on polytropic and non-barotropic neutron stars with toroidal, poloidal and mixed fields testing their stability in a Newtonian framework. Furthermore, we provide the analytical scheme for dropping the Cowling approximation in an axisymmetric system and investigate its impact. Our results confirm the instability of simple magnetized neutron star models as well as a stabilization tendency in the case of mixed fields and stratification. These findings agree with analytical studies whose spectrum of model systems we extend by lifting former simplifications.

An Efficient Algorithm for Perturbed Orbit Integration Combining Analytical Continuation and Modified Chebyshev Picard Iteration

NASA Astrophysics Data System (ADS)

Elgohary, T.; Kim, D.; Turner, J.; Junkins, J.

2014-09-01

Several methods exist for integrating the motion in high order gravity fields. Some recent methods use an approximate starting orbit, and an efficient method is needed for generating warm starts that account for specific low order gravity approximations. By introducing two scalar Lagrange-like invariants and employing Leibniz product rule, the perturbed motion is integrated by a novel recursive formulation. The Lagrange-like invariants allow exact arbitrary order time derivatives. Restricting attention to the perturbations due to the zonal harmonics J2 through J6, we illustrate an idea. The recursively generated vector-valued time derivatives for the trajectory are used to develop a continuation series-based solution for propagating position and velocity. Numerical comparisons indicate performance improvements of ~ 70X over existing explicit Runge-Kutta methods while maintaining mm accuracy for the orbit predictions. The Modified Chebyshev Picard Iteration (MCPI) is an iterative path approximation method to solve nonlinear ordinary differential equations. The MCPI utilizes Picard iteration with orthogonal Chebyshev polynomial basis functions to recursively update the states. The key advantages of the MCPI are as follows: 1) Large segments of a trajectory can be approximated by evaluating the forcing function at multiple nodes along the current approximation during each iteration. 2) It can readily handle general gravity perturbations as well as non-conservative forces. 3) Parallel applications are possible. The Picard sequence converges to the solution over large time intervals when the forces are continuous and differentiable. According to the accuracy of the starting solutions, however, the MCPI may require significant number of iterations and function evaluations compared to other integrators. In this work, we provide an efficient methodology to establish good starting solutions from the continuation series method; this warm start improves the performance of the MCPI significantly and will likely be useful for other applications where efficiently computed approximate orbit solutions are needed.

Pair and triplet approximation of a spatial lattice population model with multiscale dispersal using Markov chains for estimating spatial autocorrelation.

PubMed

Hiebeler, David E; Millett, Nicholas E

2011-06-21

We investigate a spatial lattice model of a population employing dispersal to nearest and second-nearest neighbors, as well as long-distance dispersal across the landscape. The model is studied via stochastic spatial simulations, ordinary pair approximation, and triplet approximation. The latter method, which uses the probabilities of state configurations of contiguous blocks of three sites as its state variables, is demonstrated to be greatly superior to pair approximations for estimating spatial correlation information at various scales. Correlations between pairs of sites separated by arbitrary distances are estimated by constructing spatial Markov processes using the information from both approximations. These correlations demonstrate why pair approximation misses basic qualitative features of the model, such as decreasing population density as a large proportion of offspring are dropped on second-nearest neighbors, and why triplet approximation is able to include them. Analytical and numerical results show that, excluding long-distance dispersal, the initial growth rate of an invading population is maximized and the equilibrium population density is also roughly maximized when the population spreads its offspring evenly over nearest and second-nearest neighboring sites. Copyright © 2011 Elsevier Ltd. All rights reserved.

Travelling-wave solutions of a weakly nonlinear two-dimensional higher-order Kadomtsev-Petviashvili dynamical equation for dispersive shallow-water waves

NASA Astrophysics Data System (ADS)

Seadawy, Aly R.

2017-01-01

The propagation of three-dimensional nonlinear irrotational flow of an inviscid and incompressible fluid of the long waves in dispersive shallow-water approximation is analyzed. The problem formulation of the long waves in dispersive shallow-water approximation lead to fifth-order Kadomtsev-Petviashvili (KP) dynamical equation by applying the reductive perturbation theory. By using an extended auxiliary equation method, the solitary travelling-wave solutions of the two-dimensional nonlinear fifth-order KP dynamical equation are derived. An analytical as well as a numerical solution of the two-dimensional nonlinear KP equation are obtained and analyzed with the effects of external pressure flow.

Analysis of the dynamic response of a supersonic inlet to flow-field perturbations upstream of the normal shock

NASA Technical Reports Server (NTRS)

Cole, G. L.; Willoh, R. G.

1975-01-01

A linearized mathematical analysis is presented for determining the response of normal shock position and subsonic duct pressures to flow-field perturbations upstream of the normal shock in mixed-compression supersonic inlets. The inlet duct cross-sectional area variation is approximated by constant-area sections; this approximation results in one-dimensional wave equations. A movable normal shock separates the supersonic and subsonic flow regions, and a choked exit is assumed for the inlet exit condition. The analysis leads to a closed-form matrix solution for the shock position and pressure transfer functions. Analytical frequency response results are compared with experimental data and a method of characteristics solution.

Impact desolvation of electrosprayed microdroplets--a new ionization method for mass spectrometry of large biomolecules.

PubMed

Aksyonov, S A; Williams, P

2001-01-01

Impact desolvation of electrosprayed microdroplets (IDEM) is a new method for producing gas-phase ions of large biomolecules. Analytes are dissolved in an electrolyte solution which is electrosprayed in vacuum, producing highly charged micron and sub-micron sized droplets (microdroplets). These microdroplets are accelerated through potential differences approximately 5 - 10 kV to velocities of several km/s and allowed to impact a target surface. The energetic impacts vaporize the droplets and release desolvated gas-phase ions of the analyte molecules. Oligonucleotides (2- to 12-mer) and peptides (bradykinin, neurotensin) yield singly and doubly charged molecular ions with no detectable fragmentation. Because the extent of multiple charging is significantly less than in atmospheric pressure electrospray ionization, and the method produces ions largely free of adducts from solutions of high ionic strength, IDEM has some promise as a method for coupling to liquid chromatographic techniques and for mixture analysis. Ions are produced in vacuum at a flat equipotential surface, potentially allowing efficient ion extraction. Copyright 2001 John Wiley & Sons, Ltd.

Infinite product expansion of the Fokker-Planck equation with steady-state solution.

PubMed

Martin, R J; Craster, R V; Kearney, M J

2015-07-08

We present an analytical technique for solving Fokker-Planck equations that have a steady-state solution by representing the solution as an infinite product rather than, as usual, an infinite sum. This method has many advantages: automatically ensuring positivity of the resulting approximation, and by design exactly matching both the short- and long-term behaviour. The efficacy of the technique is demonstrated via comparisons with computations of typical examples.

Infinite product expansion of the Fokker–Planck equation with steady-state solution

PubMed Central

Martin, R. J.; Craster, R. V.; Kearney, M. J.

2015-01-01

We present an analytical technique for solving Fokker–Planck equations that have a steady-state solution by representing the solution as an infinite product rather than, as usual, an infinite sum. This method has many advantages: automatically ensuring positivity of the resulting approximation, and by design exactly matching both the short- and long-term behaviour. The efficacy of the technique is demonstrated via comparisons with computations of typical examples. PMID:26346100

Analytical approximations to the dynamics of an array of coupled DC SQUIDs

NASA Astrophysics Data System (ADS)

Berggren, Susan; Palacios, Antonio

2014-04-01

Coupled dynamical systems that operate near the onset of a bifurcation can lead, under certain conditions, to strong signal amplification effects. Over the past years we have studied this generic feature on a wide range of systems, including: magnetic and electric fields sensors, gyroscopic devices, and arrays of loops of superconducting quantum interference devices, also known as SQUIDs. In this work, we consider an array of SQUID loops connected in series as a case study to derive asymptotic analytical approximations to the exact solutions through perturbation analysis. Two approaches are considered. First, a straightforward expansion in which the non-linear parameter related to the inductance of the DC SQUID is treated as the small perturbation parameter. Second, a more accurate procedure that considers the SQUID phase dynamics as non-uniform motion on a circle. This second procedure is readily extended to the series array and it could serve as a mathematical framework to find approximate solutions to related complex systems with high-dimensionality. To the best of our knowledge, an approximate analytical solutions to an array of SQUIDs has not been reported yet in the literature.

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

Rossi, Mariana; Manolopoulos, David E.; Ceriotti, Michele

Two of the most successful methods that are presently available for simulating the quantum dynamics of condensed phase systems are centroid molecular dynamics (CMD) and ring polymer molecular dynamics (RPMD). Despite their conceptual differences, practical implementations of these methods differ in just two respects: the choice of the Parrinello-Rahman mass matrix and whether or not a thermostat is applied to the internal modes of the ring polymer during the dynamics. Here, we explore a method which is halfway between the two approximations: we keep the path integral bead masses equal to the physical particle masses but attach a Langevin thermostatmore » to the internal modes of the ring polymer during the dynamics. We justify this by showing analytically that the inclusion of an internal mode thermostat does not affect any of the established features of RPMD: thermostatted RPMD is equally valid with respect to everything that has actually been proven about the method as RPMD itself. In particular, because of the choice of bead masses, the resulting method is still optimum in the short-time limit, and the transition state approximation to its reaction rate theory remains closely related to the semiclassical instanton approximation in the deep quantum tunneling regime. In effect, there is a continuous family of methods with these properties, parameterised by the strength of the Langevin friction. Here, we explore numerically how the approximation to quantum dynamics depends on this friction, with a particular emphasis on vibrational spectroscopy. We find that a broad range of frictions approaching optimal damping give similar results, and that these results are immune to both the resonance problem of RPMD and the curvature problem of CMD.« less

Response to ``Comment on `Scalings for radiation from plasma bubbles' '' [Phys. Plasmas 18, 034701 (2011)

NASA Astrophysics Data System (ADS)

Thomas, A. G. R.

2011-03-01

In the preceding Comment, Corde, Stordeur, and Malka claim that the trapping threshold derived in my recent paper is incorrect. Their principal argument is that the elliptical orbits I used are not exact solutions of the equation of motion in the fields of the bubble. The original paper never claimed this—rather I claimed that the use of elliptical orbits was a reasonable approximation, which I based on observations from particle-in-cell simulations. Integration of the equation of motion for analytical expressions for idealized bubble fields (either analytically [I. Kostyukov, E. Nerush, A. Pukhov, and V. Seredov, Phys. Rev. Lett. 103, 175003 (2009)] or numerically [S. Corde, A. Stordeur, and V. Malka, "Comment on `Scalings for radiation from plasma bubbles,' " Phys. Plasmas 18, 034701 (2011)]) produces a trapping threshold wholly inconsistent with experiments and full particle-in-cell (PIC) simulations (e.g., requiring an estimated laser intensity of a0˜30 for ne˜1019 cm-3). The inconsistency in the particle trajectories between PIC and the numeric model used by the comment authors arises due to the fact that the analytical fields are only approximately true for "real" plasma bubbles, and lack certain key features of the field structure. Two possible methods of resolution to this inconsistency are either to find ever more complicated but accurate models for the bubble fields or to find approximate solutions to the equations of motion that capture the essential features of the self-consistent electron trajectories. The latter, heuristic approach used in my recent paper produced a threshold that is better matched to experimental observations. In this reply, I will also revisit the problem and examine the relationship between bubble radius and electron momentum at the point of trapping without reference to a particular trajectory.

Super-sample covariance approximations and partial sky coverage

NASA Astrophysics Data System (ADS)

Lacasa, Fabien; Lima, Marcos; Aguena, Michel

2018-04-01

Super-sample covariance (SSC) is the dominant source of statistical error on large scale structure (LSS) observables for both current and future galaxy surveys. In this work, we concentrate on the SSC of cluster counts, also known as sample variance, which is particularly useful for the self-calibration of the cluster observable-mass relation; our approach can similarly be applied to other observables, such as galaxy clustering and lensing shear. We first examined the accuracy of two analytical approximations proposed in the literature for the flat sky limit, finding that they are accurate at the 15% and 30-35% level, respectively, for covariances of counts in the same redshift bin. We then developed a harmonic expansion formalism that allows for the prediction of SSC in an arbitrary survey mask geometry, such as large sky areas of current and future surveys. We show analytically and numerically that this formalism recovers the full sky and flat sky limits present in the literature. We then present an efficient numerical implementation of the formalism, which allows fast and easy runs of covariance predictions when the survey mask is modified. We applied our method to a mask that is broadly similar to the Dark Energy Survey footprint, finding a non-negligible negative cross-z covariance, i.e. redshift bins are anti-correlated. We also examined the case of data removal from holes due to, for example bright stars, quality cuts, or systematic removals, and find that this does not have noticeable effects on the structure of the SSC matrix, only rescaling its amplitude by the effective survey area. These advances enable analytical covariances of LSS observables to be computed for current and future galaxy surveys, which cover large areas of the sky where the flat sky approximation fails.

Model selection on solid ground: Rigorous comparison of nine ways to evaluate Bayesian model evidence

PubMed Central

Schöniger, Anneli; Wöhling, Thomas; Samaniego, Luis; Nowak, Wolfgang

2014-01-01

Bayesian model selection or averaging objectively ranks a number of plausible, competing conceptual models based on Bayes' theorem. It implicitly performs an optimal trade-off between performance in fitting available data and minimum model complexity. The procedure requires determining Bayesian model evidence (BME), which is the likelihood of the observed data integrated over each model's parameter space. The computation of this integral is highly challenging because it is as high-dimensional as the number of model parameters. Three classes of techniques to compute BME are available, each with its own challenges and limitations: (1) Exact and fast analytical solutions are limited by strong assumptions. (2) Numerical evaluation quickly becomes unfeasible for expensive models. (3) Approximations known as information criteria (ICs) such as the AIC, BIC, or KIC (Akaike, Bayesian, or Kashyap information criterion, respectively) yield contradicting results with regard to model ranking. Our study features a theory-based intercomparison of these techniques. We further assess their accuracy in a simplistic synthetic example where for some scenarios an exact analytical solution exists. In more challenging scenarios, we use a brute-force Monte Carlo integration method as reference. We continue this analysis with a real-world application of hydrological model selection. This is a first-time benchmarking of the various methods for BME evaluation against true solutions. Results show that BME values from ICs are often heavily biased and that the choice of approximation method substantially influences the accuracy of model ranking. For reliable model selection, bias-free numerical methods should be preferred over ICs whenever computationally feasible. PMID:25745272

Simultaneous quantification of the boar-taint compounds skatole and androstenone by surface-enhanced Raman scattering (SERS) and multivariate data analysis.

PubMed

Sørensen, Klavs M; Westley, Chloe; Goodacre, Royston; Engelsen, Søren Balling

2015-10-01

This study investigates the feasibility of using surface-enhanced Raman scattering (SERS) for the quantification of absolute levels of the boar-taint compounds skatole and androstenone in porcine fat. By investigation of different types of nanoparticles, pH and aggregating agents, an optimized environment that promotes SERS of the analytes was developed and tested with different multivariate spectral pre-processing techniques, and this was combined with variable selection on a series of analytical standards. The resulting method exhibited prediction errors (root mean square error of cross validation, RMSECV) of 2.4 × 10(-6) M skatole and 1.2 × 10(-7) M androstenone, with a limit of detection corresponding to approximately 2.1 × 10(-11) M for skatole and approximately 1.8 × 10(-10) for androstenone. The method was subsequently tested on porcine fat extract, leading to prediction errors (RMSECV) of 0.17 μg/g for skatole and 1.5 μg/g for androstenone. It is clear that this optimized SERS method, when combined with multivariate analysis, shows great potential for optimization into an on-line application, which will be the first of its kind, and opens up possibilities for simultaneous detection of other meat-quality metabolites or pathogen markers. Graphical abstract Artistic rendering of a laser-illuminated gold colloid sphere with skatole and androstenone adsorbed on the surface.

An orientation measurement method based on Hall-effect sensors for permanent magnet spherical actuators with 3D magnet array.

PubMed

Yan, Liang; Zhu, Bo; Jiao, Zongxia; Chen, Chin-Yin; Chen, I-Ming

2014-10-24

An orientation measurement method based on Hall-effect sensors is proposed for permanent magnet (PM) spherical actuators with three-dimensional (3D) magnet array. As there is no contact between the measurement system and the rotor, this method could effectively avoid friction torque and additional inertial moment existing in conventional approaches. Curved surface fitting method based on exponential approximation is proposed to formulate the magnetic field distribution in 3D space. The comparison with conventional modeling method shows that it helps to improve the model accuracy. The Hall-effect sensors are distributed around the rotor with PM poles to detect the flux density at different points, and thus the rotor orientation can be computed from the measured results and analytical models. Experiments have been conducted on the developed research prototype of the spherical actuator to validate the accuracy of the analytical equations relating the rotor orientation and the value of magnetic flux density. The experimental results show that the proposed method can measure the rotor orientation precisely, and the measurement accuracy could be improved by the novel 3D magnet array. The study result could be used for real-time motion control of PM spherical actuators.

HOST turbine heat transfer program summary

NASA Technical Reports Server (NTRS)

Gladden, Herbert J.; Simoneau, Robert J.

1988-01-01

The objectives of the HOST Turbine Heat Transfer subproject were to obtain a better understanding of the physics of the aerothermodynamic phenomena and to assess and improve the analytical methods used to predict the flow and heat transfer in high temperature gas turbines. At the time the HOST project was initiated, an across-the-board improvement in turbine design technology was needed. A building-block approach was utilized and the research ranged from the study of fundamental phenomena and modeling to experiments in simulated real engine environments. Experimental research accounted for approximately 75 percent of the funding with the remainder going to analytical efforts. A healthy government/industry/university partnership, with industry providing almost half of the research, was created to advance the turbine heat transfer design technology base.

Interacting steps with finite-range interactions: Analytical approximation and numerical results

NASA Astrophysics Data System (ADS)

Jaramillo, Diego Felipe; Téllez, Gabriel; González, Diego Luis; Einstein, T. L.

2013-05-01

We calculate an analytical expression for the terrace-width distribution P(s) for an interacting step system with nearest- and next-nearest-neighbor interactions. Our model is derived by mapping the step system onto a statistically equivalent one-dimensional system of classical particles. The validity of the model is tested with several numerical simulations and experimental results. We explore the effect of the range of interactions q on the functional form of the terrace-width distribution and pair correlation functions. For physically plausible interactions, we find modest changes when next-nearest neighbor 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.

A Simple Method for Deriving the Confidence Regions for the Penalized Cox’s Model via the Minimand Perturbation†

PubMed Central

Lin, Chen-Yen; Halabi, Susan

2017-01-01

We propose a minimand perturbation method to derive the confidence regions for the regularized estimators for the Cox’s proportional hazards model. Although the regularized estimation procedure produces a more stable point estimate, it remains challenging to provide an interval estimator or an analytic variance estimator for the associated point estimate. Based on the sandwich formula, the current variance estimator provides a simple approximation, but its finite sample performance is not entirely satisfactory. Besides, the sandwich formula can only provide variance estimates for the non-zero coefficients. In this article, we present a generic description for the perturbation method and then introduce a computation algorithm using the adaptive least absolute shrinkage and selection operator (LASSO) penalty. Through simulation studies, we demonstrate that our method can better approximate the limiting distribution of the adaptive LASSO estimator and produces more accurate inference compared with the sandwich formula. The simulation results also indicate the possibility of extending the applications to the adaptive elastic-net penalty. We further demonstrate our method using data from a phase III clinical trial in prostate cancer. PMID:29326496

A Simple Method for Deriving the Confidence Regions for the Penalized Cox's Model via the Minimand Perturbation.

PubMed

Lin, Chen-Yen; Halabi, Susan

2017-01-01

We propose a minimand perturbation method to derive the confidence regions for the regularized estimators for the Cox's proportional hazards model. Although the regularized estimation procedure produces a more stable point estimate, it remains challenging to provide an interval estimator or an analytic variance estimator for the associated point estimate. Based on the sandwich formula, the current variance estimator provides a simple approximation, but its finite sample performance is not entirely satisfactory. Besides, the sandwich formula can only provide variance estimates for the non-zero coefficients. In this article, we present a generic description for the perturbation method and then introduce a computation algorithm using the adaptive least absolute shrinkage and selection operator (LASSO) penalty. Through simulation studies, we demonstrate that our method can better approximate the limiting distribution of the adaptive LASSO estimator and produces more accurate inference compared with the sandwich formula. The simulation results also indicate the possibility of extending the applications to the adaptive elastic-net penalty. We further demonstrate our method using data from a phase III clinical trial in prostate cancer.

Nonlinear Schroedinger Approximations for Partial Differential Equations with Quadratic and Quasilinear Terms

NASA Astrophysics Data System (ADS)

Cummings, Patrick

We consider the approximation of solutions of two complicated, physical systems via the nonlinear Schrodinger equation (NLS). In particular, we discuss the evolution of wave packets and long waves in two physical models. Due to the complicated nature of the equations governing many physical systems and the in-depth knowledge we have for solutions of the nonlinear Schrodinger equation, it is advantageous to use approximation results of this kind to model these physical systems. The approximations are simple enough that we can use them to understand the qualitative and quantitative behavior of the solutions, and by justifying them we can show that the behavior of the approximation captures the behavior of solutions to the original equation, at least for long, but finite time. We first consider a model of the water wave equations which can be approximated by wave packets using the NLS equation. We discuss a new proof that both simplifies and strengthens previous justification results of Schneider and Wayne. Rather than using analytic norms, as was done by Schneider and Wayne, we construct a modified energy functional so that the approximation holds for the full interval of existence of the approximate NLS solution as opposed to a subinterval (as is seen in the analytic case). Furthermore, the proof avoids problems associated with inverting the normal form transform by working with a modified energy functional motivated by Craig and Hunter et al. We then consider the Klein-Gordon-Zakharov system and prove a long wave approximation result. In this case there is a non-trivial resonance that cannot be eliminated via a normal form transform. By combining the normal form transform for small Fourier modes and using analytic norms elsewhere, we can get a justification result on the order 1 over epsilon squared time scale.

A modular finite-element model (MODFE) for areal and axisymmetric ground-water-flow problems, Part 2: Derivation of finite-element equations and comparisons with analytical solutions

USGS Publications Warehouse

Cooley, Richard L.

1992-01-01

MODFE, a modular finite-element model for simulating steady- or unsteady-state, area1 or axisymmetric flow of ground water in a heterogeneous anisotropic aquifer is documented in a three-part series of reports. In this report, part 2, the finite-element equations are derived by minimizing a functional of the difference between the true and approximate hydraulic head, which produces equations that are equivalent to those obtained by either classical variational or Galerkin techniques. Spatial finite elements are triangular with linear basis functions, and temporal finite elements are one dimensional with linear basis functions. Physical processes that can be represented by the model include (1) confined flow, unconfined flow (using the Dupuit approximation), or a combination of both; (2) leakage through either rigid or elastic confining units; (3) specified recharge or discharge at points, along lines, or areally; (4) flow across specified-flow, specified-head, or head-dependent boundaries; (5) decrease of aquifer thickness to zero under extreme water-table decline and increase of aquifer thickness from zero as the water table rises; and (6) head-dependent fluxes from springs, drainage wells, leakage across riverbeds or confining units combined with aquifer dewatering, and evapotranspiration. The matrix equations produced by the finite-element method are solved by the direct symmetric-Doolittle method or the iterative modified incomplete-Cholesky conjugate-gradient method. The direct method can be efficient for small- to medium-sized problems (less than about 500 nodes), and the iterative method is generally more efficient for larger-sized problems. Comparison of finite-element solutions with analytical solutions for five example problems demonstrates that the finite-element model can yield accurate solutions to ground-water flow problems.

Ultra-performance liquid chromatography/tandem mass spectrometric quantification of structurally diverse drug mixtures using an ESI-APCI multimode ionization source.

PubMed

Yu, Kate; Di, Li; Kerns, Edward; Li, Susan Q; Alden, Peter; Plumb, Robert S

2007-01-01

We report in this paper an ultra-performance liquid chromatography/tandem mass spectrometric (UPLC(R)/MS/MS) method utilizing an ESI-APCI multimode ionization source to quantify structurally diverse analytes. Eight commercial drugs were used as test compounds. Each LC injection was completed in 1 min using a UPLC system coupled with MS/MS multiple reaction monitoring (MRM) detection. Results from three separate sets of experiments are reported. In the first set of experiments, the eight test compounds were analyzed as a single mixture. The mass spectrometer was switching rapidly among four ionization modes (ESI+, ESI-, APCI-, and APCI+) during an LC run. Approximately 8-10 data points were collected across each LC peak. This was insufficient for a quantitative analysis. In the second set of experiments, four compounds were analyzed as a single mixture. The mass spectrometer was switching rapidly among four ionization modes during an LC run. Approximately 15 data points were obtained for each LC peak. Quantification results were obtained with a limit of detection (LOD) as low as 0.01 ng/mL. For the third set of experiments, the eight test compounds were analyzed as a batch. During each LC injection, a single compound was analyzed. The mass spectrometer was detecting at a particular ionization mode during each LC injection. More than 20 data points were obtained for each LC peak. Quantification results were also obtained. This single-compound analytical method was applied to a microsomal stability test. Compared with a typical HPLC method currently used for the microsomal stability test, the injection-to-injection cycle time was reduced to 1.5 min (UPLC method) from 3.5 min (HPLC method). The microsome stability results were comparable with those obtained by traditional HPLC/MS/MS.

Optimal guidance law development for an advanced launch system

NASA Technical Reports Server (NTRS)

Calise, Anthony J.; Leung, Martin S. K.

1995-01-01

The objective of this research effort was to develop a real-time guidance approach for launch vehicles ascent to orbit injection. Various analytical approaches combined with a variety of model order and model complexity reduction have been investigated. Singular perturbation methods were first attempted and found to be unsatisfactory. The second approach based on regular perturbation analysis was subsequently investigated. It also fails because the aerodynamic effects (ignored in the zero order solution) are too large to be treated as perturbations. Therefore, the study demonstrates that perturbation methods alone (both regular and singular perturbations) are inadequate for use in developing a guidance algorithm for the atmospheric flight phase of a launch vehicle. During a second phase of the research effort, a hybrid analytic/numerical approach was developed and evaluated. The approach combines the numerical methods of collocation and the analytical method of regular perturbations. The concept of choosing intelligent interpolating functions is also introduced. Regular perturbation analysis allows the use of a crude representation for the collocation solution, and intelligent interpolating functions further reduce the number of elements without sacrificing the approximation accuracy. As a result, the combined method forms a powerful tool for solving real-time optimal control problems. Details of the approach are illustrated in a fourth order nonlinear example. The hybrid approach is then applied to the launch vehicle problem. The collocation solution is derived from a bilinear tangent steering law, and results in a guidance solution for the entire flight regime that includes both atmospheric and exoatmospheric flight phases.

Dynamic properties and damping predictions for laminated plates: High order theories - Timoshenko beam

NASA Astrophysics Data System (ADS)

Diveyev, Bohdan; Konyk, Solomija; Crocker, Malcolm J.

2018-01-01

The main aim of this study is to predict the elastic and damping properties of composite laminated plates. This problem has an exact elasticity solution for simple uniform bending and transverse loading conditions. This paper presents a new stress analysis method for the accurate determination of the detailed stress distributions in laminated plates subjected to cylindrical bending. Some approximate methods for the stress state predictions for laminated plates are presented here. The present method is adaptive and does not rely on strong assumptions about the model of the plate. The theoretical model described here incorporates deformations of each sheet of the lamina, which account for the effects of transverse shear deformation, transverse normal strain-stress and nonlinear variation of displacements with respect to the thickness coordinate. Predictions of the dynamic and damping values of laminated plates for various geometrical, mechanical and fastening properties are presented. Comparison with the Timoshenko beam theory is systematically made for analytical and approximation variants.

Fully nonlocal inelastic scattering computations for spectroscopical transmission electron microscopy methods

NASA Astrophysics Data System (ADS)

Rusz, Ján; Lubk, Axel; Spiegelberg, Jakob; Tyutyunnikov, Dmitry

2017-12-01

The complex interplay of elastic and inelastic scattering amenable to different levels of approximation constitutes the major challenge for the computation and hence interpretation of TEM-based spectroscopical methods. The two major approaches to calculate inelastic scattering cross sections of fast electrons on crystals—Yoshioka-equations-based forward propagation and the reciprocal wave method—are founded in two conceptually differing schemes—a numerical forward integration of each inelastically scattered wave function, yielding the exit density matrix, and a computation of inelastic scattering matrix elements using elastically scattered initial and final states (double channeling). Here, we compare both approaches and show that the latter is computationally competitive to the former by exploiting analytical integration schemes over multiple excited states. Moreover, we show how to include full nonlocality of the inelastic scattering event, neglected in the forward propagation approaches, at no additional computing costs in the reciprocal wave method. Detailed simulations show in some cases significant errors due to the z -locality approximation and hence pitfalls in the interpretation of spectroscopical TEM results.

Use of the thin sheath approximation for obtaining ion temperatures from the ISEE 1 limited aperture RPA. [for magnetosphere

NASA Technical Reports Server (NTRS)

Comfort, R. H.; Baugher, C. R.; Chappell, C. R.

1982-01-01

A procedure for analyzing low-energy (less than approximately 100 eV) ion data from the plasma composition experiment on ISEE 1 is set forth. The method is based on a derived analytic expression for particle flux to a limited aperture retarding potential analyzer (RPA) in the thin sheath approximation, which makes allowance for some effects of a charged spacecraft on plasma particle trajectories. Calculations using simulated data are employed in testing the efficacy and accuracy of the technique. On the basis of an analysis of these calculation results and the mathematical model, the method is seen as being able to provide accurate ion temperatures from all good plasmaspheric RPA data. It is noted that corresponding densities and spacecraft potentials should be accurate when spacecraft potentials are negative but that they are subject to error for positive spacecraft potentials, particularly when ion Mach numbers are much less than 1. An analysis of data from a representative ISEE 1 pass produces a plasmasphere temperature profile that is consistent in overall structure with previous observations.

A pertinent approach to solve nonlinear fuzzy integro-differential equations.

PubMed

Narayanamoorthy, S; Sathiyapriya, S P

2016-01-01

Fuzzy integro-differential equations is one of the important parts of fuzzy analysis theory that holds theoretical as well as applicable values in analytical dynamics and so an appropriate computational algorithm to solve them is in essence. In this article, we use parametric forms of fuzzy numbers and suggest an applicable approach for solving nonlinear fuzzy integro-differential equations using homotopy perturbation method. A clear and detailed description of the proposed method is provided. Our main objective is to illustrate that the construction of appropriate convex homotopy in a proper way leads to highly accurate solutions with less computational work. The efficiency of the approximation technique is expressed via stability and convergence analysis so as to guarantee the efficiency and performance of the methodology. Numerical examples are demonstrated to verify the convergence and it reveals the validity of the presented numerical technique. Numerical results are tabulated and examined by comparing the obtained approximate solutions with the known exact solutions. Graphical representations of the exact and acquired approximate fuzzy solutions clarify the accuracy of the approach.

Deriving the exact nonadiabatic quantum propagator in the mapping variable representation.

PubMed

Hele, Timothy J H; Ananth, Nandini

2016-12-22

We derive an exact quantum propagator for nonadiabatic dynamics in multi-state systems using the mapping variable representation, where classical-like Cartesian variables are used to represent both continuous nuclear degrees of freedom and discrete electronic states. The resulting Liouvillian is a Moyal series that, when suitably approximated, can allow for the use of classical dynamics to efficiently model large systems. We demonstrate that different truncations of the exact Liouvillian lead to existing approximate semiclassical and mixed quantum-classical methods and we derive an associated error term for each method. Furthermore, by combining the imaginary-time path-integral representation of the Boltzmann operator with the exact Liouvillian, we obtain an analytic expression for thermal quantum real-time correlation functions. These results provide a rigorous theoretical foundation for the development of accurate and efficient classical-like dynamics to compute observables such as electron transfer reaction rates in complex quantized systems.

Stochastic sampling of quadrature grids for the evaluation of vibrational expectation values

NASA Astrophysics Data System (ADS)

López Ríos, Pablo; Monserrat, Bartomeu; Needs, Richard J.

2018-02-01

The thermal lines method for the evaluation of vibrational expectation values of electronic observables [B. Monserrat, Phys. Rev. B 93, 014302 (2016), 10.1103/PhysRevB.93.014302] was recently proposed as a physically motivated approximation offering balance between the accuracy of direct Monte Carlo integration and the low computational cost of using local quadratic approximations. In this paper we reformulate thermal lines as a stochastic implementation of quadrature-grid integration, analyze the analytical form of its bias, and extend the method to multiple-point quadrature grids applicable to any factorizable harmonic or anharmonic nuclear wave function. The bias incurred by thermal lines is found to depend on the local form of the expectation value, and we demonstrate that the use of finer quadrature grids along selected modes can eliminate this bias, while still offering an ˜30 % lower computational cost than direct Monte Carlo integration in our tests.

Investigation of factors affecting the heater wire method of calibrating fine wire thermocouples

NASA Technical Reports Server (NTRS)

Keshock, E. G.

1972-01-01

An analytical investigation was made of a transient method of calibrating fine wire thermocouples. The system consisted of a 10 mil diameter standard thermocouple (Pt, Pt-13% Rh) and an 0.8 mil diameter chromel-alumel thermocouple attached to a 20 mil diameter electrically heated platinum wire. The calibration procedure consisted of electrically heating the wire to approximately 2500 F within about a seven-second period in an environment approximating atmospheric conditions at 120,000 feet. Rapid periodic readout of the standard and fine wire thermocouple signals permitted a comparison of the two temperature indications. An analysis was performed which indicated that the temperature distortion at the heater wire produced by the thermocouple junctions appears to be of negligible magnitude. Consequently, the calibration technique appears to be basically sound, although several practical changes which appear desirable are presented and discussed. Additional investigation is warranted to evaluate radiation effects and transient response characteristics.

A compatible high-order meshless method for the Stokes equations with applications to suspension flows

NASA Astrophysics Data System (ADS)

Trask, Nathaniel; Maxey, Martin; Hu, Xiaozhe

2018-02-01

A stable numerical solution of the steady Stokes problem requires compatibility between the choice of velocity and pressure approximation that has traditionally proven problematic for meshless methods. In this work, we present a discretization that couples a staggered scheme for pressure approximation with a divergence-free velocity reconstruction to obtain an adaptive, high-order, finite difference-like discretization that can be efficiently solved with conventional algebraic multigrid techniques. We use analytic benchmarks to demonstrate equal-order convergence for both velocity and pressure when solving problems with curvilinear geometries. In order to study problems in dense suspensions, we couple the solution for the flow to the equations of motion for freely suspended particles in an implicit monolithic scheme. The combination of high-order accuracy with fully-implicit schemes allows the accurate resolution of stiff lubrication forces directly from the solution of the Stokes problem without the need to introduce sub-grid lubrication models.

Analytical mass formula and nuclear surface properties in the ETF approximation. Part I: symmetric nuclei

NASA Astrophysics Data System (ADS)

Aymard, François; Gulminelli, Francesca; Margueron, Jérôme

2016-08-01

The problem of determination of nuclear surface energy is addressed within the framework of the extended Thomas Fermi (ETF) approximation using Skyrme functionals. We propose an analytical model for the density profiles with variationally determined diffuseness parameters. In this first paper, we consider the case of symmetric nuclei. In this situation, the ETF functional can be exactly integrated, leading to an analytical formula expressing the surface energy as a function of the couplings of the energy functional. The importance of non-local terms is stressed and it is shown that they cannot be deduced simply from the local part of the functional, as it was suggested in previous works.

Analytical assessment of some characteristic ratios for s-wave superconductors

NASA Astrophysics Data System (ADS)

Gonczarek, Ryszard; Krzyzosiak, Mateusz; Gonczarek, Adam; Jacak, Lucjan

2018-04-01

We evaluate some thermodynamic quantities and characteristic ratios that describe low- and high-temperature s-wave superconducting systems. Based on a set of fundamental equations derived within the conformal transformation method, a simple model is proposed and studied analytically. After including a one-parameter class of fluctuations in the density of states, the mathematical structure of the s-wave superconducting gap, the free energy difference, and the specific heat difference is found and discussed in an analytic manner. Both the zero-temperature limit T = 0 and the subcritical temperature range T ≲ T c are discussed using the method of successive approximations. The equation for the ratio R 1, relating the zero-temperature energy gap and the critical temperature, is formulated and solved numerically for various values of the model parameter. Other thermodynamic quantities are analyzed, including a characteristic ratio R 2, quantifying the dynamics of the specific heat jump at the critical temperature. It is shown that the obtained model results coincide with experimental data for low- T c superconductors. The prospect of application of the presented model in studies of high- T c superconductors and other superconducting systems of the new generation is also discussed.

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.

Analytical solutions for tomato peeling with combined heat flux and convective boundary conditions

NASA Astrophysics Data System (ADS)

Cuccurullo, G.; Giordano, L.; Metallo, A.

2017-11-01

Peeling of tomatoes by radiative heating is a valid alternative to steam or lye, which are expensive and pollutant methods. Suitable energy densities are required in order to realize short time operations, thus involving only a thin layer under the tomato surface. This paper aims to predict the temperature field in rotating tomatoes exposed to the source irradiation. Therefore, a 1D unsteady analytical model is presented, which involves a semi-infinite slab subjected to time dependent heating while convective heat transfer takes place on the exposed surface. In order to account for the tomato rotation, the heat source is described as the positive half-wave of a sinusoidal function. The problem being linear, the solution is derived following the Laplace Transform Method. In addition, an easy-to-handle solution for the problem at hand is presented, which assumes a differentiable function for approximating the source while neglecting convective cooling, the latter contribution turning out to be negligible for the context at hand. A satisfying agreement between the two analytical solutions is found, therefore, an easy procedure for a proper design of the dry heating system can be set up avoiding the use of numerical simulations.

Comparison of bipolar vs. tripolar concentric ring electrode Laplacian estimates.

PubMed

Besio, W; Aakula, R; Dai, W

2004-01-01

Potentials on the body surface from the heart are of a spatial and temporal function. The 12-lead electrocardiogram (ECG) provides useful global temporal assessment, but it yields limited spatial information due to the smoothing effect caused by the volume conductor. The smoothing complicates identification of multiple simultaneous bioelectrical events. In an attempt to circumvent the smoothing problem, some researchers used a five-point method (FPM) to numerically estimate the analytical solution of the Laplacian with an array of monopolar electrodes. The FPM is generalized to develop a bi-polar concentric ring electrode system. We have developed a new Laplacian ECG sensor, a trielectrode sensor, based on a nine-point method (NPM) numerical approximation of the analytical Laplacian. For a comparison, the NPM, FPM and compact NPM were calculated over a 400 x 400 mesh with 1/400 spacing. Tri and bi-electrode sensors were also simulated and their Laplacian estimates were compared against the analytical Laplacian. We found that tri-electrode sensors have a much-improved accuracy with significantly less relative and maximum errors in estimating the Laplacian operator. Apart from the higher accuracy, our new electrode configuration will allow better localization of the electrical activity of the heart than bi-electrode configurations.

Exact analytical solution of shear-induced flexural vibration of functionally graded piezoelectric beam

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

Sharma, Pankaj, E-mail: psharma@rtu.ac.in; Parashar, Sandeep Kumar, E-mail: parashar2@yahoo.com

The priority of this paper is to obtain the exact analytical solution for free flexural vibration of FGPM beam actuated using the d{sub 15} effect. In piezoelectric actuators, the potential use of d{sub 15} effect has been of particular interest for engineering applications since shear piezoelectric coefficient d15 is much higher than the other piezoelectric coupling constants d{sub 31} and d{sub 33}. The applications of shear actuators are to induce and control the flexural vibrations of beams and plates. In this study, a modified Timoshenko beam theory is used where electric potential is assumed to vary sinusoidaly along the thicknessmore » direction. The material properties are assumed to be graded across the thickness in accordance with power law distribution. Hamilton's principle is employed to obtain the equations of motion along with the associated boundary conditions for FGPM beams. Exact analytical solution is derived thus obtained equations of motion. Results for clamped-clamped and clamped-free boundary conditions are presented. The presented result and method shell serve as benchmark for comparing the results obtained from the other approximate methods.« less

Analytic studies of the hard dumbell fluid

NASA Astrophysics Data System (ADS)

Morriss, G. P.; Cummings, P. T.

A closed form analytic theory for the structure of the hard dumbell fluid is introduced and evaluated. It is found to be comparable in accuracy to the reference interaction site approximation (RISA) of Chandler and Andersen.

Aerodynamic design optimization using sensitivity analysis and computational fluid dynamics

NASA Technical Reports Server (NTRS)

Baysal, Oktay; Eleshaky, Mohamed E.

1991-01-01

A new and efficient method is presented for aerodynamic design optimization, which is based on a computational fluid dynamics (CFD)-sensitivity analysis algorithm. The method is applied to design a scramjet-afterbody configuration for an optimized axial thrust. The Euler equations are solved for the inviscid analysis of the flow, which in turn provides the objective function and the constraints. The CFD analysis is then coupled with the optimization procedure that uses a constrained minimization method. The sensitivity coefficients, i.e. gradients of the objective function and the constraints, needed for the optimization are obtained using a quasi-analytical method rather than the traditional brute force method of finite difference approximations. During the one-dimensional search of the optimization procedure, an approximate flow analysis (predicted flow) based on a first-order Taylor series expansion is used to reduce the computational cost. Finally, the sensitivity of the optimum objective function to various design parameters, which are kept constant during the optimization, is computed to predict new optimum solutions. The flow analysis of the demonstrative example are compared with the experimental data. It is shown that the method is more efficient than the traditional methods.

Tank Vapor Characterization Project: Tank 241-S-102 fourth temporal study: Headspace gas and vapor characterization results from samples collected on December 19, 1996

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

Pool, K.H.; Evans, J.C.; Olsen, K.B.

1997-08-01

This report presents the results from analyses of samples taken from the headspace of waste storage tank 241-S-102 (Tank S-102) at the Hanford Site in Washington State. Tank headspace samples collected by SGN Eurisys Service Corporation (SESC) were analyzed by Pacific Northwest National Laboratory (PNNL) to determine headspace concentrations of selected non-radioactive analytes. Analyses were performed by the Vapor Analytical Laboratory (VAL) at PNNL. Vapor concentrations from sorbent trap samples are based on measured sample volumes provided by SESC. Ammonia was determined to be above the immediate notification limit of 150 ppm as specified by the sampling and analysis planmore » (SAP). Hydrogen was the principal flammable constituent of the Tank S-102 headspace, determined to be present at approximately 2.410% of its lower flammability limit (LFL). Total headspace flammability was estimated to be <2.973% of its lower flammability limit (LFL). Total headspace flammability was estimated to be <2.973% of the LFL. Average measured concentrations of targeted gases, inorganic vapors, and selected organic vapors are provided in Table S.1. A summary of experimental methods, including sampling methodology, analytical procedures, and quality assurance and control methods are presented in Section 2.0. Detailed descriptions of the analytical results are provided in Section 3.0.« less

Use of dirichlet distributions and orthogonal projection techniques for the fluctuation analysis of steady-state multivariate birth-death systems

NASA Astrophysics Data System (ADS)

Palombi, Filippo; Toti, Simona

2015-05-01

Approximate weak solutions of the Fokker-Planck equation represent a useful tool to analyze the equilibrium fluctuations of birth-death systems, as they provide a quantitative knowledge lying in between numerical simulations and exact analytic arguments. In this paper, we adapt the general mathematical formalism known as the Ritz-Galerkin method for partial differential equations to the Fokker-Planck equation with time-independent polynomial drift and diffusion coefficients on the simplex. Then, we show how the method works in two examples, namely the binary and multi-state voter models with zealots.

Diagnostics for insufficiencies of posterior calculations in Bayesian signal inference.

PubMed

Dorn, Sebastian; Oppermann, Niels; Ensslin, Torsten A

2013-11-01

We present an error-diagnostic validation method for posterior distributions in Bayesian signal inference, an advancement of a previous work. It transfers deviations from the correct posterior into characteristic deviations from a uniform distribution of a quantity constructed for this purpose. We show that this method is able to reveal and discriminate several kinds of numerical and approximation errors, as well as their impact on the posterior distribution. For this we present four typical analytical examples of posteriors with incorrect variance, skewness, position of the maximum, or normalization. We show further how this test can be applied to multidimensional signals.

A new analytical solar radiation pressure model for current BeiDou satellites: IGGBSPM

PubMed Central

Tan, Bingfeng; Yuan, Yunbin; Zhang, Baocheng; Hsu, Hou Ze; Ou, Jikun

2016-01-01

An analytical solar radiation pressure (SRP) model, IGGBSPM (an abbreviation for Institute of Geodesy and Geophysics BeiDou Solar Pressure Model), has been developed for three BeiDou satellite types, namely, geostationary orbit (GEO), inclined geosynchronous orbit (IGSO) and medium earth orbit (MEO), based on a ray-tracing method. The performance of IGGBSPM was assessed based on numerical integration, SLR residuals and analyses of empirical SRP parameters (except overlap computations). The numerical results show that the integrated orbit resulting from IGGBSPM differs from the precise ephemerides by approximately 5 m and 2 m for GEO and non-GEO satellites, respectively. Moreover, when IGGBSPM is used as an a priori model to enhance the ECOM (5-parameter) model with stochastic pulses, named ECOM + APR, for precise orbit determination, the SLR RMS residual improves by approximately 20–25 percent over the ECOM-only solution during the yaw-steering period and by approximately 40 percent during the yaw-fixed period. For the BeiDou GEO01 satellite, improvements of 18 and 32 percent can be achieved during the out-of-eclipse season and during the eclipse season, respectively. An investigation of the estimated ECOM D0 parameters indicated that the β-angle dependence that is evident in the ECOM-only solution is no longer present in the ECOM + APR solution. PMID:27595795