Assessment of Optimal Flexibility in Ensemble of Frequency Responsive Loads

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

Kundu, Soumya; Hansen, Jacob; Lian, Jianming

2018-04-19

Potential of electrical loads in providing grid ancillary services is often limited due to the uncertainties associated with the load behavior. A knowledge of the expected uncertainties with a load control program would invariably yield to better informed control policies, opening up the possibility of extracting the maximal load control potential without affecting grid operations. In the context of frequency responsive load control, a probabilistic uncertainty analysis framework is presented to quantify the expected error between the target and actual load response, under uncertainties in the load dynamics. A closed-form expression of an optimal demand flexibility, minimizing the expected errormore » in actual and committed flexibility, is provided. Analytical results are validated through Monte Carlo simulations of ensembles of electric water heaters.« less

Unsteady seepage flow over sloping beds in response to multiple localized recharge

NASA Astrophysics Data System (ADS)

Bansal, Rajeev K.

2017-05-01

New generalized solutions of linearized Boussinesq equation are derived to approximate the dynamic behavior of subsurface seepage flow induced by multiple localized time-varying recharges over sloping ditch-drain aquifer system. The mathematical model is based on extended Dupuit-Forchheimer assumption and treats the spatial location of recharge basins as additional parameter. Closed form analytic expressions for spatio-temporal variations in water head distribution and discharge rate into the drains are obtained by solving the governing flow equation using eigenvalue-eigenfunction method. Downward and zero-sloping aquifers are treated as special cases of main results. A numerical example is used for illustration of combined effects of various parameters such as spatial coordinates of the recharge basin, aquifer's bed slope, and recharge rate on the dynamic profiles of phreatic surface.

Model reduction by weighted Component Cost Analysis

NASA Technical Reports Server (NTRS)

Kim, Jae H.; Skelton, Robert E.

1990-01-01

Component Cost Analysis considers any given system driven by a white noise process as an interconnection of different components, and assigns a metric called 'component cost' to each component. These component costs measure the contribution of each component to a predefined quadratic cost function. A reduced-order model of the given system may be obtained by deleting those components that have the smallest component costs. The theory of Component Cost Analysis is extended to include finite-bandwidth colored noises. The results also apply when actuators have dynamics of their own. Closed-form analytical expressions of component costs are also derived for a mechanical system described by its modal data. This is very useful to compute the modal costs of very high order systems. A numerical example for MINIMAST system is presented.

Multi-layer membrane model for mass transport in a direct ethanol fuel cell using an alkaline anion exchange membrane

NASA Astrophysics Data System (ADS)

Bahrami, Hafez; Faghri, Amir

2012-11-01

A one-dimensional, isothermal, single-phase model is presented to investigate the mass transport in a direct ethanol fuel cell incorporating an alkaline anion exchange membrane. The electrochemistry is analytically solved and the closed-form solution is provided for two limiting cases assuming Tafel expressions for both oxygen reduction and ethanol oxidation. A multi-layer membrane model is proposed to properly account for the diffusive and electroosmotic transport of ethanol through the membrane. The fundamental differences in fuel crossover for positive and negative electroosmotic drag coefficients are discussed. It is found that ethanol crossover is significantly reduced upon using an alkaline anion exchange membrane instead of a proton exchange membrane, especially at current densities higher than 500 A m

Stochastic modelling of the hydrologic operation of rainwater harvesting systems

NASA Astrophysics Data System (ADS)

Guo, Rui; Guo, Yiping

2018-07-01

Rainwater harvesting (RWH) systems are an effective low impact development practice that provides both water supply and runoff reduction benefits. A stochastic modelling approach is proposed in this paper to quantify the water supply reliability and stormwater capture efficiency of RWH systems. The input rainfall series is represented as a marked Poisson process and two typical water use patterns are analytically described. The stochastic mass balance equation is solved analytically, and based on this, explicit expressions relating system performance to system characteristics are derived. The performances of a wide variety of RWH systems located in five representative climatic regions of the United States are examined using the newly derived analytical equations. Close agreements between analytical and continuous simulation results are shown for all the compared cases. In addition, an analytical equation is obtained expressing the required storage size as a function of the desired water supply reliability, average water use rate, as well as rainfall and catchment characteristics. The equations developed herein constitute a convenient and effective tool for sizing RWH systems and evaluating their performances.

QCD triple Pomeron coupling from string amplitudes

NASA Astrophysics Data System (ADS)

Bialas, A.; Navelet, H.; Peschanski, R.

1998-06-01

Using the recent solution of the triple Pomeron coupling in the QCD dipole picture as a closed string amplitude with six legs, its analytical form in terms of hypergeometric functions and numerical value are derived.

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

NASA Technical Reports Server (NTRS)

Goldberg, Louis F.

1992-01-01

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

An analytical model with flexible accuracy for deep submicron DCVSL cells

NASA Astrophysics Data System (ADS)

Valiollahi, Sepideh; Ardeshir, Gholamreza

2018-07-01

Differential cascoded voltage switch logic (DCVSL) cells are among the best candidates of circuit designers for a wide range of applications due to advantages such as low input capacitance, high switching speed, small area and noise-immunity; nevertheless, a proper model has not yet been developed to analyse them. This paper analyses deep submicron DCVSL cells based on a flexible accuracy-simplicity trade-off including the following key features: (1) the model is capable of producing closed-form expressions with an acceptable accuracy; (2) model equations can be solved numerically to offer higher accuracy; (3) the short-circuit currents occurring in high-low/low-high transitions are accounted in analysis and (4) the changes in the operating modes of transistors during transitions together with an efficient submicron I-V model, which incorporates the most important non-ideal short-channel effects, are considered. The accuracy of the proposed model is validated in IBM 0.13 µm CMOS technology through comparisons with the accurate physically based BSIM3 model. The maximum error caused by analytical solutions is below 10%, while this amount is below 7% for numerical solutions.

Modeling and experimental characterization of a new piezoelectric sensor for low-amplitude vibration measurement

NASA Astrophysics Data System (ADS)

Hou, X. Y.; Koh, C. G.; Kuang, K. S. C.; Lee, W. H.

2017-07-01

This paper investigates the capability of a novel piezoelectric sensor for low-frequency and low-amplitude vibration measurement. The proposed design effectively amplifies the input acceleration via two amplifying mechanisms and thus eliminates the use of the external charge amplifier or conditioning amplifier typically employed for measurement system. The sensor is also self-powered, i.e. no external power unit is required. Consequently, wiring and electrical insulation for on-site measurement are considerably simpler. In addition, the design also greatly reduces the interference from rotational motion which often accompanies the translational acceleration to be measured. An analytical model is developed based on a set of piezoelectric constitutive equations and beam theory. Closed-form expression is derived to correlate sensor geometry and material properties with its dynamic performance. Experimental calibration is then carried out to validate the analytical model. After calibration, experiments are carried out to check the feasibility of the new sensor in structural vibration detection. From experimental results, it is concluded that the proposed sensor is suitable for measuring low-frequency and low-amplitude vibrations.

Analytic integrable systems: Analytic normalization and embedding flows

NASA Astrophysics Data System (ADS)

Zhang, Xiang

In this paper we mainly study the existence of analytic normalization and the normal form of finite dimensional complete analytic integrable dynamical systems. More details, we will prove that any complete analytic integrable diffeomorphism F(x)=Bx+f(x) in (Cn,0) with B having eigenvalues not modulus 1 and f(x)=O(|) is locally analytically conjugate to its normal form. Meanwhile, we also prove that any complete analytic integrable differential system x˙=Ax+f(x) in (Cn,0) with A having nonzero eigenvalues and f(x)=O(|) is locally analytically conjugate to its normal form. Furthermore we will prove that any complete analytic integrable diffeomorphism defined on an analytic manifold can be embedded in a complete analytic integrable flow. We note that parts of our results are the improvement of Moser's one in J. Moser, The analytic invariants of an area-preserving mapping near a hyperbolic fixed point, Comm. Pure Appl. Math. 9 (1956) 673-692 and of Poincaré's one in H. Poincaré, Sur l'intégration des équations différentielles du premier order et du premier degré, II, Rend. Circ. Mat. Palermo 11 (1897) 193-239. These results also improve the ones in Xiang Zhang, Analytic normalization of analytic integrable systems and the embedding flows, J. Differential Equations 244 (2008) 1080-1092 in the sense that the linear part of the systems can be nonhyperbolic, and the one in N.T. Zung, Convergence versus integrability in Poincaré-Dulac normal form, Math. Res. Lett. 9 (2002) 217-228 in the way that our paper presents the concrete expression of the normal form in a restricted case.

Generalised quasiprobability distribution for Hermite polynomial squeezed states

NASA Astrophysics Data System (ADS)

Datta, Sunil; D'Souza, Richard

1996-02-01

Generalized quasiprobability distributions (QPD) for Hermite polynomial states are presented. These states are solutions of an eigenvalue equation which is quadratic in creation and annihilation operators. Analytical expressions for the QPD are presented for some special cases of the eigenvalues. For large squeezing these analytical expressions for the QPD take the form of a finite series in even Hermite functions. These expressions very transparently exhibit the transition between, P, Q and W functions corresponding to the change of the s-parameter of the QPD. Further, they clearly show the two-photon nature of the processes involved in the generation of these states.

Asymptotic co- and post-seismic displacements in a homogeneous Maxwell sphere

NASA Astrophysics Data System (ADS)

Tang, He; Sun, Wenke

2018-07-01

The deformations of the Earth caused by internal and external forces are usually expressed through Green's functions or the superposition of normal modes, that is, via numerical methods, which are applicable for computing both co- and post-seismic deformations. It is difficult to express these deformations in an analytical form, even for a uniform viscoelastic sphere. In this study, we present a set of asymptotic solutions for computing co- and post-seismic displacements; these solutions can be further applied to solving co- and post-seismic geoid, gravity and strain changes. Expressions are derived for a uniform Maxwell Earth by combining the reciprocity theorem, which links earthquake, tidal, shear and loading deformations, with the asymptotic solutions of these three external forces (tidal, shear and loading) and analytical inverse Laplace transformation formulae. Since the asymptotic solutions are given in a purely analytical form without series summations or extra convergence skills, they can be practically applied in an efficient way, especially when computing post-seismic deformations and glacial isotactic adjustments of the Earth over long timescales.

Asymptotic Co- and Post-seismic displacements in a homogeneous Maxwell sphere

NASA Astrophysics Data System (ADS)

Tang, He; Sun, Wenke

2018-05-01

The deformations of the Earth caused by internal and external forces are usually expressed through Green's functions or the superposition of normal modes, i.e. via numerical methods, which are applicable for computing both co- and post-seismic deformations. It is difficult to express these deformations in an analytical form, even for a uniform viscoelastic sphere. In this study, we present a set of asymptotic solutions for computing co- and post-seismic displacements; these solutions can be further applied to solving co- and post-seismic geoid, gravity, and strain changes. Expressions are derived for a uniform Maxwell Earth by combining the reciprocity theorem, which links earthquake, tidal, shear and loading deformations, with the asymptotic solutions of these three external forces (tidal, shear and loading) and analytical inverse Laplace transformation formulae. Since the asymptotic solutions are given in a purely analytical form without series summations or extra convergence skills, they can be practically applied in an efficient way, especially when computing post-seismic deformations and glacial isotactic adjustments of the Earth over long timescales.

Diagnostics of seeded RF plasmas: An experimental study related to the gaseous core reactor

NASA Technical Reports Server (NTRS)

Thompson, S. D.; Clement, J. D.; Williams, J. R.

1974-01-01

Measurements of the temperature profiles in an RF argon plasma were made over magnetic field intensities ranging from 20 amp turns/cm to 80 amp turns/cm. The results were compared with a one-dimensional numerical treatment of the governing equations and with an approximate closed form analytical solution that neglected radiation losses. The average measured temperatures in the plasma compared well with the numerical treatment, though the experimental profile showed less of an off center temperature peak than predicted by theory. This may be a result of the complex turbulent flow pattern present in the experimental torch and not modeled in the numerical treatment. The radiation term cannot be neglected for argon at the power levels investigated. The closed form analytical approximation that neglected radiation led to temperature predictions on the order of 1000 K to 2000 K higher than measured or predicted by the numerical treatment which considered radiation losses.

An Analytical Study of Prostate-Specific Antigen Dynamics.

PubMed

Esteban, Ernesto P; Deliz, Giovanni; Rivera-Rodriguez, Jaileen; Laureano, Stephanie M

2016-01-01

The purpose of this research is to carry out a quantitative study of prostate-specific antigen dynamics for patients with prostatic diseases, such as benign prostatic hyperplasia (BPH) and localized prostate cancer (LPC). The proposed PSA mathematical model was implemented using clinical data of 218 Japanese patients with histological proven BPH and 147 Japanese patients with LPC (stages T2a and T2b). For prostatic diseases (BPH and LPC) a nonlinear equation was obtained and solved in a close form to predict PSA progression with patients' age. The general solution describes PSA dynamics for patients with both diseases LPC and BPH. Particular solutions allow studying PSA dynamics for patients with BPH or LPC. Analytical solutions have been obtained and solved in a close form to develop nomograms for a better understanding of PSA dynamics in patients with BPH and LPC. This study may be useful to improve the diagnostic and prognosis of prostatic diseases.

A forecast for extinction debt in the presence of speciation.

PubMed

Sgardeli, Vasiliki; Iwasa, Yoh; Varvoglis, Harry; Halley, John M

2017-02-21

Predicting biodiversity relaxation following a disturbance is of great importance to conservation biology. Recently-developed models of stochastic community assembly allow us to predict the evolution of communities on the basis of mechanistic processes at the level of individuals. The neutral model of biodiversity, in particular, has provided closed-form solutions for the relaxation of biodiversity in isolated communities (no immigration or speciation). Here, we extend these results by deriving a relaxation curve for a neutral community in which new species are introduced through the mechanism of random fission speciation (RFS). The solution provides simple closed-form expressions for the equilibrium species richness, the relaxation time and the species-individual curve, which are good approximation to the more complicated formulas existing for the same model. The derivation of the relaxation curve is based on the assumption of a broken-stick species-abundance distribution (SAD) as an initial community configuration; yet for commonly observed SADs, the maximum deviation from the curve does not exceed 10%. Importantly, the solution confirms theoretical results and observations showing that the relaxation time increases with community size and thus habitat area. Such simple and analytically tractable models can help crystallize our ideas on the leading factors affecting biodiversity loss. Copyright © 2016 Elsevier Ltd. All rights reserved.

Power-efficient distributed resource allocation under goodput QoS constraints for heterogeneous networks

NASA Astrophysics Data System (ADS)

Andreotti, Riccardo; Del Fiorentino, Paolo; Giannetti, Filippo; Lottici, Vincenzo

2016-12-01

This work proposes a distributed resource allocation (RA) algorithm for packet bit-interleaved coded OFDM transmissions in the uplink of heterogeneous networks (HetNets), characterized by small cells deployed over a macrocell area and sharing the same band. Every user allocates its transmission resources, i.e., bits per active subcarrier, coding rate, and power per subcarrier, to minimize the power consumption while both guaranteeing a target quality of service (QoS) and accounting for the interference inflicted by other users transmitting over the same band. The QoS consists of the number of information bits delivered in error-free packets per unit of time, or goodput (GP), estimated at the transmitter by resorting to an efficient effective SNR mapping technique. First, the RA problem is solved in the point-to-point case, thus deriving an approximate yet accurate closed-form expression for the power allocation (PA). Then, the interference-limited HetNet case is examined, where the RA problem is described as a non-cooperative game, providing a solution in terms of generalized Nash equilibrium. Thanks to the closed-form of the PA, the solution analysis is based on the best response concept. Hence, sufficient conditions for existence and uniqueness of the solution are analytically derived, along with a distributed algorithm capable of reaching the game equilibrium.

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

Performance Analysis for Joint Target Parameter Estimation in UMTS-Based Passive Multistatic Radar with Antenna Arrays Using Modified Cramér-Rao Lower Bounds.

PubMed

Shi, Chenguang; Wang, Fei; Salous, Sana; Zhou, Jianjiang

2017-10-18

In this study, the modified Cramér-Rao lower bounds (MCRLBs) on the joint estimation of target position and velocity is investigated for a universal mobile telecommunication system (UMTS)-based passive multistatic radar system with antenna arrays. First, we analyze the log-likelihood redfunction of the received signal for a complex Gaussian extended target. Then, due to the non-deterministic transmitted data symbols, the analytically closed-form expressions of the MCRLBs on the Cartesian coordinates of target position and velocity are derived for a multistatic radar system with N t UMTS-based transmit station of L t antenna elements and N r receive stations of L r antenna elements. With the aid of numerical simulations, it is shown that increasing the number of receiving elements in each receive station can reduce the estimation errors. In addition, it is demonstrated that the MCRLB is not only a function of signal-to-noise ratio (SNR), the number of receiving antenna elements and the properties of the transmitted UMTS signals, but also a function of the relative geometric configuration between the target and the multistatic radar system.The analytical expressions for MCRLB will open up a new dimension for passive multistatic radar system by aiding the optimal placement of receive stations to improve the target parameter estimation performance.

Performance Analysis for Joint Target Parameter Estimation in UMTS-Based Passive Multistatic Radar with Antenna Arrays Using Modified Cramér-Rao Lower Bounds

PubMed Central

Wang, Fei; Salous, Sana; Zhou, Jianjiang

2017-01-01

In this study, the modified Cramér-Rao lower bounds (MCRLBs) on the joint estimation of target position and velocity is investigated for a universal mobile telecommunication system (UMTS)-based passive multistatic radar system with antenna arrays. First, we analyze the log-likelihood redfunction of the received signal for a complex Gaussian extended target. Then, due to the non-deterministic transmitted data symbols, the analytically closed-form expressions of the MCRLBs on the Cartesian coordinates of target position and velocity are derived for a multistatic radar system with Nt UMTS-based transmit station of Lt antenna elements and Nr receive stations of Lr antenna elements. With the aid of numerical simulations, it is shown that increasing the number of receiving elements in each receive station can reduce the estimation errors. In addition, it is demonstrated that the MCRLB is not only a function of signal-to-noise ratio (SNR), the number of receiving antenna elements and the properties of the transmitted UMTS signals, but also a function of the relative geometric configuration between the target and the multistatic radar system.The analytical expressions for MCRLB will open up a new dimension for passive multistatic radar system by aiding the optimal placement of receive stations to improve the target parameter estimation performance. PMID:29057805

Advances in the theory of box integrals

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

Bailey, David H.; Borwein, J.M.; Crandall, R.E.

2009-06-25

Box integrals - expectations <|{rvec r}|{sup s}> or <|{rvec r}-{rvec q}|{sup s}> over the unit n-cube (or n-box) - have over three decades been occasionally given closed forms for isolated n,s. By employing experimental mathematics together with a new, global analytic strategy, we prove that for n {le} 4 dimensions the box integrals are for any integer s hypergeometrically closed in a sense we clarify herein. For n = 5 dimensions, we show that a single unresolved integral we call K{sub 5} stands in the way of such hyperclosure proofs. We supply a compendium of exemplary closed forms that naturallymore » arise algorithmically from this theory.« less

Analytical and finite element simulation of a three-bar torsion spring

NASA Astrophysics Data System (ADS)

Rădoi, M.; Cicone, T.

2016-08-01

The present study is dedicated to the innovative 3-bar torsion spring used as suspension solution for the first time at Lunokhod-1, the first autonomous vehicle sent for the exploration of the Moon in the early 70-ies by the former USSR. The paper describes a simple analytical model for calculation of spring static characteristics, taking into account both torsion and bending effects. Closed form solutions of this model allows quick and elegant parametric analysis. A comparison with a single torsion bar with the same stiffness reveal an increase of the maximum stress with more than 50%. A 3D finite element (FE) simulation is proposed to evaluate the accuracy of the analytical model. The model was meshed in an automated pattern (sweep for hubs and tetrahedrons for bars) with mesh morphing. Very close results between analytical and numerical solutions have been found, concluding that the analytical model is accurate. The 3-D finite element simulation was used to evaluate the effects of design details like fillet radius of the bars or contact stresses in the hex hub.

Analytical approximations of the firing rate of an adaptive exponential integrate-and-fire neuron in the presence of synaptic noise.

PubMed

Hertäg, Loreen; Durstewitz, Daniel; Brunel, Nicolas

2014-01-01

Computational models offer a unique tool for understanding the network-dynamical mechanisms which mediate between physiological and biophysical properties, and behavioral function. A traditional challenge in computational neuroscience is, however, that simple neuronal models which can be studied analytically fail to reproduce the diversity of electrophysiological behaviors seen in real neurons, while detailed neuronal models which do reproduce such diversity are intractable analytically and computationally expensive. A number of intermediate models have been proposed whose aim is to capture the diversity of firing behaviors and spike times of real neurons while entailing the simplest possible mathematical description. One such model is the exponential integrate-and-fire neuron with spike rate adaptation (aEIF) which consists of two differential equations for the membrane potential (V) and an adaptation current (w). Despite its simplicity, it can reproduce a wide variety of physiologically observed spiking patterns, can be fit to physiological recordings quantitatively, and, once done so, is able to predict spike times on traces not used for model fitting. Here we compute the steady-state firing rate of aEIF in the presence of Gaussian synaptic noise, using two approaches. The first approach is based on the 2-dimensional Fokker-Planck equation that describes the (V,w)-probability distribution, which is solved using an expansion in the ratio between the time constants of the two variables. The second is based on the firing rate of the EIF model, which is averaged over the distribution of the w variable. These analytically derived closed-form expressions were tested on simulations from a large variety of model cells quantitatively fitted to in vitro electrophysiological recordings from pyramidal cells and interneurons. Theoretical predictions closely agreed with the firing rate of the simulated cells fed with in-vivo-like synaptic noise.

Binegativity of two qubits under noise

NASA Astrophysics Data System (ADS)

Sazim, Sk; Awasthi, Natasha

2018-07-01

Recently, it was argued that the binegativity might be a good quantifier of entanglement for two-qubit states. Like the concurrence and the negativity, the binegativity is also analytically computable quantifier for all two qubits. Based on numerical evidence, it was conjectured that it is a PPT (positive partial transposition) monotone and thus fulfills the criterion to be a good measure of entanglement. In this work, we investigate its behavior under noisy channels which indicate that the binegativity is decreasing monotonically with respect to increasing noise. We also find that the binegativity is closely connected to the negativity and has closed analytical form for arbitrary two qubits. Our study supports the conjecture that the binegativity is a monotone.

Exact solutions for laminated composite cylindrical shells in cylindrical bending

NASA Technical Reports Server (NTRS)

Yuan, F. G.

1992-01-01

Analytic elasticity solutions for laminated composite cylindrical shells under cylindrical bending are presented. The material of the shell is assumed to be general cylindrically anisotropic. Based on the theory of cylindrical anisotropic elasticity, coupled governing partial differential equations are developed. The general expressions for the stresses and displacements in the laminated composite cylinders are discussed. The closed form solutions based on Classical Shell Theory (CST) and Donnell's (1933) theory are also derived for comparison purposes. Three examples illustrate the effect of radius-to-thickness ratio, coupling and stacking sequence. The results show that, in general, CST yields poor stress and displacement distributions for thick-section composite shells, but converges to the exact elasticity solution as the radius-to-thickness ratio increases. It is also shown that Donnell's theory significantly underestimates the stress and displacement response.

A simple method to calculate first-passage time densities with arbitrary initial conditions

NASA Astrophysics Data System (ADS)

Nyberg, Markus; Ambjörnsson, Tobias; Lizana, Ludvig

2016-06-01

Numerous applications all the way from biology and physics to economics depend on the density of first crossings over a boundary. Motivated by the lack of general purpose analytical tools for computing first-passage time densities (FPTDs) for complex problems, we propose a new simple method based on the independent interval approximation (IIA). We generalise previous formulations of the IIA to include arbitrary initial conditions as well as to deal with discrete time and non-smooth continuous time processes. We derive a closed form expression for the FPTD in z and Laplace-transform space to a boundary in one dimension. Two classes of problems are analysed in detail: discrete time symmetric random walks (Markovian) and continuous time Gaussian stationary processes (Markovian and non-Markovian). Our results are in good agreement with Langevin dynamics simulations.

Combined effects of suction/injection and wall surface curvature on natural convection flow in a vertical micro-porous annulus

NASA Astrophysics Data System (ADS)

Jha, B. K.; Aina, B.; Muhammad, S. A.

2015-03-01

This study investigates analytically the hydrodynamic and thermal behaviour of a fully developed natural convection flow in a vertical micro-porous-annulus (MPA) taking into account the velocity slip and temperature jump at the outer surface of inner porous cylinder and inner surface of outer porous cylinder. A closed — form solution is presented for velocity, temperature, volume flow rate, skin friction and rate of heat transfer expressed as a Nusselt number. The influence of each governing parameter on hydrodynamic and thermal behaviour is discussed with the aid of graphs. During the course of investigation, it is found that as suction/injection on the cylinder walls increases, the fluid velocity and temperature is enhanced. In addition, it is observed that wall surface curvature has a significant effect on flow and thermal characteristics.

Axially grooved heat pipe study

NASA Technical Reports Server (NTRS)

1977-01-01

A technology evaluation study on axially grooved heat pipes is presented. The state-of-the-art is reviewed and present and future requirements are identified. Analytical models, the Groove Analysis Program (GAP) and a closed form solution, were developed to facilitate parametric performance evaluations. GAP provides a numerical solution of the differential equations which govern the hydrodynamic flow. The model accounts for liquid recession, liquid/vapor shear interaction, puddle flow as well as laminar and turbulent vapor flow conditions. The closed form solution was developed to reduce computation time and complexity in parametric evaluations. It is applicable to laminar and ideal charge conditions, liquid/vapor shear interaction, and an empirical liquid flow factor which accounts for groove geometry and liquid recession effects. The validity of the closed form solution is verified by comparison with GAP predictions and measured data.

Derivation of three closed loop kinematic velocity models using normalized quaternion feedback for an autonomous redundant manipulator with application to inverse kinematics

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

Unseren, M.A.

1993-04-01

The report discusses the orientation tracking control problem for a kinematically redundant, autonomous manipulator moving in a three dimensional workspace. The orientation error is derived using the normalized quaternion error method of Ickes, the Luh, Walker, and Paul error method, and a method suggested here utilizing the Rodrigues parameters, all of which are expressed in terms of normalized quaternions. The analytical time derivatives of the orientation errors are determined. The latter, along with the translational velocity error, form a dosed loop kinematic velocity model of the manipulator using normalized quaternion and translational position feedback. An analysis of the singularities associatedmore » with expressing the models in a form suitable for solving the inverse kinematics problem is given. Two redundancy resolution algorithms originally developed using an open loop kinematic velocity model of the manipulator are extended to properly take into account the orientation tracking control problem. This report furnishes the necessary mathematical framework required prior to experimental implementation of the orientation tracking control schemes on the seven axis CESARm research manipulator or on the seven-axis Robotics Research K1207i dexterous manipulator, the latter of which is to be delivered to the Oak Ridge National Laboratory in 1993.« less

Neutron Multiplicity: LANL W Covariance Matrix for Curve Fitting

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

Wendelberger, James G.

2016-12-08

In neutron multiplicity counting one may fit a curve by minimizing an objective function, χmore » $$2\\atop{n}$$. The objective function includes the inverse of an n by n matrix of covariances, W. The inverse of the W matrix has a closed form solution. In addition W -1 is a tri-diagonal matrix. The closed form and tridiagonal nature allows for a simpler expression of the objective function χ$$2\\atop{n}$$. Minimization of this simpler expression will provide the optimal parameters for the fitted curve.« less

Analytical model for investigation of interior noise characteristics in aircraft with multiple propellers including synchrophasing

NASA Technical Reports Server (NTRS)

Fuller, C. R.

1986-01-01

A simplified analytical model of transmission of noise into the interior of propeller-driven aircraft has been developed. The analysis includes directivity and relative phase effects of the propeller noise sources, and leads to a closed form solution for the coupled motion between the interior and exterior fields via the shell (fuselage) vibrational response. Various situations commonly encountered in considering sound transmission into aircraft fuselages are investigated analytically and the results obtained are compared to measurements in real aircraft. In general the model has proved successful in identifying basic mechanisms behind noise transmission phenomena.

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.

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

NASA Astrophysics Data System (ADS)

Chen, Wenting; Du, Meiyu; Xu, Xiang

2017-01-01

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

Converging Towards the Optimal Path to Extinction

DTIC Science & Technology

2011-01-01

the reproductive rate R0 should be greater than but very close to 1. However, most real diseases have R0 larger than 1.5, which translates into a...can analytically find an expression for the action along the optimal path. The expression for the action is a function of k and the reproductive number...the optimal path for a range of values of the reproductive number R0. In contrast to the prior two examples, here the action must be computed

Analysis of social optimum for staggered shifts in a single-entry traffic corridor with no late arrivals

NASA Astrophysics Data System (ADS)

Li, Chuan-Yao; Huang, Hai-Jun; Tang, Tie-Qiao

2017-03-01

This paper investigates the traffic flow dynamics under the social optimum (SO) principle in a single-entry traffic corridor with staggered shifts from the analytical and numerical perspectives. The LWR (Lighthill-Whitham and Richards) model and the Greenshield's velocity-density function are utilized to describe the dynamic properties of traffic flow. The closed-form SO solution is analytically derived and some numerical examples are used to further testify the analytical solution. The optimum proportion of the numbers of commuters with different desired arrival times is further discussed, where the analytical and numerical results both indicate that the cumulative outflow curve under the SO principle is piecewise smooth.

A closed-form analytical model for predicting 3D boundary layer displacement thickness for the validation of viscous flow solvers

NASA Astrophysics Data System (ADS)

Kumar, V. R. Sanal; Sankar, Vigneshwaran; Chandrasekaran, Nichith; Saravanan, Vignesh; Natarajan, Vishnu; Padmanabhan, Sathyan; Sukumaran, Ajith; Mani, Sivabalan; Rameshkumar, Tharikaa; Nagaraju Doddi, Hema Sai; Vysaprasad, Krithika; Sharan, Sharad; Murugesh, Pavithra; Shankar, S. Ganesh; Nejaamtheen, Mohammed Niyasdeen; Baskaran, Roshan Vignesh; Rahman Mohamed Rafic, Sulthan Ariff; Harisrinivasan, Ukeshkumar; Srinivasan, Vivek

2018-02-01

A closed-form analytical model is developed for estimating the 3D boundary-layer-displacement thickness of an internal flow system at the Sanal flow choking condition for adiabatic flows obeying the physics of compressible viscous fluids. At this unique condition the boundary-layer blockage induced fluid-throat choking and the adiabatic wall-friction persuaded flow choking occur at a single sonic-fluid-throat location. The beauty and novelty of this model is that without missing the flow physics we could predict the exact boundary-layer blockage of both 2D and 3D cases at the sonic-fluid-throat from the known values of the inlet Mach number, the adiabatic index of the gas and the inlet port diameter of the internal flow system. We found that the 3D blockage factor is 47.33 % lower than the 2D blockage factor with air as the working fluid. We concluded that the exact prediction of the boundary-layer-displacement thickness at the sonic-fluid-throat provides a means to correctly pinpoint the causes of errors of the viscous flow solvers. The methodology presented herein with state-of-the-art will play pivotal roles in future physical and biological sciences for a credible verification, calibration and validation of various viscous flow solvers for high-fidelity 2D/3D numerical simulations of real-world flows. Furthermore, our closed-form analytical model will be useful for the solid and hybrid rocket designers for the grain-port-geometry optimization of new generation single-stage-to-orbit dual-thrust-motors with the highest promising propellant loading density within the given envelope without manifestation of the Sanal flow choking leading to possible shock waves causing catastrophic failures.

Scaling and percolation in the small-world network model

NASA Astrophysics Data System (ADS)

Newman, M. E. J.; Watts, D. J.

1999-12-01

In this paper we study the small-world network model of Watts and Strogatz, which mimics some aspects of the structure of networks of social interactions. We argue that there is one nontrivial length-scale in the model, analogous to the correlation length in other systems, which is well-defined in the limit of infinite system size and which diverges continuously as the randomness in the network tends to zero, giving a normal critical point in this limit. This length-scale governs the crossover from large- to small-world behavior in the model, as well as the number of vertices in a neighborhood of given radius on the network. We derive the value of the single critical exponent controlling behavior in the critical region and the finite size scaling form for the average vertex-vertex distance on the network, and, using series expansion and Padé approximants, find an approximate analytic form for the scaling function. We calculate the effective dimension of small-world graphs and show that this dimension varies as a function of the length-scale on which it is measured, in a manner reminiscent of multifractals. We also study the problem of site percolation on small-world networks as a simple model of disease propagation, and derive an approximate expression for the percolation probability at which a giant component of connected vertices first forms (in epidemiological terms, the point at which an epidemic occurs). The typical cluster radius satisfies the expected finite size scaling form with a cluster size exponent close to that for a random graph. All our analytic results are confirmed by extensive numerical simulations of the model.

Generating survival times to simulate Cox proportional hazards models with time-varying covariates.

PubMed

Austin, Peter C

2012-12-20

Simulations and Monte Carlo methods serve an important role in modern statistical research. They allow for an examination of the performance of statistical procedures in settings in which analytic and mathematical derivations may not be feasible. A key element in any statistical simulation is the existence of an appropriate data-generating process: one must be able to simulate data from a specified statistical model. We describe data-generating processes for the Cox proportional hazards model with time-varying covariates when event times follow an exponential, Weibull, or Gompertz distribution. We consider three types of time-varying covariates: first, a dichotomous time-varying covariate that can change at most once from untreated to treated (e.g., organ transplant); second, a continuous time-varying covariate such as cumulative exposure at a constant dose to radiation or to a pharmaceutical agent used for a chronic condition; third, a dichotomous time-varying covariate with a subject being able to move repeatedly between treatment states (e.g., current compliance or use of a medication). In each setting, we derive closed-form expressions that allow one to simulate survival times so that survival times are related to a vector of fixed or time-invariant covariates and to a single time-varying covariate. We illustrate the utility of our closed-form expressions for simulating event times by using Monte Carlo simulations to estimate the statistical power to detect as statistically significant the effect of different types of binary time-varying covariates. This is compared with the statistical power to detect as statistically significant a binary time-invariant covariate. Copyright © 2012 John Wiley & Sons, Ltd.

Comprehensive Reactive Receiver Modeling for Diffusive Molecular Communication Systems: Reversible Binding, Molecule Degradation, and Finite Number of Receptors.

PubMed

Ahmadzadeh, Arman; Arjmandi, Hamidreza; Burkovski, Andreas; Schober, Robert

2016-10-01

This paper studies the problem of receiver modeling in molecular communication systems. We consider the diffusive molecular communication channel between a transmitter nano-machine and a receiver nano-machine in a fluid environment. The information molecules released by the transmitter nano-machine into the environment can degrade in the channel via a first-order degradation reaction and those that reach the receiver nano-machine can participate in a reversible bimolecular reaction with receiver receptor proteins. Thereby, we distinguish between two scenarios. In the first scenario, we assume that the entire surface of the receiver is covered by receptor molecules. We derive a closed-form analytical expression for the expected received signal at the receiver, i.e., the expected number of activated receptors on the surface of the receiver. Then, in the second scenario, we consider the case where the number of receptor molecules is finite and the uniformly distributed receptor molecules cover the receiver surface only partially. We show that the expected received signal for this scenario can be accurately approximated by the expected received signal for the first scenario after appropriately modifying the forward reaction rate constant. The accuracy of the derived analytical results is verified by Brownian motion particle-based simulations of the considered environment, where we also show the impact of the effect of receptor occupancy on the derived analytical results.

Heterogeneous network epidemics: real-time growth, variance and extinction of infection.

PubMed

Ball, Frank; House, Thomas

2017-09-01

Recent years have seen a large amount of interest in epidemics on networks as a way of representing the complex structure of contacts capable of spreading infections through the modern human population. The configuration model is a popular choice in theoretical studies since it combines the ability to specify the distribution of the number of contacts (degree) with analytical tractability. Here we consider the early real-time behaviour of the Markovian SIR epidemic model on a configuration model network using a multitype branching process. We find closed-form analytic expressions for the mean and variance of the number of infectious individuals as a function of time and the degree of the initially infected individual(s), and write down a system of differential equations for the probability of extinction by time t that are numerically fast compared to Monte Carlo simulation. We show that these quantities are all sensitive to the degree distribution-in particular we confirm that the mean prevalence of infection depends on the first two moments of the degree distribution and the variance in prevalence depends on the first three moments of the degree distribution. In contrast to most existing analytic approaches, the accuracy of these results does not depend on having a large number of infectious individuals, meaning that in the large population limit they would be asymptotically exact even for one initial infectious individual.

Erotic Language as Dramatic Action in Plays by Lyly and Shakespeare

ERIC Educational Resources Information Center

Knoll, Gillian

2012-01-01

This study closely examines the language of desire in the dramatic works of John Lyly and William Shakespeare, and argues that contemplative and analytical speeches about desire function as modes of action in their plays. Erotic speeches do more than express desire in a purely descriptive or perlocutionary capacity distinct from the action of the…

Analytical model of the optical vortex microscope.

PubMed

Płocinniczak, Łukasz; Popiołek-Masajada, Agnieszka; Masajada, Jan; Szatkowski, Mateusz

2016-04-20

This paper presents an analytical model of the optical vortex scanning microscope. In this microscope the Gaussian beam with an embedded optical vortex is focused into the sample plane. Additionally, the optical vortex can be moved inside the beam, which allows fine scanning of the sample. We provide an analytical solution of the whole path of the beam in the system (within paraxial approximation)-from the vortex lens to the observation plane situated on the CCD camera. The calculations are performed step by step from one optical element to the next. We show that at each step, the expression for light complex amplitude has the same form with only four coefficients modified. We also derive a simple expression for the vortex trajectory of small vortex displacements.

Oscillatory flow past a slip cylindrical inclusion embedded in a Brinkman medium

NASA Astrophysics Data System (ADS)

Palaniappan, D.

2015-11-01

Transient flow past a circular cylinder embedded in a porous medium is studied based on Brinkman model with Navier slip conditions. Closed form analytic solution for the stream-function describing slow oscillatory flow around a solid cylindrical inclusion is obtained in the limit of low-Reynolds-number. The key parameters such as the frequency of oscillation λ, the permeability constant δ, and the slip coefficient ξ dictate the flow fields and physical quantities in the entire flow domain. Asymptotic steady-state analysis when δ --> 0 reveals the paradoxical behavior detected by Stokes. Local streamlines for small times demonstrate interesting flow patterns. Rapid transitions including flow separations and eddies are observed far away from the solid inclusion. Analytic expressions for the wall shear stress and the force acting on the cylinder are computed and compared with existing results. It is noted that the slip parameter in the range 0 <= ξ <= 0 . 5 has a significant effect in reducing the stress and force. In the limit of large permeability, Darcy (potential) flow is recovered outside a boundary layer. The results are of some interest in predicting maximum wall stress and pressure drop associated with biological models in fibrous media.

Capillary Flow in Containers of Polygonal Section: Theory and Experiment

NASA Technical Reports Server (NTRS)

Weislogel, Mark M.; Rame, Enrique (Technical Monitor)

2001-01-01

An improved understanding of the large-length-scale capillary flows arising in a low-gravity environment is critical to that engineering community concerned with the design and analysis of spacecraft fluids management systems. Because a significant portion of liquid behavior in spacecraft is capillary dominated it is natural to consider designs that best exploit the spontaneous character of such flows. In the present work, a recently verified asymptotic analysis is extended to approximate spontaneous capillary flows in a large class of cylindrical containers of irregular polygonal section experiencing a step reduction in gravitational acceleration. Drop tower tests are conducted using partially-filled irregular triangular containers for comparison with the theoretical predictions. The degree to which the experimental data agree with the theory is a testament to the robustness of the basic analytical assumption of predominantly parallel flow. As a result, the closed form analytical expressions presented serve as simple, accurate tools for predicting bulk flow characteristics essential to practical low-g system design and analysis. Equations for predicting corner wetting rates, total container flow rates, and transient surfaces shapes are provided that are relevant also to terrestrial applications such as capillary flow in porous media.

Capillary Flow in an Interior Corner

NASA Technical Reports Server (NTRS)

Weislogel, Mark Milton

1996-01-01

The design of fluids management processes in the low-gravity environment of space requires an accurate model and description of capillarity-controlled flow in containers of irregular geometry. Here we consider the capillary rise of a fluid along an interior corner of a container following a rapid reduction in gravity. The analytical portion of the work presents an asymptotic formulation in the limit of a slender fluid column, slight surface curvature along the corner, small inertia, and low gravity. New similarity solutions are found and a list of closed form expressions is provided for flow rate and column length. In particular, it is found that the flow is proportional to t(exp 1/2) for a constant height boundary condition, t(exp 2/5) for a spreading drop, and t(exp 3/5) for constant flow. In the experimental portion of the work, measurements from a 2.2s drop tower are reported. An extensive data set, collected over a previously unexplored range of flow parameters, includes estimates of repeatability and accuracy, the role of inertia and column slenderness, and the effects of corner angle, container geometry, and fluid properties. Comprehensive comparisons are made which illustrate the applicability of the analytic results to low-g fluid systems design.

Analytical Solution for the Aeroelastic Response of a Two-Dimensional Elastic Plate in Axial Flow

NASA Astrophysics Data System (ADS)

Medina, Cory; Kang, Chang-Kwon

2017-11-01

The aeroelastic response of an elastic plate in an unsteady flow describes many engineering problems from bio-locomotion, deforming airfoils, to energy harvesting. However, the analysis is challenging because the shape of the plate is a priori unknown. This study presents an analytical model that can predict the two-way tightly coupled aeroelastic response of a two-dimensional elastic plate including the effects of plate curvature along the flow direction. The plate deforms due to the dynamic balance of wing inertia, elastic restoring force, and aerodynamic force. The coupled model utilizes the linearized Euler-Bernoulli beam theory for the structural model and thin airfoil theory as presented by Theodorsen, which assumes incompressible potential flow, for the aerodynamic model. The coupled equations of motion are solved via Galerkin's method, where closed form solutions for the plate deformation are obtained by deriving the unsteady aerodynamic pressure with respect to the plate normal functions, expressed in a Chebyshev polynomial expansion. Stability analysis is performed for a range of mass ratios obtaining the flutter velocities and corresponding frequencies and the results agree well with the results reported in the literature.

Trapped-mode-induced Fano resonance and acoustical transparency in a one-dimensional solid-fluid phononic crystal

NASA Astrophysics Data System (ADS)

Quotane, Ilyasse; El Boudouti, El Houssaine; Djafari-Rouhani, Bahram

2018-01-01

We investigate theoretically and numerically the possibility of existence of Fano and acoustic-induced transparency (AIT) resonances in a simple though realistic one-dimensional acoustic structure made of solid-fluid layers inserted between two fluids. These resonances are obtained by combining appropriately the zeros of transmission (antiresonance) induced by the solid layers and the local resonances induced by the solid or combined solid-fluid layers with surface free boundary conditions. In particular, we show the possibility of trapped modes, also called bound states in continuum, which have recently found a high renewal interest. These modes appear as resonances with zero width in the transmission spectra as well as in the density of states (DOS). We consider three different structures: (i) a single solid layer inserted between two fluids. This simple structure shows the possibility of existence of trapped modes, which are discrete modes of the solid layer that lie in the continuum modes of the surrounding fluids. We give explicit analytical expressions of the dispersion relation of these eigenmodes of the solid layer which are found independent of the nature of the surrounding fluids. By slightly detuning the angle of incidence from that associated to the trapped mode, we get a well-defined Fano resonance characterized by an asymmetric Fano profile in the transmission spectra. (ii) The second structure consists of a solid-fluid-solid triple layer embedded between two fluids. This structure is found more appropriate to show both Fano and acoustic-induced transparency resonances. We provide detailed analytical expressions for the transmission and reflection coefficients that enable us to deduce a closed-form expression of the dispersion relation giving the trapped modes. Two situations can be distinguished in the triple-layer system: in the case of a symmetric structure (i.e., the same solid layers) we show, by detuning the incidence angle θ , the possibility of existence of Fano resonances that can be fitted following a Fano-type expression. The variation of the Fano parameter that describes the asymmetry of such resonances as well as their width versus θ is studied in detail. In the case of an asymmetric structure (i.e., different solid layers), we show the existence of an incidence angle that enables to squeeze a resonance between two transmission zeros induced by the two solid layers. This resonance behaves like an AIT resonance, its position and width depend on the nature of the fluid and solid layers as well as on the difference between the thicknesses of the solid layers. (iii) In the case of a periodic structure (phononic crystal), we show that trapped modes and Fano resonances give rise, respectively, to dispersionless flat bands with zero group velocity and nearly flat bands with negative or positive group velocities. The analytical results presented here are obtained by means of the Green's function method which enables to deduce in closed form: dispersion curves, transmission and reflection coefficients, DOS, as well as the displacement fields. The proposed solid-fluid layered structures should have important applications for designing acoustic mirrors and acoustic filters as well as supersonic and subsonic materials.

Spline-based high-accuracy piecewise-polynomial phase-to-sinusoid amplitude converters.

PubMed

Petrinović, Davor; Brezović, Marko

2011-04-01

We propose a method for direct digital frequency synthesis (DDS) using a cubic spline piecewise-polynomial model for a phase-to-sinusoid amplitude converter (PSAC). This method offers maximum smoothness of the output signal. Closed-form expressions for the cubic polynomial coefficients are derived in the spectral domain and the performance analysis of the model is given in the time and frequency domains. We derive the closed-form performance bounds of such DDS using conventional metrics: rms and maximum absolute errors (MAE) and maximum spurious free dynamic range (SFDR) measured in the discrete time domain. The main advantages of the proposed PSAC are its simplicity, analytical tractability, and inherent numerical stability for high table resolutions. Detailed guidelines for a fixed-point implementation are given, based on the algebraic analysis of all quantization effects. The results are verified on 81 PSAC configurations with the output resolutions from 5 to 41 bits by using a bit-exact simulation. The VHDL implementation of a high-accuracy DDS based on the proposed PSAC with 28-bit input phase word and 32-bit output value achieves SFDR of its digital output signal between 180 and 207 dB, with a signal-to-noise ratio of 192 dB. Its implementation requires only one 18 kB block RAM and three 18-bit embedded multipliers in a typical field-programmable gate array (FPGA) device. © 2011 IEEE

Nonlinear pressure-flow relationships for passive microfluidic valves.

PubMed

Seker, Erkin; Leslie, Daniel C; Haj-Hariri, Hossein; Landers, James P; Utz, Marcel; Begley, Matthew R

2009-09-21

An analytical solution is presented for the nonlinear pressure-flow relationship of deformable passive valves, which are formed by bonding a deformable film over etched channels separated by a weir. A fluidic pathway connecting the channels is opened when the upstream pressure creates a tunnel along a predefined narrow strip where the film is not bonded to the weir. When the width of the strip is comparable to the inlet channel width, the predicted closed-form pressure-flow rate relationship is in excellent agreement with experiments, which determine pressures by measuring film deflections for prescribed flow rates. The validated closed-form models involve no fitting parameters, and provide the foundation to design passive diodes with specific nonlinear pressure-flow characteristics.

Instability of a planar expansion wave.

PubMed

Velikovich, A L; Zalesak, S T; Metzler, N; Wouchuk, J G

2005-10-01

An expansion wave is produced when an incident shock wave interacts with a surface separating a fluid from a vacuum. Such an interaction starts the feedout process that transfers perturbations from the rippled inner (rear) to the outer (front) surface of a target in inertial confinement fusion. Being essentially a standing sonic wave superimposed on a centered expansion wave, a rippled expansion wave in an ideal gas, like a rippled shock wave, typically produces decaying oscillations of all fluid variables. Its behavior, however, is different at large and small values of the adiabatic exponent gamma. At gamma > 3, the mass modulation amplitude delta(m) in a rippled expansion wave exhibits a power-law growth with time alpha(t)beta, where beta = (gamma - 3)/(gamma - 1). This is the only example of a hydrodynamic instability whose law of growth, dependent on the equation of state, is expressed in a closed analytical form. The growth is shown to be driven by a physical mechanism similar to that of a classical Richtmyer-Meshkov instability. In the opposite extreme gamma - 1 < 1, delta(m) exhibits oscillatory growth, approximately linear with time, until it reaches its peak value approximately (gamma - 1)(-1/2), and then starts to decrease. The mechanism driving the growth is the same as that of Vishniac's instability of a blast wave in a gas with low . Exact analytical expressions for the growth rates are derived for both cases and favorably compared to hydrodynamic simulation results.

Preliminary numerical analysis of improved gas chromatograph model

NASA Technical Reports Server (NTRS)

Woodrow, P. T.

1973-01-01

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

Optical Peaking Enhancement in High-Speed Ring Modulators

PubMed Central

Müller, J.; Merget, F.; Azadeh, S. Sharif; Hauck, J.; García, S. Romero; Shen, B.; Witzens, J.

2014-01-01

Ring resonator modulators (RRM) combine extreme compactness, low power consumption and wavelength division multiplexing functionality, making them a frontrunner for addressing the scalability requirements of short distance optical links. To extend data rates beyond the classically assumed bandwidth capability, we derive and experimentally verify closed form equations of the electro-optic response and asymmetric side band generation resulting from inherent transient time dynamics and leverage these to significantly improve device performance. An equivalent circuit description with a commonly used peaking amplifier model allows straightforward assessment of the effect on existing communication system architectures. A small signal analytical expression of peaking in the electro-optic response of RRMs is derived and used to extend the electro-optic bandwidth of the device above 40 GHz as well as to open eye diagrams penalized by intersymbol interference at 32, 40 and 44 Gbps. Predicted peaking and asymmetric side band generation are in excellent agreement with experiments. PMID:25209255

Enzyme efficiency: An open reaction system perspective

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

Banerjee, Kinshuk, E-mail: kb36@rice.edu; Bhattacharyya, Kamal, E-mail: pchemkb@gmail.com

2015-12-21

A measure of enzyme efficiency is proposed for an open reaction network that, in suitable form, applies to closed systems as well. The idea originates from the description of classical enzyme kinetics in terms of cycles. We derive analytical expressions for the efficiency measure by treating the network not only deterministically but also stochastically. The latter accounts for any significant amount of noise that can be present in biological systems and hence reveals its impact on efficiency. Numerical verification of the results is also performed. It is found that the deterministic equation overestimates the efficiency, the more so for verymore » small system sizes. Roles of various kinetics parameters and system sizes on the efficiency are thoroughly explored and compared with the standard definition k{sub 2}/K{sub M}. Study of substrate fluctuation also indicates an interesting efficiency-accuracy balance.« less

Performance prediction of a synchronization link for distributed aerospace wireless systems.

PubMed

Wang, Wen-Qin; Shao, Huaizong

2013-01-01

For reasons of stealth and other operational advantages, distributed aerospace wireless systems have received much attention in recent years. In a distributed aerospace wireless system, since the transmitter and receiver placed on separated platforms which use independent master oscillators, there is no cancellation of low-frequency phase noise as in the monostatic cases. Thus, high accurate time and frequency synchronization techniques are required for distributed wireless systems. The use of a dedicated synchronization link to quantify and compensate oscillator frequency instability is investigated in this paper. With the mathematical statistical models of phase noise, closed-form analytic expressions for the synchronization link performance are derived. The possible error contributions including oscillator, phase-locked loop, and receiver noise are quantified. The link synchronization performance is predicted by utilizing the knowledge of the statistical models, system error contributions, and sampling considerations. Simulation results show that effective synchronization error compensation can be achieved by using this dedicated synchronization link.

Particle–hole ring diagrams for fermions in two dimensions

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

Kaiser, N., E-mail: nkaiser@ph.tum.de

2014-11-15

The set of particle–hole ring diagrams for a many-fermion system in two dimensions is studied. The complex-valued polarization function is derived in detail and shown to be expressible in terms of square-root functions. For a contact-interaction the perturbative contributions to the energy per particle Ē(k{sub f}) are calculated in a closed analytical form from third up to twelfth order. The resummation of the particle–hole ring diagrams to all orders is studied and a pronounced dependence on the dimensionless coupling parameter α is found. There is a substantial difference between the complete ring-sum with all exchange-type diagrams included and the standardmore » resummation of the leading n-ring diagrams only. The spin factor S{sub n}(g) associated to the nth order ring diagrams is derived for arbitrary spin-degeneracy g.« less

Time-frequency dynamics of superluminal pulse transition to the subluminal regime.

PubMed

Dorrah, Ahmed H; Ramakrishnan, Abhinav; Mojahedi, Mo

2015-03-01

Spectral reshaping and nonuniform phase delay associated with an electromagnetic pulse propagating in a temporally dispersive medium may lead to interesting observations in which the group velocity becomes superluminal or even negative. In such cases, the finite bandwidth of the superluminal region implies the inevitable existence of a cutoff distance beyond which a superluminal pulse becomes subluminal. In this paper, we derive a closed-form analytic expression to estimate this cutoff distance in abnormal dispersive media with gain. Moreover, the method of steepest descent is used to track the time-frequency dynamics associated with the evolution of the center of mass of a superluminal pulse to the subluminal regime. This evolution takes place at longer propagation depths as a result of the subluminal components affecting the behavior of the pulse. Finally, the analysis presents the fundamental limitations of superluminal propagation in light of factors such as the medium depth, pulse width, and the medium dispersion strength.

Current fluctuations in quantum absorption refrigerators

NASA Astrophysics Data System (ADS)

Segal, Dvira

2018-05-01

Absorption refrigerators transfer thermal energy from a cold bath to a hot bath without input power by utilizing heat from an additional "work" reservoir. Particularly interesting is a three-level design for a quantum absorption refrigerator, which can be optimized to reach the maximal (Carnot) cooling efficiency. Previous studies of three-level chillers focused on the behavior of the averaged cooling current. Here, we go beyond that and study the full counting statistics of heat exchange in a three-level chiller model. We explain how to obtain the complete cumulant generating function of the refrigerator in a steady state, then derive a partial cumulant generating function, which yields closed-form expressions for both the averaged cooling current and its noise. Our analytical results and simulations are beneficial for the design of nanoscale engines and cooling systems far from equilibrium, with their performance optimized according to different criteria, efficiency, power, fluctuations, and dissipation.

Power Spectrum of a Noisy System Close to a Heteroclinic Orbit

NASA Astrophysics Data System (ADS)

Giner-Baldó, Jordi; Thomas, Peter J.; Lindner, Benjamin

2017-07-01

We consider a two-dimensional dynamical system that possesses a heteroclinic orbit connecting four saddle points. This system is not able to show self-sustained oscillations on its own. If endowed with white Gaussian noise it displays stochastic oscillations, the frequency and quality factor of which are controlled by the noise intensity. This stochastic oscillation of a nonlinear system with noise is conveniently characterized by the power spectrum of suitable observables. In this paper we explore different analytical and semianalytical ways to compute such power spectra. Besides a number of explicit expressions for the power spectrum, we find scaling relations for the frequency, spectral width, and quality factor of the stochastic heteroclinic oscillator in the limit of weak noise. In particular, the quality factor shows a slow logarithmic increase with decreasing noise of the form Q˜ [ln (1/D)]^2. Our results are compared to numerical simulations of the respective Langevin equations.

A Spherical Harmonic Analysis of the Ooty Wide Field Array (OWFA) Visibility Signal

NASA Astrophysics Data System (ADS)

Chatterjee, Suman; Bharadwaj, Somnath

2018-04-01

Considering redshifted 21-cm intensity mapping with the upcoming OWFA whose field of view subtends ˜57° in the N-S direction, we present a formalism which relates the measured visibilities to the spherical harmonic coefficients of the sky signal. We use this to calculate window functions which relate the two-visibility correlations i.e. the correlation between the visibilities measured at two baselines and two frequencies, to different multipoles of the multi-frequency angular power spectrum Cℓ(ν1, ν2). The formalism here is validated using simulations. We also present approximate closed form analytical expressions which can be used to calculate the window functions. Comparing the widely adopted flat sky approximation, we find that its predictions match those of our spherical harmonic formalism to within 16% across the entire OWFA baseline range. The match improves at large baselines where we have <5% deviations.

Effect of Longitudinal Magnetic Field on Vibration Characteristics of Single-Walled Carbon Nanotubes in a Viscoelastic Medium

NASA Astrophysics Data System (ADS)

Zhang, D. P.; Lei, Y.; Shen, Z. B.

2017-12-01

The effect of longitudinal magnetic field on vibration response of a sing-walled carbon nanotube (SWCNT) embedded in viscoelastic medium is investigated. Based on nonlocal Euler-Bernoulli beam theory, Maxwell's relations, and Kelvin viscoelastic foundation model, the governing equations of motion for vibration analysis are established. The complex natural frequencies and corresponding mode shapes in closed form for the embedded SWCNT with arbitrary boundary conditions are obtained using transfer function method (TFM). The new analytical expressions for the complex natural frequencies are also derived for certain typical boundary conditions and Kelvin-Voigt model. Numerical results from the model are presented to show the effects of nonlocal parameter, viscoelastic parameter, boundary conditions, aspect ratio, and strength of the magnetic field on vibration characteristics for the embedded SWCNT in longitudinal magnetic field. The results demonstrate the efficiency of the proposed methods for vibration analysis of embedded SWCNTs under magnetic field.

Modeling of the reactant conversion rate in a turbulent shear flow

NASA Technical Reports Server (NTRS)

Frankel, S. H.; Madnia, C. K.; Givi, P.

1992-01-01

Results are presented of direct numerical simulations (DNS) of spatially developing shear flows under the influence of infinitely fast chemical reactions of the type A + B yields Products. The simulation results are used to construct the compositional structure of the scalar field in a statistical manner. The results of this statistical analysis indicate that the use of a Beta density for the probability density function (PDF) of an appropriate Shvab-Zeldovich mixture fraction provides a very good estimate of the limiting bounds of the reactant conversion rate within the shear layer. This provides a strong justification for the implementation of this density in practical modeling of non-homogeneous turbulent reacting flows. However, the validity of the model cannot be generalized for predictions of higher order statistical quantities. A closed form analytical expression is presented for predicting the maximum rate of reactant conversion in non-homogeneous reacting turbulence.

Micromechanics analysis of space simulated thermal deformations and stresses in continuous fiber reinforced composites

NASA Technical Reports Server (NTRS)

Bowles, David E.

1990-01-01

Space simulated thermally induced deformations and stresses in continuous fiber reinforced composites were investigated with a micromechanics analysis. The investigation focused on two primary areas. First, available explicit expressions for predicting the effective coefficients of thermal expansion (CTEs) for a composite were compared with each other, and with a finite element (FE) analysis, developed specifically for this study. Analytical comparisons were made for a wide range of fiber/matrix systems, and predicted values were compared with experimental data. The second area of investigation focused on the determination of thermally induced stress fields in the individual constituents. Stresses predicted from the FE analysis were compared to those predicted from a closed-form solution to the composite cylinder (CC) model, for two carbon fiber/epoxy composites. A global-local formulation, combining laminated plate theory and FE analysis, was used to determine the stresses in multidirectional laminates. Thermally induced damage initiation predictions were also made.

Analysis of dual coupler nested coupled cavities.

PubMed

Adib, George A; Sabry, Yasser M; Khalil, Diaa

2017-12-01

Coupled ring resonators are now forming the basic building blocks in several optical systems serving different applications. In many of these applications, a small full width at half maximum is required, along with a large free spectral range. In this work, a configuration of passive coupled cavities constituting dual coupler nested cavities is proposed. A theoretical study of the configuration is presented allowing us to obtain analytical expressions of its different spectral characteristics. The transfer function of the configuration is also used to generate design curves while comparing these results with analytical expressions. Finally, the configuration is compared with other coupled cavity configurations.

Closed-form nonlinear frequency of flexoelectric nanobeams with surface and nonlocal effects under closed circuit electric field

NASA Astrophysics Data System (ADS)

Barati, Mohammad Reza

2018-02-01

Nonlocal and surface effects on nonlinear vibration characteristics of a flexoelectric nanobeams under magnetic field are examined. Eringen’s nonlocal elasticity as well as surface elasticity theories are employed to describe the size-dependency of the flexoelectric nanobeam. Also, flexoelectricity is an important size-dependent phenomena for piezoelectric structures at nanoscale, related to the strain gradient-electric polarization coupling. After the derivation of governing equation via Hamilton’s principle, Galerkin method is employed to satisfy boundary conditions. Also, analytical procedures are implemented to obtain the closed-form nonlinear frequency of flexoelectric nanobeam. It is showed that magnetic field intensity, flexoelectric parameter, nonlocal parameter, elastic foundation and applied voltage on the top surface of the nanobeam have great influences on nonlinear vibration frequency.

Frequency-dependent laminar electroosmotic flow in a closed-end rectangular microchannel.

PubMed

Marcos; Yang, C; Ooi, K T; Wong, T N; Masliyah, J H

2004-07-15

This article presents an analysis of the frequency- and time-dependent electroosmotic flow in a closed-end rectangular microchannel. An exact solution to the modified Navier-Stokes equation governing the ac electroosmotic flow field is obtained by using the Green's function formulation in combination with a complex variable approach. An analytical expression for the induced backpressure gradient is derived. With the Debye-Hückel approximation, the electrical double-layer potential distribution in the channel is obtained by analytically solving the linearized two-dimensional Poisson-Boltzmann equation. Since the counterparts of the flow rate and the electrical current are shown to be linearly proportional to the applied electric field and the pressure gradient, Onsager's principle of reciprocity is demonstrated for transient and ac electroosmotic flows. The time evolution of the electroosmotic flow and the effect of a frequency-dependent ac electric field on the oscillating electroosmotic flow in a closed-end rectangular microchannel are examined. Specifically, the induced pressure gradient is analyzed under effects of the channel dimension and the frequency of electric field. In addition, based on the Stokes second problem, the solution of the slip velocity approximation is presented for comparison with the results obtained from the analytical scheme developed in this study. Copyright 2004 Elsevier Inc.

Analytical solutions to optimal underactuated spacecraft formation reconfiguration

NASA Astrophysics Data System (ADS)

Huang, Xu; Yan, Ye; Zhou, Yang

2015-11-01

Underactuated systems can generally be defined as systems with fewer number of control inputs than that of the degrees of freedom to be controlled. In this paper, analytical solutions to optimal underactuated spacecraft formation reconfiguration without either the radial or the in-track control are derived. By using a linear dynamical model of underactuated spacecraft formation in circular orbits, controllability analysis is conducted for either underactuated case. Indirect optimization methods based on the minimum principle are then introduced to generate analytical solutions to optimal open-loop underactuated reconfiguration problems. Both fixed and free final conditions constraints are considered for either underactuated case and comparisons between these two final conditions indicate that the optimal control strategies with free final conditions require less control efforts than those with the fixed ones. Meanwhile, closed-loop adaptive sliding mode controllers for both underactuated cases are designed to guarantee optimal trajectory tracking in the presence of unmatched external perturbations, linearization errors, and system uncertainties. The adaptation laws are designed via a Lyapunov-based method to ensure the overall stability of the closed-loop system. The explicit expressions of the terminal convergent regions of each system states have also been obtained. Numerical simulations demonstrate the validity and feasibility of the proposed open-loop and closed-loop control schemes for optimal underactuated spacecraft formation reconfiguration in circular orbits.

An improved shear beam method for the characterization of bonded composite joints

NASA Technical Reports Server (NTRS)

Hiel, Clem C.; Brinson, Hal F.

1989-01-01

Closed-form analytical solutions, which govern the displacements and stresses in an adhesive shear beam, are discussed. The remarkable precision with which the shear stresses in the adhesive can be predicted forms the basis of the proposed characterization procedure. The shear modulus of the adhesive is obtained by means of a parameter estimation procedure which requires a symbiosis of theoretical and experimental stress analysis.

Size analysis of polyglutamine protein aggregates using fluorescence detection in an analytical ultracentrifuge.

PubMed

Polling, Saskia; Hatters, Danny M; Mok, Yee-Foong

2013-01-01

Defining the aggregation process of proteins formed by poly-amino acid repeats in cells remains a challenging task due to a lack of robust techniques for their isolation and quantitation. Sedimentation velocity methodology using fluorescence detected analytical ultracentrifugation is one approach that can offer significant insight into aggregation formation and kinetics. While this technique has traditionally been used with purified proteins, it is now possible for substantial information to be collected with studies using cell lysates expressing a GFP-tagged protein of interest. In this chapter, we describe protocols for sample preparation and setting up the fluorescence detection system in an analytical ultracentrifuge to perform sedimentation velocity experiments on cell lysates containing aggregates formed by poly-amino acid repeat proteins.

Modeling convection-diffusion-reaction systems for microfluidic molecular communications with surface-based receivers in Internet of Bio-Nano Things

PubMed Central

Akan, Ozgur B.

2018-01-01

We consider a microfluidic molecular communication (MC) system, where the concentration-encoded molecular messages are transported via fluid flow-induced convection and diffusion, and detected by a surface-based MC receiver with ligand receptors placed at the bottom of the microfluidic channel. The overall system is a convection-diffusion-reaction system that can only be solved by numerical methods, e.g., finite element analysis (FEA). However, analytical models are key for the information and communication technology (ICT), as they enable an optimisation framework to develop advanced communication techniques, such as optimum detection methods and reliable transmission schemes. In this direction, we develop an analytical model to approximate the expected time course of bound receptor concentration, i.e., the received signal used to decode the transmitted messages. The model obviates the need for computationally expensive numerical methods by capturing the nonlinearities caused by laminar flow resulting in parabolic velocity profile, and finite number of ligand receptors leading to receiver saturation. The model also captures the effects of reactive surface depletion layer resulting from the mass transport limitations and moving reaction boundary originated from the passage of finite-duration molecular concentration pulse over the receiver surface. Based on the proposed model, we derive closed form analytical expressions that approximate the received pulse width, pulse delay and pulse amplitude, which can be used to optimize the system from an ICT perspective. We evaluate the accuracy of the proposed model by comparing model-based analytical results to the numerical results obtained by solving the exact system model with COMSOL Multiphysics. PMID:29415019

The {sech}( {\\hat{ξ }} ) -Type Profiles: A Swiss-Army Knife for Exact Analytical Modeling of Thermal Diffusion and Wave Propagation in Graded Media

NASA Astrophysics Data System (ADS)

Krapez, J.-C.

2018-07-01

This work deals with the exact analytical modeling of transfer phenomena in heterogeneous materials exhibiting one-dimensional continuous variations of their properties. Regarding heat transfer, it has recently been shown that by applying a Liouville transformation and multiple Darboux transformations, infinite sequences of solvable profiles of thermal effusivity can be constructed together with the associated temperature (exact) solutions, all in closed-form expressions (vs. the diffusion-time variable and with a growing number of parameters). In addition, a particular class of profiles, the so-called {sech}( {\\hat{ξ }} ) -type profiles, exhibit high agility and at the same time parsimony. In this paper we delve further into the description of these solvable profiles and their properties. Most importantly, their quadrupole formulation is provided, enabling smooth synthetic profiles of effusivity of arbitrary complexity to be built, and allowing the corresponding temperature dynamic response to be obtained very easily thereafter. Examples are given with increasing variability of the effusivity and an increasing number of elementary profiles. These highly flexible profiles are equally relevant to providing an exact analytical solution to wave propagation problems in 1D graded media (i.e., Maxwell's equations, the acoustic equation, the telegraph equation, etc.). From now on, whether it be for diffusion-like or wave-like problems, when the leading properties present (possibly piecewise-) continuously heterogeneous profiles, the classical staircase model can be advantageously replaced by a "high-level" quadrupole model consisting of one or more {sech}( {\\hat{ξ }} ) -type profiles, which makes the latter a true Swiss-Army knife for analytical modeling.

Modeling convection-diffusion-reaction systems for microfluidic molecular communications with surface-based receivers in Internet of Bio-Nano Things.

PubMed

Kuscu, Murat; Akan, Ozgur B

2018-01-01

We consider a microfluidic molecular communication (MC) system, where the concentration-encoded molecular messages are transported via fluid flow-induced convection and diffusion, and detected by a surface-based MC receiver with ligand receptors placed at the bottom of the microfluidic channel. The overall system is a convection-diffusion-reaction system that can only be solved by numerical methods, e.g., finite element analysis (FEA). However, analytical models are key for the information and communication technology (ICT), as they enable an optimisation framework to develop advanced communication techniques, such as optimum detection methods and reliable transmission schemes. In this direction, we develop an analytical model to approximate the expected time course of bound receptor concentration, i.e., the received signal used to decode the transmitted messages. The model obviates the need for computationally expensive numerical methods by capturing the nonlinearities caused by laminar flow resulting in parabolic velocity profile, and finite number of ligand receptors leading to receiver saturation. The model also captures the effects of reactive surface depletion layer resulting from the mass transport limitations and moving reaction boundary originated from the passage of finite-duration molecular concentration pulse over the receiver surface. Based on the proposed model, we derive closed form analytical expressions that approximate the received pulse width, pulse delay and pulse amplitude, which can be used to optimize the system from an ICT perspective. We evaluate the accuracy of the proposed model by comparing model-based analytical results to the numerical results obtained by solving the exact system model with COMSOL Multiphysics.

Anomaly formulas for the complex-valued analytic torsion on compact bordisms

PubMed Central

Maldonado Molina, Osmar

2013-01-01

We extend the complex-valued analytic torsion, introduced by Burghelea and Haller on closed manifolds, to compact Riemannian bordisms. We do so by considering a flat complex vector bundle over a compact Riemannian manifold, endowed with a fiberwise nondegenerate symmetric bilinear form. The Riemmanian metric and the bilinear form are used to define non-selfadjoint Laplacians acting on vector-valued smooth forms under absolute and relative boundary conditions. In order to define the complex-valued analytic torsion in this situation, we study spectral properties of these generalized Laplacians. Then, as main results, we obtain so-called anomaly formulas for this torsion. Our reasoning takes into account that the coefficients in the heat trace asymptotic expansion associated to the boundary value problem under consideration, are locally computable. The anomaly formulas for the complex-valued Ray–Singer torsion are derived first by using the corresponding ones for the Ray–Singer metric, obtained by Brüning and Ma on manifolds with boundary, and then an argument of analytic continuation. In odd dimensions, our anomaly formulas are in accord with the corresponding results of Su, without requiring the variations of the Riemannian metric and bilinear structures to be supported in the interior of the manifold. PMID:27087744

Polarization effects in recoil-induced resonances

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

Lazebnyi, D. B., E-mail: becks.ddf@gmail.com; Brazhnikov, D. V.; Taichenachev, A. V.

2017-01-15

The effect of the field polarization on the amplitude of recoil-induced resonances (RIRs) is considered for laser-cooled free atoms and for atoms in a working magneto-optical trap (MOT). For all closed dipole transitions, explicit analytical expressions are obtained for the polarization dependence of the resonance amplitudes within a perturbation theory. Optimal polarization conditions are found for the observation of resonances.

What Makes a Senior Thesis Good?

ERIC Educational Resources Information Center

Trosset, Carol; Weisler, Steven

2018-01-01

Kuh (2008) describes the capstone as a "culminating experience" students undertake close to graduation often involving "a project of some sort that integrates and applies what they've learned" (p. 11). The senior thesis is one form of the capstone in which students write an analytic paper under faculty supervision, typically as…

Roy-Steiner equations for pion-nucleon scattering

NASA Astrophysics Data System (ADS)

Ditsche, C.; Hoferichter, M.; Kubis, B.; Meißner, U.-G.

2012-06-01

Starting from hyperbolic dispersion relations, we derive a closed system of Roy-Steiner equations for pion-nucleon scattering that respects analyticity, unitarity, and crossing symmetry. We work out analytically all kernel functions and unitarity relations required for the lowest partial waves. In order to suppress the dependence on the high energy regime we also consider once- and twice-subtracted versions of the equations, where we identify the subtraction constants with subthreshold parameters. Assuming Mandelstam analyticity we determine the maximal range of validity of these equations. As a first step towards the solution of the full system we cast the equations for the π π to overline N N partial waves into the form of a Muskhelishvili-Omnès problem with finite matching point, which we solve numerically in the single-channel approximation. We investigate in detail the role of individual contributions to our solutions and discuss some consequences for the spectral functions of the nucleon electromagnetic form factors.

Order of wetting transitions in electrolyte solutions.

PubMed

Ibagon, Ingrid; Bier, Markus; Dietrich, S

2014-05-07

For wetting films in dilute electrolyte solutions close to charged walls we present analytic expressions for their effective interface potentials. The analysis of these expressions renders the conditions under which corresponding wetting transitions can be first- or second-order. Within mean field theory we consider two models, one with short- and one with long-ranged solvent-solvent and solvent-wall interactions. The analytic results reveal in a transparent way that wetting transitions in electrolyte solutions, which occur far away from their critical point (i.e., the bulk correlation length is less than half of the Debye length) are always first-order if the solvent-solvent and solvent-wall interactions are short-ranged. In contrast, wetting transitions close to the bulk critical point of the solvent (i.e., the bulk correlation length is larger than the Debye length) exhibit the same wetting behavior as the pure, i.e., salt-free, solvent. If the salt-free solvent is governed by long-ranged solvent-solvent as well as long-ranged solvent-wall interactions and exhibits critical wetting, adding salt can cause the occurrence of an ion-induced first-order thin-thick transition which precedes the subsequent continuous wetting as for the salt-free solvent.

Perturbations of a close-earth satellite due to sunlight diffusely reflected from the earth. I - Uniform albedo

NASA Technical Reports Server (NTRS)

Lautman, D. A.

1977-01-01

A semianalytic method has been developed to calculate the radiation-pressure perturbations of a close-earth satellite due to sunlight reflected from the earth. The assumptions made are that the satellite is spherically symmetric and that the solar radiation is reflected from the earth according to Lambert's Law with uniform albedo. By using expressions for the components of the radiation-pressure force due to Lochry, the expressions for the perturbations of the elements were developed into series in the true anomaly. The perturbations within a given revolution can be obtained analytically by integrating with respect to v while holding all slowly varying quantities constant. The long-range perturbations are then obtained by accumulating the net perturbations at the end of each revolution.

Part-2: Analytical Expressions of Concentrations of Glucose, Oxygen, and Gluconic Acid in a Composite Membrane for Closed-Loop Insulin Delivery for the Non-steady State Conditions.

PubMed

Mehala, N; Rajendran, L; Meena, V

2017-02-01

A mathematical model developed by Abdekhodaie and Wu (J Membr Sci 335:21-31, 2009), which describes a dynamic process involving an enzymatic reaction and diffusion of reactants and product inside glucose-sensitive composite membrane has been discussed. This theoretical model depicts a system of non-linear non-steady state reaction diffusion equations. These equations have been solved using new approach of homotopy perturbation method and analytical solutions pertaining to the concentrations of glucose, oxygen, and gluconic acid are derived. These analytical results are compared with the numerical results, and limiting case results for steady state conditions and a good agreement is observed. The influence of various kinetic parameters involved in the model has been presented graphically. Theoretical evaluation of the kinetic parameters like the maximal reaction velocity (V max ) and Michaelis-Menten constants for glucose and oxygen (K g and K ox ) is also reported. This predicted model is very much useful for designing the glucose-responsive composite membranes for closed-loop insulin delivery.

Theoretical study of production of unique glasses in space

NASA Technical Reports Server (NTRS)

Larsen, D. C.

1974-01-01

Analytical functional relationships describing homogeneous nucleation and crystallization in various supercooled liquids were developed. The time and temperature dependent relationships of nucleation and crystallization (intrinsic properties) are being used to relate glass forming tendency to extrinsic parameters such as cooling rate through computer simulation. Single oxide systems are being studied initially to aid in developing workable kinetic models and to indicate the primary materials parameters affecting glass formation. The theory and analytical expressions developed for simple systems is then extended to complex oxide systems. A thorough understanding of nucleation and crystallization kinetics of glass forming systems provides a priori knowledge of the ability of a given system to form a glass.

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

The HVT technique and the 'uncertainty' relation for central potentials

NASA Astrophysics Data System (ADS)

Grypeos, M. E.; Koutroulos, C. G.; Oyewumi, K. J.; Petridou, Th

2004-08-01

The quantum mechanical hypervirial theorems (HVT) technique is used to treat the so-called 'uncertainty' relation for quite a general class of central potential wells, including the (reduced) Poeschl-Teller and the Gaussian one. It is shown that this technique is quite suitable in deriving an approximate analytic expression in the form of a truncated power series expansion for the dimensionless product Pnl equiv langr2rangnllangp2rangnl/planck2, for every (deeply) bound state of a particle moving non-relativistically in the well, provided that a (dimensionless) parameter s is sufficiently small. Attention is also paid to a number of cases, among the limited existing ones, in which exact analytic or semi-analytic expressions for Pnl can be derived. Finally, numerical results are given and discussed.

Fusion Molecules of Heat Shock Protein HSPX with Other Antigens of Mycobacterium tuberculosis Show High Potential in Serodiagnosis of Tuberculosis.

PubMed

Khalid, Ruqyya; Afzal, Madeeha; Khurshid, Sana; Paracha, Rehan Zafar; Khan, Imran H; Akhtar, Muhammad Waheed

Variable individual response against the antigens of Mycobacterium tuberculosis necessitates detection of multiple antibodies for enhancing reliability of serodiagnosis of tuberculosis. Fusion molecules consisting of two or more antigens showing high sensitivity would be helpful in achieving this objective. Antigens of M. tuberculosis HSPX and PE35 were expressed in a soluble form whereas tnPstS1 and FbpC1 were expressed as inclusion bodies at 37°C. Heat shock protein HSPX when attached to the N-termini of the antigens PE35, tnPstS1 and FbpC1, all the fusion molecules were expressed at high levels in E. coli in a soluble form. ELISA analysis of the plasma samples of TB patients against HSPX-tnPstS1 showed 57.7% sensitivity which is nearly the same as the expected combined value obtained after deducting the number of plasma samples (32) containing the antibodies against both the individual antigens. Likewise, the 54.4% sensitivity of HSPX-PE35 was nearly the same as that expected from the combined values of the contributing antigens. Structural analysis of all the fusion molecules by CD spectroscopy showed that α-helical and β-sheet contents were found close to those obtained through molecular modeling. Molecular modeling studies of HSPX-tnPstS1 and HSPX-PE35 support the analytical results as most of the epitopes of the contributing antigens were found to be available for binding to the corresponding antibodies. Using these fusion molecules in combination with other antigenic molecules should reduce the number of antigenic proteins required for a more reliable and economical serodiagnosis of tuberculosis. Also, HSPX seems to have potential application in soluble expression of heterologous proteins in E. coli.

Exact partition functions for deformed N=2 theories with N_f=4 flavours

NASA Astrophysics Data System (ADS)

Beccaria, Matteo; Fachechi, Alberto; Macorini, Guido; Martina, Luigi

2016-12-01

We consider the Ω-deformed N=2 SU(2) gauge theory in four dimensions with N f = 4 massive fundamental hypermultiplets. The low energy effective action depends on the deformation parameters ɛ 1 , ɛ 2, the scalar field expectation value a, and the hypermultiplet masses m = ( m 1 , m 2 , m 3 , m 4). Motivated by recent findings in the N={2}^{*} theory, we explore the theories that are characterized by special fixed ratios ɛ 2 /ɛ 1 and m /ɛ 1 and propose a simple condition on the structure of the multi-instanton contributions to the prepotential determining the effective action. This condition determines a finite set Π N of special points such that the prepotential has N poles at fixed positions independent on the instanton number. In analogy with what happens in the N={2}^{*} gauge theory, the full prepotential of the Π N theories may be given in closed form as an explicit function of a and the modular parameter q appearing in special combinations of Eisenstein series and Jacobi theta functions with well defined modular properties. The resulting finite pole partition functions are related by AGT correspondence to special 4-point spherical conformal blocks of the Virasoro algebra. We examine in full details special cases where the closed expression of the block is known and confirms our Ansatz. We systematically study the special features of Zamolodchikov's recursion for the Π N conformal blocks. As a result, we provide a novel effective recursion relation that can be exactly solved and allows to prove the conjectured closed expressions analytically in the case of the Π1 and Π2 conformal blocks.

Tables for determining lead, uranium, and thorium isotope ages

NASA Technical Reports Server (NTRS)

Schonfeld, E.

1974-01-01

Tables for determining lead, uranium, and thorium isotope ages are presented in the form of computer printouts. Decay constants, analytical expressions for the functions evaluated, and the precision of the calculations are briefly discussed.

Closed-form Capacity Expressions for the α-μ Fading Channel with SC Diversity under Different Adaptive Transmission Strategies

NASA Astrophysics Data System (ADS)

Mohamed, Refaat; Ismail, Mahmoud H.; Newagy, Fatma; Mourad, Heba M.

2013-03-01

Stemming from the fact that the α-μ fading distribution is one of the very general fading models used in the literature to describe the small scale fading phenomenon, in this paper, closed-form expressions for the Shannon capacity of the α-μ fading channel operating under four main adaptive transmission strategies are derived assuming integer values for μ. These expressions are derived for the case of no diversity as well as for selection combining diversity with independent and identically distributed branches. The obtained expressions reduce to those previously derived in the literature for the Weibull as well as the Rayleigh fading cases, which are both special cases of the α-μ channel. Numerical results are presented for the capacity under the four adaptive transmission strategies and the effect of the fading parameter as well as the number of diversity branches is studied.

Analytical approach to peel stresses in bonded composite stiffened panels

NASA Technical Reports Server (NTRS)

Barkey, Derek A.; Madan, Ram C.; Sutton, Jason O.

1987-01-01

A closed-form solution was obtained for the stresses and displacements of two bonded beams. A system of two fourth-order and two second-order differential equations with the associated boundary equations was determined using a variational work approach. A FORTRAN computer program was devised to solve for the eigenvalues and eigenvectors of this system and to calculate the coefficients from the boundary conditions. The results were then compared with NASTRAN finite-element solutions and shown to agree closely.

Morse oscillator propagator in the high temperature limit I: Theory

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

Toutounji, Mohamad, E-mail: Mtoutounji@uaeu.ac.ae

2017-02-15

In an earlier work of the author the time evolution of Morse oscillator was studied analytically and exactly at low temperatures whereupon optical correlation functions were calculated using Morse oscillator coherent states were employed. Morse oscillator propagator in the high temperature limit is derived and a closed form of its corresponding canonical partition function is obtained. Both diagonal and off-diagonal forms of Morse oscillator propagator are derived in the high temperature limit. Partition functions of diatomic molecules are calculated. - Highlights: • Derives the quantum propagator of Morse oscillator in the high temperature limit. • Uses the resulting diagonal propagatormore » to derive a closed form of Morse oscillator partition function. • Provides a more sophisticated formula of the quantum propagator to test the accuracy of the herein results.« less

The Effective Resistance of the -Cycle Graph with Four Nearest Neighbors

NASA Astrophysics Data System (ADS)

Chair, Noureddine

2014-02-01

The exact expression for the effective resistance between any two vertices of the -cycle graph with four nearest neighbors , is given. It turns out that this expression is written in terms of the effective resistance of the -cycle graph , the square of the Fibonacci numbers, and the bisected Fibonacci numbers. As a consequence closed form formulas for the total effective resistance, the first passage time, and the mean first passage time for the simple random walk on the the -cycle graph with four nearest neighbors are obtained. Finally, a closed form formula for the effective resistance of with all first neighbors removed is obtained.

A plasma membrane sucrose-binding protein that mediates sucrose uptake shares structural and sequence similarity with seed storage proteins but remains functionally distinct.

PubMed

Overvoorde, P J; Chao, W S; Grimes, H D

1997-06-20

Photoaffinity labeling of a soybean cotyledon membrane fraction identified a sucrose-binding protein (SBP). Subsequent studies have shown that the SBP is a unique plasma membrane protein that mediates the linear uptake of sucrose in the presence of up to 30 mM external sucrose when ectopically expressed in yeast. Analysis of the SBP-deduced amino acid sequence indicates it lacks sequence similarity with other known transport proteins. Data presented here, however, indicate that the SBP shares significant sequence and structural homology with the vicilin-like seed storage proteins that organize into homotrimers. These similarities include a repeated sequence that forms the basis of the reiterated domain structure characteristic of the vicilin-like protein family. In addition, analytical ultracentrifugation and nonreducing SDS-polyacrylamide gel electrophoresis demonstrate that the SBP appears to be organized into oligomeric complexes with a Mr indicative of the existence of SBP homotrimers and homodimers. The structural similarity shared by the SBP and vicilin-like proteins provides a novel framework to explore the mechanistic basis of SBP-mediated sucrose uptake. Expression of the maize Glb protein (a vicilin-like protein closely related to the SBP) in yeast demonstrates that a closely related vicilin-like protein is unable to mediate sucrose uptake. Thus, despite sequence and structural similarities shared by the SBP and the vicilin-like protein family, the SBP is functionally divergent from other members of this group.

"Living turned inside out": the musical expression of psychotic and schizoid experience in talking heads' Remain in light.

PubMed

Brog, Michael A

2002-06-01

The rock album form, advantageously suited for the artistic expression of intrapsychic experience, has been virtually ignored by analytic writers. Remain in Light by Talking Heads is presented as an example of an album that effectively utilizes the potentialities of this form to give powerful and disconcerning musical expression to a variety of psychotic and schizoid experiences, bringing to life the formulations for these phenomena of Bion, Winnicott, Guntrip, Ogden, Grotstein and others. The album suggests a variety of mechanisms by which sound and music may serve both defensive and compensatory functions in relation to these phenomena.

A comparison between observed and analytical velocity dispersion profiles of 20 nearby galaxy clusters

NASA Astrophysics Data System (ADS)

Khan, Mohammad S.; Abdullah, Mohamed H.; Ali, Gamal B.

2014-05-01

We derive analytical expression for the velocity dispersion of galaxy clusters, using the statistical mechanical approach. We compare the observed velocity dispersion profiles for 20 nearby ( z≤0.1) galaxy clusters with the analytical ones. It is interesting to find that the analytical results closely match with the observed velocity dispersion profiles only if the presence of the diffuse matter in clusters is taken into consideration. This takes us to introduce a new approach to detect the ratio of diffuse mass, M diff , within a galaxy cluster. For the present sample, the ratio f= M diff / M, where M the cluster's total mass is found to has an average value of 45±12 %. This leads us to the result that nearly 45 % of the cluster mass is impeded outside the galaxies, while around 55 % of the cluster mass is settled in the galaxies.

Analytic solution of the Spencer-Lewis angular-spatial moments equations

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

Filippone, W.L.

A closed-form solution for the angular-spatial moments of the Spencer-Lewis equation is presented that is valid for infinite homogeneous media. From the moments, the electron density distribution as a function of position and path length (energy) is reconstructed for several sample problems involving plane isotropic sources of electrons in aluminium. The results are in excellent agreement with those determined numerically using the streaming ray method. The primary use of the closed form solution will most likely be to generate accurate electron transport benchmark solutions. In principle, the electron density as a function of space, path length, and direction can bemore » determined for planar sources of arbitrary angular distribution.« less

Axisymmetric deformations and stresses of unsymmetrically laminated composite cylinders in axial compression with thermally-induced preloading effects

NASA Technical Reports Server (NTRS)

Paraska, Peter J.

1993-01-01

This report documents an analytical study of the response of unsymmetrically laminated cylinders subjected to thermally-induced preloading effects and compressive axial load. Closed-form solutions are obtained for the displacements and intralaminar stresses and recursive relations for the interlaminar shear stress were obtained using the closed-form intralaminar stress solutions. For the cylinder geometries and stacking sequence examples analyzed, several important and as yet undocumented effects of including thermally-induced preloading in the analysis are observed. It should be noted that this work is easily extended to include uniform internal and/or external pressure loadings and the application of strain and stress failure theories.

Generalized bipartite quantum state discrimination problems with sequential measurements

NASA Astrophysics Data System (ADS)

Nakahira, Kenji; Kato, Kentaro; Usuda, Tsuyoshi Sasaki

2018-02-01

We investigate an optimization problem of finding quantum sequential measurements, which forms a wide class of state discrimination problems with the restriction that only local operations and one-way classical communication are allowed. Sequential measurements from Alice to Bob on a bipartite system are considered. Using the fact that the optimization problem can be formulated as a problem with only Alice's measurement and is convex programming, we derive its dual problem and necessary and sufficient conditions for an optimal solution. Our results are applicable to various practical optimization criteria, including the Bayes criterion, the Neyman-Pearson criterion, and the minimax criterion. In the setting of the problem of finding an optimal global measurement, its dual problem and necessary and sufficient conditions for an optimal solution have been widely used to obtain analytical and numerical expressions for optimal solutions. Similarly, our results are useful to obtain analytical and numerical expressions for optimal sequential measurements. Examples in which our results can be used to obtain an analytical expression for an optimal sequential measurement are provided.

Phase transitions in single macromolecules: Loop-stretch transition versus loop adsorption transition in end-grafted polymer chains

NASA Astrophysics Data System (ADS)

Zhang, Shuangshuang; Qi, Shuanhu; Klushin, Leonid I.; Skvortsov, Alexander M.; Yan, Dadong; Schmid, Friederike

2018-01-01

We use Brownian dynamics simulations and analytical theory to compare two prominent types of single molecule transitions. One is the adsorption transition of a loop (a chain with two ends bound to an attractive substrate) driven by an attraction parameter ɛ and the other is the loop-stretch transition in a chain with one end attached to a repulsive substrate, driven by an external end-force F applied to the free end. Specifically, we compare the behavior of the respective order parameters of the transitions, i.e., the mean number of surface contacts in the case of the adsorption transition and the mean position of the chain end in the case of the loop-stretch transition. Close to the transition points, both the static behavior and the dynamic behavior of chains with different length N are very well described by a scaling ansatz with the scaling parameters (ɛ - ɛ*)Nϕ (adsorption transition) and (F - F*)Nν (loop-stretch transition), respectively, where ϕ is the crossover exponent of the adsorption transition and ν is the Flory exponent. We show that both the loop-stretch and the loop adsorption transitions provide an exceptional opportunity to construct explicit analytical expressions for the crossover functions which perfectly describe all simulation results on static properties in the finite-size scaling regime. Explicit crossover functions are based on the ansatz for the analytical form of the order parameter distributions at the respective transition points. In contrast to the close similarity in equilibrium static behavior, the dynamic relaxation at the two transitions shows qualitative differences, especially in the strongly ordered regimes. This is attributed to the fact that the surface contact dynamics in a strongly adsorbed chain is governed by local processes, whereas the end height relaxation of a strongly stretched chain involves the full spectrum of Rouse modes.

Exact solutions to model surface and volume charge distributions

NASA Astrophysics Data System (ADS)

Mukhopadhyay, S.; Majumdar, N.; Bhattacharya, P.; Jash, A.; Bhattacharya, D. S.

2016-10-01

Many important problems in several branches of science and technology deal with charges distributed along a line, over a surface and within a volume. Recently, we have made use of new exact analytic solutions of surface charge distributions to develop the nearly exact Boundary Element Method (neBEM) toolkit. This 3D solver has been successful in removing some of the major drawbacks of the otherwise elegant Green's function approach and has been found to be very accurate throughout the computational domain, including near- and far-field regions. Use of truly distributed singularities (in contrast to nodally concentrated ones) on rectangular and right-triangular elements used for discretizing any three-dimensional geometry has essentially removed many of the numerical and physical singularities associated with the conventional BEM. In this work, we will present this toolkit and the development of several numerical models of space charge based on exact closed-form expressions. In one of the models, Particles on Surface (ParSur), the space charge inside a small elemental volume of any arbitrary shape is represented as being smeared on several surfaces representing the volume. From the studies, it can be concluded that the ParSur model is successful in getting the estimates close to those obtained using the first-principles, especially close to and within the cell. In the paper, we will show initial applications of ParSur and other models in problems related to high energy physics.

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.

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.

Electromagnetic pulses, localized and causal

NASA Astrophysics Data System (ADS)

Lekner, John

2018-01-01

We show that pulse solutions of the wave equation can be expressed as time Fourier superpositions of scalar monochromatic beam wave functions (solutions of the Helmholtz equation). This formulation is shown to be equivalent to Bateman's integral expression for solutions of the wave equation, for axially symmetric solutions. A closed-form one-parameter solution of the wave equation, containing no backward-propagating parts, is constructed from a beam which is the tight-focus limit of two families of beams. Application is made to transverse electric and transverse magnetic pulses, with evaluation of the energy, momentum and angular momentum for a pulse based on the general localized and causal form. Such pulses can be represented as superpositions of photons. Explicit total energy and total momentum values are given for the one-parameter closed-form pulse.

Thermodynamics of Quantum Gases for the Entire Range of Temperature

ERIC Educational Resources Information Center

Biswas, Shyamal; Jana, Debnarayan

2012-01-01

We have analytically explored the thermodynamics of free Bose and Fermi gases for the entire range of temperature, and have extended the same for harmonically trapped cases. We have obtained approximate chemical potentials for the quantum gases in closed forms of temperature so that the thermodynamic properties of the quantum gases become…

Cylindrical and spherical solitary waves in an electron-acoustic plasma with vortex electron distribution

NASA Astrophysics Data System (ADS)

Demiray, Hilmi; El-Zahar, Essam R.

2018-04-01

We consider the nonlinear propagation of electron-acoustic waves in a plasma composed of a cold electron fluid, hot electrons obeying a trapped/vortex-like distribution, and stationary ions. The basic nonlinear equations of the above described plasma are re-examined in the cylindrical (spherical) coordinates by employing the reductive perturbation technique. The modified cylindrical (spherical) KdV equation with fractional power nonlinearity is obtained as the evolution equation. Due to the nature of nonlinearity, this evolution equation cannot be reduced to the conventional KdV equation. A new family of closed form analytical approximate solution to the evolution equation and a comparison with numerical solution are presented and the results are depicted in some 2D and 3D figures. The results reveal that both solutions are in good agreement and the method can be used to obtain a new progressive wave solution for such evolution equations. Moreover, the resulting closed form analytical solution allows us to carry out a parametric study to investigate the effect of the physical parameters on the solution behavior of the modified cylindrical (spherical) KdV equation.

Polar decomposition for attitude determination from vector observations

NASA Technical Reports Server (NTRS)

Bar-Itzhack, Itzhack Y.

1993-01-01

This work treats the problem of weighted least squares fitting of a 3D Euclidean-coordinate transformation matrix to a set of unit vectors measured in the reference and transformed coordinates. A closed-form analytic solution to the problem is re-derived. The fact that the solution is the closest orthogonal matrix to some matrix defined on the measured vectors and their weights is clearly demonstrated. Several known algorithms for computing the analytic closed form solution are considered. An algorithm is discussed which is based on the polar decomposition of matrices into the closest unitary matrix to the decomposed matrix and a Hermitian matrix. A somewhat longer improved algorithm is suggested too. A comparison of several algorithms is carried out using simulated data as well as real data from the Upper Atmosphere Research Satellite. The comparison is based on accuracy and time consumption. It is concluded that the algorithms based on polar decomposition yield a simple although somewhat less accurate solution. The precision of the latter algorithms increase with the number of the measured vectors and with the accuracy of their measurement.

An analytical study of the endoreversible Curzon-Ahlborn cycle for a non-linear heat transfer law

NASA Astrophysics Data System (ADS)

Páez-Hernández, Ricardo T.; Portillo-Díaz, Pedro; Ladino-Luna, Delfino; Ramírez-Rojas, Alejandro; Pacheco-Paez, Juan C.

2016-01-01

In the present article, an endoreversible Curzon-Ahlborn engine is studied by considering a non-linear heat transfer law, particularly the Dulong-Petit heat transfer law, using the `componendo and dividendo' rule as well as a simple differentiation to obtain the Curzon-Ahlborn efficiency as proposed by Agrawal in 2009. This rule is actually a change of variable that simplifies a two-variable problem to a one-variable problem. From elemental calculus, we obtain an analytical expression of efficiency and the power output. The efficiency is given only in terms of the temperatures of the reservoirs, such as both Carnot and Curzon-Ahlborn cycles. We make a comparison between efficiencies measured in real power plants and theoretical values from analytical expressions obtained in this article and others found in literature from several other authors. This comparison shows that the theoretical values of efficiency are close to real efficiency, and in some cases, they are exactly the same. Therefore, we can say that the Agrawal method is good in calculating thermal engine efficiencies approximately.

An analytic resolution to the competition between Lyman-Werner radiation and metal winds in direct collapse black hole hosts

NASA Astrophysics Data System (ADS)

Agarwal, Bhaskar; Regan, John; Klessen, Ralf S.; Downes, Turlough P.; Zackrisson, Erik

2017-10-01

A near pristine atomic cooling halo close to a star forming galaxy offers a natural pathway for forming massive direct collapse black hole (DCBH) seeds, which could be the progenitors of the z > 6 redshift quasars. The close proximity of the haloes enables a sufficient Lyman-Werner flux to effectively dissociate H2 in the core of the atomic cooling halo. A mild background may also be required to delay star formation in the atomic cooling halo, often attributed to distant background galaxies. In this paper, we investigate the impact of metal pollution from both the background galaxies and the close star forming galaxy under extremely unfavourable conditions such as instantaneous metal mixing. We find that within the time window of DCBH formation, the level of pollution never exceeds the critical threshold (Zcr ˜ 1 × 10-5 Z⊙) and attains a maximum metallicity of Z ˜ 2 × 10- 6 Z⊙. As the system evolves, the metallicity eventually exceeds the critical threshold, long after the DCBH has formed.

Capillary Driven Flows Along Differentially Wetted Interior Corners

NASA Technical Reports Server (NTRS)

Golliher, Eric L. (Technical Monitor); Nardin, C. L.; Weislogel, M. M.

2005-01-01

Closed-form analytic solutions useful for the design of capillary flows in a variety of containers possessing interior corners were recently collected and reviewed. Low-g drop tower and aircraft experiments performed at NASA to date show excellent agreement between theory and experiment for perfectly wetting fluids. The analytical expressions are general in terms of contact angle, but do not account for variations in contact angle between the various surfaces within the system. Such conditions may be desirable for capillary containment or to compute the behavior of capillary corner flows in containers consisting of different materials with widely varying wetting characteristics. A simple coordinate rotation is employed to recast the governing system of equations for flows in containers with interior corners with differing contact angles on the faces of the corner. The result is that a large number of capillary driven corner flows may be predicted with only slightly modified geometric functions dependent on corner angle and the two (or more) contact angles of the system. A numerical solution is employed to verify the new problem formulation. The benchmarked computations support the use of the existing theoretical approach to geometries with variable wettability. Simple experiments to confirm the theoretical findings are recommended. Favorable agreement between such experiments and the present theory may argue well for the extension of the analytic results to predict fluid performance in future large length scale capillary fluid systems for spacecraft as well as for small scale capillary systems on Earth.

An Analytical Dynamics Approach to the Control of Mechanical Systems

NASA Astrophysics Data System (ADS)

Mylapilli, Harshavardhan

A new and novel approach to the control of nonlinear mechanical systems is presented in this study. The approach is inspired by recent results in analytical dynamics that deal with the theory of constrained motion. The control requirements on the dynamical system are viewed from an analytical dynamics perspective and the theory of constrained motion is used to recast these control requirements as constraints on the dynamical system. Explicit closed form expressions for the generalized nonlinear control forces are obtained by using the fundamental equation of mechanics. The control so obtained is optimal at each instant of time and causes the constraints to be exactly satisfied. No linearizations and/or approximations of the nonlinear dynamical system are made, and no a priori structure is imposed on the nature of nonlinear controller. Three examples dealing with highly nonlinear complex dynamical systems that are chosen from diverse areas of discrete and continuum mechanics are presented to demonstrate the control approach. The first example deals with the energy control of underactuated inhomogeneous nonlinear lattices (or chains), the second example deals with the synchronization of the motion of multiple coupled slave gyros with that of a master gyro, and the final example deals with the control of incompressible hyperelastic rubber-like thin cantilever beams. Numerical simulations accompanying these examples show the ease, simplicity and the efficacy with which the control methodology can be applied and the accuracy with which the desired control objectives can be met.

Acoustic attraction, repulsion and radiation force cancellation on a pair of rigid particles with arbitrary cross-sections in 2D: Circular cylinders example

NASA Astrophysics Data System (ADS)

Mitri, F. G.

2017-11-01

The acoustic radiation forces arising on a pair of sound impenetrable cylindrical particles of arbitrary cross-sections are derived. Plane progressive, standing or quasi-standing waves with an arbitrary incidence angle are considered. Multiple scattering effects are described using the multipole expansion formalism and the addition theorem of cylindrical wave functions. An effective incident acoustic field on a particular object is determined, and used with the scattered field to derive closed-form analytical expressions for the radiation force vector components. The mathematical expressions for the radiation force components are exact, and have been formulated in partial-wave series expansions in cylindrical coordinates involving the angle of incidence, the reflection coefficient forming the progressive or the (quasi)standing wave field, the addition theorem, and the expansion coefficients. Numerical examples illustrate the analysis for two rigid circular cross-sections immersed in a non-viscous fluid. Computations for the dimensionless radiation force functions are performed with emphasis on varying the angle of incidence, the interparticle distance, the sizes of the particles as well as the characteristics of the incident field. Depending on the interparticle distance and angle of incidence, one of the particles yields neutrality; it experiences no force and becomes unresponsive (i.e., ;invisible;) to the linear momentum transfer of the effective incident field due to multiple scattering cancellation effects. Moreover, attractive or repulsive forces between the two particles may arise depending on the interparticle distance, the angle of incidence and size parameters of the particles. This study provides a complete analytical method and computations for the axial and transverse radiation force components in multiple acoustic scattering encompassing the cases of plane progressive, standing or quasi-standing waves of arbitrary incidence by a pair of scatterers. Potential applications concern the prediction of the forces used in acoustically-engineered metamaterials with reconfigurable periodicities, cloaking devices, and liquid crystals to name a few examples.

Comparison of three methods for wind turbine capacity factor estimation.

PubMed

Ditkovich, Y; Kuperman, A

2014-01-01

Three approaches to calculating capacity factor of fixed speed wind turbines are reviewed and compared using a case study. The first "quasiexact" approach utilizes discrete wind raw data (in the histogram form) and manufacturer-provided turbine power curve (also in discrete form) to numerically calculate the capacity factor. On the other hand, the second "analytic" approach employs a continuous probability distribution function, fitted to the wind data as well as continuous turbine power curve, resulting from double polynomial fitting of manufacturer-provided power curve data. The latter approach, while being an approximation, can be solved analytically thus providing a valuable insight into aspects, affecting the capacity factor. Moreover, several other merits of wind turbine performance may be derived based on the analytical approach. The third "approximate" approach, valid in case of Rayleigh winds only, employs a nonlinear approximation of the capacity factor versus average wind speed curve, only requiring rated power and rotor diameter of the turbine. It is shown that the results obtained by employing the three approaches are very close, enforcing the validity of the analytically derived approximations, which may be used for wind turbine performance evaluation.

An approach to the determination of aircraft handling qualities using pilot transfer functions

NASA Technical Reports Server (NTRS)

Adams, J. J.; Hatch, H. G., Jr.

1978-01-01

It was shown that a correlation exists between pilot-aircraft system closed-loop characteristics, determined by using analytical expressions for pilot response along with the analytical expression for the aircraft response, and pilot ratings obtained in many previous flight and simulation studies. Two different levels of preferred pilot response were used. These levels were: (1) a static gain and a second-order lag function with a lag time constant of 0.2 second; and (2) a static gain, a lead time constant of 1 second, and a 0.2-second lag time constant. If a system response with a pitch-angle time constant of 2.6 seconds and a stable oscillatory mode of motion with a period of 2.5 seconds could be achieved with the first-level pilot model, it was shown that the pilot rating will be satisfactory for that vehicle.

A robust power spectrum split cancellation-based spectrum sensing method for cognitive radio systems

NASA Astrophysics Data System (ADS)

Qi, Pei-Han; Li, Zan; Si, Jiang-Bo; Gao, Rui

2014-12-01

Spectrum sensing is an essential component to realize the cognitive radio, and the requirement for real-time spectrum sensing in the case of lacking prior information, fading channel, and noise uncertainty, indeed poses a major challenge to the classical spectrum sensing algorithms. Based on the stochastic properties of scalar transformation of power spectral density (PSD), a novel spectrum sensing algorithm, referred to as the power spectral density split cancellation method (PSC), is proposed in this paper. The PSC makes use of a scalar value as a test statistic, which is the ratio of each subband power to the full band power. Besides, by exploiting the asymptotic normality and independence of Fourier transform, the distribution of the ratio and the mathematical expressions for the probabilities of false alarm and detection in different channel models are derived. Further, the exact closed-form expression of decision threshold is calculated in accordance with Neyman—Pearson criterion. Analytical and simulation results show that the PSC is invulnerable to noise uncertainty, and can achive excellent detection performance without prior knowledge in additive white Gaussian noise and flat slow fading channels. In addition, the PSC benefits from a low computational cost, which can be completed in microseconds.

Lightweight Expression of Granular Objects (LEGO) Content Modeling Using the SNOMED CT Observables Model to Represent Nursing Assessment Data.

PubMed

Johnson, Christie

2016-01-01

This poster presentation presents a content modeling strategy using the SNOMED CT Observable Model to represent large amounts of detailed clinical data in a consistent and computable manner that can support multiple use cases. Lightweight Expression of Granular Objects (LEGOs) represent question/answer pairs on clinical data collection forms, where a question is modeled by a (usually) post-coordinated SNOMED CT expression. LEGOs transform electronic patient data into a normalized consumable, which means that the expressions can be treated as extensions of the SNOMED CT hierarchies for the purpose of performing subsumption queries and other analytics. Utilizing the LEGO approach for modeling clinical data obtained from a nursing admission assessment provides a foundation for data exchange across disparate information systems and software applications. Clinical data exchange of computable LEGO patient information enables the development of more refined data analytics, data storage and clinical decision support.

Elastic strain field due to an inclusion of a polyhedral shape with a non-uniform lattice misfit

NASA Astrophysics Data System (ADS)

Nenashev, A. V.; Dvurechenskii, A. V.

2017-03-01

An analytical solution in a closed form is obtained for the three-dimensional elastic strain distribution in an unlimited medium containing an inclusion with a coordinate-dependent lattice mismatch (an eigenstrain). Quantum dots consisting of a solid solution with a spatially varying composition are examples of such inclusions. It is assumed that both the inclusion and the surrounding medium (the matrix) are elastically isotropic and have the same Young's modulus and Poisson ratio. The inclusion shape is supposed to be an arbitrary polyhedron, and the coordinate dependence of the lattice misfit, with respect to the matrix, is assumed to be a polynomial of any degree. It is shown that, both inside and outside the inclusion, the strain tensor is expressed as a sum of contributions of all faces, edges, and vertices of the inclusion. Each of these contributions, as a function of the observation point's coordinates, is a product of some polynomial and a simple analytical function, which is the solid angle subtended by the face from the observation point (for a contribution of a face), or the potential of the uniformly charged edge (for a contribution of an edge), or the distance from the vertex to the observation point (for a contribution of a vertex). The method of constructing the relevant polynomial functions is suggested. We also found out that similar expressions describe an electrostatic or gravitational potential, as well as its first and second derivatives, of a polyhedral body with a charge/mass density that depends on coordinates polynomially.

Optimal CCD readout by digital correlated double sampling

NASA Astrophysics Data System (ADS)

Alessandri, C.; Abusleme, A.; Guzman, D.; Passalacqua, I.; Alvarez-Fontecilla, E.; Guarini, M.

2016-01-01

Digital correlated double sampling (DCDS), a readout technique for charge-coupled devices (CCD), is gaining popularity in astronomical applications. By using an oversampling ADC and a digital filter, a DCDS system can achieve a better performance than traditional analogue readout techniques at the expense of a more complex system analysis. Several attempts to analyse and optimize a DCDS system have been reported, but most of the work presented in the literature has been experimental. Some approximate analytical tools have been presented for independent parameters of the system, but the overall performance and trade-offs have not been yet modelled. Furthermore, there is disagreement among experimental results that cannot be explained by the analytical tools available. In this work, a theoretical analysis of a generic DCDS readout system is presented, including key aspects such as the signal conditioning stage, the ADC resolution, the sampling frequency and the digital filter implementation. By using a time-domain noise model, the effect of the digital filter is properly modelled as a discrete-time process, thus avoiding the imprecision of continuous-time approximations that have been used so far. As a result, an accurate, closed-form expression for the signal-to-noise ratio at the output of the readout system is reached. This expression can be easily optimized in order to meet a set of specifications for a given CCD, thus providing a systematic design methodology for an optimal readout system. Simulated results are presented to validate the theory, obtained with both time- and frequency-domain noise generation models for completeness.

Symmetric tops in combined electric fields: Conditional quasisolvability via the quantum Hamilton-Jacobi theory

NASA Astrophysics Data System (ADS)

Schatz, Konrad; Friedrich, Bretislav; Becker, Simon; Schmidt, Burkhard

2018-05-01

We make use of the quantum Hamilton-Jacobi (QHJ) theory to investigate conditional quasisolvability of the quantum symmetric top subject to combined electric fields (symmetric top pendulum). We derive the conditions of quasisolvability of the time-independent Schrödinger equation as well as the corresponding finite sets of exact analytic solutions. We do so for this prototypical trigonometric system as well as for its anti-isospectral hyperbolic counterpart. An examination of the algebraic and numerical spectra of these two systems reveals mutually closely related patterns. The QHJ approach allows us to retrieve the closed-form solutions for the spherical and planar pendula and the Razavy system that had been obtained in our earlier work via supersymmetric quantum mechanics as well as to find a cornucopia of additional exact analytic solutions.

Fluctuation-induced forces in confined ideal and imperfect Bose gases

NASA Astrophysics Data System (ADS)

Diehl, H. W.; Rutkevich, Sergei B.

2017-06-01

Fluctuation-induced ("Casimir") forces caused by thermal and quantum fluctuations are investigated for ideal and imperfect Bose gases confined to d -dimensional films of size ∞d -1×D under periodic (P), antiperiodic (A), Dirichlet-Dirichlet (DD), Neumann-Neumann (NN), and Robin (R) boundary conditions (BCs). The full scaling functions ΥdBC(xλ=D /λth ,xξ=D /ξ ) of the residual reduced grand potential per area φres,dBC(T ,μ ,D ) =D-(d -1 )ΥdBC(xλ,xξ) are determined for the ideal gas case with these BCs, where λth and ξ are the thermal de Broglie wavelength and the bulk correlation length, respectively. The associated limiting scaling functions ΘdBC(xξ) ≡ΥdBC(∞ ,xξ) describing the critical behavior at the bulk condensation transition are shown to agree with those previously determined from a massive free O (2 ) theory for BC=P,A,DD,DN,NN . For d =3 , they are expressed in closed analytical form in terms of polylogarithms. The analogous scaling functions ΥdBC(xλ,xξ,c1D ,c2D ) and ΘdR(xξ,c1D ,c2D ) under the RBCs (∂z-c1) ϕ |z=0=(∂z+c2) ϕ | z =D=0 with c1≥0 and c2≥0 are also determined. The corresponding scaling functions Υ∞,d P(xλ,xξ) and Θ∞,d P(xξ) for the imperfect Bose gas are shown to agree with those of the interacting Bose gas with n internal degrees of freedom in the limit n →∞ . Hence, for d =3 , Θ∞,d P(xξ) is known exactly in closed analytic form. To account for the breakdown of translation invariance in the direction perpendicular to the boundary planes implied by free BCs such as DDBCs, a modified imperfect Bose gas model is introduced that corresponds to the limit n →∞ of this interacting Bose gas. Numerically and analytically exact results for the scaling function Θ∞,3 DD(xξ) therefore follow from those of the O (2 n ) ϕ4 model for n →∞ .

Fluctuation-induced forces in confined ideal and imperfect Bose gases.

PubMed

Diehl, H W; Rutkevich, Sergei B

2017-06-01

Fluctuation-induced ("Casimir") forces caused by thermal and quantum fluctuations are investigated for ideal and imperfect Bose gases confined to d-dimensional films of size ∞^{d-1}×D under periodic (P), antiperiodic (A), Dirichlet-Dirichlet (DD), Neumann-Neumann (NN), and Robin (R) boundary conditions (BCs). The full scaling functions Υ_{d}^{BC}(x_{λ}=D/λ_{th},x_{ξ}=D/ξ) of the residual reduced grand potential per area φ_{res,d}^{BC}(T,μ,D)=D^{-(d-1)}Υ_{d}^{BC}(x_{λ},x_{ξ}) are determined for the ideal gas case with these BCs, where λ_{th} and ξ are the thermal de Broglie wavelength and the bulk correlation length, respectively. The associated limiting scaling functions Θ_{d}^{BC}(x_{ξ})≡Υ_{d}^{BC}(∞,x_{ξ}) describing the critical behavior at the bulk condensation transition are shown to agree with those previously determined from a massive free O(2) theory for BC=P,A,DD,DN,NN. For d=3, they are expressed in closed analytical form in terms of polylogarithms. The analogous scaling functions Υ_{d}^{BC}(x_{λ},x_{ξ},c_{1}D,c_{2}D) and Θ_{d}^{R}(x_{ξ},c_{1}D,c_{2}D) under the RBCs (∂_{z}-c_{1})ϕ|_{z=0}=(∂_{z}+c_{2})ϕ|_{z=D}=0 with c_{1}≥0 and c_{2}≥0 are also determined. The corresponding scaling functions Υ_{∞,d}^{P}(x_{λ},x_{ξ}) and Θ_{∞,d}^{P}(x_{ξ}) for the imperfect Bose gas are shown to agree with those of the interacting Bose gas with n internal degrees of freedom in the limit n→∞. Hence, for d=3, Θ_{∞,d}^{P}(x_{ξ}) is known exactly in closed analytic form. To account for the breakdown of translation invariance in the direction perpendicular to the boundary planes implied by free BCs such as DDBCs, a modified imperfect Bose gas model is introduced that corresponds to the limit n→∞ of this interacting Bose gas. Numerically and analytically exact results for the scaling function Θ_{∞,3}^{DD}(x_{ξ}) therefore follow from those of the O(2n)ϕ^{4} model for n→∞.

A semi-analytical method for near-trapped mode and fictitious frequencies of multiple scattering by an array of elliptical cylinders in water waves

NASA Astrophysics Data System (ADS)

Chen, Jeng-Tzong; Lee, Jia-Wei

2013-09-01

In this paper, we focus on the water wave scattering by an array of four elliptical cylinders. The null-field boundary integral equation method (BIEM) is used in conjunction with degenerate kernels and eigenfunctions expansion. The closed-form fundamental solution is expressed in terms of the degenerate kernel containing the Mathieu and the modified Mathieu functions in the elliptical coordinates. Boundary densities are represented by using the eigenfunction expansion. To avoid using the addition theorem to translate the Mathieu functions, the present approach can solve the water wave problem containing multiple elliptical cylinders in a semi-analytical manner by introducing the adaptive observer system. Regarding water wave problems, the phenomena of numerical instability of fictitious frequencies may appear when the BIEM/boundary element method (BEM) is used. Besides, the near-trapped mode for an array of four identical elliptical cylinders is observed in a special layout. Both physical (near-trapped mode) and mathematical (fictitious frequency) resonances simultaneously appear in the present paper for a water wave problem by an array of four identical elliptical cylinders. Two regularization techniques, the combined Helmholtz interior integral equation formulation (CHIEF) method and the Burton and Miller approach, are adopted to alleviate the numerical resonance due to fictitious frequency.

Capture zones for simple aquifers

USGS Publications Warehouse

McElwee, Carl D.

1991-01-01

Capture zones showing the area influenced by a well within a certain time are useful for both aquifer protection and cleanup. If hydrodynamic dispersion is neglected, a deterministic curve defines the capture zone. Analytical expressions for the capture zones can be derived for simple aquifers. However, the capture zone equations are transcendental and cannot be explicitly solved for the coordinates of the capture zone boundary. Fortunately, an iterative scheme allows the solution to proceed quickly and efficiently even on a modest personal computer. Three forms of the analytical solution must be used in an iterative scheme to cover the entire region of interest, after the extreme values of the x coordinate are determined by an iterative solution. The resulting solution is a discrete one, and usually 100-1000 intervals along the x-axis are necessary for a smooth definition of the capture zone. The presented program is written in FORTRAN and has been used in a variety of computing environments. No graphics capability is included with the program; it is assumed the user has access to a commercial package. The superposition of capture zones for multiple wells is expected to be satisfactory if the spacing is not too close. Because this program deals with simple aquifers, the results rarely will be the final word in a real application.

Lateral Torsional Buckling of Anisotropic Laminated Composite Beams Subjected to Various Loading and Boundary Conditions

NASA Astrophysics Data System (ADS)

Ahmadi, Habiburrahman

Thin-walled structures are major components in many engineering applications. When a thin-walled slender beam is subjected to lateral loads, causing moments, the beam may buckle by a combined lateral bending and twisting of cross-section, which is called lateral-torsional buckling. A generalized analytical approach for lateral-torsional buckling of anisotropic laminated, thin-walled, rectangular cross-section composite beams under various loading conditions (namely, pure bending and concentrated load) and boundary conditions (namely, simply supported and cantilever) was developed using the classical laminated plate theory (CLPT), with all considered assumptions, as a basis for the constitutive equations. Buckling of such type of members has not been addressed in the literature. Closed form buckling expressions were derived in terms of the lateral, torsional and coupling stiffness coefficients of the overall composite. These coefficients were obtained through dimensional reduction by static condensation of the 6x6 constitutive matrix mapped into an effective 2x2 coupled weak axis bending-twisting relationship. The stability of the beam under different geometric and material parameters, like length/height ratio, ply thickness, and ply orientation, was investigated. The analytical formulas were verified against finite element buckling solutions using ABAQUS for different lamination orientations showing excellent accuracy.

Viscoelastic shear lag model to predict the micromechanical behavior of tendon under dynamic tensile loading.

PubMed

Wu, Jiayu; Yuan, Hong; Li, Longyuan; Fan, Kunjie; Qian, Shanguang; Li, Bing

2018-01-21

Owing to its viscoelastic nature, tendon exhibits stress rate-dependent breaking and stiffness function. A Kelvin-Voigt viscoelastic shear lag model is proposed to illustrate the micromechanical behavior of the tendon under dynamic tensile conditions. Theoretical closed-form expressions are derived to predict the deformation and stress transfer between fibrils and interfibrillar matrix while tendon is dynamically stretched. The results from the analytical solutions demonstrate that how the fibril overlap length and fibril volume fraction affect the stress transfer and mechanical properties of tendon. We find that the viscoelastic property of interfibrillar matrix mainly results in collagen fibril failure under fast loading rate or creep rupture of tendon. However, discontinuous fibril model and hierarchical structure of tendon ensure relative sliding under slow loading rate, helping dissipate energy and protecting fibril from damage, which may be a key reason why regularly staggering alignment microstructure is widely selected in nature. According to the growth, injury, healing and healed process of tendon observed by many researchers, the conclusions presented in this paper agrees well with the experimental findings. Additionally, the emphasis of this paper is on micromechanical behavior of tendon, whereas this analytical viscoelastic shear lag model can be equally applicable to other soft or hard tissues, owning the similar microstructure. Copyright © 2017 Elsevier Ltd. All rights reserved.

Propagation of sound waves through a linear shear layer: A closed form solution

NASA Technical Reports Server (NTRS)

Scott, J. N.

1978-01-01

Closed form solutions are presented for sound propagation from a line source in or near a shear layer. The analysis was exact for all frequencies and was developed assuming a linear velocity profile in the shear layer. This assumption allowed the solution to be expressed in terms of parabolic cyclinder functions. The solution is presented for a line monopole source first embedded in the uniform flow and then in the shear layer. Solutions are also discussed for certain types of dipole and quadrupole sources. Asymptotic expansions of the exact solutions for small and large values of Strouhal number gave expressions which correspond to solutions previously obtained for these limiting cases.

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

NASA Astrophysics Data System (ADS)

Akyurek, Bengu Ozge

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

Two important limitations relating to the spiking of environmental samples with contaminants of emerging concern: How close to the real analyte concentrations are the reported recovered values?

PubMed

Michael, Costas; Bayona, Josep Maria; Lambropoulou, Dimitra; Agüera, Ana; Fatta-Kassinos, Despo

2017-06-01

Occurrence and effects of contaminants of emerging concern pose a special challenge to environmental scientists. The investigation of these effects requires reliable, valid, and comparable analytical data. To this effect, two critical aspects are raised herein, concerning the limitations of the produced analytical data. The first relates to the inherent difficulty that exists in the analysis of environmental samples, which is related to the lack of knowledge (information), in many cases, of the form(s) of the contaminant in which is present in the sample. Thus, the produced analytical data can only refer to the amount of the free contaminant ignoring the amount in which it may be present in other forms; e.g., as in chelated and conjugated form. The other important aspect refers to the way with which the spiking procedure is generally performed to determine the recovery of the analytical method. Spiking environmental samples, in particular solid samples, with standard solution followed by immediate extraction, as is the common practice, can lead to an overestimation of the recovery. This is so, because no time is given to the system to establish possible equilibria between the solid matter-inorganic and/or organic-and the contaminant. Therefore, the spiking procedure need to be reconsidered by including a study of the extractable amount of the contaminant versus the time elapsed between spiking and the extraction of the sample. This study can become an element of the validation package of the method.

Controlling the spectral shape of nonlinear Thomson scattering with proper laser chirping

DOE PAGES

Rykovanov, S. G.; Geddes, C. G. R.; Schroeder, C. B.; ...

2016-03-18

Effects of nonlinearity in Thomson scattering of a high intensity laser pulse from electrons are analyzed. Analytic expressions for laser pulse shaping in frequency (chirping) are obtained which control spectrum broadening for high laser pulse intensities. These analytic solutions allow prediction of the spectral form and required laser parameters to avoid broadening. Results of analytical and numerical calculations agree well. The control over the scattered radiation bandwidth allows narrow bandwidth sources to be produced using high scattering intensities, which in turn greatly improves scattering yield for future x- and gamma-ray sources.

Randomization-Based Inference about Latent Variables from Complex Samples: The Case of Two-Stage Sampling

ERIC Educational Resources Information Center

Li, Tiandong

2012-01-01

In large-scale assessments, such as the National Assessment of Educational Progress (NAEP), plausible values based on Multiple Imputations (MI) have been used to estimate population characteristics for latent constructs under complex sample designs. Mislevy (1991) derived a closed-form analytic solution for a fixed-effect model in creating…

Analytic saddlepoint approximation for ionization energy loss distributions

DOE PAGES

Sjue, Sky K. L.; George, Jr., Richard Neal; Mathews, David Gregory

2017-07-27

Here, we present a saddlepoint approximation for ionization energy loss distributions, valid for arbitrary relativistic velocities of the incident particle 0 < v/c < 1, provided that ionizing collisions are still the dominant energy loss mechanism. We derive a closed form solution closely related to Moyal’s distribution. This distribution is intended for use in simulations with relatively low computational overhead. The approximation generally reproduces the Vavilov most probable energy loss and full width at half maximum to better than 1% and 10%, respectively, with significantly better agreement as Vavilov’s κ approaches 1.

Analytic saddlepoint approximation for ionization energy loss distributions

NASA Astrophysics Data System (ADS)

Sjue, S. K. L.; George, R. N.; Mathews, D. G.

2017-09-01

We present a saddlepoint approximation for ionization energy loss distributions, valid for arbitrary relativistic velocities of the incident particle 0 < v / c < 1 , provided that ionizing collisions are still the dominant energy loss mechanism. We derive a closed form solution closely related to Moyal's distribution. This distribution is intended for use in simulations with relatively low computational overhead. The approximation generally reproduces the Vavilov most probable energy loss and full width at half maximum to better than 1% and 10%, respectively, with significantly better agreement as Vavilov's κ approaches 1.

Analytic saddlepoint approximation for ionization energy loss distributions

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

Sjue, Sky K. L.; George, Jr., Richard Neal; Mathews, David Gregory

Here, we present a saddlepoint approximation for ionization energy loss distributions, valid for arbitrary relativistic velocities of the incident particle 0 < v/c < 1, provided that ionizing collisions are still the dominant energy loss mechanism. We derive a closed form solution closely related to Moyal’s distribution. This distribution is intended for use in simulations with relatively low computational overhead. The approximation generally reproduces the Vavilov most probable energy loss and full width at half maximum to better than 1% and 10%, respectively, with significantly better agreement as Vavilov’s κ approaches 1.

Analytical solutions for avalanche-breakdown voltages of single-diffused Gaussian junctions

NASA Astrophysics Data System (ADS)

Shenai, K.; Lin, H. C.

1983-03-01

Closed-form solutions of the potential difference between the two edges of the depletion layer of a single diffused Gaussian p-n junction are obtained by integrating Poisson's equation and equating the magnitudes of the positive and negative charges in the depletion layer. By using the closed form solution of the static Poisson's equation and Fulop's average ionization coefficient, the ionization integral in the depletion layer is computed, which yields the correct values of avalanche breakdown voltage, depletion layer thickness at breakdown, and the peak electric field as a function of junction depth. Newton's method is used for rapid convergence. A flowchart to perform the calculations with a programmable hand-held calculator, such as the TI-59, is shown.

Conservation-form equations of unsteady open-channel flow

USGS Publications Warehouse

Lai, C.; Baltzer, R.A.; Schaffranek, R.W.

2002-01-01

The unsteady open-channel flow equations are typically expressed in a variety of forms due to the imposition of differing assumptions, use of varied dependent variables, and inclusion of different source/sink terms. Questions often arise as to whether a particular equation set is expressed in a form consistent with the conservation-law definition. The concept of conservation form is developed to clarify the meaning mathematically. Six sets of unsteady-flow equations typically used in engineering practice are presented and their conservation properties are identified and discussed. Results of the theoretical development and analysis of the equations are substantiated in a set of numerical experiments conducted using alternate equation forms. Findings of these analytical and numerical efforts demonstrate that the choice of dependent variable is the fundamental factor determining the nature of the conservation properties of any particular equation form.

An analytical method for the inverse Cauchy problem of Lame equation in a rectangle

NASA Astrophysics Data System (ADS)

Grigor’ev, Yu

2018-04-01

In this paper, we present an analytical computational method for the inverse Cauchy problem of Lame equation in the elasticity theory. A rectangular domain is frequently used in engineering structures and we only consider the analytical solution in a two-dimensional rectangle, wherein a missing boundary condition is recovered from the full measurement of stresses and displacements on an accessible boundary. The essence of the method consists in solving three independent Cauchy problems for the Laplace and Poisson equations. For each of them, the Fourier series is used to formulate a first-kind Fredholm integral equation for the unknown function of data. Then, we use a Lavrentiev regularization method, and the termwise separable property of kernel function allows us to obtain a closed-form regularized solution. As a result, for the displacement components, we obtain solutions in the form of a sum of series with three regularization parameters. The uniform convergence and error estimation of the regularized solutions are proved.

Applications of computer algebra to distributed parameter systems

NASA Technical Reports Server (NTRS)

Storch, Joel A.

1993-01-01

In the analysis of vibrations of continuous elastic systems, one often encounters complicated transcendental equations with roots directly related to the system's natural frequencies. Typically, these equations contain system parameters whose values must be specified before a numerical solution can be obtained. The present paper presents a method whereby the fundamental frequency can be obtained in analytical form to any desired degree of accuracy. The method is based upon truncation of rapidly converging series involving inverse powers of the system natural frequencies. A straightforward method to developing these series and summing them in closed form is presented. It is demonstrated how Computer Algebra can be exploited to perform the intricate analytical procedures which otherwise would render the technique difficult to apply in practice. We illustrate the method by developing two analytical approximations to the fundamental frequency of a vibrating cantilever carrying a rigid tip body. The results are compared to the numerical solution of the exact (transcendental) frequency equation over a range of system parameters.

Higher-n triangular dilatonic black holes

NASA Astrophysics Data System (ADS)

Zadora, Anton; Gal'tsov, Dmitri V.; Chen, Chiang-Mei

2018-04-01

Dilaton gravity with the form fields is known to possess dyon solutions with two horizons for the discrete "triangular" values of the dilaton coupling constant a =√{ n (n + 1) / 2 }. This sequence first obtained numerically and then explained analytically as consequence of the regularity of the dilaton, should have some higher-dimensional and/or group theoretical origin. Meanwhile, this origin was explained earlier only for n = 1 , 2 in which cases the solutions were known analytically. We extend this explanation to n = 3 , 5 presenting analytical triangular solutions for the theory with different dilaton couplings a , b in electric and magnetic sectors in which case the quantization condition reads ab = n (n + 1) / 2. The solutions are derived via the Toda chains for B2 and G2 Lie algebras. They are found in the closed form in general D space-time dimensions. Solutions satisfy the entropy product rules indicating on the microscopic origin of their entropy and have negative binding energy in the extremal case.

Creation of quasi-Dirac points in the Floquet band structure of bilayer graphene.

PubMed

Cheung, W M; Chan, K S

2017-06-01

We study the Floquet quasi-energy band structure of bilayer graphene when it is illuminated by two laser lights with frequencies [Formula: see text] and [Formula: see text] using Floquet theory. We focus on the dynamical gap formed by the conduction band with Floquet index  =  -1 and the valence band with Floquet index  =  +1 to understand how Dirac points can be formed. It is found that the dynamical gap does not have rotation symmetry in the momentum space, and quasi-Dirac points, where the conduction and valence bands almost touch, can be created when the dynamical gap closes along some directions with suitably chosen radiation parameters. We derive analytical expressions for the direction dependence of the dynamical gaps using Lowdin perturbation theory to gain a better understanding of the formation of quasi-Dirac points. When both radiations are circularly polarized, the gap can be exactly zero along some directions, when only the first and second order perturbations are considered. Higher order perturbations can open a very small gap in this case. When both radiations are linearly polarized, the gap can be exactly zero up to the fourth order perturbation and more than one quasi-Dirac point is formed. We also study the electron velocity around a dynamical gap and show that the magnitude of the velocity drops to values close to zero when the k vector is near to the gap minimum. The direction of the velocity also changes around the gap minimum, and when the gap is larger in value the change in the velocity direction is more gradual. The warping effect does not affect the formation of a Dirac point along the k x axis, while it prevents its formation when there is phase shift between the two radiations.

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.

Heterodyne efficiency of a coherent free-space optical communication model through atmospheric turbulence.

PubMed

Ren, Yongxiong; Dang, Anhong; Liu, Ling; Guo, Hong

2012-10-20

The heterodyne efficiency of a coherent free-space optical (FSO) communication model under the effects of atmospheric turbulence and misalignment is studied in this paper. To be more general, both the transmitted beam and local oscillator beam are assumed to be partially coherent based on the Gaussian Schell model (GSM). By using the derived analytical form of the cross-spectral function of a GSM beam propagating through atmospheric turbulence, a closed-form expression of heterodyne efficiency is derived, assuming that the propagation directions for the transmitted and local oscillator beams are slightly different. Then the impacts of atmospheric turbulence, configuration of the two beams (namely, beam radius and spatial coherence width), detector radius, and misalignment angle over heterodyne efficiency are examined. Numerical results suggest that the beam radius of the two overlapping beams can be optimized to achieve a maximum heterodyne efficiency according to the turbulence conditions and the detector radius. It is also found that atmospheric turbulence conditions will significantly degrade the efficiency of heterodyne detection, and compared to fully coherent beams, partially coherent beams are less sensitive to the changes in turbulence conditions and more robust against misalignment at the receiver.

Distribution of shortest path lengths in a class of node duplication network models

NASA Astrophysics Data System (ADS)

Steinbock, Chanania; Biham, Ofer; Katzav, Eytan

2017-09-01

We present analytical results for the distribution of shortest path lengths (DSPL) in a network growth model which evolves by node duplication (ND). The model captures essential properties of the structure and growth dynamics of social networks, acquaintance networks, and scientific citation networks, where duplication mechanisms play a major role. Starting from an initial seed network, at each time step a random node, referred to as a mother node, is selected for duplication. Its daughter node is added to the network, forming a link to the mother node, and with probability p to each one of its neighbors. The degree distribution of the resulting network turns out to follow a power-law distribution, thus the ND network is a scale-free network. To calculate the DSPL we derive a master equation for the time evolution of the probability Pt(L =ℓ ) , ℓ =1 ,2 ,⋯ , where L is the distance between a pair of nodes and t is the time. Finding an exact analytical solution of the master equation, we obtain a closed form expression for Pt(L =ℓ ) . The mean distance 〈L〉 t and the diameter Δt are found to scale like lnt , namely, the ND network is a small-world network. The variance of the DSPL is also found to scale like lnt . Interestingly, the mean distance and the diameter exhibit properties of a small-world network, rather than the ultrasmall-world network behavior observed in other scale-free networks, in which 〈L〉 t˜lnlnt .

Effective description of domain wall strings

NASA Astrophysics Data System (ADS)

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

2018-04-01

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

Experimental investigations of stability of static liquid fillets and liquid-gas interface in capillary passages for gas-free liquid acquisition in zero gravity

NASA Astrophysics Data System (ADS)

Purohit, Ghanshyam Purshottamdas

Experimental investigations of static liquid fillets formed between small gaps of a cylindrical surface and a flat surface are carried out. The minimum volume of liquid required to form a stable fillet and the maximum liquid content the fillet can hold before becoming unstable are studied. Fillet shapes are captured in photographs obtained by a high speed image system. Experiments were conducted using water, UPA and PF 5060 on two surfaces-stand-blasted titanium and polished copper for different surface inclinations. Experimental data are generalized using appropriate non-dimensional groups. Analytical model are developed to describe the fillet curvature. Fillet curvature data are compared against model predictions and are found to be in close agreement. Bubble point experiments were carried out to measure the capillary pressure difference across the liquid-gas interface in the channels of photo-chemically etched disk stacks. Experiments were conducted using titanium stacks of five different geometrical configurations. Both well wetting liquids (IPA and PF5060) and partially wetting liquid (water) were used during experiments. Test results are found to be in close agreement with analytical predictions. Experiments were carried out to measure the frictional pressure drop across the stack as a function of liquid flow rate using two different liquids (water and IPA) and five stacks of different geometrical configurations. A channel pressure drop model is developed by treating the flow within stack channels as fully developed laminar flow between parallel plates and solving the one-dimensional Navier Stokes equation. An alternate model is developed by treating the flow in channels as flow within porous media. Expressions are developed for effective porosity and permeability for the stacks and the pressure drop is related to these parameters. Pressure drop test results are found to be in close agreement with model predictions. As a specific application of this work, a surface tension propellant management device (PMD) that uses photo-chemically etched disk stacks as capillary elements is examined. These PMDs are used in gas pressurized liquid propellant tanks to supply gas-free propellant to rocket engines in near zero-gravity environment. The experimentally validated models are integrated to perform key analyses for predicting PMD performance in zero gravity.

Geometry of Conservation Laws for a Class of Parabolic Partial Differential Equations

NASA Astrophysics Data System (ADS)

Clelland, Jeanne Nielsen

1996-08-01

I consider the problem of computing the space of conservation laws for a second-order, parabolic partial differential equation for one function of three independent variables. The PDE is formulated as an exterior differential system {cal I} on a 12 -manifold M, and its conservation laws are identified with the vector space of closed 3-forms in the infinite prolongation of {cal I} modulo the so -called "trivial" conservation laws. I use the tools of exterior differential systems and Cartan's method of equivalence to study the structure of the space of conservation laws. My main result is:. Theorem. Any conservation law for a second-order, parabolic PDE for one function of three independent variables can be represented by a closed 3-form in the differential ideal {cal I} on the original 12-manifold M. I show that if a nontrivial conservation law exists, then {cal I} has a deprolongation to an equivalent system {cal J} on a 7-manifold N, and any conservation law for {cal I} can be expressed as a closed 3-form on N which lies in {cal J}. Furthermore, any such system in the real analytic category is locally equivalent to a system generated by a (parabolic) equation of the formA(u _{xx}u_{yy}-u_sp {xy}{2}) + B_1u_{xx }+2B_2u_{xy} +B_3u_ {yy}+C=0crwhere A, B_{i}, C are functions of x, y, t, u, u_{x}, u _{y}, u_{t}. I compute the space of conservation laws for several examples, and I begin the process of analyzing the general case using Cartan's method of equivalence. I show that the non-linearizable equation u_{t} = {1over2}e ^{-u}(u_{xx}+u_ {yy})has an infinite-dimensional space of conservation laws. This stands in contrast to the two-variable case, for which Bryant and Griffiths showed that any equation whose space of conservation laws has dimension 4 or more is locally equivalent to a linear equation, i.e., is linearizable.

Comparative Kinetic Analysis of Closed-Ended and Open-Ended Porous Sensors

NASA Astrophysics Data System (ADS)

Zhao, Yiliang; Gaur, Girija; Mernaugh, Raymond L.; Laibinis, Paul E.; Weiss, Sharon M.

2016-09-01

Efficient mass transport through porous networks is essential for achieving rapid response times in sensing applications utilizing porous materials. In this work, we show that open-ended porous membranes can overcome diffusion challenges experienced by closed-ended porous materials in a microfluidic environment. A theoretical model including both transport and reaction kinetics is employed to study the influence of flow velocity, bulk analyte concentration, analyte diffusivity, and adsorption rate on the performance of open-ended and closed-ended porous sensors integrated with flow cells. The analysis shows that open-ended pores enable analyte flow through the pores and greatly reduce the response time and analyte consumption for detecting large molecules with slow diffusivities compared with closed-ended pores for which analytes largely flow over the pores. Experimental confirmation of the results was carried out with open- and closed-ended porous silicon (PSi) microcavities fabricated in flow-through and flow-over sensor configurations, respectively. The adsorption behavior of small analytes onto the inner surfaces of closed-ended and open-ended PSi membrane microcavities was similar. However, for large analytes, PSi membranes in a flow-through scheme showed significant improvement in response times due to more efficient convective transport of analytes. The experimental results and theoretical analysis provide quantitative estimates of the benefits offered by open-ended porous membranes for different analyte systems.

Purely cubic action for string field theory

NASA Technical Reports Server (NTRS)

Horowitz, G. T.; Lykken, J.; Rohm, R.; Strominger, A.

1986-01-01

It is shown that Witten's (1986) open-bosonic-string field-theory action and a closed-string analog can be written as a purely cubic interaction term. The conventional form of the action arises by expansion around particular solutions of the classical equations of motion. The explicit background dependence of the conventional action via the Becchi-Rouet-Stora-Tyutin operator is eliminated in the cubic formulation. A closed-form expression is found for the full nonlinear gauge-transformation law.

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

Price, Morgan N.; Arkin, Adam P.; Alm, Eric J.

Operons are a major feature of all prokaryotic genomes, but how and why operon structures vary is not well understood. To elucidate the life-cycle of operons, we compared gene order between Escherichia coli K12 and its relatives and identified the recently formed and destroyed operons in E. coli. This allowed us to determine how operons form, how they become closely spaced, and how they die. Our findings suggest that operon evolution is driven by selection on gene expression patterns. First, both operon creation and operon destruction lead to large changes in gene expression patterns. For example, the removal of lysAmore » and ruvA from ancestral operons that contained essential genes allowed their expression to respond to lysine levels and DNA damage, respectively. Second, some operons have undergone accelerated evolution, with multiple new genes being added during a brief period. Third, although most operons are closely spaced because of a neutral bias towards deletion and because of selection against large overlaps, highly expressed operons tend to be widely spaced because of regulatory fine-tuning by intervening sequences. Although operon evolution seems to be adaptive, it need not be optimal: new operons often comprise functionally unrelated genes that were already in proximity before the operon formed.« less

Exact Closed-form Solutions for Lamb's Problem

NASA Astrophysics Data System (ADS)

Feng, Xi; Zhang, Haiming

2018-04-01

In this article, we report on an exact closed-form solution for the displacement at the surface of an elastic half-space elicited by a buried point source that acts at some point underneath that surface. This is commonly referred to as the 3-D Lamb's problem, for which previous solutions were restricted to sources and receivers placed at the free surface. By means of the reciprocity theorem, our solution should also be valid as a means to obtain the displacements at interior points when the source is placed at the free surface. We manage to obtain explicit results by expressing the solution in terms of elementary algebraic expression as well as elliptic integrals. We anchor our developments on Poisson's ratio 0.25 starting from Johnson's (1974) integral solutions which must be computed numerically. In the end, our closed-form results agree perfectly with the numerical results of Johnson (1974), which strongly confirms the correctness of our explicit formulas. It is hoped that in due time, these formulas may constitute a valuable canonical solution that will serve as a yardstick against which other numerical solutions can be compared and measured.

Exact closed-form solutions for Lamb's problem

NASA Astrophysics Data System (ADS)

Feng, Xi; Zhang, Haiming

2018-07-01

In this paper, we report on an exact closed-form solution for the displacement at the surface of an elastic half-space elicited by a buried point source that acts at some point underneath that surface. This is commonly referred to as the 3-D Lamb's problem for which previous solutions were restricted to sources and receivers placed at the free surface. By means of the reciprocity theorem, our solution should also be valid as a means to obtain the displacements at interior points when the source is placed at the free surface. We manage to obtain explicit results by expressing the solution in terms of elementary algebraic expression as well as elliptic integrals. We anchor our developments on Poisson's ratio 0.25 starting from Johnson's integral solutions which must be computed numerically. In the end, our closed-form results agree perfectly with the numerical results of Johnson, which strongly confirms the correctness of our explicit formulae. It is hoped that in due time, these formulae may constitute a valuable canonical solution that will serve as a yardstick against which other numerical solutions can be compared and measured.

Spectral Characteristics of Wake Vortex Sound During Roll-Up

NASA Technical Reports Server (NTRS)

Booth, Earl R., Jr. (Technical Monitor); Zhang, Yan; Wang, Frank Y.; Hardin, Jay C.

2003-01-01

This report presents an analysis of the sound spectra generated by a trailing aircraft vortex during its rolling-up process. The study demonstrates that a rolling-up vortex could produce low frequency (less than 100 Hz) sound with very high intensity (60 dB above threshold of human hearing) at a distance of 200 ft from the vortex core. The spectrum then drops o rapidly thereafter. A rigorous analytical approach has been adopted in this report to derive the spectrum of vortex sound. First, the sound pressure was solved from an alternative treatment of the Lighthill s acoustic analogy approach [1]. After the application of Green s function for free space, a tensor analysis was applied to permit the removal of the source term singularity of the wave equation in the far field. Consequently, the sound pressure is expressed in terms of the retarded time that indicates the time history and spacial distribution of the sound source. The Fourier transformation is then applied to the sound pressure to compute its spectrum. As a result, the Fourier transformation greatly simplifies the expression of the vortex sound pressure involving the retarded time, so that the numerical computation is applicable with ease for axisymmetric line vortices during the rolling-up process. The vortex model assumes that the vortex circulation is proportional to the time and the core radius is a constant. In addition, the velocity profile is assumed to be self-similar along the aircraft flight path, so that a benchmark vortex velocity profile can be devised to obtain a closed form solution, which is then used to validate the numerical calculations for other more realistic vortex profiles for which no closed form solutions are available. The study suggests that acoustic sensors operating at low frequency band could be profitably deployed for detecting the vortex sound during the rolling-up process.

Light-induced cross transport phenomena in a single-component gas

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

Chermyaninov, I. V.; Chernyak, V. G., E-mail: Vladimir.Chernyak@usu.ru

2013-07-15

The cross transport processes that occur in a single-component gas in a capillary and are caused by resonance laser radiation and pressure and temperature gradients are studied. An expression for entropy production is derived using a system of kinetic Boltzmann equations in a linear approximation. The kinetic coefficients that determine the transport processes are shown to satisfy the Onsager reciprocal relations at any Knudsen numbers and any character of the elastic interaction of gas particles with the capillary surface. The light-induced baro- and thermoeffects that take place in a closed heat-insulated system in the field of resonance laser radiation aremore » considered. Analytical expressions are obtained for the Onsager coefficients in an almost free-molecular regime. The light-induced pressure and temperature gradients that appear in a closed heat-insulated capillary under typical experimental conditions are numerically estimated.« less

Spectral simplicity of apparent complexity. II. Exact complexities and complexity spectra

NASA Astrophysics Data System (ADS)

Riechers, Paul M.; Crutchfield, James P.

2018-03-01

The meromorphic functional calculus developed in Part I overcomes the nondiagonalizability of linear operators that arises often in the temporal evolution of complex systems and is generic to the metadynamics of predicting their behavior. Using the resulting spectral decomposition, we derive closed-form expressions for correlation functions, finite-length Shannon entropy-rate approximates, asymptotic entropy rate, excess entropy, transient information, transient and asymptotic state uncertainties, and synchronization information of stochastic processes generated by finite-state hidden Markov models. This introduces analytical tractability to investigating information processing in discrete-event stochastic processes, symbolic dynamics, and chaotic dynamical systems. Comparisons reveal mathematical similarities between complexity measures originally thought to capture distinct informational and computational properties. We also introduce a new kind of spectral analysis via coronal spectrograms and the frequency-dependent spectra of past-future mutual information. We analyze a number of examples to illustrate the methods, emphasizing processes with multivariate dependencies beyond pairwise correlation. This includes spectral decomposition calculations for one representative example in full detail.

Utterance selection model of language change

NASA Astrophysics Data System (ADS)

Baxter, G. J.; Blythe, R. A.; Croft, W.; McKane, A. J.

2006-04-01

We present a mathematical formulation of a theory of language change. The theory is evolutionary in nature and has close analogies with theories of population genetics. The mathematical structure we construct similarly has correspondences with the Fisher-Wright model of population genetics, but there are significant differences. The continuous time formulation of the model is expressed in terms of a Fokker-Planck equation. This equation is exactly soluble in the case of a single speaker and can be investigated analytically in the case of multiple speakers who communicate equally with all other speakers and give their utterances equal weight. Whilst the stationary properties of this system have much in common with the single-speaker case, time-dependent properties are richer. In the particular case where linguistic forms can become extinct, we find that the presence of many speakers causes a two-stage relaxation, the first being a common marginal distribution that persists for a long time as a consequence of ultimate extinction being due to rare fluctuations.

Shear wave in a pre-stressed poroelastic medium diffracted by a rigid strip

NASA Astrophysics Data System (ADS)

Singh, Abhishek Kumar; Yadav, Ram Prasad; Kumar, Santan; Chattopadhyay, Amares

2017-10-01

The investigated work analytically addresses the diffraction of horizontally polarised shear wave by a rigid strip in a pre-stressed transversely isotropic poroelastic infinite medium. The far field solution for the diffracted displacement of shear wave has been established in closed form. The diffraction patterns for displacement in the said medium have been computed numerically and its dependence on wave number has been depicted graphically. Further, the study also delineates the pronounced influence of various affecting parameters viz. anisotropy parameter, porosity parameter, speed of the shear wave, and incident angle on the diffracted displacement of the propagating wave. The effects of horizontal as well as vertical compressive and tensile pre-stresses on diffracted displacement of propagating wave have been examined meticulously in a comparative manner. It can be remarkably quoted that porosity prevailing in the medium disfavors the diffracted displacement of the propagating wave. In addition, some special cases have been deduced from the determined expression of the diffracted displacement of shear wave at a large distance from the strip.

Fundamental limits in 3D landmark localization.

PubMed

Rohr, Karl

2005-01-01

This work analyses the accuracy of estimating the location of 3D landmarks and characteristic image structures. Based on nonlinear estimation theory we study the minimal stochastic errors of the position estimate caused by noisy data. Given analytic models of the image intensities we derive closed-form expressions for the Cramér-Rao bound for different 3D structures such as 3D edges, 3D ridges, 3D lines, and 3D blobs. It turns out, that the precision of localization depends on the noise level, the size of the region-of-interest, the width of the intensity transitions, as well as on other parameters describing the considered image structure. The derived lower bounds can serve as benchmarks and the performance of existing algorithms can be compared with them. To give an impression of the achievable accuracy numeric examples are presented. Moreover, by experimental investigations we demonstrate that the derived lower bounds can be achieved by fitting parametric intensity models directly to the image data.

Performance Prediction of a Synchronization Link for Distributed Aerospace Wireless Systems

PubMed Central

Shao, Huaizong

2013-01-01

For reasons of stealth and other operational advantages, distributed aerospace wireless systems have received much attention in recent years. In a distributed aerospace wireless system, since the transmitter and receiver placed on separated platforms which use independent master oscillators, there is no cancellation of low-frequency phase noise as in the monostatic cases. Thus, high accurate time and frequency synchronization techniques are required for distributed wireless systems. The use of a dedicated synchronization link to quantify and compensate oscillator frequency instability is investigated in this paper. With the mathematical statistical models of phase noise, closed-form analytic expressions for the synchronization link performance are derived. The possible error contributions including oscillator, phase-locked loop, and receiver noise are quantified. The link synchronization performance is predicted by utilizing the knowledge of the statistical models, system error contributions, and sampling considerations. Simulation results show that effective synchronization error compensation can be achieved by using this dedicated synchronization link. PMID:23970828

Three-dimensional couette flow of dusty fluid with heat transfer in the presence of magnetic field

NASA Astrophysics Data System (ADS)

Gayathri, R.; Govindarajan, A.; Sasikala, R.

2018-04-01

This paper is focused on the mathematical modelling of three-dimensional couette flow and heat transfer of a dusty fluid between two infinite horizontal parallel porous flat plates in the presence of an induced magnetic field. The problem is formulated using a continuum two-phase model and the resulting equations are solved analytically. The lower plate is stationary while the upper plate is undergoing uniform motion in its plane. These plates are, respectively subjected to transverse exponential injection and its corresponding removal by constant suction. Due to this type of injection velocity, the flow becomes three dimensional. The closed-form expressions for velocity and temperature fields of both the fluid and dust phase are obtained by solving the governing partial differentiation equations using the perturbation method. A selective set of graphical results is presented and discussed to show interesting features of the problem. It is found that the velocity profiles of both fluid and dust particles decrease due to the increase of (magnetic parameter) Hartmann number.

Analytical approximations for the collapse of an empty spherical bubble.

PubMed

Obreschkow, D; Bruderer, M; Farhat, M

2012-06-01

The Rayleigh equation 3/2R+RR+pρ(-1)=0 with initial conditions R(0)=R(0), R(0)=0 models the collapse of an empty spherical bubble of radius R(T) in an ideal, infinite liquid with far-field pressure p and density ρ. The solution for r≡R/R(0) as a function of time t≡T/T(c), where R(T(c))≡0, is independent of R(0), p, and ρ. While no closed-form expression for r(t) is known, we find that r(0)(t)=(1-t(2))(2/5) approximates r(t) with an error below 1%. A systematic development in orders of t(2) further yields the 0.001% approximation r(*)(t)=r(0)(t)[1-a(1)Li(2.21)(t(2))], where a(1)≈-0.01832099 is a constant and Li is the polylogarithm. The usefulness of these approximations is demonstrated by comparison to high-precision cavitation data obtained in microgravity.

Packet error rate analysis of decode-and-forward free-space optical cooperative networks in the presence of random link blockage

NASA Astrophysics Data System (ADS)

Zdravković, Nemanja; Cvetkovic, Aleksandra; Milic, Dejan; Djordjevic, Goran T.

2017-09-01

This paper analyses end-to-end packet error rate (PER) of a free-space optical decode-and-forward cooperative network over a gamma-gamma atmospheric turbulence channel in the presence of temporary random link blockage. Closed-form analytical expressions for PER are derived for the cases with and without transmission links being prone to blockage. Two cooperation protocols (denoted as 'selfish' and 'pilot-adaptive') are presented and compared, where the latter accounts for the presence of blockage and adapts transmission power. The influence of scintillation, link distance, average transmitted signal power, network topology and probability of an uplink and/or internode link being blocked are discussed when the destination applies equal gain combining. The results show that link blockage caused by obstacles can degrade system performance, causing an unavoidable PER floor. The implementation of the pilot-adaptive protocol improves performance when compared to the selfish protocol, diminishing internode link blockage and lowering the PER floor, especially for larger networks.

A Model for Developing Clinical Analytics Capacity: Closing the Loops on Outcomes to Optimize Quality.

PubMed

Eggert, Corinne; Moselle, Kenneth; Protti, Denis; Sanders, Dale

2017-01-01

Closed Loop Analytics© is receiving growing interest in healthcare as a term referring to information technology, local data and clinical analytics working together to generate evidence for improvement. The Closed Loop Analytics model consists of three loops corresponding to the decision-making levels of an organization and the associated data within each loop - Patients, Protocols, and Populations. The authors propose that each of these levels should utilize the same ecosystem of electronic health record (EHR) and enterprise data warehouse (EDW) enabled data, in a closed-loop fashion, with that data being repackaged and delivered to suit the analytic and decision support needs of each level, in support of better outcomes.

Dynamic imaging model and parameter optimization for a star tracker.

PubMed

Yan, Jinyun; Jiang, Jie; Zhang, Guangjun

2016-03-21

Under dynamic conditions, star spots move across the image plane of a star tracker and form a smeared star image. This smearing effect increases errors in star position estimation and degrades attitude accuracy. First, an analytical energy distribution model of a smeared star spot is established based on a line segment spread function because the dynamic imaging process of a star tracker is equivalent to the static imaging process of linear light sources. The proposed model, which has a clear physical meaning, explicitly reflects the key parameters of the imaging process, including incident flux, exposure time, velocity of a star spot in an image plane, and Gaussian radius. Furthermore, an analytical expression of the centroiding error of the smeared star spot is derived using the proposed model. An accurate and comprehensive evaluation of centroiding accuracy is obtained based on the expression. Moreover, analytical solutions of the optimal parameters are derived to achieve the best performance in centroid estimation. Finally, we perform numerical simulations and a night sky experiment to validate the correctness of the dynamic imaging model, the centroiding error expression, and the optimal parameters.

A combinatorial model of malware diffusion via bluetooth connections.

PubMed

Merler, Stefano; Jurman, Giuseppe

2013-01-01

We outline here the mathematical expression of a diffusion model for cellphones malware transmitted through Bluetooth channels. In particular, we provide the deterministic formula underlying the proposed infection model, in its equivalent recursive (simple but computationally heavy) and closed form (more complex but efficiently computable) expression.

Technical analysis of the Slosson Written Expression Test.

PubMed

Erford, Bradley T; Hofler, Donald B

2004-06-01

The Slosson Written Expression Test was designed to assess students ages 8-17 years at risk for difficulties in written expression. Scores from three independent samples were used to evaluate the test's reliability and validity for measuring students' written expression. Test-retest reliability of the SWET subscales ranged from .80 to .94 (n = 151), and .95 for the Written Expression Total Standard Scores. The median alternate-form reliability for students' Written Expression Total Standard Scores was .81 across the three forms. Scores on the Slosson test yielded concurrent validity coefficients (n = 143) of .60 with scores from the Woodcock-Johnson: Tests of Achievement-Third Edition Broad Written Language Domain and .49 with scores on the Test of Written Language-Third Edition Spontaneous Writing Quotient. Exploratory factor analytic procedures suggested the Slosson test is comprised of two dimensions, Writing Mechanics and Writing Maturity (47.1% and 20.1% variance accounted for, respectively). In general, the Slosson Written Expression Test presents with sufficient technical characteristics to be considered a useful written expression screening test.

Analytical solution for the transient wave propagation of a buried cylindrical P-wave line source in a semi-infinite elastic medium with a fluid surface layer

NASA Astrophysics Data System (ADS)

Shan, Zhendong; Ling, Daosheng

2018-02-01

This article develops an analytical solution for the transient wave propagation of a cylindrical P-wave line source in a semi-infinite elastic solid with a fluid layer. The analytical solution is presented in a simple closed form in which each term represents a transient physical wave. The Scholte equation is derived, through which the Scholte wave velocity can be determined. The Scholte wave is the wave that propagates along the interface between the fluid and solid. To develop the analytical solution, the wave fields in the fluid and solid are defined, their analytical solutions in the Laplace domain are derived using the boundary and interface conditions, and the solutions are then decomposed into series form according to the power series expansion method. Each item of the series solution has a clear physical meaning and represents a transient wave path. Finally, by applying Cagniard's method and the convolution theorem, the analytical solutions are transformed into the time domain. Numerical examples are provided to illustrate some interesting features in the fluid layer, the interface and the semi-infinite solid. When the P-wave velocity in the fluid is higher than that in the solid, two head waves in the solid, one head wave in the fluid and a Scholte wave at the interface are observed for the cylindrical P-wave line source.

Simulation of Aluminum Micro-mirrors for Space Applications at Cryogenic Temperatures

NASA Technical Reports Server (NTRS)

Kuhn, J. L.; Dutta, S. B.; Greenhouse, M. A.; Mott, D. B.

2000-01-01

Closed form and finite element models are developed to predict the device response of aluminum electrostatic torsion micro-mirrors fabricated on silicon substrate for space applications at operating temperatures of 30K. Initially, closed form expressions for electrostatic pressure arid mechanical restoring torque are used to predict the pull-in and release voltages at room temperature. Subsequently, a detailed mechanical finite element model is developed to predict stresses and vertical beam deflection induced by the electrostatic and thermal loads. An incremental and iterative solution method is used in conjunction with the nonlinear finite element model and closed form electrostatic equations to solve. the coupled electro-thermo-mechanical problem. The simulation results are compared with experimental measurements at room temperature of fabricated micro-mirror devices.

A simple closed-form solution for assessing concentration uncertainty

NASA Astrophysics Data System (ADS)

de Barros, F. P. J.; Fiori, Aldo; Bellin, Alberto

2011-12-01

We propose closed-form approximate solutions for the moments of a nonreactive tracer that can be used in applications, such as risk analysis. This is in line with the tenet that analytical solutions provide useful information, with minimum cost, during initial site characterization efforts and can serve as a preliminary screening tool when used with prior knowledge. We show that with the help of a few assumptions, the first-order solutions of the concentration moments proposed by Fiori and Dagan (2000) can be further simplified to assume a form similar to well-known deterministic solutions, therefore facilitating their use in applications. A highly anisotropic formation is assumed, and we neglect the transverse components of the two-particle correlation trajectory. The proposed solution compares well with the work of Fiori and Dagan while presenting the same simplicity of use of existing solutions for homogeneous porous media.

Angles-only navigation for autonomous orbital rendezvous

NASA Astrophysics Data System (ADS)

Woffinden, David C.

The proposed thesis of this dissertation has both a practical element and theoretical component which aim to answer key questions related to the use of angles-only navigation for autonomous orbital rendezvous. The first and fundamental principle to this work argues that an angles-only navigation filter can determine the relative position and orientation (pose) between two spacecraft to perform the necessary maneuvers and close proximity operations for autonomous orbital rendezvous. Second, the implementation of angles-only navigation for on-orbit applications is looked upon with skeptical eyes because of its perceived limitation of determining the relative range between two vehicles. This assumed, yet little understood subtlety can be formally characterized with a closed-form analytical observability criteria which specifies the necessary and sufficient conditions for determining the relative position and velocity with only angular measurements. With a mathematical expression of the observability criteria, it can be used to (1) identify the orbital rendezvous trajectories and maneuvers that ensure the relative position and velocity are observable for angles-only navigation, (2) quantify the degree or level of observability and (3) compute optimal maneuvers that maximize observability. In summary, the objective of this dissertation is to provide both a practical and theoretical foundation for the advancement of autonomous orbital rendezvous through the use of angles-only navigation.

Analytical solution for boundary heat fluxes from a radiating rectangular medium

NASA Technical Reports Server (NTRS)

Siegel, R.

1991-01-01

Reference is made to the work of Shah (1979) which demonstrated the possibility of partially integrating the radiative equations analytically to obtain an 'exact' solution. Shah's solution was given as a double integration of the modified Bessel function of order zero. Here, it is shown that the 'exact' solution for a rectangular region radiating to cold black walls can be conveniently derived, and expressed in simple form, by using an integral function, Sn, analogous to the exponential integral function appearing in plane-layer solutions.

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.

Too hot to handle? Analytic solutions for massive neutrino or warm dark matter cosmologies

NASA Astrophysics Data System (ADS)

Slepian, Zachary; Portillo, Stephen K. N.

2018-05-01

We obtain novel closed-form solutions to the Friedmann equation for cosmological models containing a component whose equation of state is that of radiation (w = 1/3) at early times and that of cold pressureless matter (w = 0) at late times. The equation of state smoothly transitions from the early to late-time behavior and exactly describes the evolution of a species with a Dirac Delta function distribution in momentum magnitudes |p_0| (i.e. all particles have the same |p_0|). Such a component, here termed "hot matter", is an approximate model for both neutrinos and warm dark matter. We consider it alone and in combination with cold matter and with radiation, also obtaining closed-form solutions for the growth of super-horizon perturbations in each case. The idealized model recovers t(a) to better than 1.5% accuracy for all a relative to a Fermi-Dirac distribution (as describes neutrinos). We conclude by adding the second moment of the distribution to our exact solution and then generalizing to include all moments of an arbitrary momentum distribution in a closed-form solution.

Too hot to handle? Analytic solutions for massive neutrino or warm dark matter cosmologies

NASA Astrophysics Data System (ADS)

Slepian, Zachary; Portillo, Stephen K. N.

2018-07-01

We obtain novel closed-form solutions to the Friedmann equation for cosmological models containing a component whose equation of state is that of radiation (w = 1/3) at early times and that of cold pressureless matter (w= 0) at late times. The equation of state smoothly transitions from the early- to late-time behaviour and exactly describes the evolution of a species with a Dirac delta function distribution in momentum magnitudes |{p}_0| (i.e. all particles have the same |{p}_0|). Such a component, here termed `hot matter', is an approximate model for both neutrinos and warm dark matter. We consider it alone and in combination with cold matter and with radiation, also obtaining closed-form solutions for the growth of superhorizon perturbations in each case. The idealized model recovers t(a) to better than 1.5 per cent accuracy for all a relative to a Fermi-Dirac distribution (as describes neutrinos). We conclude by adding the second moment of the distribution to our exact solution and then generalizing to include all moments of an arbitrary momentum distribution in a closed-form solution.

Closed-form analytical solutions for assessing the consequences of sea-level rise on unconfined sloping island aquifers

NASA Astrophysics Data System (ADS)

Chesnaux, R.

2016-04-01

Closed-form analytical solutions for assessing the consequences of sea-level rise on fresh groundwater oceanic island lenses are provided for the cases of both strip and circular islands. Solutions are proposed for directly calculating the change in the thickness of the lens, the changes in volume and the changes in travel time of fresh groundwater within island aquifers. The solutions apply for homogenous aquifers recharged by surface infiltration and discharged by a down-gradient, fixed-head boundary. They also take into account the inland shift of the ocean due to land surface inundation, this shift being determined by the coastal slope of inland aquifers. The solutions are given for two simple island geometries: circular islands and strip islands. Base case examples are presented to illustrate, on one hand, the amplitude of the change of the fresh groundwater lens thickness and the volume depletion of the lens in oceanic island with sea-level rise, and on the other hand, the shortening of time required for groundwater to discharge into the ocean. These consequences can now be quantified and may help decision-makers to anticipate the effects of sea-level rise on fresh groundwater availability in oceanic island aquifers.

Fast and Analytical EAP Approximation from a 4th-Order Tensor.

PubMed

Ghosh, Aurobrata; Deriche, Rachid

2012-01-01

Generalized diffusion tensor imaging (GDTI) was developed to model complex apparent diffusivity coefficient (ADC) using higher-order tensors (HOTs) and to overcome the inherent single-peak shortcoming of DTI. However, the geometry of a complex ADC profile does not correspond to the underlying structure of fibers. This tissue geometry can be inferred from the shape of the ensemble average propagator (EAP). Though interesting methods for estimating a positive ADC using 4th-order diffusion tensors were developed, GDTI in general was overtaken by other approaches, for example, the orientation distribution function (ODF), since it is considerably difficult to recuperate the EAP from a HOT model of the ADC in GDTI. In this paper, we present a novel closed-form approximation of the EAP using Hermite polynomials from a modified HOT model of the original GDTI-ADC. Since the solution is analytical, it is fast, differentiable, and the approximation converges well to the true EAP. This method also makes the effort of computing a positive ADC worthwhile, since now both the ADC and the EAP can be used and have closed forms. We demonstrate our approach with 4th-order tensors on synthetic data and in vivo human data.

Fast and Analytical EAP Approximation from a 4th-Order Tensor

PubMed Central

Ghosh, Aurobrata; Deriche, Rachid

2012-01-01

Generalized diffusion tensor imaging (GDTI) was developed to model complex apparent diffusivity coefficient (ADC) using higher-order tensors (HOTs) and to overcome the inherent single-peak shortcoming of DTI. However, the geometry of a complex ADC profile does not correspond to the underlying structure of fibers. This tissue geometry can be inferred from the shape of the ensemble average propagator (EAP). Though interesting methods for estimating a positive ADC using 4th-order diffusion tensors were developed, GDTI in general was overtaken by other approaches, for example, the orientation distribution function (ODF), since it is considerably difficult to recuperate the EAP from a HOT model of the ADC in GDTI. In this paper, we present a novel closed-form approximation of the EAP using Hermite polynomials from a modified HOT model of the original GDTI-ADC. Since the solution is analytical, it is fast, differentiable, and the approximation converges well to the true EAP. This method also makes the effort of computing a positive ADC worthwhile, since now both the ADC and the EAP can be used and have closed forms. We demonstrate our approach with 4th-order tensors on synthetic data and in vivo human data. PMID:23365552

A Combinatorial Model of Malware Diffusion via Bluetooth Connections

PubMed Central

Merler, Stefano; Jurman, Giuseppe

2013-01-01

We outline here the mathematical expression of a diffusion model for cellphones malware transmitted through Bluetooth channels. In particular, we provide the deterministic formula underlying the proposed infection model, in its equivalent recursive (simple but computationally heavy) and closed form (more complex but efficiently computable) expression. PMID:23555677

Customized Steady-State Constraints for Parameter Estimation in Non-Linear Ordinary Differential Equation Models

PubMed Central

Rosenblatt, Marcus; Timmer, Jens; Kaschek, Daniel

2016-01-01

Ordinary differential equation models have become a wide-spread approach to analyze dynamical systems and understand underlying mechanisms. Model parameters are often unknown and have to be estimated from experimental data, e.g., by maximum-likelihood estimation. In particular, models of biological systems contain a large number of parameters. To reduce the dimensionality of the parameter space, steady-state information is incorporated in the parameter estimation process. For non-linear models, analytical steady-state calculation typically leads to higher-order polynomial equations for which no closed-form solutions can be obtained. This can be circumvented by solving the steady-state equations for kinetic parameters, which results in a linear equation system with comparatively simple solutions. At the same time multiplicity of steady-state solutions is avoided, which otherwise is problematic for optimization. When solved for kinetic parameters, however, steady-state constraints tend to become negative for particular model specifications, thus, generating new types of optimization problems. Here, we present an algorithm based on graph theory that derives non-negative, analytical steady-state expressions by stepwise removal of cyclic dependencies between dynamical variables. The algorithm avoids multiple steady-state solutions by construction. We show that our method is applicable to most common classes of biochemical reaction networks containing inhibition terms, mass-action and Hill-type kinetic equations. Comparing the performance of parameter estimation for different analytical and numerical methods of incorporating steady-state information, we show that our approach is especially well-tailored to guarantee a high success rate of optimization. PMID:27243005

Customized Steady-State Constraints for Parameter Estimation in Non-Linear Ordinary Differential Equation Models.

PubMed

Rosenblatt, Marcus; Timmer, Jens; Kaschek, Daniel

2016-01-01

Ordinary differential equation models have become a wide-spread approach to analyze dynamical systems and understand underlying mechanisms. Model parameters are often unknown and have to be estimated from experimental data, e.g., by maximum-likelihood estimation. In particular, models of biological systems contain a large number of parameters. To reduce the dimensionality of the parameter space, steady-state information is incorporated in the parameter estimation process. For non-linear models, analytical steady-state calculation typically leads to higher-order polynomial equations for which no closed-form solutions can be obtained. This can be circumvented by solving the steady-state equations for kinetic parameters, which results in a linear equation system with comparatively simple solutions. At the same time multiplicity of steady-state solutions is avoided, which otherwise is problematic for optimization. When solved for kinetic parameters, however, steady-state constraints tend to become negative for particular model specifications, thus, generating new types of optimization problems. Here, we present an algorithm based on graph theory that derives non-negative, analytical steady-state expressions by stepwise removal of cyclic dependencies between dynamical variables. The algorithm avoids multiple steady-state solutions by construction. We show that our method is applicable to most common classes of biochemical reaction networks containing inhibition terms, mass-action and Hill-type kinetic equations. Comparing the performance of parameter estimation for different analytical and numerical methods of incorporating steady-state information, we show that our approach is especially well-tailored to guarantee a high success rate of optimization.

A simple analytical aerodynamic model of Langley Winged-Cone Aerospace Plane concept

NASA Technical Reports Server (NTRS)

Pamadi, Bandu N.

1994-01-01

A simple three DOF analytical aerodynamic model of the Langley Winged-Coned Aerospace Plane concept is presented in a form suitable for simulation, trajectory optimization, and guidance and control studies. The analytical model is especially suitable for methods based on variational calculus. Analytical expressions are presented for lift, drag, and pitching moment coefficients from subsonic to hypersonic Mach numbers and angles of attack up to +/- 20 deg. This analytical model has break points at Mach numbers of 1.0, 1.4, 4.0, and 6.0. Across these Mach number break points, the lift, drag, and pitching moment coefficients are made continuous but their derivatives are not. There are no break points in angle of attack. The effect of control surface deflection is not considered. The present analytical model compares well with the APAS calculations and wind tunnel test data for most angles of attack and Mach numbers.

Concise expression of a classical radiation spectrum

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

Wang, C.

1993-06-01

In this paper we present a concise expression of the classical electromagnetic radiation spectrum of a moving charge. It is shown to be equivalent to the often used and much more complicated form derived from the Lienard-Wiechert potentials when the observation distance [ital R] satisfies the condition [ital R][much gt][gamma][lambda]. The expression reveals a relationship between the radiation spectrum and the motion of the radiation source. It also forms the basis of an efficient computing approach, which is of practical value in numerical calculations of the spectral output of accelerated charges. The advantages of this approach for analytical and numericalmore » applications are discussed and the bending-magnet synchrotron radiation spectrum is calculated according to the approach.« less

Development and Applications of Liquid Sample Desorption Electrospray Ionization Mass Spectrometry

NASA Astrophysics Data System (ADS)

Zheng, Qiuling; Chen, Hao

2016-06-01

Desorption electrospray ionization mass spectrometry (DESI-MS) is a recent advance in the field of analytical chemistry. This review surveys the development of liquid sample DESI-MS (LS-DESI-MS), a variant form of DESI-MS that focuses on fast analysis of liquid samples, and its novel analy-tical applications in bioanalysis, proteomics, and reaction kinetics. Due to the capability of directly ionizing liquid samples, liquid sample DESI (LS-DESI) has been successfully used to couple MS with various analytical techniques, such as microfluidics, microextraction, electrochemistry, and chromatography. This review also covers these hyphenated techniques. In addition, several closely related ionization methods, including transmission mode DESI, thermally assisted DESI, and continuous flow-extractive DESI, are briefly discussed. The capabilities of LS-DESI extend and/or complement the utilities of traditional DESI and electrospray ionization and will find extensive and valuable analytical application in the future.

On the nonlinear dynamics of trolling-mode AFM: Analytical solution using multiple time scales method

NASA Astrophysics Data System (ADS)

Sajjadi, Mohammadreza; Pishkenari, Hossein Nejat; Vossoughi, Gholamreza

2018-06-01

Trolling mode atomic force microscopy (TR-AFM) has resolved many imaging problems by a considerable reduction of the liquid-resonator interaction forces in liquid environments. The present study develops a nonlinear model of the meniscus force exerted to the nanoneedle of TR-AFM and presents an analytical solution to the distributed-parameter model of TR-AFM resonator utilizing multiple time scales (MTS) method. Based on the developed analytical solution, the frequency-response curves of the resonator operation in air and liquid (for different penetration length of the nanoneedle) are obtained. The closed-form analytical solution and the frequency-response curves are validated by the comparison with both the finite element solution of the main partial differential equations and the experimental observations. The effect of excitation angle of the resonator on horizontal oscillation of the probe tip and the effect of different parameters on the frequency-response of the system are investigated.

Optimal Mortgage Refinancing: A Closed Form Solution.

PubMed

Agarwal, Sumit; Driscoll, John C; Laibson, David I

2013-06-01

We derive the first closed-form optimal refinancing rule: Refinance when the current mortgage interest rate falls below the original rate by at least [Formula: see text] In this formula W (.) is the Lambert W -function, [Formula: see text] ρ is the real discount rate, λ is the expected real rate of exogenous mortgage repayment, σ is the standard deviation of the mortgage rate, κ/M is the ratio of the tax-adjusted refinancing cost and the remaining mortgage value, and τ is the marginal tax rate. This expression is derived by solving a tractable class of refinancing problems. Our quantitative results closely match those reported by researchers using numerical methods.

Analysis of delamination related fracture processes in composites

NASA Technical Reports Server (NTRS)

Armanios, Erian A.

1992-01-01

An anisotropic thin walled closed section beam theory was developed based on an asymptotical analysis of the shell energy functional. The displacement field is not assumed a priori and emerges as a result of the analysis. In addition to the classical out-of-plane torsional warping, two new contributions are identified namely, axial strain and bending warping. A comparison of the derived governing equations confirms the theory developed by Reissner and Tsai. Also, explicit closed form expressions for the beam stiffness coefficients, the stress and displacement fields are provided. The predictions of the present theory were validated by comparison with finite element simulation, other closed form analyses and test data.

Semi-analytical Karhunen-Loeve representation of irregular waves based on the prolate spheroidal wave functions

NASA Astrophysics Data System (ADS)

Lee, Gibbeum; Cho, Yeunwoo

2018-01-01

A new semi-analytical approach is presented to solving the matrix eigenvalue problem or the integral equation in Karhunen-Loeve (K-L) representation of random data such as irregular ocean waves. Instead of direct numerical approach to this matrix eigenvalue problem, which may suffer from the computational inaccuracy for big data, a pair of integral and differential equations are considered, which are related to the so-called prolate spheroidal wave functions (PSWF). First, the PSWF is expressed as a summation of a small number of the analytical Legendre functions. After substituting them into the PSWF differential equation, a much smaller size matrix eigenvalue problem is obtained than the direct numerical K-L matrix eigenvalue problem. By solving this with a minimal numerical effort, the PSWF and the associated eigenvalue of the PSWF differential equation are obtained. Then, the eigenvalue of the PSWF integral equation is analytically expressed by the functional values of the PSWF and the eigenvalues obtained in the PSWF differential equation. Finally, the analytically expressed PSWFs and the eigenvalues in the PWSF integral equation are used to form the kernel matrix in the K-L integral equation for the representation of exemplary wave data such as ordinary irregular waves. It is found that, with the same accuracy, the required memory size of the present method is smaller than that of the direct numerical K-L representation and the computation time of the present method is shorter than that of the semi-analytical method based on the sinusoidal functions.

Analytic evaluation of the weighting functions for remote sensing of blackbody planetary atmospheres : the case of limb viewing geometry

NASA Technical Reports Server (NTRS)

Ustinov, Eugene A.

2006-01-01

In a recent publication (Ustinov, 2002), we proposed an analytic approach to evaluation of radiative and geophysical weighting functions for remote sensing of a blackbody planetary atmosphere, based on general linearization approach applied to the case of nadir viewing geometry. In this presentation, the general linearization approach is applied to the limb viewing geometry. The expressions, similar to those obtained in (Ustinov, 2002), are obtained for weighting functions with respect to the distance along the line of sight. Further on, these expressions are converted to the expressions for weighting functions with respect to the vertical coordinate in the atmosphere. Finally, the numerical representation of weighting functions in the form of matrices of partial derivatives of grid limb radiances with respect to the grid values of atmospheric parameters is used for a convolution with the finite field of view of the instrument.

Self-diffusion coefficients and shear viscosity of inverse power fluids: from hard- to soft-spheres.

PubMed

Heyes, D M; Brańka, A C

2008-07-21

Molecular dynamics computer simulation has been used to compute the self-diffusion coefficient, D, and shear viscosity, eta(s), of soft-sphere fluids, in which the particles interact through the soft-sphere or inverse power pair potential, phi(r) = epsilon(sigma/r)(n), where n measures the steepness or stiffness of the potential, and epsilon and sigma are a characteristic energy and distance, respectively. The simulations were carried out on monodisperse systems for a range of n values from the hard-sphere (n --> infinity) limit down to n = 4, and up to densities in excess of the fluid-solid co-existence value. A new analytical procedure is proposed which reproduces the transport coefficients at high densities, and can be used to extrapolate the data to densities higher than accurately accessible by simulation or experiment, and tending to the glass transition. This formula, DX(c-1) proportional, variant A/X + B, where c is an adjustable parameter, and X is either the packing fraction or the pressure, is a development of one proposed by Dymond. In the expression, -A/B is the value of X at the ideal glass transition (i.e., where D and eta(s)(-1) --> 0). Estimated values are presented for the packing fraction and the pressure at the glass transition for n values between the hard and soft particle limits. The above expression is also shown to reproduce the high density viscosity data of supercritical argon, krypton and nitrogen. Fits to the soft-sphere simulation transport coefficients close to solid-fluid co-existence are also made using the analytic form, ln(D) = alpha(X)X, and n-dependence of the alpha(X) is presented (X is either the packing fraction or the pressure).

Broadening of Analyte Streams due to a Transverse Pressure Gradient in Free-Flow Isoelectric Focusing

PubMed Central

Dutta, Debashis

2017-01-01

Pressure-driven cross-flows can arise in free-flow isoelectric focusing systems (FFIEF) due to a non-uniform electroosmotic flow velocity along the channel width induced by the pH gradient in this direction. In addition, variations in the channel cross-section as well as unwanted differences in hydrostatic heads at the buffer/sample inlet ports can also lead to such pressure-gradients which besides altering the equilibrium position of the sample zones have a tendency to substantially broaden their widths deteriorating the separations. In this situation, a thorough assessment of stream broadening due to transverse pressure-gradients in FFIEF devices is necessary in order to establish accurate design rules for the assay. The present article describes a mathematical framework to estimate the noted zone dispersion in FFIEF separations based on the method-of-moments approach under laminar flow conditions. A closed-form expression has been derived for the spatial variance of the analyte streams at their equilibrium positions as a function of the various operating parameters governing the assay performance. This expression predicts the normalized stream variance under the chosen conditions to be determined by two dimensionless Péclet numbers evaluated based on the transverse pressure-driven and electrophoretic solute velocities in the separation chamber, respectively. Moreover, the analysis shows that while the stream width can be expected to increase with an increase in the value of the first Péclet number, the opposite trend will be followed with respect to the latter. The noted results have been validated using Monte Carlo simulations that also establish a time/length scale over which the predicted equilibrium stream width is attained in the system. PMID:28081900

Analytical expressions for the closure probability of a stiff wormlike chain for finite capture radius.

PubMed

Guérin, T

2017-08-01

Estimating the probability that two monomers of the same polymer chain are close together is a key ingredient to characterize intramolecular reactions and polymer looping. In the case of stiff wormlike polymers (rigid fluctuating elastic rods), for which end-to-end encounters are rare events, we derive an explicit analytical formula for the probability η(r_{c}) that the distance between the chain extremities is smaller than some capture radius r_{c}. The formula is asymptotically exact in the limit of stiff chains, and it leads to the identification of two distinct scaling regimes for the closure factor, originating from a strong variation of the fluctuations of the chain orientation at closure. Our theory is compatible with existing analytical results from the literature that cover the cases of a vanishing capture radius and of nearly fully extended chains.

Excitation of low-frequency residual currents at combination frequencies of an ionising two-colour laser pulse

NASA Astrophysics Data System (ADS)

Vvedenskii, N. V.; Kostin, V. A.; Laryushin, I. D.; Silaev, A. A.

2016-05-01

We have studied the processes of excitation of low-frequency residual currents in a plasma produced through ionisation of gases by two-colour laser pulses in laser-plasma schemes for THz generation. We have developed an analytical approach that allows one to find residual currents in the case when one of the components of a two-colour pulse is weak enough. The derived analytical expressions show that the effective generation of the residual current (and hence the effective THz generation) is possible if the ratio of the frequencies in the two-colour laser pulse is close to a rational fraction with a not very big odd sum of the numerator and denominator. The results of numerical calculations (including those based on the solution of the three-dimensional time-dependent Schrödinger equation) agree well with the analytical results.

Elliptic Double-Box Integrals: Massless Scattering Amplitudes beyond Polylogarithms

NASA Astrophysics Data System (ADS)

Bourjaily, Jacob L.; McLeod, Andrew J.; Spradlin, Marcus; von Hippel, Matt; Wilhelm, Matthias

2018-03-01

We derive an analytic representation of the ten-particle, two-loop double-box integral as an elliptic integral over weight-three polylogarithms. To obtain this form, we first derive a fourfold, rational (Feynman-)parametric representation for the integral, expressed directly in terms of dual-conformally invariant cross ratios; from this, the desired form is easily obtained. The essential features of this integral are illustrated by means of a simplified toy model, and we attach the relevant expressions for both integrals in ancillary files. We propose a normalization for such integrals that renders all of their polylogarithmic degenerations pure, and we discuss the need for a new "symbology" of mixed iterated elliptic and polylogarithmic integrals in order to bring them to a more canonical form.

Analytic Scattering and Refraction Models for Exoplanet Transit Spectra

NASA Astrophysics Data System (ADS)

Robinson, Tyler D.; Fortney, Jonathan J.; Hubbard, William B.

2017-12-01

Observations of exoplanet transit spectra are essential to understanding the physics and chemistry of distant worlds. The effects of opacity sources and many physical processes combine to set the shape of a transit spectrum. Two such key processes—refraction and cloud and/or haze forward-scattering—have seen substantial recent study. However, models of these processes are typically complex, which prevents their incorporation into observational analyses and standard transit spectrum tools. In this work, we develop analytic expressions that allow for the efficient parameterization of forward-scattering and refraction effects in transit spectra. We derive an effective slant optical depth that includes a correction for forward-scattered light, and present an analytic form of this correction. We validate our correction against a full-physics transit spectrum model that includes scattering, and we explore the extent to which the omission of forward-scattering effects may bias models. Also, we verify a common analytic expression for the location of a refractive boundary, which we express in terms of the maximum pressure probed in a transit spectrum. This expression is designed to be easily incorporated into existing tools, and we discuss how the detection of a refractive boundary could help indicate the background atmospheric composition by constraining the bulk refractivity of the atmosphere. Finally, we show that opacity from Rayleigh scattering and collision-induced absorption will outweigh the effects of refraction for Jupiter-like atmospheres whose equilibrium temperatures are above 400-500 K.

Diffusion Influenced Adsorption Kinetics.

PubMed

Miura, Toshiaki; Seki, Kazuhiko

2015-08-27

When the kinetics of adsorption is influenced by the diffusive flow of solutes, the solute concentration at the surface is influenced by the surface coverage of solutes, which is given by the Langmuir-Hinshelwood adsorption equation. The diffusion equation with the boundary condition given by the Langmuir-Hinshelwood adsorption equation leads to the nonlinear integro-differential equation for the surface coverage. In this paper, we solved the nonlinear integro-differential equation using the Grünwald-Letnikov formula developed to solve fractional kinetics. Guided by the numerical results, analytical expressions for the upper and lower bounds of the exact numerical results were obtained. The upper and lower bounds were close to the exact numerical results in the diffusion- and reaction-controlled limits, respectively. We examined the validity of the two simple analytical expressions obtained in the diffusion-controlled limit. The results were generalized to include the effect of dispersive diffusion. We also investigated the effect of molecular rearrangement of anisotropic molecules on surface coverage.

POTENTIOLOGY (noun): study focusing on the development of new interatomic pair potential forms; sometimes pursued in an obsessive compulsive manner [The New Yorel Dictionary (2002, unpublished)].} IN SPECTROSCOPY: IT MATTERS

NASA Astrophysics Data System (ADS)

Le Roy, Robert J.

2009-06-01

Spectroscopists have long attempted to summarize what they know about small molecules in terms of a knowledge of potential energy curves or surfaces. For most of the past century, this involved deducing polynomial-expansion force-field coefficients from energy level expressions fitted to experimental data, or for diatomic molecules, by generating tables of many-digit RKR turning points from such expressions. In recent years, however, it has become increasingly common either to use high-level ab initio calculations to compute the desired potentials, or to determine parametrized global analytic potential functions from direct fits to spectroscopic data. In the former case, this invoked a need for robust, flexible, compact, and `portable' analytic potentials for summarizing the information contained in the (sometimes very large numbers of) ab initio points, and making them `user friendly'. In the latter case, the same properties are required for potentials used in the least-squares fitting procedure. In both cases, there is also a cardinal need for potential function forms that extrapolate sensibly, beyond the range of the experimental data or ab initio points. This talk will describe some recent developments in this area, and make a case for what is arguably the `best' general-purpose analytic potential function form now available. Applications to both diatomic molecules and simple polyatomic molecules will be discussed. footnote

Cascaded analysis of signal and noise propagation through a heterogeneous breast model.

PubMed

Mainprize, James G; Yaffe, Martin J

2010-10-01

The detectability of lesions in radiographic images can be impaired by patterns caused by the surrounding anatomic structures. The presence of such patterns is often referred to as anatomic noise. Others have previously extended signal and noise propagation theory to include variable background structure as an additional noise term and used in simulations for analysis by human and ideal observers. Here, the analytic forms of the signal and noise transfer are derived to obtain an exact expression for any input random distribution and the "power law" filter used to generate the texture of the tissue distribution. A cascaded analysis of propagation through a heterogeneous model is derived for x-ray projection through simulated heterogeneous backgrounds. This is achieved by considering transmission through the breast as a correlated amplification point process. The analytic forms of the cascaded analysis were compared to monoenergetic Monte Carlo simulations of x-ray propagation through power law structured backgrounds. As expected, it was found that although the quantum noise power component scales linearly with the x-ray signal, the anatomic noise will scale with the square of the x-ray signal. There was a good agreement between results obtained using analytic expressions for the noise power and those from Monte Carlo simulations for different background textures, random input functions, and x-ray fluence. Analytic equations for the signal and noise properties of heterogeneous backgrounds were derived. These may be used in direct analysis or as a tool to validate simulations in evaluating detectability.

Ring dark and antidark solitons in nonlocal media.

PubMed

Horikis, Theodoros P; Frantzeskakis, Dimitrios J

2016-02-01

Ring dark and antidark solitons in nonlocal media are found. These structures have, respectively, the form of annular dips or humps on top of a stable CW background, and exist in a weak or strong nonlocality regime, defined by the sign of a characteristic parameter. It is demonstrated analytically that these solitons satisfy an effective cylindrical Kadomtsev-Petviashvili (aka Johnson's) equation and, as such, can be written explicitly in closed form. Numerical simulations show that they propagate undistorted and undergo quasi-elastic collisions, attesting to their stability properties.

The asymptotic form of non-global logarithms, black disc saturation, and gluonic deserts

NASA Astrophysics Data System (ADS)

Neill, Duff

2017-01-01

We develop an asymptotic perturbation theory for the large logarithmic behavior of the non-linear integro-differential equation describing the soft correlations of QCD jet measurements, the Banfi-Marchesini-Smye (BMS) equation. This equation captures the late-time evolution of radiating color dipoles after a hard collision. This allows us to prove that at large values of the control variable (the non-global logarithm, a function of the infra-red energy scales associated with distinct hard jets in an event), the distribution has a gaussian tail. We compute the decay width analytically, giving a closed form expression, and find it to be jet geometry independent, up to the number of legs of the dipole in the active jet. Enabling the asymptotic expansion is the correct perturbative seed, where we perturb around an anzats encoding formally no real emissions, an intuition motivated by the buffer region found in jet dynamics. This must be supplemented with the correct application of the BFKL approximation to the BMS equation in collinear limits. Comparing to the asymptotics of the conformally related evolution equation encountered in small-x physics, the Balitisky-Kovchegov (BK) equation, we find that the asymptotic form of the non-global logarithms directly maps to the black-disc unitarity limit of the BK equation, despite the contrasting physical pictures. Indeed, we recover the equations of saturation physics in the final state dynamics of QCD.

Segmented strings and the McMillan map

DOE PAGES

Gubser, Steven S.; Parikh, Sarthak; Witaszczyk, Przemek

2016-07-25

We present new exact solutions describing motions of closed segmented strings in AdS 3 in terms of elliptic functions. The existence of analytic expressions is due to the integrability of the classical equations of motion, which in our examples reduce to instances of the McMillan map. Here, we also obtain a discrete evolution rule for the motion in AdS 3 of arbitrary bound states of fundamental strings and D1-branes in the test approximation.

Two kinds of Airy-related beams

NASA Astrophysics Data System (ADS)

Xu, Yiqing; Zhou, Guoquan; Zhang, Lijun; Ru, Guoyun

2015-08-01

Two kinds of Airy-related beams are introduced in this manuscript. The normalized intensity distribution in the x-direction of the two kinds of Airy-related beams is close to that of the Gaussian beam. The normalized intensity distribution in the y-direction of the two kinds of Airy-related beams is close to that of the second-order and the third-order elegant Hermite-Gaussian beams, respectively. Analytical expressions of the two kinds of Airy-related beams passing through an ABCD paraxial optical system are derived. The beam propagation factors for the two kinds of Airy-related beams are 1.933 and 2.125, respectively. Analytical expressions of the beam half widths and the kurtosis parameters of the two kinds of Airy-related beams passing through an ABCD paraxial optical system are also presented. As a numerical example, the propagation properties of the two kinds of Airy-related beams are demonstrated in free space. Moreover, the comparison between the two kinds of Airy-related beams and their corresponding elegant Hermite-Gaussian beams along the two transverse directions are performed in detail. Upon propagation, the former kind of Airy-related beam will evolve from the central bright beam into the dark hollow beam. Contrarily, the latter kind of Airy-related beam will evolve from the dark hollow beam into the central bright beam. These two kinds of Airy-related beams can be used to describe specially distributed beams.

Thermostatistical description of gas mixtures from space partitions

NASA Astrophysics Data System (ADS)

Rohrmann, R. D.; Zorec, J.

2006-10-01

The new mathematical framework based on the free energy of pure classical fluids presented by Rohrmann [Physica A 347, 221 (2005)] is extended to multicomponent systems to determine thermodynamic and structural properties of chemically complex fluids. Presently, the theory focuses on D -dimensional mixtures in the low-density limit (packing factor η<0.01 ). The formalism combines the free-energy minimization technique with space partitions that assign an available volume v to each particle. v is related to the closeness of the nearest neighbor and provides a useful tool to evaluate the perturbations experimented by particles in a fluid. The theory shows a close relationship between statistical geometry and statistical mechanics. New, unconventional thermodynamic variables and mathematical identities are derived as a result of the space division. Thermodynamic potentials μil , conjugate variable of the populations Nil of particles class i with the nearest neighbors of class l are defined and their relationships with the usual chemical potentials μi are established. Systems of hard spheres are treated as illustrative examples and their thermodynamics functions are derived analytically. The low-density expressions obtained agree nicely with those of scaled-particle theory and Percus-Yevick approximation. Several pair distribution functions are introduced and evaluated. Analytical expressions are also presented for hard spheres with attractive forces due to Kac-tails and square-well potentials. Finally, we derive general chemical equilibrium conditions.

Analytic solutions for single and multiple cylinders of gravitating polytropes in magnetostatic equilibrium

NASA Technical Reports Server (NTRS)

Lerche, I.; Low, B. C.

1980-01-01

Exact analytic solutions for the static equilibrium of a gravitating plasma polytrope in the presence of magnetic fields are presented. The means of generating various equilibrium configurations to illustrate directly the complex physical relationships between pressure, magnetic fields, and gravity in self-gravitating systems is demonstrated. One of the solutions is used to model interstellar clouds suspended by magnetic fields against the galactic gravity such as may be formed by the Parker (1966) instability. It is concluded that the pinching effect of closed loops of magnetic fields in the clouds may be a dominant agent in further collapsing the clouds following their formation.

Combined effects of heat and mass transfer to magneto hydrodynamics oscillatory dusty fluid flow in a porous channel

NASA Astrophysics Data System (ADS)

Govindarajan, A.; Vijayalakshmi, R.; Ramamurthy, V.

2018-04-01

The main aim of this article is to study the combined effects of heat and mass transfer to radiative Magneto Hydro Dynamics (MHD) oscillatory optically thin dusty fluid in a saturated porous medium channel. Based on certain assumptions, the momentum, energy, concentration equations are obtained.The governing equations are non-dimensionalised, simplified and solved analytically. The closed analytical form solutions for velocity, temperature, concentration profiles are obtained. Numerical computations are presented graphically to show the salient features of various physical parameters. The shear stress, the rate of heat transfer and the rate of mass transfer are also presented graphically.

Analytical skin friction and heat transfer formula for compressible internal flows

NASA Technical Reports Server (NTRS)

Dechant, Lawrence J.; Tattar, Marc J.

1994-01-01

An analytic, closed-form friction formula for turbulent, internal, compressible, fully developed flow was derived by extending the incompressible law-of-the-wall relation to compressible cases. The model is capable of analyzing heat transfer as a function of constant surface temperatures and surface roughness as well as analyzing adiabatic conditions. The formula reduces to Prandtl's law of friction for adiabatic, smooth, axisymmetric flow. In addition, the formula reduces to the Colebrook equation for incompressible, adiabatic, axisymmetric flow with various roughnesses. Comparisons with available experiments show that the model averages roughly 12.5 percent error for adiabatic flow and 18.5 percent error for flow involving heat transfer.

Quantifying and managing uncertainty in operational modal analysis

NASA Astrophysics Data System (ADS)

Au, Siu-Kui; Brownjohn, James M. W.; Mottershead, John E.

2018-03-01

Operational modal analysis aims at identifying the modal properties (natural frequency, damping, etc.) of a structure using only the (output) vibration response measured under ambient conditions. Highly economical and feasible, it is becoming a common practice in full-scale vibration testing. In the absence of (input) loading information, however, the modal properties have significantly higher uncertainty than their counterparts identified from free or forced vibration (known input) tests. Mastering the relationship between identification uncertainty and test configuration is of great interest to both scientists and engineers, e.g., for achievable precision limits and test planning/budgeting. Addressing this challenge beyond the current state-of-the-art that are mostly concerned with identification algorithms, this work obtains closed form analytical expressions for the identification uncertainty (variance) of modal parameters that fundamentally explains the effect of test configuration. Collectively referred as 'uncertainty laws', these expressions are asymptotically correct for well-separated modes, small damping and long data; and are applicable under non-asymptotic situations. They provide a scientific basis for planning and standardization of ambient vibration tests, where factors such as channel noise, sensor number and location can be quantitatively accounted for. The work is reported comprehensively with verification through synthetic and experimental data (laboratory and field), scientific implications and practical guidelines for planning ambient vibration tests.

A formal power series expansion-regularization approach for Lévy stable distributions: the symmetric case with \\alpha =2/M (M positive integer)

NASA Astrophysics Data System (ADS)

Crisanto-Neto, J. C.; da Luz, M. G. E.; Raposo, E. P.; Viswanathan, G. M.

2016-09-01

In practice, the Lévy α-stable distribution is usually expressed in terms of the Fourier integral of its characteristic function. Indeed, known closed form expressions are relatively scarce given the huge parameters space: 0\\lt α ≤slant 2 ({{L\\'{e}vy}} {{index}}), -1≤slant β ≤slant 1 ({{skewness}}),σ \\gt 0 ({{scale}}), and -∞ \\lt μ \\lt ∞ ({{shift}}). Hence, systematic efforts have been made towards the development of proper methods for analytically solving the mentioned integral. As a further contribution in this direction, here we propose a new way to tackle the problem. We consider an approach in which one first solves the Fourier integral through a formal (thus not necessarily convergent) series representation. Then, one uses (if necessary) a pertinent sum-regularization procedure to the resulting divergent series, so as to obtain an exact formula for the distribution, which is amenable to direct numerical calculations. As a concrete study, we address the centered, symmetric, unshifted and unscaled distribution (β =0, μ =0, σ =1), with α ={α }M=2/M, M=1,2,3\\ldots . Conceivably, the present protocol could be applied to other sets of parameter values.

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.

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

NASA Astrophysics Data System (ADS)

Sakhr, Jamal; Nieminen, John M.

2018-03-01

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

Thermodynamics of Thomas-Fermi screened Coulomb systems

NASA Technical Reports Server (NTRS)

Firey, B.; Ashcroft, N. W.

1977-01-01

We obtain in closed analytic form, estimates for the thermodynamic properties of classical fluids with pair potentials of Yukawa type, with special reference to dense fully ionized plasmas with Thomas-Fermi or Debye-Hueckel screening. We further generalize the hard-sphere perturbative approach used for similarly screened two-component mixtures, and demonstrate phase separation in this simple model of a liquid mixture of metallic helium and hydrogen.

Combinatoric analysis of heterogeneous stochastic self-assembly.

PubMed

D'Orsogna, Maria R; Zhao, Bingyu; Berenji, Bijan; Chou, Tom

2013-09-28

We analyze a fully stochastic model of heterogeneous nucleation and self-assembly in a closed system with a fixed total particle number M, and a fixed number of seeds Ns. Each seed can bind a maximum of N particles. A discrete master equation for the probability distribution of the cluster sizes is derived and the corresponding cluster concentrations are found using kinetic Monte-Carlo simulations in terms of the density of seeds, the total mass, and the maximum cluster size. In the limit of slow detachment, we also find new analytic expressions and recursion relations for the cluster densities at intermediate times and at equilibrium. Our analytic and numerical findings are compared with those obtained from classical mass-action equations and the discrepancies between the two approaches analyzed.

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.

Role of partial miscibility on pressure buildup due to constant rate injection of CO2 into closed and open brine aquifers

NASA Astrophysics Data System (ADS)

Mathias, Simon A.; Gluyas, Jon G.; GonzáLez MartíNez de Miguel, Gerardo J.; Hosseini, Seyyed A.

2011-12-01

This work extends an existing analytical solution for pressure buildup because of CO2 injection in brine aquifers by incorporating effects associated with partial miscibility. These include evaporation of water into the CO2 rich phase and dissolution of CO2 into brine and salt precipitation. The resulting equations are closed-form, including the locations of the associated leading and trailing shock fronts. Derivation of the analytical solution involves making a number of simplifying assumptions including: vertical pressure equilibrium, negligible capillary pressure, and constant fluid properties. The analytical solution is compared to results from TOUGH2 and found to accurately approximate the extent of the dry-out zone around the well, the resulting permeability enhancement due to residual brine evaporation, the volumetric saturation of precipitated salt, and the vertically averaged pressure distribution in both space and time for the four scenarios studied. While brine evaporation is found to have a considerable effect on pressure, the effect of CO2 dissolution is found to be small. The resulting equations remain simple to evaluate in spreadsheet software and represent a significant improvement on current methods for estimating pressure-limited CO2 storage capacity.

ASVCP guidelines: quality assurance for point-of-care testing in veterinary medicine.

PubMed

Flatland, Bente; Freeman, Kathleen P; Vap, Linda M; Harr, Kendal E

2013-12-01

Point-of-care testing (POCT) refers to any laboratory testing performed outside the conventional reference laboratory and implies close proximity to patients. Instrumental POCT systems consist of small, handheld or benchtop analyzers. These have potential utility in many veterinary settings, including private clinics, academic veterinary medical centers, the community (eg, remote area veterinary medical teams), and for research applications in academia, government, and industry. Concern about the quality of veterinary in-clinic testing has been expressed in published veterinary literature; however, little guidance focusing on POCT is available. Recognizing this void, the ASVCP formed a subcommittee in 2009 charged with developing quality assurance (QA) guidelines for veterinary POCT. Guidelines were developed through literature review and a consensus process. Major recommendations include (1) taking a formalized approach to POCT within the facility, (2) use of written policies, standard operating procedures, forms, and logs, (3) operator training, including periodic assessment of skills, (4) assessment of instrument analytical performance and use of both statistical quality control and external quality assessment programs, (5) use of properly established or validated reference intervals, (6) and ensuring accurate patient results reporting. Where possible, given instrument analytical performance, use of a validated 13s control rule for interpretation of control data is recommended. These guidelines are aimed at veterinarians and veterinary technicians seeking to improve management of POCT in their clinical or research setting, and address QA of small chemistry and hematology instruments. These guidelines are not intended to be all-inclusive; rather, they provide a minimum standard for maintenance of POCT instruments in the veterinary setting. © 2013 American Society for Veterinary Clinical Pathology and European Society for Veterinary Clinical Pathology.

Dynamics of protein aggregation and oligomer formation governed by secondary nucleation

NASA Astrophysics Data System (ADS)

Michaels, Thomas C. T.; Lazell, Hamish W.; Arosio, Paolo; Knowles, Tuomas P. J.

2015-08-01

The formation of aggregates in many protein systems can be significantly accelerated by secondary nucleation, a process where existing assemblies catalyse the nucleation of new species. In particular, secondary nucleation has emerged as a central process controlling the proliferation of many filamentous protein structures, including molecular species related to diseases such as sickle cell anemia and a range of neurodegenerative conditions. Increasing evidence suggests that the physical size of protein filaments plays a key role in determining their potential for deleterious interactions with living cells, with smaller aggregates of misfolded proteins, oligomers, being particularly toxic. It is thus crucial to progress towards an understanding of the factors that control the sizes of protein aggregates. However, the influence of secondary nucleation on the time evolution of aggregate size distributions has been challenging to quantify. This difficulty originates in large part from the fact that secondary nucleation couples the dynamics of species distant in size space. Here, we approach this problem by presenting an analytical treatment of the master equation describing the growth kinetics of linear protein structures proliferating through secondary nucleation and provide closed-form expressions for the temporal evolution of the resulting aggregate size distribution. We show how the availability of analytical solutions for the full filament distribution allows us to identify the key physical parameters that control the sizes of growing protein filaments. Furthermore, we use these results to probe the dynamics of the populations of small oligomeric species as they are formed through secondary nucleation and discuss the implications of our work for understanding the factors that promote or curtail the production of these species with a potentially high deleterious biological activity.

Dynamics of protein aggregation and oligomer formation governed by secondary nucleation

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

Michaels, Thomas C. T., E-mail: tctm3@cam.ac.uk; Lazell, Hamish W.; Arosio, Paolo

2015-08-07

The formation of aggregates in many protein systems can be significantly accelerated by secondary nucleation, a process where existing assemblies catalyse the nucleation of new species. In particular, secondary nucleation has emerged as a central process controlling the proliferation of many filamentous protein structures, including molecular species related to diseases such as sickle cell anemia and a range of neurodegenerative conditions. Increasing evidence suggests that the physical size of protein filaments plays a key role in determining their potential for deleterious interactions with living cells, with smaller aggregates of misfolded proteins, oligomers, being particularly toxic. It is thus crucial tomore » progress towards an understanding of the factors that control the sizes of protein aggregates. However, the influence of secondary nucleation on the time evolution of aggregate size distributions has been challenging to quantify. This difficulty originates in large part from the fact that secondary nucleation couples the dynamics of species distant in size space. Here, we approach this problem by presenting an analytical treatment of the master equation describing the growth kinetics of linear protein structures proliferating through secondary nucleation and provide closed-form expressions for the temporal evolution of the resulting aggregate size distribution. We show how the availability of analytical solutions for the full filament distribution allows us to identify the key physical parameters that control the sizes of growing protein filaments. Furthermore, we use these results to probe the dynamics of the populations of small oligomeric species as they are formed through secondary nucleation and discuss the implications of our work for understanding the factors that promote or curtail the production of these species with a potentially high deleterious biological activity.« less

Theoretical and Numerical Investigations on Shallow Tunnelling in Unsaturated Soils

NASA Astrophysics Data System (ADS)

Soranzo, Enrico; Wu, Wei

2013-04-01

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

Frequency response of a supported thermocouple wire: Effects of axial conduction

NASA Technical Reports Server (NTRS)

Forney, L. J.; Meeks, E. L.; Fralick, G. C.

1991-01-01

Theoretical expressions are derived for the steady-state frequency response of a supported thermocouple wire. In particular, the effects of axial heat conduction are demonstrated for both a supported one material wire and a two material wire with unequal material properties across the junction. For the case of a one material supported wire, an exact solution is derived which compares favorably with an approximate expression that only matches temperatures at the support junction. Moreover, for the case of a two material supported wire, an analytical expression is derived that closely correlates numerical results. Experimental data were taken with a type K supported thermocouple. The test thermocouple was constructed with dimensions to demonstrate the effects of axial heat conduction assuming constant physical properties across the junction.

Effect of damage on elastically tailored composite laminates

NASA Technical Reports Server (NTRS)

Armanios, Erian; Badir, Ashraf; Berdichevsky, Victor

1991-01-01

A variationally consistent theory is derived in order to predict the response of anisotropic thin-walled closed sections subjected to axial load, torsion and bending. The theory is valid for arbitrary cross-sections made of laminated composite materials with variable thickness and stiffness. Closed form expressions for the stiffness coefficients are provided as integrals in terms of lay-ups parameters and cross-sectional geometry. A comparison of stiffness coefficients and response with finite element predictions and a closed form solution is performed. The theory is applied to the investigation of the effect of damage on the extension-twist coupling in a thin-walled closed section beam. The damage is simulated as a progressive ply-by-ply failure. Results show that damage can have a significant effect on the extension-twist coupling.

Relaxation of structural parameters and potential coefficients of nonrigid molecules. General symmetry properties and application to ab initio study of 1,2-difluoroethane

NASA Astrophysics Data System (ADS)

Ha, T.-K.; Günthard, H. H.

1989-07-01

Structural parameters like bond length, bond angles, etc. and harmonic and anharmonic potential coefficients of molecules with internal rotation, inversion or puckering modes are generally assumed to vary with the large amplitude internal coordinates in a concerted manner (relaxation). Taking the coordinate vectors of the nuclear configuration of semirigid molecules with relaxation (SRMRs) as functions of relaxing structural parameters and finite amplitude internal coordinate, the isometric group of SRMRs is discussed and the irreducible representations of the latter are shown to classify into engendered and nonengendered ones. On this basis a concept of equivalent sets of nuclei SRMRs is introduced and an analytical expression is derived which defines the most general functional form of relaxation increments of all common types of structural parameters compatible with isometric symmetry. This formula is shown to be a close analog of an analytical expression defining the transformations induced by the isometric group of infinitesimal internal coordinates associated with typical structural parameters. Furthermore analogous formulae are given for the most general form of the relaxation of harmonic potential coefficients as a function of finite internal coordinates. The general relations are illustrated by ab initio calculations for 1,2-difluoroethane at the MP4/DZP//HF/4-31G* level for twelve values of the dihedral angle including complete structure optimization. The potential to internal rotation is found to be in essential agreement with experimentally derived data. For a complete set of ab initio structural parameters the associated relaxation increments are represented as Fourier series, which are shown to confirm the form predicted by the general formula and the isometric group of 1,2-difluoroethane. Depending on type of the structural parameters (bond length, bond angles, etc.), the associated relaxation increments appear to follow some simple rules. Similarly a complete set of harmonic potential coefficients derived from the ab initio calculations will be analyzed in terms of Fourier series and shown to conform to the symmetry requirements of the symmetry group. Relaxation of potential coefficients is found to amount to up to ≈5% for some types of diagonal and nondiagonal terms and to reflect certain "topological" rules similar to regularities of harmonic potential constants of quasi-rigid molecules found in empirical determinations of valence force fields.

Fractional Fourier transform of truncated elliptical Gaussian beams.

PubMed

Du, Xinyue; Zhao, Daomu

2006-12-20

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

Conditions and Linear Stability Analysis at the Transition to Synchronization of Three Coupled Phase Oscillators in a Ring

NASA Astrophysics Data System (ADS)

El-Nashar, Hassan F.

2017-06-01

We consider a system of three nonidentical coupled phase oscillators in a ring topology. We explore the conditions that must be satisfied in order to obtain the phases at the transition to a synchrony state. These conditions lead to the correct mathematical expressions of phases that aid to find a simple analytic formula for critical coupling when the oscillators transit to a synchronization state having a common frequency value. The finding of a simple expression for the critical coupling allows us to perform a linear stability analysis at the transition to the synchronization stage. The obtained analytic forms of the eigenvalues show that the three coupled phase oscillators with periodic boundary conditions transit to a synchrony state when a saddle-node bifurcation occurs.

Erratum: A Comparison of Closures for Stochastic Advection-Diffusion Equations

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

Jarman, Kenneth D.; Tartakovsky, Alexandre M.

2015-01-01

This note corrects an error in the authors' article [SIAM/ASA J. Uncertain. Quantif., 1 (2013), pp. 319 347] in which the cited work [Neuman, Water Resour. Res., 29(3) (1993), pp. 633 645] was incorrectly represented and attributed. Concentration covariance equations presented in our article as new were in fact previously derived in the latter work. In the original abstract, the phrase " . . .we propose a closed-form approximation to two-point covariance as a measure of uncertainty. . ." should be replaced by the phrase " . . .we study a closed-form approximation to two-point covariance, previously derived in [Neumanmore » 1993], as a measure of uncertainty." The primary results in our article--the analytical and numerical comparison of existing closure methods for specific example problems are not changed by this correction.« less

A general method for the derivation of the functional forms of the effective energy terms in coarse-grained energy functions of polymers. I. Backbone potentials of coarse-grained polypeptide chains

NASA Astrophysics Data System (ADS)

Sieradzan, Adam K.; Makowski, Mariusz; Augustynowicz, Antoni; Liwo, Adam

2017-03-01

A general and systematic method for the derivation of the functional expressions for the effective energy terms in coarse-grained force fields of polymer chains is proposed. The method is based on the expansion of the potential of mean force of the system studied in the cluster-cumulant series and expanding the all-atom energy in the Taylor series in the squares of interatomic distances about the squares of the distances between coarse-grained centers, to obtain approximate analytical expressions for the cluster cumulants. The primary degrees of freedom to average about are the angles for collective rotation of the atoms contained in the coarse-grained interaction sites about the respective virtual-bond axes. The approach has been applied to the revision of the virtual-bond-angle, virtual-bond-torsional, and backbone-local-and-electrostatic correlation potentials for the UNited RESidue (UNRES) model of polypeptide chains, demonstrating the strong dependence of the torsional and correlation potentials on virtual-bond angles, not considered in the current UNRES. The theoretical considerations are illustrated with the potentials calculated from the ab initio potential-energy surface of terminally blocked alanine by numerical integration and with the statistical potentials derived from known protein structures. The revised torsional potentials correctly indicate that virtual-bond angles close to 90° result in the preference for the turn and helical structures, while large virtual-bond angles result in the preference for polyproline II and extended backbone geometry. The revised correlation potentials correctly reproduce the preference for the formation of β-sheet structures for large values of virtual-bond angles and for the formation of α-helical structures for virtual-bond angles close to 90°.

Analytical solutions of the planar cyclic voltammetry process for two soluble species with equal diffusivities and fast electron transfer using the method of eigenfunction expansions

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

Samin, Adib; Lahti, Erik; Zhang, Jinsuo, E-mail: zhang.3558@osu.edu

Cyclic voltammetry is a powerful tool that is used for characterizing electrochemical processes. Models of cyclic voltammetry take into account the mass transport of species and the kinetics at the electrode surface. Analytical solutions of these models are not well-known due to the complexity of the boundary conditions. In this study we present closed form analytical solutions of the planar voltammetry model for two soluble species with fast electron transfer and equal diffusivities using the eigenfunction expansion method. Our solution methodology does not incorporate Laplace transforms and yields good agreement with the numerical solution. This solution method can be extendedmore » to cases that are more general and may be useful for benchmarking purposes.« less

Solving the optimal attention allocation problem in manual control

NASA Technical Reports Server (NTRS)

Kleinman, D. L.

1976-01-01

Within the context of the optimal control model of human response, analytic expressions for the gradients of closed-loop performance metrics with respect to human operator attention allocation are derived. These derivatives serve as the basis for a gradient algorithm that determines the optimal attention that a human should allocate among several display indicators in a steady-state manual control task. Application of the human modeling techniques are made to study the hover control task for a CH-46 VTOL flight tested by NASA.

Exploring the nonlinear cloud and rain equation

NASA Astrophysics Data System (ADS)

Koren, Ilan; Tziperman, Eli; Feingold, Graham

2017-01-01

Marine stratocumulus cloud decks are regarded as the reflectors of the climate system, returning back to space a significant part of the income solar radiation, thus cooling the atmosphere. Such clouds can exist in two stable modes, open and closed cells, for a wide range of environmental conditions. This emergent behavior of the system, and its sensitivity to aerosol and environmental properties, is captured by a set of nonlinear equations. Here, using linear stability analysis, we express the transition from steady to a limit-cycle state analytically, showing how it depends on the model parameters. We show that the control of the droplet concentration (N), the environmental carrying-capacity (H0), and the cloud recovery parameter (τ) can be linked by a single nondimensional parameter (μ=√{N }/(ατH0) ) , suggesting that for deeper clouds the transition from open (oscillating) to closed (stable fixed point) cells will occur for higher droplet concentration (i.e., higher aerosol loading). The analytical calculations of the possible states, and how they are affected by changes in aerosol and the environmental variables, provide an enhanced understanding of the complex interactions of clouds and rain.

Explicit robust schemes for implementation of general principal value-based constitutive models

NASA Technical Reports Server (NTRS)

Arnold, S. M.; Saleeb, A. F.; Tan, H. Q.; Zhang, Y.

1993-01-01

The issue of developing effective and robust schemes to implement general hyperelastic constitutive models is addressed. To this end, special purpose functions are used to symbolically derive, evaluate, and automatically generate the associated FORTRAN code for the explicit forms of the corresponding stress function and material tangent stiffness tensors. These explicit forms are valid for the entire deformation range. The analytical form of these explicit expressions is given here for the case in which the strain-energy potential is taken as a nonseparable polynomial function of the principle stretches.

A continued fraction resummation form of bath relaxation effect in the spin-boson model

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

Gong, Zhihao; Tang, Zhoufei; Wu, Jianlan, E-mail: jianlanwu@zju.edu.cn

2015-02-28

In the spin-boson model, a continued fraction form is proposed to systematically resum high-order quantum kinetic expansion (QKE) rate kernels, accounting for the bath relaxation effect beyond the second-order perturbation. In particular, the analytical expression of the sixth-order QKE rate kernel is derived for resummation. With higher-order correction terms systematically extracted from higher-order rate kernels, the resummed quantum kinetic expansion approach in the continued fraction form extends the Pade approximation and can fully recover the exact quantum dynamics as the expansion order increases.

Microscopic theory of Brownian motion revisited: The Rayleigh model

NASA Astrophysics Data System (ADS)

Kim, Changho; Karniadakis, George Em

2013-03-01

We investigate three force autocorrelation functions , , and and the friction coefficient γ for the Rayleigh model (a massive particle in an ideal gas) by analytic methods and molecular-dynamics (MD) simulations. Here, F and F+ are the total force and the Mori fluctuating force, respectively, whereas F0 is the force on the Brownian particle under the frozen dynamics, where the Brownian particle is held fixed and the solvent particles move under the external potential due to the presence of the Brownian particle. By using ensemble averaging and the ray representation approach, we obtain two expressions for in terms of the one-particle trajectory and corresponding expressions for γ by the time integration of these expressions. Performing MD simulations of the near-Brownian-limit (NBL) regime, we investigate the convergence of and and compare them with . We show that for a purely repulsive potential between the Brownian particle and a solvent particle, both expressions for produce in the NBL regime. On the other hand, for a potential containing an attractive component, the ray representation expression produces only the contribution of the nontrapped solvent particles. However, we show that the net contribution of the trapped particles to γ disappears, and hence we confirm that both the ensemble-averaged expression and the ray representation expression for γ are valid even if the potential contains an attractive component. We also obtain a closed-form expression of γ for the square-well potential. Finally, we discuss theoretical and practical aspects for the evaluation of and γ.

Analytic family of post-merger template waveforms

NASA Astrophysics Data System (ADS)

Del Pozzo, Walter; Nagar, Alessandro

2017-06-01

Building on the analytical description of the post-merger (ringdown) waveform of coalescing, nonprecessing, spinning binary black holes introduced by Damour and Nagar [Phys. Rev. D 90, 024054 (2014), 10.1103/PhysRevD.90.024054], we propose an analytic, closed form, time-domain, representation of the ℓ=m =2 gravitational radiation mode emitted after merger. This expression is given as a function of the component masses and dimensionless spins (m1 ,2,χ1 ,2) of the two inspiraling objects, as well as of the mass MBH and (complex) frequency σ1 of the fundamental quasinormal mode of the remnant black hole. Our proposed template is obtained by fitting the post-merger waveform part of several publicly available numerical relativity simulations from the Simulating eXtreme Spacetimes (SXS) catalog and then suitably interpolating over (symmetric) mass ratio and spins. We show that this analytic expression accurately reproduces (˜0.01 rad ) the phasing of the post-merger data of other data sets not used in its construction. This is notably the case of the spin-aligned run SXS:BBH:0305, whose intrinsic parameters are consistent with the 90% credible intervals reported in the parameter-estimation followup of GW150914 by B.P. Abbott et al. [Phys. Rev. Lett. 116, 241102 (2016), 10.1103/PhysRevLett.116.241102]. Using SXS waveforms as "experimental" data, we further show that our template could be used on the actual GW150914 data to perform a new measure of the complex frequency of the fundamental quasinormal mode so as to exploit the complete (high signal-to-noise-ratio) post-merger waveform. We assess the usefulness of our proposed template by analyzing, in a realistic setting, SXS full inspiral-merger-ringdown waveforms and constructing posterior probability distribution functions for the central frequency damping time of the first overtone of the fundamental quasinormal mode as well as for the physical parameters of the systems. We also briefly explore the possibility opened by our waveform model to test the second law of black hole dynamics. Our model will help improve current tests of general relativity, in particular the general-relativistic no-hair theorem, and allow for novel tests, such as that of the area theorem.

Asymptotic Effect of Misspecification in the Random Part of the Multilevel Model

ERIC Educational Resources Information Center

Berkhof, Johannes; Kampen, Jarl Kennard

2004-01-01

The authors examine the asymptotic effect of omitting a random coefficient in the multilevel model and derive expressions for the change in (a) the variance components estimator and (b) the estimated variance of the fixed effects estimator. They apply the method of moments, which yields a closed form expression for the omission effect. In…

Pricing timer options and variance derivatives with closed-form partial transform under the 3/2 model

PubMed Central

Zheng, Wendong; Zeng, Pingping

2016-01-01

ABSTRACT Most of the empirical studies on stochastic volatility dynamics favour the 3/2 specification over the square-root (CIR) process in the Heston model. In the context of option pricing, the 3/2 stochastic volatility model (SVM) is reported to be able to capture the volatility skew evolution better than the Heston model. In this article, we make a thorough investigation on the analytic tractability of the 3/2 SVM by proposing a closed-form formula for the partial transform of the triple joint transition density which stand for the log asset price, the quadratic variation (continuous realized variance) and the instantaneous variance, respectively. Two distinct formulations are provided for deriving the main result. The closed-form partial transform enables us to deduce a variety of marginal partial transforms and characteristic functions and plays a crucial role in pricing discretely sampled variance derivatives and exotic options that depend on both the asset price and quadratic variation. Various applications and numerical examples on pricing moment swaps and timer options with discrete monitoring feature are given to demonstrate the versatility of the partial transform under the 3/2 model. PMID:28706460

iCOSSY: An Online Tool for Context-Specific Subnetwork Discovery from Gene Expression Data

PubMed Central

Saha, Ashis; Jeon, Minji; Tan, Aik Choon; Kang, Jaewoo

2015-01-01

Pathway analyses help reveal underlying molecular mechanisms of complex biological phenotypes. Biologists tend to perform multiple pathway analyses on the same dataset, as there is no single answer. It is often inefficient for them to implement and/or install all the algorithms by themselves. Online tools can help the community in this regard. Here we present an online gene expression analytical tool called iCOSSY which implements a novel pathway-based COntext-specific Subnetwork discoverY (COSSY) algorithm. iCOSSY also includes a few modifications of COSSY to increase its reliability and interpretability. Users can upload their gene expression datasets, and discover important subnetworks of closely interacting molecules to differentiate between two phenotypes (context). They can also interactively visualize the resulting subnetworks. iCOSSY is a web server that finds subnetworks that are differentially expressed in two phenotypes. Users can visualize the subnetworks to understand the biology of the difference. PMID:26147457

Misfit stresses in a composite core-shell nanowire with an eccentric parallelepipedal core subjected to one-dimensional cross dilatation eigenstrain

NASA Astrophysics Data System (ADS)

Krasnitckii, S. A.; Kolomoetc, D. R.; Smirnov, A. M.; Gutkin, M. Yu

2017-05-01

The boundary-value problem in the classical theory of elasticity for a core-shell nanowire with an eccentric parallelepipedal core of an arbitrary rectangular cross section is solved. The core is subjected to one-dimensional cross dilatation eigenstrain. The misfit stresses are given in a closed analytical form suitable for theoretical modeling of misfit accommodation in relevant heterostructures.

Combined structures-controls optimization of lattice trusses

NASA Technical Reports Server (NTRS)

Balakrishnan, A. V.

1991-01-01

The role that distributed parameter model can play in CSI is demonstrated, in particular in combined structures controls optimization problems of importance in preliminary design. Closed form solutions can be obtained for performance criteria such as rms attitude error, making possible analytical solutions of the optimization problem. This is in contrast to the need for numerical computer solution involving the inversion of large matrices in traditional finite element model (FEM) use. Another advantage of the analytic solution is that it can provide much needed insight into phenomena that can otherwise be obscured or difficult to discern from numerical computer results. As a compromise in level of complexity between a toy lab model and a real space structure, the lattice truss used in the EPS (Earth Pointing Satellite) was chosen. The optimization problem chosen is a generic one: of minimizing the structure mass subject to a specified stability margin and to a specified upper bond on the rms attitude error, using a co-located controller and sensors. Standard FEM treating each bar as a truss element is used, while the continuum model is anisotropic Timoshenko beam model. Performance criteria are derived for each model, except that for the distributed parameter model, explicit closed form solutions was obtained. Numerical results obtained by the two model show complete agreement.

A closed form slug test theory for high permeability aquifers.

PubMed

Ostendorf, David W; DeGroot, Don J; Dunaj, Philip J; Jakubowski, Joseph

2005-01-01

We incorporate a linear estimate of casing friction into the analytical slug test theory of Springer and Gelhar (1991) for high permeability aquifers. The modified theory elucidates the influence of inertia and casing friction on consistent, closed form equations for the free surface, pressure, and velocity fluctuations for overdamped and underdamped conditions. A consistent, but small, correction for kinetic energy is included as well. A characteristic velocity linearizes the turbulent casing shear stress so that an analytical solution for attenuated, phase shifted pressure fluctuations fits a single parameter (damping frequency) to transducer data from any depth in the casing. Underdamped slug tests of 0.3, 0.6, and 1 m amplitudes at five transducer depths in a 5.1 cm diameter PVC well 21 m deep in the Plymouth-Carver Aquifer yield a consistent hydraulic conductivity of 1.5 x 10(-3) m/s. The Springer and Gelhar (1991) model underestimates the hydraulic conductivity for these tests by as much as 25% by improperly ascribing smooth turbulent casing friction to the aquifer. The match point normalization of Butler (1998) agrees with our fitted hydraulic conductivity, however, when friction is included in the damping frequency. Zurbuchen et al. (2002) use a numerical model to establish a similar sensitivity of hydraulic conductivity to nonlinear casing friction.

Creep and stress relaxation induced by interface diffusion in metal matrix composites

NASA Astrophysics Data System (ADS)

Li, Yinfeng; Li, Zhonghua

2013-03-01

An analytical solution is developed to predict the creep rate induced by interface diffusion in unidirectional fiber-reinforced and particle reinforced composites. The driving force for the interface diffusion is the normal stress acting on the interface, which is obtained from rigorous Eshelby inclusion theory. The closed-form solution is an explicit function of the applied stress, volume fraction and radius of the fiber, as well as the modulus ratio between the fiber and the matrix. It is interesting that the solution is formally similar to that of Coble creep in polycrystalline materials. For the application of the present solution in the realistic composites, the scale effect is taken into account by finite element analysis based on a unit cell. Based on the solution, a closed-form solution is also given as a description of stress relaxation induced by interfacial diffusion under constant strain. In addition, the analytical solution for the interface stress presented in this study gives some insight into the relationship between the interface diffusion and interface slip. This work was supported by the financial support from the Nature Science Foundation of China (No. 10932007), the National Basic Research Program of China (No. 2010CB631003/5), and the Doctoral Program of Higher Education of China (No. 20100073110006).

A New Approach to the Numerical Evaluation of the Inverse Radon Transform with Discrete, Noisy Data.

DTIC Science & Technology

1980-07-01

spline form. The resulting analytic expression for the inner integral in the inverse transform is then readily evaluated, and the outer (periodic...integral is replaced by a sum. The work involved to obtain the inverse transform appears to be within the capability of existing computing equipment for

Quantifying wall turbulence via a symmetry approach: A Lie group theory

NASA Astrophysics Data System (ADS)

She, Zhen-Su; Chen, Xi; Hussain, Fazle

2017-11-01

We present a symmetry-based approach which yields analytic expressions for the mean velocity and kinetic energy profiles from a Lie-group analysis. After verifying the dilation-group invariance of the Reynolds averaged Navier-Stokes equation in the presence of a wall, we select a stress and energy length function as similarity variables which are assumed to have a simple dilation-invariant form. Three kinds of (local) invariant forms of the length functions are postulated, a combination of which yields a multi-layer formula giving its distribution in the entire flow region normal to the wall. The mean velocity profile is then predicted using the mean momentum equation, which yields, in particular, analytic expressions for the (universal) wall function and separate wake functions for pipe and channel - which are validated by data from direct numerical simulations (DNS). Future applications to a variety of wall flows such as flows around flat plate or airfoil, in a Rayleigh-Benard cell or Taylor-Couette system, etc., are discussed, for which the dilation group invariance is valid in the wall-normal direction.

Analytic and numeric Green's functions for a two-dimensional electron gas in an orthogonal magnetic field

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

Cresti, Alessandro; Grosso, Giuseppe; Parravicini, Giuseppe Pastori

2006-05-15

We have derived closed analytic expressions for the Green's function of an electron in a two-dimensional electron gas threaded by a uniform perpendicular magnetic field, also in the presence of a uniform electric field and of a parabolic spatial confinement. A workable and powerful numerical procedure for the calculation of the Green's functions for a large infinitely extended quantum wire is considered exploiting a lattice model for the wire, the tight-binding representation for the corresponding matrix Green's function, and the Peierls phase factor in the Hamiltonian hopping matrix element to account for the magnetic field. The numerical evaluation of themore » Green's function has been performed by means of the decimation-renormalization method, and quite satisfactorily compared with the analytic results worked out in this paper. As an example of the versatility of the numerical and analytic tools here presented, the peculiar semilocal character of the magnetic Green's function is studied in detail because of its basic importance in determining magneto-transport properties in mesoscopic systems.« less

Point spread functions and deconvolution of ultrasonic images.

PubMed

Dalitz, Christoph; Pohle-Fröhlich, Regina; Michalk, Thorsten

2015-03-01

This article investigates the restoration of ultrasonic pulse-echo C-scan images by means of deconvolution with a point spread function (PSF). The deconvolution concept from linear system theory (LST) is linked to the wave equation formulation of the imaging process, and an analytic formula for the PSF of planar transducers is derived. For this analytic expression, different numerical and analytic approximation schemes for evaluating the PSF are presented. By comparing simulated images with measured C-scan images, we demonstrate that the assumptions of LST in combination with our formula for the PSF are a good model for the pulse-echo imaging process. To reconstruct the object from a C-scan image, we compare different deconvolution schemes: the Wiener filter, the ForWaRD algorithm, and the Richardson-Lucy algorithm. The best results are obtained with the Richardson-Lucy algorithm with total variation regularization. For distances greater or equal twice the near field distance, our experiments show that the numerically computed PSF can be replaced with a simple closed analytic term based on a far field approximation.

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

Guzatov, D. V., E-mail: dm-guzatov@mail.ru

Analytic expressions for the radiative and nonradiative decay rates for an electric quadrupole source (atom, molecule) in the vicinity of a spherical particle (dielectric, metal) have been derived and analyzed within the classical electrodynamics. It has been shown that the highest increase in the decay rates appears in the quasi-static case, when the wavelength of the transition in question is much larger than the characteristic size of the system formed by the particle and the quadrupole. Asymptotic expressions for the decay rates have been derived for this case.

Performance of Orbital Neutron Instruments for Spatially Resolved Hydrogen Measurements of Airless Planetary Bodies

PubMed Central

Elphic, Richard C.; Feldman, William C.; Funsten, Herbert O.; Prettyman, Thomas H.

2010-01-01

Abstract Orbital neutron spectroscopy has become a standard technique for measuring planetary surface compositions from orbit. While this technique has led to important discoveries, such as the deposits of hydrogen at the Moon and Mars, a limitation is its poor spatial resolution. For omni-directional neutron sensors, spatial resolutions are 1–1.5 times the spacecraft's altitude above the planetary surface (or 40–600 km for typical orbital altitudes). Neutron sensors with enhanced spatial resolution have been proposed, and one with a collimated field of view is scheduled to fly on a mission to measure lunar polar hydrogen. No quantitative studies or analyses have been published that evaluate in detail the detection and sensitivity limits of spatially resolved neutron measurements. Here, we describe two complementary techniques for evaluating the hydrogen sensitivity of spatially resolved neutron sensors: an analytic, closed-form expression that has been validated with Lunar Prospector neutron data, and a three-dimensional modeling technique. The analytic technique, called the Spatially resolved Neutron Analytic Sensitivity Approximation (SNASA), provides a straightforward method to evaluate spatially resolved neutron data from existing instruments as well as to plan for future mission scenarios. We conclude that the existing detector—the Lunar Exploration Neutron Detector (LEND)—scheduled to launch on the Lunar Reconnaissance Orbiter will have hydrogen sensitivities that are over an order of magnitude poorer than previously estimated. We further conclude that a sensor with a geometric factor of ∼ 100 cm2 Sr (compared to the LEND geometric factor of ∼ 10.9 cm2 Sr) could make substantially improved measurements of the lunar polar hydrogen spatial distribution. Key Words: Planetary instrumentation—Planetary science—Moon—Spacecraft experiments—Hydrogen. Astrobiology 10, 183–200. PMID:20298147

Delay in polling systems in heavy traffic

NASA Astrophysics Data System (ADS)

van der Mei, Robert D.

1998-10-01

We study the delay in asymmetric cyclic polling systems with general mixtures of gated and exhaustive service, with generally distributed service times and switch-over times, in heavy traffic. We obtain closed-form expressions for all moments of the delay incurred at each of the queues. The expressions are strikingly simple and can even be expressed as finite products of known factors. The results provide new insights into the heavy-traffic behavior of polling systems.

Moving from Within The Maternal: The Choreography of Analytic Eroticism.

PubMed

Elise, Dianne

2017-02-01

With Kristeva's concept of maternal eroticism (2014) as starting point, the "multiverse" of mother/child erotic sensibilities-the dance of the semiotic chora-is explored and a parallel engagement proposed within the analytic dyad. The dance of psychoanalysis is not the creative product of the patient's mind alone. Clinical work invites, requires, a choreographic engagement by the clinician in interplay with the patient. The clinician's analytic activity is thus akin to choreography: the structuring of a dance, or of a session, expresses an inner impulse brought into narrative form. The embodied art of dance parallels the clinician's creative vitality in contributing to the shaping of the movement of a session. Through formulation of an analytic eroticism, the terrain of what traditionally has been viewed as erotic transference and countertransference can be expanded to clinical benefit.

Analytic Results for e+e- --> tt-bar and gammagamma --> tt-bar Observables near the Threshold up to the Next-to-Next-to-Leading Order of NRQCD

NASA Astrophysics Data System (ADS)

Penin, A. A.; Pivovarov, A. A.

2001-02-01

We present an analytical description of top-antitop pair production near the threshold in $e^+e^-$ annihilation and $\\g\\g$ collisions. A set of basic observables considered includes the total cross sections, forward-backward asymmetry and top quark polarization. The threshold effects relevant for the basic observables are described by three universal functions related to S wave production, P wave production and S-P interference. These functions are computed analytically up to the next-to-next-to-leading order of NRQCD. The total $e^+e^-\\to t\\bar t$ cross section near the threshold is obtained in the next-to-next-to-leading order in the closed form including the contribution due to the axial coupling of top quark and mediated by the Z-boson. The effects of the running of the strong coupling constant and of the finite top quark width are taken into account analytically for the P wave production and S-P wave interference.

Analytical saturated domain orientation textures and electromechanical properties of ferroelectric ceramics due to electric/mechanical poling

NASA Astrophysics Data System (ADS)

Li, F. X.; Rajapakse, R. K. N. D.

2007-03-01

Saturated domain orientation textures of three types of pseudocubic (tetragonal, rhombohedral, and orthorhombic) ferroelectric ceramics after complete electric and uniaxial tension (compression) poling is studied analytically in this paper. A one-dimensional orientation distribution function (ODF) of the domain polar vectors is explicitly derived from the uniform inverse pole figures of the poling field axes on a stereographic projection with respect to the fixed crystallite coordinates. The analytical ODF is used to obtain the analytical solutions of saturated polarization and strain after electric/mechanical poling. Based on the closed form solution of the saturated domain orientation textures, the resultant intrinsic electromechanical properties of ferroelectric ceramics, which depend only on the ODF and properties of the corresponding single crystals, are obtained. The results show how the macroscopic symmetries of ferroelectric crystals change from 4mm (tetragonal), 3m (rhombohedral), and mm2 (orthorhombic) single crystals to a ∞mm (transversely isotropic) completely poled ceramic.

Analytical Model for Thermal Elastoplastic Stresses of Functionally Graded Materials

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

Zhai, P. C.; Chen, G.; Liu, L. S.

2008-02-15

A modification analytical model is presented for the thermal elastoplastic stresses of functionally graded materials subjected to thermal loading. The presented model follows the analytical scheme presented by Y. L. Shen and S. Suresh [6]. In the present model, the functionally graded materials are considered as multilayered materials. Each layer consists of metal and ceramic with different volume fraction. The ceramic layer and the FGM interlayers are considered as elastic brittle materials. The metal layer is considered as elastic-perfectly plastic ductile materials. Closed-form solutions for different characteristic temperature for thermal loading are presented as a function of the structure geometriesmore » and the thermomechanical properties of the materials. A main advance of the present model is that the possibility of the initial and spread of plasticity from the two sides of the ductile layers taken into account. Comparing the analytical results with the results from the finite element analysis, the thermal stresses and deformation from the present model are in good agreement with the numerical ones.« less

Current-sheet formation in two-dimensional coronal fields

NASA Astrophysics Data System (ADS)

Billinghurst, M. N.; Craig, I. J. D.; Sneyd, A. D.

1993-11-01

The formation of current sheets by shearing motions in line-tied twin-lobed fields is examined. A general analytic argument shows that current sheets form along the fieldline bounding the two lobes in the case of both symmetric and asymmetric footpoint motions. In the case of strictly antisymmetric motions however no current sheets can form. These findings are reinforced by magnetic relaxation experiments involving sheared two-lobed fields represented by Clebsh variables. It is pointed out that, although current singularites cannot be expected to form when the line-tying assumption is relaxed, the two-lobed geometry is still consistent with the formation of highly localised currents - and strong resistive dissipation - along field lines close to the bounding fieldline.

Harmonics analysis of the photonic time stretch system.

PubMed

Mei, Yuan; Xu, Boyu; Chi, Hao; Jin, Tao; Zheng, Shilie; Jin, Xiaofeng; Zhang, Xianmin

2016-09-10

Photonic time stretch (PTS) has been intensively investigated in recent decades due to its potential application to ultra-wideband analog-to-digital conversion. A high-speed analog signal can be captured by an electronic analog-to-digital converter (ADC) with the help of the PTS technique, which slows down the speed of signal in the photonic domain. Unfortunately, the process of the time stretch is not linear due to the nonlinear modulation of the electro-optic intensity modulator in the PTS system, which means the undesired harmonics distortion. In this paper, we present an exact analytical model to fully characterize the harmonics generation in the PTS systems for the first time, to the best of our knowledge. We obtain concise and closed-form expressions for all harmonics of the PTS system with either a single-arm Mach-Zehnder modulator (MZM) or a push-pull MZM. The presented model can largely simplify the PTS system design and the system parameters estimation, such as system bandwidth, harmonics power, time-bandwidth product, and dynamic range. The correctness of the mathematic model is verified by the numerical and experimental results.

Multi-diversity combining and selection for relay-assisted mixed RF/FSO system

NASA Astrophysics Data System (ADS)

Chen, Li; Wang, Weidong

2017-12-01

We propose and analyze multi-diversity combining and selection to enhance the performance of relay-assisted mixed radio frequency/free-space optics (RF/FSO) system. We focus on a practical scenario for cellular network where a single-antenna source is communicating to a multi-apertures destination through a relay equipped with multiple receive antennas and multiple transmit apertures. The RF single input multiple output (SIMO) links employ either maximal-ratio combining (MRC) or receive antenna selection (RAS), and the FSO multiple input multiple output (MIMO) links adopt either repetition coding (RC) or transmit laser selection (TLS). The performance is evaluated via an outage probability analysis over Rayleigh fading RF links and Gamma-Gamma atmospheric turbulence FSO links with pointing errors where channel state information (CSI) assisted amplify-and-forward (AF) scheme is considered. Asymptotic closed-form expressions at high signal-to-noise ratio (SNR) are also derived. Coding gain and diversity order for different combining and selection schemes are further discussed. Numerical results are provided to verify and illustrate the analytical results.

Analytical model and figures of merit for filtered Microwave Photonic Links.

PubMed

Gasulla, Ivana; Capmany, José

2011-09-26

The concept of filtered Microwave Photonic Links is proposed in order to provide the most general and versatile description of complex analog photonic systems. We develop a field propagation model where a global optical filter, characterized by its optical transfer function, embraces all the intermediate optical components in a linear link. We assume a non-monochromatic light source characterized by an arbitrary spectral distribution which has a finite linewidth spectrum and consider both intensity modulation and phase modulation with balanced and single detection. Expressions leading to the computation of the main figures of merit concerning the link gain, noise and intermodulation distortion are provided which, to our knowledge, are not available in the literature. The usefulness of this derivation resides in the capability to directly provide performance criteria results for complex links just by substituting in the overall closed-form formulas the numerical or measured optical transfer function characterizing the link. This theory is presented thus as a potential tool for a wide range of relevant microwave photonic application cases which is extendable to multiport radio over fiber systems. © 2011 Optical Society of America

Nonlinear Schrödinger equations with single power nonlinearity and harmonic potential

NASA Astrophysics Data System (ADS)

Cipolatti, R.; de Macedo Lira, Y.; Trallero-Giner, C.

2018-03-01

We consider a generalized nonlinear Schrödinger equation (GNLS) with a single power nonlinearity of the form λ ≤ft\\vert \\varphi \\right\\vert p , with p  >  0 and λ\\in{R} , in the presence of a harmonic confinement. We report the conditions that p and λ must fulfill for the existence and uniqueness of ground states of the GNLS. We discuss the Cauchy problem and summarize which conditions are required for the nonlinear term λ ≤ft\\vert \\varphi \\right\\vert p to render the ground state solutions orbitally stable. Based on a new variational method we provide exact formulæ for the minimum energy for each index p and the changing range of values of the nonlinear parameter λ. Also, we report an approximate close analytical expression for the ground state energy, performing a comparative analysis of the present variational calculations with those obtained by a generalized Thomas-Fermi approach, and soliton solutions for the respective ranges of p and λ where these solutions can be implemented to describe the minimum energy.

Regularized inversion of controlled source audio-frequency magnetotelluric data in horizontally layered transversely isotropic media

NASA Astrophysics Data System (ADS)

Zhou, Jianmei; Wang, Jianxun; Shang, Qinglong; Wang, Hongnian; Yin, Changchun

2014-04-01

We present an algorithm for inverting controlled source audio-frequency magnetotelluric (CSAMT) data in horizontally layered transversely isotropic (TI) media. The popular inversion method parameterizes the media into a large number of layers which have fixed thickness and only reconstruct the conductivities (e.g. Occam's inversion), which does not enable the recovery of the sharp interfaces between layers. In this paper, we simultaneously reconstruct all the model parameters, including both the horizontal and vertical conductivities and layer depths. Applying the perturbation principle and the dyadic Green's function in TI media, we derive the analytic expression of Fréchet derivatives of CSAMT responses with respect to all the model parameters in the form of Sommerfeld integrals. A regularized iterative inversion method is established to simultaneously reconstruct all the model parameters. Numerical results show that the inverse algorithm, including the depths of the layer interfaces, can significantly improve the inverse results. It can not only reconstruct the sharp interfaces between layers, but also can obtain conductivities close to the true value.

Theoretical limits of localizing 3-D landmarks and features.

PubMed

Rohr, Karl

2007-09-01

In this paper, we analyze the accuracy of estimating the location of 3-D landmarks and characteristic image structures. Based on nonlinear estimation theory, we study the minimal stochastic errors of the position estimate caused by noisy data. Given analytic models of the image intensities, we derive closed-form expressions of the Cramér-Rao bound for different 3-D structures such as 3-D edges, 3-D ridges, 3-D lines, 3-D boxes, and 3-D blobs. It turns out that the precision of localization depends on the noise level, the size of the region-of-interest, the image contrast, the width of the intensity transitions, as well as on other parameters describing the considered image structure. The derived lower bounds can serve as benchmarks and the performance of existing algorithms can be compared with them. To give an impression of the achievable accuracy, numeric examples are presented. Moreover, by experimental investigations, we demonstrate that the derived lower bounds can be achieved by fitting parametric intensity models directly to the image data.

Semidefinite Relaxation-Based Optimization of Multiple-Input Wireless Power Transfer Systems

NASA Astrophysics Data System (ADS)

Lang, Hans-Dieter; Sarris, Costas D.

2017-11-01

An optimization procedure for multi-transmitter (MISO) wireless power transfer (WPT) systems based on tight semidefinite relaxation (SDR) is presented. This method ensures physical realizability of MISO WPT systems designed via convex optimization -- a robust, semi-analytical and intuitive route to optimizing such systems. To that end, the nonconvex constraints requiring that power is fed into rather than drawn from the system via all transmitter ports are incorporated in a convex semidefinite relaxation, which is efficiently and reliably solvable by dedicated algorithms. A test of the solution then confirms that this modified problem is equivalent (tight relaxation) to the original (nonconvex) one and that the true global optimum has been found. This is a clear advantage over global optimization methods (e.g. genetic algorithms), where convergence to the true global optimum cannot be ensured or tested. Discussions of numerical results yielded by both the closed-form expressions and the refined technique illustrate the importance and practicability of the new method. It, is shown that this technique offers a rigorous optimization framework for a broad range of current and emerging WPT applications.

Modelling mass diffusion for a multi-layer sphere immersed in a semi-infinite medium: application to drug delivery.

PubMed

Carr, Elliot J; Pontrelli, Giuseppe

2018-04-12

We present a general mechanistic model of mass diffusion for a composite sphere placed in a large ambient medium. The multi-layer problem is described by a system of diffusion equations coupled via interlayer boundary conditions such as those imposing a finite mass resistance at the external surface of the sphere. While the work is applicable to the generic problem of heat or mass transfer in a multi-layer sphere, the analysis and results are presented in the context of drug kinetics for desorbing and absorbing spherical microcapsules. We derive an analytical solution for the concentration in the sphere and in the surrounding medium that avoids any artificial truncation at a finite distance. The closed-form solution in each concentric layer is expressed in terms of a suitably-defined inverse Laplace transform that can be evaluated numerically. Concentration profiles and drug mass curves in the spherical layers and in the external environment are presented and the dependency of the solution on the mass transfer coefficient at the surface of the sphere analyzed. Copyright © 2018 Elsevier Inc. All rights reserved.

Effects of film injection angle on turbine vane cooling

NASA Technical Reports Server (NTRS)

Gauntner, J. W.

1977-01-01

Film ejection from discrete holes in the suction surface of a turbine vane was studied for hole axes (1) slanted 30 deg to the surface in the streamwise direction and (2) slanted 30 deg to the surface and 45 deg from the streamwise direction toward the hub. The holes were near the throat area in a five-row staggered array with 8-diameter spacing. Mass flux ratios were as high as 1.2. The data were obtained in an annular sector cascade at conditions where both the ratio of the boundary layer momentum thickness-to-hole diameter and the momentum thickness Reynolds number were typical of an advanced turbofan engine at both takeoff and cruise. Wall temperatures were measured downstream of each of the rows of holes. Results of this study are expressed as a comparison of cooling effectiveness between the in-line angle injection and the compound-angle injection as a function of mass flux ratio. These heat transfer results are also compared with the results of a referenced flow visualization study. Also included is a closed-form analytical solution for temperature within the film cooled wall.

Large Angle Unsteady Aerodynamic Theory of a Flat Plate

NASA Astrophysics Data System (ADS)

Manar, Field; Jones, Anya

2016-11-01

A purely analytical approach is taken for the evaluation of the unsteady loads on a flat plate. This allows for an extremely low cost theoretical prediction of the plate loads in the style of Wagner and Theodorsen, without making the assumption of small angle of attack or small disturbance flow. The forces and moments are evaluated using the time rate of change of fluid momentum, expressed as an integral of the vorticity field. The flow is taken as inviscid and incompressible with isolated vorticity bound to the plate and in the shed wake. The bound vorticity distribution on the plate is solved exactly using conformal mapping of the plate to a cylinder. In keeping with the original assumption of Wagner, the wake vorticity is assumed to remain stationary in an inertial reference frame and convection is disregarded. Formulation in this manner allows for a closed form solution of Wagner's problem valid at all angles of attack. Separation from the leading edge of the plate can also be included to further increase the fidelity of the model at high angles.

Effects of shear flow on phase nucleation and crystallization.

PubMed

Mura, Federica; Zaccone, Alessio

2016-04-01

Classical nucleation theory offers a good framework for understanding the common features of new phase formation processes in metastable homogeneous media at rest. However, nucleation processes in liquids are ubiquitously affected by hydrodynamic flow, and there is no satisfactory understanding of whether shear promotes or slows down the nucleation process. We developed a classical nucleation theory for sheared systems starting from the molecular level of the Becker-Doering master kinetic equation and we analytically derived a closed-form expression for the nucleation rate. The theory accounts for the effect of flow-mediated transport of molecules to the nucleus of the new phase, as well as for the mechanical deformation imparted to the nucleus by the flow field. The competition between flow-induced molecular transport, which accelerates nucleation, and flow-induced nucleus straining, which lowers the nucleation rate by increasing the nucleation energy barrier, gives rise to a marked nonmonotonic dependence of the nucleation rate on the shear rate. The theory predicts an optimal shear rate at which the nucleation rate is one order of magnitude larger than in the absence of flow.

The frequency-dependent response of single aerosol particles to vapour phase oscillations and its application in measuring diffusion coefficients

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

Preston, Thomas C.; Davies, James F.; Wilson, Kevin R.

A new method for measuring diffusion in the condensed phase of single aerosol particles is proposed and demonstrated. The technique is based on the frequency-dependent response of a binary particle to oscillations in the vapour phase of one of its chemical components. Here, we discuss how this physical situation allows for what would typically be a non-linear boundary value problem to be approximately reduced to a linear boundary value problem. For the case of aqueous aerosol particles, we investigate the accuracy of the closed-form analytical solution to this linear problem through a comparison with the numerical solution of the fullmore » problem. Then, using experimentally measured whispering gallery modes to track the frequency-dependent response of aqueous particles to relative humidity oscillations, we determine diffusion coefficients as a function of water activity. The measured diffusion coefficients are compared to previously reported values found using the two common experiments: (i) the analysis of the sorption/desorption of water from a particle after a step-wise change to the surrounding relative humidity and (ii) the isotopic exchange of water between a particle and the vapour phase. The technique presented here has two main strengths: first, when compared to the sorption/desorption experiment, it does not require the numerical evaluation of a boundary value problem during the fitting process as a closed-form expression is available. Second, when compared to the isotope exchange experiment, it does not require the use of labeled molecules. Therefore, the frequency-dependent experiment retains the advantages of these two commonly used methods but does not suffer from their drawbacks.« less

The frequency-dependent response of single aerosol particles to vapour phase oscillations and its application in measuring diffusion coefficients

DOE PAGES

Preston, Thomas C.; Davies, James F.; Wilson, Kevin R.

2017-01-13

A new method for measuring diffusion in the condensed phase of single aerosol particles is proposed and demonstrated. The technique is based on the frequency-dependent response of a binary particle to oscillations in the vapour phase of one of its chemical components. Here, we discuss how this physical situation allows for what would typically be a non-linear boundary value problem to be approximately reduced to a linear boundary value problem. For the case of aqueous aerosol particles, we investigate the accuracy of the closed-form analytical solution to this linear problem through a comparison with the numerical solution of the fullmore » problem. Then, using experimentally measured whispering gallery modes to track the frequency-dependent response of aqueous particles to relative humidity oscillations, we determine diffusion coefficients as a function of water activity. The measured diffusion coefficients are compared to previously reported values found using the two common experiments: (i) the analysis of the sorption/desorption of water from a particle after a step-wise change to the surrounding relative humidity and (ii) the isotopic exchange of water between a particle and the vapour phase. The technique presented here has two main strengths: first, when compared to the sorption/desorption experiment, it does not require the numerical evaluation of a boundary value problem during the fitting process as a closed-form expression is available. Second, when compared to the isotope exchange experiment, it does not require the use of labeled molecules. Therefore, the frequency-dependent experiment retains the advantages of these two commonly used methods but does not suffer from their drawbacks.« less

Modeling and control of non-square MIMO system using relay feedback.

PubMed

Kalpana, D; Thyagarajan, T; Gokulraj, N

2015-11-01

This paper proposes a systematic approach for the modeling and control of non-square MIMO systems in time domain using relay feedback. Conventionally, modeling, selection of the control configuration and controller design of non-square MIMO systems are performed using input/output information of direct loop, while the output of undesired responses that bears valuable information on interaction among the loops are not considered. However, in this paper, the undesired response obtained from relay feedback test is also taken into consideration to extract the information about the interaction between the loops. The studies are performed on an Air Path Scheme of Turbocharged Diesel Engine (APSTDE) model, which is a typical non-square MIMO system, with input and output variables being 3 and 2 respectively. From the relay test response, the generalized analytical expressions are derived and these analytical expressions are used to estimate unknown system parameters and also to evaluate interaction measures. The interaction is analyzed by using Block Relative Gain (BRG) method. The model thus identified is later used to design appropriate controller to carry out closed loop studies. Closed loop simulation studies were performed for both servo and regulatory operations. Integral of Squared Error (ISE) performance criterion is employed to quantitatively evaluate performance of the proposed scheme. The usefulness of the proposed method is demonstrated on a lab-scale Two-Tank Cylindrical Interacting System (TTCIS), which is configured as a non-square system. Copyright © 2015 ISA. Published by Elsevier Ltd. All rights reserved.

Closed Form Solutions for Unsteady Free Convection Flow of a Second Grade Fluid over an Oscillating Vertical Plate

PubMed Central

Ali, Farhad; Khan, Ilyas; Shafie, Sharidan

2014-01-01

Closed form solutions for unsteady free convection flows of a second grade fluid near an isothermal vertical plate oscillating in its plane using the Laplace transform technique are established. Expressions for velocity and temperature are obtained and displayed graphically for different values of Prandtl number Pr, thermal Grashof number Gr, viscoelastic parameter α, phase angle ωτ and time τ. Numerical values of skin friction τ 0 and Nusselt number Nu are shown in tables. Some well-known solutions in literature are reduced as the limiting cases of the present solutions. PMID:24551033

On the closed form mechanistic modeling of milling: Specific cutting energy, torque, and power

NASA Astrophysics Data System (ADS)

Bayoumi, A. E.; Yücesan, G.; Hutton, D. V.

1994-02-01

Specific energy in metal cutting, defined as the energy expended in removing a unit volume of workpiece material, is formulated and determined using a previously developed closed form mechanistic force model for milling operations. Cutting power is computed from the cutting torque, cutting force, kinematics of the cutter, and the volumetric material removal rate. Closed form expressions for specific cutting energy were formulated and found to be functions of the process parameters: pressure and friction for both rake and flank surfaces and chip flow angle at the rake face of the tool. Friction is found to play a very important role in cutting torque and power. Experiments were carried out to determine the effects of feedrate, cutting speed, workpiece material, and flank wear land width on specific cutting energy. It was found that the specific cutting energy increases with a decrease in the chip thickness and with an increase in flank wear land.

Cumulant-based expressions for the multibody terms for the correlation between local and electrostatic interactions in the united-residue force field

NASA Astrophysics Data System (ADS)

Liwo, Adam; Czaplewski, Cezary; Pillardy, Jarosław; Scheraga, Harold A.

2001-08-01

A general method to derive site-site or united-residue potentials is presented. The basic principle of the method is the separation of the degrees of freedom of a system into the primary and secondary ones. The primary degrees of freedom describe the basic features of the system, while the secondary ones are averaged over when calculating the potential of mean force, which is hereafter referred to as the restricted free energy (RFE) function. The RFE can be factored into one-, two-, and multibody terms, using the cluster-cumulant expansion of Kubo. These factors can be assigned the functional forms of the corresponding lowest-order nonzero generalized cumulants, which can, in most cases, be evaluated analytically, after making some simplifying assumptions. This procedure to derive coarse-grain force fields is very valuable when applied to multibody terms, whose functional forms are hard to deduce in another way (e.g., from structural databases). After the functional forms have been derived, they can be parametrized based on the RFE surfaces of model systems obtained from all-atom models or on the statistics derived from structural databases. The approach has been applied to our united-residue force field for proteins. Analytical expressions were derived for the multibody terms pertaining to the correlation between local and electrostatic interactions within the polypeptide backbone; these expressions correspond to up to sixth-order terms in the cumulant expansion of the RFE. These expressions were subsequently parametrized by fitting to the RFEs of selected peptide fragments, calculated with the empirical conformational energy program for peptides force field. The new multibody terms enable not only the heretofore predictable α-helical segments, but also regular β-sheets, to form as the lowest-energy structures, as assessed by test calculations on a model helical protein A, as well as a model 20-residue polypeptide (betanova); the latter was not possible without introducing these new terms.

Closed-Form Jensen-Renyi Divergence for Mixture of Gaussians and Applications to Group-Wise Shape Registration*

PubMed Central

Wang, Fei; Syeda-Mahmood, Tanveer; Vemuri, Baba C.; Beymer, David; Rangarajan, Anand

2010-01-01

In this paper, we propose a generalized group-wise non-rigid registration strategy for multiple unlabeled point-sets of unequal cardinality, with no bias toward any of the given point-sets. To quantify the divergence between the probability distributions – specifically Mixture of Gaussians – estimated from the given point sets, we use a recently developed information-theoretic measure called Jensen-Renyi (JR) divergence. We evaluate a closed-form JR divergence between multiple probabilistic representations for the general case where the mixture models differ in variance and the number of components. We derive the analytic gradient of the divergence measure with respect to the non-rigid registration parameters, and apply it to numerical optimization of the group-wise registration, leading to a computationally efficient and accurate algorithm. We validate our approach on synthetic data, and evaluate it on 3D cardiac shapes. PMID:20426043

Closed-form Jensen-Renyi divergence for mixture of Gaussians and applications to group-wise shape registration.

PubMed

Wang, Fei; Syeda-Mahmood, Tanveer; Vemuri, Baba C; Beymer, David; Rangarajan, Anand

2009-01-01

In this paper, we propose a generalized group-wise non-rigid registration strategy for multiple unlabeled point-sets of unequal cardinality, with no bias toward any of the given point-sets. To quantify the divergence between the probability distributions--specifically Mixture of Gaussians--estimated from the given point sets, we use a recently developed information-theoretic measure called Jensen-Renyi (JR) divergence. We evaluate a closed-form JR divergence between multiple probabilistic representations for the general case where the mixture models differ in variance and the number of components. We derive the analytic gradient of the divergence measure with respect to the non-rigid registration parameters, and apply it to numerical optimization of the group-wise registration, leading to a computationally efficient and accurate algorithm. We validate our approach on synthetic data, and evaluate it on 3D cardiac shapes.

Analysis of Tank PMD Rewetting Following Thrust Resettling

NASA Astrophysics Data System (ADS)

Weislogel, M. M.; Sala, M. A.; Collicott, S. H.

2002-10-01

Recent investigations have successfully demonstrated closed-form analytical solutions of spontaneous capillary flows in idealized cylindrical containers with interior corners. In this report, the theory is extended and applied to complex containers modeling spacecraft fuel tanks employing propellant management devices (PMDs). The specific problem investigated is one of spontaneous rewetting of a typical partially filled liquid fuel/cryogen tank with PMD after thrust resettling. The transients of this flow impact the logistics of orbital maneuvers and potentially tank thermal control. The general procedure to compute the initial condition (mean radius of curvature for the interface) for the closed-form transient flows is first outlined then solved for several 'complex' cylindrical tanks exhibiting symmetry. The utility and limitations of the technique as a design tool are discussed in a summary, which also highlights comparisons with NASA flight data of a model propellant tank with PMD.

Nanoscale phase transition behavior of shape memory alloys — closed form solution of 1D effective modelling

NASA Astrophysics Data System (ADS)

Li, M. P.; Sun, Q. P.

2018-01-01

We investigate the roles of grain size (lg) and grain boundary thickness (lb) on the stress-induced phase transition (PT) behaviors of nanocrystalline shape memory alloys (SMAs) by using a Core-shell type "crystallite-amorphous composite" model. A non-dimensionalized length scale lbarg(=lg /lb) is identified as the governing parameter which is indicative of the energy competition between the crystallite and the grain boundary. Closed form analytical solutions of a reduced effective 1D model with embedded microstructure length scales of lg and lb are presented in this paper. It is shown that, with lbarg reduction, the energy of the elastic non-transformable grain boundary will gradually become dominant in the phase transition process, and eventually bring fundamental changes of the deformation behaviors: breakdown of two-phase coexistence and vanishing of superelastic hysteresis. The predictions are supported by experimental data of nanocrystalline NiTi SMAs.

Analysis of Tank PMD Rewetting Following Thrust Resettling

NASA Technical Reports Server (NTRS)

Weislogel, M. M.; Sala, M. A.; Collicott, S. H.; Rame, Enrique (Technical Monitor)

2002-01-01

Recent investigations have successfully demonstrated closed-form analytical solutions of spontaneous capillary flows in idealized cylindrical containers with interior corners. In this report, the theory is extended and applied to complex containers modeling spacecraft fuel tanks employing propellant management devices (PMDs). The specific problem investigated is one of spontaneous rewetting of a typical partially filled liquid fuel/cryogen tank with PMD after thrust resettling. The transients of this flow impact the logistics of orbital maneuvers and potentially tank thermal control. The general procedure to compute the initial condition (mean radius of curvature for the interface) for the closed-form transient flows is first outlined then solved for several 'complex' cylindrical tanks exhibiting symmetry. The utility and limitations of the technique as a design tool are discussed in a summary, which also highlights comparisons with NASA flight data of a model propellant tank with PMD.

Self-induced transparency of an extremely short pulse

NASA Technical Reports Server (NTRS)

Lee, C. T.

1973-01-01

An extremely short pulse propagation in a resonant medium is properly described by a closed form steady-state analytic solution. The usual slowly varying envelope approximation (SVEA) is not made. Instead, different assumptions with respect to pulse speed and pulse duration are used, and any possible nonresonant loss is ignored. This study indicates that the results obtained by the SVEA approach are much better than they have been intuitively expected to be.

Existence of small loops in a bifurcation diagram near degenerate eigenvalues

NASA Astrophysics Data System (ADS)

Hmidi, Taoufik; Renault, Coralie

2017-10-01

In this paper, we study the global structure of a bifurcation diagram for rotating doubly connected patches near a degenerate case for incompressible Euler equations. We show that branches with the same symmetry merge, forming a small loop, provided that they are close enough. This gives an analytical proof for the numerical observations conducted in the recent work by de la Hoz et al (2016 SIAM J. Math. Anal. 48 1892-928).

Radiation Pattern of Chair Armed Microstrip Antenna

NASA Astrophysics Data System (ADS)

Mishra, Rabindra Kishore; Sahu, Kumar Satyabrat

2016-12-01

This work analyzes planar antenna conformable to chair arm shaped surfaces for WLAN application. Closed form expressions for its radiation pattern are developed and validated using measurements on prototype and commercial EM code at 2.4 GHz.

A general formula for Rayleigh-Schroedinger perturbation energy utilizing a power series expansion of the quantum mechanical Hamiltonian

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

Herbert, J.M.

1997-02-01

Perturbation theory has long been utilized by quantum chemists as a method for approximating solutions to the Schroedinger equation. Perturbation treatments represent a system`s energy as a power series in which each additional term further corrects the total energy; it is therefore convenient to have an explicit formula for the nth-order energy correction term. If all perturbations are collected into a single Hamiltonian operator, such a closed-form expression for the nth-order energy correction is well known; however, use of a single perturbed Hamiltonian often leads to divergent energy series, while superior convergence behavior is obtained by expanding the perturbed Hamiltonianmore » in a power series. This report presents a closed-form expression for the nth-order energy correction obtained using Rayleigh-Schroedinger perturbation theory and a power series expansion of the Hamiltonian.« less

Closed-form analytical solutions incorporating pumping and tidal effects in various coastal aquifer systems

NASA Astrophysics Data System (ADS)

Wang, Chaoyue; Li, Hailong; Wan, Li; Wang, Xusheng; Jiang, Xiaowei

2014-07-01

Pumping wells are common in coastal aquifers affected by tides. Here we present analytical solutions of groundwater table or head variations during a constant rate pumping from a single, fully-penetrating well in coastal aquifer systems comprising an unconfined aquifer, a confined aquifer and semi-permeable layer between them. The unconfined aquifer terminates at the coastline (or river bank) and the other two layers extend under tidal water (sea or tidal river) for a certain distance L. Analytical solutions are derived for 11 reasonable combinations of different situations of the L-value (zero, finite, and infinite), of the middle layer's permeability (semi-permeable and impermeable), of the boundary condition at the aquifer's submarine terminal (Dirichlet describing direct connection with seawater and no-flow describing the existence of an impermeable capping), and of the tidal water body (sea and tidal river). Solutions are discussed with application examples in fitting field observations and parameter estimations.

Note on use of slope diffraction coefficients for aperture antennas on finite ground planes

NASA Technical Reports Server (NTRS)

Cockrell, C. R.; Beck, F. B.

1995-01-01

The use of slope diffraction coefficients along with regular diffraction coefficients for calculating the radiation patterns of aperture antennas in a finite ground plane is investigated. Explicit expressions for regular diffraction coefficients and slope diffraction coefficients are presented. The expressions for the incident magnetic field in terms of the magnetic current in the aperture are given. The slope of the incident magnetic field is calculated and closed form expressions are presented.

Production of Computer Animated Movies for Educational Purposes.

ERIC Educational Resources Information Center

Elberg, H. H.

A detailed account is given in this paper of the procedures and the equipment used in producing six computer-animated instructional movies. First, the sequence of events were described in a script, which, together with the analytical expressions that were dealt with, formed the basis of a program. Then, the program was run on a computer and the…

Emergent rogue wave structures and statistics in spontaneous modulation instability.

PubMed

Toenger, Shanti; Godin, Thomas; Billet, Cyril; Dias, Frédéric; Erkintalo, Miro; Genty, Goëry; Dudley, John M

2015-05-20

The nonlinear Schrödinger equation (NLSE) is a seminal equation of nonlinear physics describing wave packet evolution in weakly-nonlinear dispersive media. The NLSE is especially important in understanding how high amplitude "rogue waves" emerge from noise through the process of modulation instability (MI) whereby a perturbation on an initial plane wave can evolve into strongly-localised "breather" or "soliton on finite background (SFB)" structures. Although there has been much study of such structures excited under controlled conditions, there remains the open question of how closely the analytic solutions of the NLSE actually model localised structures emerging in noise-seeded MI. We address this question here using numerical simulations to compare the properties of a large ensemble of emergent peaks in noise-seeded MI with the known analytic solutions of the NLSE. Our results show that both elementary breather and higher-order SFB structures are observed in chaotic MI, with the characteristics of the noise-induced peaks clustering closely around analytic NLSE predictions. A significant conclusion of our work is to suggest that the widely-held view that the Peregrine soliton forms a rogue wave prototype must be revisited. Rather, we confirm earlier suggestions that NLSE rogue waves are most appropriately identified as collisions between elementary SFB solutions.

Emergent rogue wave structures and statistics in spontaneous modulation instability

PubMed Central

Toenger, Shanti; Godin, Thomas; Billet, Cyril; Dias, Frédéric; Erkintalo, Miro; Genty, Goëry; Dudley, John M.

2015-01-01

The nonlinear Schrödinger equation (NLSE) is a seminal equation of nonlinear physics describing wave packet evolution in weakly-nonlinear dispersive media. The NLSE is especially important in understanding how high amplitude “rogue waves” emerge from noise through the process of modulation instability (MI) whereby a perturbation on an initial plane wave can evolve into strongly-localised “breather” or “soliton on finite background (SFB)” structures. Although there has been much study of such structures excited under controlled conditions, there remains the open question of how closely the analytic solutions of the NLSE actually model localised structures emerging in noise-seeded MI. We address this question here using numerical simulations to compare the properties of a large ensemble of emergent peaks in noise-seeded MI with the known analytic solutions of the NLSE. Our results show that both elementary breather and higher-order SFB structures are observed in chaotic MI, with the characteristics of the noise-induced peaks clustering closely around analytic NLSE predictions. A significant conclusion of our work is to suggest that the widely-held view that the Peregrine soliton forms a rogue wave prototype must be revisited. Rather, we confirm earlier suggestions that NLSE rogue waves are most appropriately identified as collisions between elementary SFB solutions. PMID:25993126

The asymptotic form of non-global logarithms, black disc saturation, and gluonic deserts

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

Neill, Duff

Here, we develop an asymptotic perturbation theory for the large logarithmic behavior of the non-linear integro-differential equation describing the soft correlations of QCD jet measurements, the Banfi-Marchesini-Smye (BMS) equation. Furthermore, this equation captures the late-time evolution of radiating color dipoles after a hard collision. This allows us to prove that at large values of the control variable (the non-global logarithm, a function of the infra-red energy scales associated with distinct hard jets in an event), the distribution has a gaussian tail. We also compute the decay width analytically, giving a closed form expression, and find it to be jet geometrymore » independent, up to the number of legs of the dipole in the active jet. By enabling the asymptotic expansion we find that the perturbative seed is correct; we perturb around an anzats encoding formally no real emissions, an intuition motivated by the buffer region found in jet dynamics. This must be supplemented with the correct application of the BFKL approximation to the BMS equation in collinear limits. Comparing to the asymptotics of the conformally related evolution equation encountered in small-x physics, the Balitisky-Kovchegov (BK) equation, we find that the asymptotic form of the non-global logarithms directly maps to the black-disc unitarity limit of the BK equation, despite the contrasting physical pictures. Indeed, we recover the equations of saturation physics in the final state dynamics of QCD.« less

A classifier neural network for rotordynamic systems

NASA Astrophysics Data System (ADS)

Ganesan, R.; Jionghua, Jin; Sankar, T. S.

1995-07-01

A feedforward backpropagation neural network is formed to identify the stability characteristic of a high speed rotordynamic system. The principal focus resides in accounting for the instability due to the bearing clearance effects. The abnormal operating condition of 'normal-loose' Coulomb rub, that arises in units supported by hydrodynamic bearings or rolling element bearings, is analysed in detail. The multiple-parameter stability problem is formulated and converted to a set of three-parameter algebraic inequality equations. These three parameters map the wider range of physical parameters of commonly-used rotordynamic systems into a narrow closed region, that is used in the supervised learning of the neural network. A binary-type state of the system is expressed through these inequalities that are deduced from the analytical simulation of the rotor system. Both the hidden layer as well as functional-link networks are formed and the superiority of the functional-link network is established. Considering the real time interpretation and control of the rotordynamic system, the network reliability and the learning time are used as the evaluation criteria to assess the superiority of the functional-link network. This functional-link network is further trained using the parameter values of selected rotor systems, and the classifier network is formed. The success rate of stability status identification is obtained to assess the potentials of this classifier network. The classifier network is shown that it can also be used, for control purposes, as an 'advisory' system that suggests the optimum way of parameter adjustment.

The asymptotic form of non-global logarithms, black disc saturation, and gluonic deserts

DOE PAGES

Neill, Duff

2017-01-25

Here, we develop an asymptotic perturbation theory for the large logarithmic behavior of the non-linear integro-differential equation describing the soft correlations of QCD jet measurements, the Banfi-Marchesini-Smye (BMS) equation. Furthermore, this equation captures the late-time evolution of radiating color dipoles after a hard collision. This allows us to prove that at large values of the control variable (the non-global logarithm, a function of the infra-red energy scales associated with distinct hard jets in an event), the distribution has a gaussian tail. We also compute the decay width analytically, giving a closed form expression, and find it to be jet geometrymore » independent, up to the number of legs of the dipole in the active jet. By enabling the asymptotic expansion we find that the perturbative seed is correct; we perturb around an anzats encoding formally no real emissions, an intuition motivated by the buffer region found in jet dynamics. This must be supplemented with the correct application of the BFKL approximation to the BMS equation in collinear limits. Comparing to the asymptotics of the conformally related evolution equation encountered in small-x physics, the Balitisky-Kovchegov (BK) equation, we find that the asymptotic form of the non-global logarithms directly maps to the black-disc unitarity limit of the BK equation, despite the contrasting physical pictures. Indeed, we recover the equations of saturation physics in the final state dynamics of QCD.« less

Almost analytical Karhunen-Loeve representation of irregular waves based on the prolate spheroidal wave functions

NASA Astrophysics Data System (ADS)

Lee, Gibbeum; Cho, Yeunwoo

2017-11-01

We present an almost analytical new approach to solving the matrix eigenvalue problem or the integral equation in Karhunen-Loeve (K-L) representation of random data such as irregular ocean waves. Instead of solving this matrix eigenvalue problem purely numerically, which may suffer from the computational inaccuracy for big data, first, we consider a pair of integral and differential equations, which are related to the so-called prolate spheroidal wave functions (PSWF). For the PSWF differential equation, the pair of the eigenvectors (PSWF) and eigenvalues can be obtained from a relatively small number of analytical Legendre functions. Then, the eigenvalues in the PSWF integral equation are expressed in terms of functional values of the PSWF and the eigenvalues of the PSWF differential equation. Finally, the analytically expressed PSWFs and the eigenvalues in the PWSF integral equation are used to form the kernel matrix in the K-L integral equation for the representation of exemplary wave data; ordinary irregular waves and rogue waves. We found that the present almost analytical method is better than the conventional data-independent Fourier representation and, also, the conventional direct numerical K-L representation in terms of both accuracy and computational cost. This work was supported by the National Research Foundation of Korea (NRF). (NRF-2017R1D1A1B03028299).

Modeling and control of flexible space platforms with articulated payloads

NASA Technical Reports Server (NTRS)

Graves, Philip C.; Joshi, Suresh M.

1989-01-01

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

A series solution for horizontal infiltration in an initially dry aquifer

NASA Astrophysics Data System (ADS)

Furtak-Cole, Eden; Telyakovskiy, Aleksey S.; Cooper, Clay A.

2018-06-01

The porous medium equation (PME) is a generalization of the traditional Boussinesq equation for hydraulic conductivity as a power law function of height. We analyze the horizontal recharge of an initially dry unconfined aquifer of semi-infinite extent, as would be found in an aquifer adjacent a rising river. If the water level can be modeled as a power law function of time, similarity variables can be introduced and the original problem can be reduced to a boundary value problem for a nonlinear ordinary differential equation. The position of the advancing front is not known ahead of time and must be found in the process of solution. We present an analytical solution in the form of a power series, with the coefficients of the series given by a recurrence relation. The analytical solution compares favorably with a highly accurate numerical solution, and only a small number of terms of the series are needed to achieve high accuracy in the scenarios considered here. We also conduct a series of physical experiments in an initially dry wedged Hele-Shaw cell, where flow is modeled by a special form of the PME. Our analytical solution closely matches the hydraulic head profiles in the Hele-Shaw cell experiment.

Comparison of Three Methods for Wind Turbine Capacity Factor Estimation

PubMed Central

Ditkovich, Y.; Kuperman, A.

2014-01-01

Three approaches to calculating capacity factor of fixed speed wind turbines are reviewed and compared using a case study. The first “quasiexact” approach utilizes discrete wind raw data (in the histogram form) and manufacturer-provided turbine power curve (also in discrete form) to numerically calculate the capacity factor. On the other hand, the second “analytic” approach employs a continuous probability distribution function, fitted to the wind data as well as continuous turbine power curve, resulting from double polynomial fitting of manufacturer-provided power curve data. The latter approach, while being an approximation, can be solved analytically thus providing a valuable insight into aspects, affecting the capacity factor. Moreover, several other merits of wind turbine performance may be derived based on the analytical approach. The third “approximate” approach, valid in case of Rayleigh winds only, employs a nonlinear approximation of the capacity factor versus average wind speed curve, only requiring rated power and rotor diameter of the turbine. It is shown that the results obtained by employing the three approaches are very close, enforcing the validity of the analytically derived approximations, which may be used for wind turbine performance evaluation. PMID:24587755

Eye movement dysfunction in first-degree relatives of patients with schizophrenia: a meta-analytic evaluation of candidate endophenotypes.

PubMed

Calkins, Monica E; Iacono, William G; Ones, Deniz S

2008-12-01

Several forms of eye movement dysfunction (EMD) are regarded as promising candidate endophenotypes of schizophrenia. Discrepancies in individual study results have led to inconsistent conclusions regarding particular aspects of EMD in relatives of schizophrenia patients. To quantitatively evaluate and compare the candidacy of smooth pursuit, saccade and fixation deficits in first-degree biological relatives, we conducted a set of meta-analytic investigations. Among 18 measures of EMD, memory-guided saccade accuracy and error rate, global smooth pursuit dysfunction, intrusive saccades during fixation, antisaccade error rate and smooth pursuit closed-loop gain emerged as best differentiating relatives from controls (standardized mean differences ranged from .46 to .66), with no significant differences among these measures. Anticipatory saccades, but no other smooth pursuit component measures were also increased in relatives. Visually-guided reflexive saccades were largely normal. Moderator analyses examining design characteristics revealed few variables affecting the magnitude of the meta-analytically observed effects. Moderate effect sizes of relatives v. controls in selective aspects of EMD supports their endophenotype potential. Future work should focus on facilitating endophenotype utility through attention to heterogeneity of EMD performance, relationships among forms of EMD, and application in molecular genetics studies.

A maximum power point tracking algorithm for photovoltaic applications

NASA Astrophysics Data System (ADS)

Nelatury, Sudarshan R.; Gray, Robert

2013-05-01

The voltage and current characteristic of a photovoltaic (PV) cell is highly nonlinear and operating a PV cell for maximum power transfer has been a challenge for a long time. Several techniques have been proposed to estimate and track the maximum power point (MPP) in order to improve the overall efficiency of a PV panel. A strategic use of the mean value theorem permits obtaining an analytical expression for a point that lies in a close neighborhood of the true MPP. But hitherto, an exact solution in closed form for the MPP is not published. This problem can be formulated analytically as a constrained optimization, which can be solved using the Lagrange method. This method results in a system of simultaneous nonlinear equations. Solving them directly is quite difficult. However, we can employ a recursive algorithm to yield a reasonably good solution. In graphical terms, suppose the voltage current characteristic and the constant power contours are plotted on the same voltage current plane, the point of tangency between the device characteristic and the constant power contours is the sought for MPP. It is subject to change with the incident irradiation and temperature and hence the algorithm that attempts to maintain the MPP should be adaptive in nature and is supposed to have fast convergence and the least misadjustment. There are two parts in its implementation. First, one needs to estimate the MPP. The second task is to have a DC-DC converter to match the given load to the MPP thus obtained. Availability of power electronics circuits made it possible to design efficient converters. In this paper although we do not show the results from a real circuit, we use MATLAB to obtain the MPP and a buck-boost converter to match the load. Under varying conditions of load resistance and irradiance we demonstrate MPP tracking in case of a commercially available solar panel MSX-60. The power electronics circuit is simulated by PSIM software.

Characterizing ISI and sub-threshold membrane potential distributions: Ensemble of IF neurons with random squared-noise intensity.

PubMed

Kumar, Sanjeev; Karmeshu

2018-04-01

A theoretical investigation is presented that characterizes the emerging sub-threshold membrane potential and inter-spike interval (ISI) distributions of an ensemble of IF neurons that group together and fire together. The squared-noise intensity σ 2 of the ensemble of neurons is treated as a random variable to account for the electrophysiological variations across population of nearly identical neurons. Employing superstatistical framework, both ISI distribution and sub-threshold membrane potential distribution of neuronal ensemble are obtained in terms of generalized K-distribution. The resulting distributions exhibit asymptotic behavior akin to stretched exponential family. Extensive simulations of the underlying SDE with random σ 2 are carried out. The results are found to be in excellent agreement with the analytical results. The analysis has been extended to cover the case corresponding to independent random fluctuations in drift in addition to random squared-noise intensity. The novelty of the proposed analytical investigation for the ensemble of IF neurons is that it yields closed form expressions of probability distributions in terms of generalized K-distribution. Based on a record of spiking activity of thousands of neurons, the findings of the proposed model are validated. The squared-noise intensity σ 2 of identified neurons from the data is found to follow gamma distribution. The proposed generalized K-distribution is found to be in excellent agreement with that of empirically obtained ISI distribution of neuronal ensemble. Copyright © 2018 Elsevier B.V. All rights reserved.

A semi-analytical beam model for the vibration of railway tracks

NASA Astrophysics Data System (ADS)

Kostovasilis, D.; Thompson, D. J.; Hussein, M. F. M.

2017-04-01

The high frequency dynamic behaviour of railway tracks, in both vertical and lateral directions, strongly affects the generation of rolling noise as well as other phenomena such as rail corrugation. An improved semi-analytical model of a beam on an elastic foundation is introduced that accounts for the coupling of the vertical and lateral vibration. The model includes the effects of cross-section asymmetry, shear deformation, rotational inertia and restrained warping. Consideration is given to the fact that the loads at the rail head, as well as those exerted by the railpads at the rail foot, may not act through the centroid of the section. The response is evaluated for a harmonic load and the solution is obtained in the wavenumber domain. Results are presented as dispersion curves for free and supported rails and are validated with the aid of a Finite Element (FE) and a waveguide finite element (WFE) model. Closed form expressions are derived for the forced response, and validated against the WFE model. Track mobilities and decay rates are presented to assess the potential implications for rolling noise and the influence of the various sources of vertical-lateral coupling. Comparison is also made with measured data. Overall, the model presented performs very well, especially for the lateral vibration, although it does not contain the high frequency cross-section deformation modes. The most significant effects on the response are shown to be the inclusion of torsion and foundation eccentricity, which mainly affect the lateral response.

Statistical Properties of Real-Time Amplitude Estimate of Harmonics Affected by Frequency Instability

NASA Astrophysics Data System (ADS)

Bellan, Diego; Pignari, Sergio A.

2016-07-01

This work deals with the statistical characterization of real-time digital measurement of the amplitude of harmonics affected by frequency instability. In fact, in modern power systems both the presence of harmonics and frequency instability are well-known and widespread phenomena mainly due to nonlinear loads and distributed generation, respectively. As a result, real-time monitoring of voltage/current frequency spectra is of paramount importance as far as power quality issues are addressed. Within this framework, a key point is that in many cases real-time continuous monitoring prevents the application of sophisticated algorithms to extract all the information from the digitized waveforms because of the required computational burden. In those cases only simple evaluations such as peak search of discrete Fourier transform are implemented. It is well known, however, that a slight change in waveform frequency results in lack of sampling synchronism and uncertainty in amplitude estimate. Of course the impact of this phenomenon increases with the order of the harmonic to be measured. In this paper an approximate analytical approach is proposed in order to describe the statistical properties of the measured magnitude of harmonics affected by frequency instability. By providing a simplified description of the frequency behavior of the windows used against spectral leakage, analytical expressions for mean value, variance, cumulative distribution function, and probability density function of the measured harmonics magnitude are derived in closed form as functions of waveform frequency treated as a random variable.

Extending the Constant Power Speed Range of the Brushless DC Motor through Dual Mode Inverter Control -- Part I: Theory and Simulation

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

Lawler, J.S.

2001-10-29

An inverter topology and control scheme has been developed that can drive low-inductance, surface-mounted permanent magnet motors over the wide constant power speed range required in electric vehicle applications. This new controller is called the dual-mode inverter control (DMIC) [1]. The DMIC can drive either the Permanent Magnet Synchronous Machine (PMSM) with sinusoidal back emf, or the brushless dc machine (BDCM) with trapezoidal emf in the motoring and regenerative braking modes. In this paper we concentrate on the BDCM under high-speed motoring conditions. Simulation results show that if all motor and inverter loss mechanisms are neglected, the constant power speedmore » range of the DMIC is infinite. The simulation results are supported by closed form expressions for peak and rms motor current and average power derived from analytical solution to the differential equations governing the DMIC/BDCM drive for the lossless case. The analytical solution shows that the range of motor inductance that can be accommodated by the DMIC is more than an order of magnitude such that the DMIC is compatible with both low- and high-inductance BDCMs. Finally, method is given for integrating the classical hysteresis band current control, used for motor control below base speed, with the phase advance of DMIC that is applied above base speed. The power versus speed performance of the DMIC is then simulated across the entire speed range.« less

Semi-analytical solution for the generalized absorbing boundary condition in molecular dynamics simulations

NASA Astrophysics Data System (ADS)

Lee, Chung-Shuo; Chen, Yan-Yu; Yu, Chi-Hua; Hsu, Yu-Chuan; Chen, Chuin-Shan

2017-07-01

We present a semi-analytical solution of a time-history kernel for the generalized absorbing boundary condition in molecular dynamics (MD) simulations. To facilitate the kernel derivation, the concept of virtual atoms in real space that can conform with an arbitrary boundary in an arbitrary lattice is adopted. The generalized Langevin equation is regularized using eigenvalue decomposition and, consequently, an analytical expression of an inverse Laplace transform is obtained. With construction of dynamical matrices in the virtual domain, a semi-analytical form of the time-history kernel functions for an arbitrary boundary in an arbitrary lattice can be found. The time-history kernel functions for different crystal lattices are derived to show the generality of the proposed method. Non-equilibrium MD simulations in a triangular lattice with and without the absorbing boundary condition are conducted to demonstrate the validity of the solution.

On the performance of dual-hop mixed RF/FSO wireless communication system in urban area over aggregated exponentiated Weibull fading channels with pointing errors

NASA Astrophysics Data System (ADS)

Wang, Yue; Wang, Ping; Liu, Xiaoxia; Cao, Tian

2018-03-01

The performance of decode-and-forward dual-hop mixed radio frequency / free-space optical system in urban area is studied. The RF link is modeled by the Nakagami-m distribution and the FSO link is described by the composite exponentiated Weibull (EW) fading channels with nonzero boresight pointing errors (NBPE). For comparison, the ABER results without pointing errors (PE) and those with zero boresight pointing errors (ZBPE) are also provided. The closed-form expression for the average bit error rate (ABER) in RF link is derived with the help of hypergeometric function, and that in FSO link is obtained by Meijer's G and generalized Gauss-Laguerre quadrature functions. Then, the end-to-end ABERs with binary phase shift keying modulation are achieved on the basis of the computed ABER results of RF and FSO links. The end-to-end ABER performance is further analyzed with different Nakagami-m parameters, turbulence strengths, receiver aperture sizes and boresight displacements. The result shows that with ZBPE and NBPE considered, FSO link suffers a severe ABER degradation and becomes the dominant limitation of the mixed RF/FSO system in urban area. However, aperture averaging can bring significant ABER improvement of this system. Monte Carlo simulation is provided to confirm the validity of the analytical ABER expressions.

Closed-form expressions for flip angle variation that maximize total signal in T1-weighted rapid gradient echo MRI.

PubMed

Drobnitzky, Matthias; Klose, Uwe

2017-03-01

Magnetization-prepared rapid gradient-echo (MPRAGE) sequences are commonly employed for T1-weighted structural brain imaging. Following a contrast preparation radiofrequency (RF) pulse, the data acquisition proceeds under nonequilibrium conditions of the relaxing longitudinal magnetization. Variation of the flip angle can be used to maximize total available signal. Simulated annealing or greedy algorithms have so far been published to numerically solve this problem, with signal-to-noise ratios optimized for clinical imaging scenarios by adhering to a predefined shape of the signal evolution. We propose an unconstrained optimization of the MPRAGE experiment that employs techniques from resource allocation theory. A new dynamic programming solution is introduced that yields closed-form expressions for optimal flip angle variation. Flip angle series are proposed that maximize total transverse magnetization (Mxy) for a range of physiologic T1 values. A 3D MPRAGE sequence is modified to allow for a controlled variation of the excitation angle. Experiments employing a T1 contrast phantom are performed at 3T. 1D acquisitions without phase encoding permit measurement of the temporal development of Mxy. Image mean signal and standard deviation for reference flip angle trains are compared in 2D measurements. Signal profiles at sharp phantom edges are acquired to access image blurring related to nonuniform Mxy development. A novel closed-form expression for flip angle variation is found that constitutes the optimal policy to reach maximum total signal. It numerically equals previously published results of other authors when evaluated under their simplifying assumptions. Longitudinal magnetization (Mz) is exhaustively used without causing abrupt changes in the measured MR signal, which is a prerequisite for artifact free images. Phantom experiments at 3T verify the expected benefit for total accumulated k-space signal when compared with published flip angle series. Describing the MR signal collection in MPRAGE sequences as a Bellman problem is a new concept. By means of recursively solving a series of overlapping subproblems, this leads to an elegant solution for the problem of maximizing total available MR signal in k-space. A closed-form expression for flip angle variation avoids the complexity of numerical optimization and eases access to controlled variation in an attempt to identify potential clinical applications. © 2017 American Association of Physicists in Medicine.

Neural Network Based Modeling and Analysis of LP Control Surface Allocation

NASA Technical Reports Server (NTRS)

Langari, Reza; Krishnakumar, Kalmanje; Gundy-Burlet, Karen

2003-01-01

This paper presents an approach to interpretive modeling of LP based control allocation in intelligent flight control. The emphasis is placed on a nonlinear interpretation of the LP allocation process as a static map to support analytical study of the resulting closed loop system, albeit in approximate form. The approach makes use of a bi-layer neural network to capture the essential functioning of the LP allocation process. It is further shown via Lyapunov based analysis that under certain relatively mild conditions the resulting closed loop system is stable. Some preliminary conclusions from a study at Ames are stated and directions for further research are given at the conclusion of the paper.

SEMI-ANALYTIC CALCULATION OF THE TEMPERATURE DISTRIBUTION IN A PERFORATED CIRCLE

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

Kennedy, J.M.; Fowler, J.K.

The flow of heat in a tube-in-shell fuel element is closely related to the two-dimensional heat flow in a circular region perforated by a number of circular holes. Mathematical expressions for the two-dimensional temperature distribution were obtained in terms of sources and sinks of increasing complexity located within the holes and beyond the outer circle. A computer program, TINS, which solves the temperature problem for an array of one or two rings of holes, with or without a center hole, is also described. (auth)

PFLOTRAN Verification: Development of a Testing Suite to Ensure Software Quality

NASA Astrophysics Data System (ADS)

Hammond, G. E.; Frederick, J. M.

2016-12-01

In scientific computing, code verification ensures the reliability and numerical accuracy of a model simulation by comparing the simulation results to experimental data or known analytical solutions. The model is typically defined by a set of partial differential equations with initial and boundary conditions, and verification ensures whether the mathematical model is solved correctly by the software. Code verification is especially important if the software is used to model high-consequence systems which cannot be physically tested in a fully representative environment [Oberkampf and Trucano (2007)]. Justified confidence in a particular computational tool requires clarity in the exercised physics and transparency in its verification process with proper documentation. We present a quality assurance (QA) testing suite developed by Sandia National Laboratories that performs code verification for PFLOTRAN, an open source, massively-parallel subsurface simulator. PFLOTRAN solves systems of generally nonlinear partial differential equations describing multiphase, multicomponent and multiscale reactive flow and transport processes in porous media. PFLOTRAN's QA test suite compares the numerical solutions of benchmark problems in heat and mass transport against known, closed-form, analytical solutions, including documentation of the exercised physical process models implemented in each PFLOTRAN benchmark simulation. The QA test suite development strives to follow the recommendations given by Oberkampf and Trucano (2007), which describes four essential elements in high-quality verification benchmark construction: (1) conceptual description, (2) mathematical description, (3) accuracy assessment, and (4) additional documentation and user information. Several QA tests within the suite will be presented, including details of the benchmark problems and their closed-form analytical solutions, implementation of benchmark problems in PFLOTRAN simulations, and the criteria used to assess PFLOTRAN's performance in the code verification procedure. References Oberkampf, W. L., and T. G. Trucano (2007), Verification and Validation Benchmarks, SAND2007-0853, 67 pgs., Sandia National Laboratories, Albuquerque, NM.

Quantum kinetic theory of the filamentation instability

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

Bret, A.; Haas, F.

2011-07-15

The quantum electromagnetic dielectric tensor for a multi-species plasma is re-derived from the gauge-invariant Wigner-Maxwell system and presented under a form very similar to the classical one. The resulting expression is then applied to a quantum kinetic theory of the electromagnetic filamentation instability. Comparison is made with the quantum fluid theory including a Bohm pressure term and with the cold classical plasma result. A number of analytical expressions are derived for the cutoff wave vector, the largest growth rate, and the most unstable wave vector.

The MGGB equation-of-state for multifield applications: a numerical recipe for analytic expression of sesame EOS data

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

Kashiwa, B. A.

2010-12-01

Abstract A thermodynamically consistent and fully general equation–of– state (EOS) for multifield applications is described. EOS functions are derived from a Helmholtz free energy expressed as the sum of thermal (fluctuational) and collisional (condensed–phase) contributions; thus the free energy is of the Mie–Gr¨uneisen1 form. The phase–coexistence region is defined using a parameterized saturation curve by extending the form introduced by Guggenheim,2 which scales the curve relative to conditions at the critical point. We use the zero–temperature condensed–phase contribution developed by Barnes,3 which extends the Thomas–Fermi–Dirac equation to zero pressure. Thus, the functional form of the EOS could be called MGGBmore » (for Mie– Gr¨uneisen–Guggenheim–Barnes). Substance–specific parameters are obtained by fitting the low–density energy to data from the Sesame4 library; fitting the zero–temperature pressure to the Sesame cold curve; and fitting the saturation curve and latent heat to laboratory data,5 if available. When suitable coexistence data, or Sesame data, are not available, then we apply the Principle of Corresponding States.2 Thus MGGB can be thought of as a numerical recipe for rendering the tabular Sesame EOS data in an analytic form that includes a proper coexistence region, and which permits the accurate calculation of derivatives associated with compressibility, expansivity, Joule coefficient, and specific heat, all of which are required for multifield applications. 1« less

Distributed reacceleration of cosmic rays

NASA Technical Reports Server (NTRS)

Wandel, Amri; Eichler, David; Letaw, John R.; Silberberg, Rein; Tsao, C. H.

1985-01-01

A model is developed in which cosmic rays, in addition to their initial acceleration by a strong shock, are continuously reaccelerated while propagating through the Galaxy. The equations describing this acceleration scheme are solved analytically and numerically. Solutions for the spectra of primary and secondary cosmic rays are given in a closed analytic form, allowing a rapid search in parameter space for viable propagation models with distributed reeacceleration included. The observed boron-to-carbon ratio can be reproduced by the reacceleration theory over a range of escape parameters, some of them quite different from the standard leaky-box model. It is also shown that even a very modest amount of reacceleration by strong shocks causes the boron-to-carbon ratio to level off at sufficiently high energies.

Optimal Mortgage Refinancing: A Closed Form Solution

PubMed Central

Agarwal, Sumit; Driscoll, John C.; Laibson, David I.

2013-01-01

We derive the first closed-form optimal refinancing rule: Refinance when the current mortgage interest rate falls below the original rate by at least 1ψ[ϕ+W(−exp(−ϕ))]. In this formula W(.) is the Lambert W-function, ψ=2(ρ+λ)σ,ϕ=1+ψ(ρ+λ)κ∕M(1−τ), ρ is the real discount rate, λ is the expected real rate of exogenous mortgage repayment, σ is the standard deviation of the mortgage rate, κ/M is the ratio of the tax-adjusted refinancing cost and the remaining mortgage value, and τ is the marginal tax rate. This expression is derived by solving a tractable class of refinancing problems. Our quantitative results closely match those reported by researchers using numerical methods. PMID:25843977

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

Griffin, Brian M.; Larson, Vincent E.

Microphysical processes, such as the formation, growth, and evaporation of precipitation, interact with variability and covariances (e.g., fluxes) in moisture and heat content. For instance, evaporation of rain may produce cold pools, which in turn may trigger fresh convection and precipitation. These effects are usually omitted or else crudely parameterized at subgrid scales in weather and climate models.A more formal approach is pursued here, based on predictive, horizontally averaged equations for the variances, covariances, and fluxes of moisture and heat content. These higher-order moment equations contain microphysical source terms. The microphysics terms can be integrated analytically, given a suitably simplemore » warm-rain microphysics scheme and an approximate assumption about the multivariate distribution of cloud-related and precipitation-related variables. Performing the integrations provides exact expressions within an idealized context.A large-eddy simulation (LES) of a shallow precipitating cumulus case is performed here, and it indicates that the microphysical effects on (co)variances and fluxes can be large. In some budgets and altitude ranges, they are dominant terms. The analytic expressions for the integrals are implemented in a single-column, higher-order closure model. Interactive single-column simulations agree qualitatively with the LES. The analytic integrations form a parameterization of microphysical effects in their own right, and they also serve as benchmark solutions that can be compared to non-analytic integration methods.« less

Experiences with Probabilistic Analysis Applied to Controlled Systems

NASA Technical Reports Server (NTRS)

Kenny, Sean P.; Giesy, Daniel P.

2004-01-01

This paper presents a semi-analytic method for computing frequency dependent means, variances, and failure probabilities for arbitrarily large-order closed-loop dynamical systems possessing a single uncertain parameter or with multiple highly correlated uncertain parameters. The approach will be shown to not suffer from the same computational challenges associated with computing failure probabilities using conventional FORM/SORM techniques. The approach is demonstrated by computing the probabilistic frequency domain performance of an optimal feed-forward disturbance rejection scheme.