Li, Rui; Ye, Hongfei; Zhang, Weisheng; Ma, Guojun; Su, Yewang
2015-01-01
Spring constant calibration of the atomic force microscope (AFM) cantilever is of fundamental importance for quantifying the force between the AFM cantilever tip and the sample. The calibration within the framework of thin plate theory undoubtedly has a higher accuracy and broader scope than that within the well-established beam theory. However, thin plate theory-based accurate analytic determination of the constant has been perceived as an extremely difficult issue. In this paper, we implement the thin plate theory-based analytic modeling for the static behavior of rectangular AFM cantilevers, which reveals that the three-dimensional effect and Poisson effect play important roles in accurate determination of the spring constants. A quantitative scaling law is found that the normalized spring constant depends only on the Poisson's ratio, normalized dimension and normalized load coordinate. Both the literature and our refined finite element model validate the present results. The developed model is expected to serve as the benchmark for accurate calibration of rectangular AFM cantilevers. PMID:26510769
NASA Astrophysics Data System (ADS)
Li, Rui; Ye, Hongfei; Zhang, Weisheng; Ma, Guojun; Su, Yewang
2015-10-01
Spring constant calibration of the atomic force microscope (AFM) cantilever is of fundamental importance for quantifying the force between the AFM cantilever tip and the sample. The calibration within the framework of thin plate theory undoubtedly has a higher accuracy and broader scope than that within the well-established beam theory. However, thin plate theory-based accurate analytic determination of the constant has been perceived as an extremely difficult issue. In this paper, we implement the thin plate theory-based analytic modeling for the static behavior of rectangular AFM cantilevers, which reveals that the three-dimensional effect and Poisson effect play important roles in accurate determination of the spring constants. A quantitative scaling law is found that the normalized spring constant depends only on the Poisson’s ratio, normalized dimension and normalized load coordinate. Both the literature and our refined finite element model validate the present results. The developed model is expected to serve as the benchmark for accurate calibration of rectangular AFM cantilevers.
Li, Rui; Ye, Hongfei; Zhang, Weisheng; Ma, Guojun; Su, Yewang
2015-01-01
Spring constant calibration of the atomic force microscope (AFM) cantilever is of fundamental importance for quantifying the force between the AFM cantilever tip and the sample. The calibration within the framework of thin plate theory undoubtedly has a higher accuracy and broader scope than that within the well-established beam theory. However, thin plate theory-based accurate analytic determination of the constant has been perceived as an extremely difficult issue. In this paper, we implement the thin plate theory-based analytic modeling for the static behavior of rectangular AFM cantilevers, which reveals that the three-dimensional effect and Poisson effect play important roles in accurate determination of the spring constants. A quantitative scaling law is found that the normalized spring constant depends only on the Poisson’s ratio, normalized dimension and normalized load coordinate. Both the literature and our refined finite element model validate the present results. The developed model is expected to serve as the benchmark for accurate calibration of rectangular AFM cantilevers. PMID:26510769
Accurate Analytic Results for the Steady State Distribution of the Eigen Model
NASA Astrophysics Data System (ADS)
Huang, Guan-Rong; Saakian, David B.; Hu, Chin-Kun
2016-04-01
Eigen model of molecular evolution is popular in studying complex biological and biomedical systems. Using the Hamilton-Jacobi equation method, we have calculated analytic equations for the steady state distribution of the Eigen model with a relative accuracy of O(1/N), where N is the length of genome. Our results can be applied for the case of small genome length N, as well as the cases where the direct numerics can not give accurate result, e.g., the tail of distribution.
Development and application of accurate analytical models for single active electron potentials
NASA Astrophysics Data System (ADS)
Miller, Michelle; Jaron-Becker, Agnieszka; Becker, Andreas
2015-05-01
The single active electron (SAE) approximation is a theoretical model frequently employed to study scenarios in which inner-shell electrons may productively be treated as frozen spectators to a physical process of interest, and accurate analytical approximations for these potentials are sought as a useful simulation tool. Density function theory is often used to construct a SAE potential, requiring that a further approximation for the exchange correlation functional be enacted. In this study, we employ the Krieger, Li, and Iafrate (KLI) modification to the optimized-effective-potential (OEP) method to reduce the complexity of the problem to the straightforward solution of a system of linear equations through simple arguments regarding the behavior of the exchange-correlation potential in regions where a single orbital dominates. We employ this method for the solution of atomic and molecular potentials, and use the resultant curve to devise a systematic construction for highly accurate and useful analytical approximations for several systems. Supported by the U.S. Department of Energy (Grant No. DE-FG02-09ER16103), and the U.S. National Science Foundation (Graduate Research Fellowship, Grants No. PHY-1125844 and No. PHY-1068706).
NASA Astrophysics Data System (ADS)
Walter, Johannes; Thajudeen, Thaseem; Süß, Sebastian; Segets, Doris; Peukert, Wolfgang
2015-04-01
Analytical centrifugation (AC) is a powerful technique for the characterisation of nanoparticles in colloidal systems. As a direct and absolute technique it requires no calibration or measurements of standards. Moreover, it offers simple experimental design and handling, high sample throughput as well as moderate investment costs. However, the full potential of AC for nanoparticle size analysis requires the development of powerful data analysis techniques. In this study we show how the application of direct boundary models to AC data opens up new possibilities in particle characterisation. An accurate analysis method, successfully applied to sedimentation data obtained by analytical ultracentrifugation (AUC) in the past, was used for the first time in analysing AC data. Unlike traditional data evaluation routines for AC using a designated number of radial positions or scans, direct boundary models consider the complete sedimentation boundary, which results in significantly better statistics. We demonstrate that meniscus fitting, as well as the correction of radius and time invariant noise significantly improves the signal-to-noise ratio and prevents the occurrence of false positives due to optical artefacts. Moreover, hydrodynamic non-ideality can be assessed by the residuals obtained from the analysis. The sedimentation coefficient distributions obtained by AC are in excellent agreement with the results from AUC. Brownian dynamics simulations were used to generate numerical sedimentation data to study the influence of diffusion on the obtained distributions. Our approach is further validated using polystyrene and silica nanoparticles. In particular, we demonstrate the strength of AC for analysing multimodal distributions by means of gold nanoparticles.
Fast and accurate analytical model to solve inverse problem in SHM using Lamb wave propagation
NASA Astrophysics Data System (ADS)
Poddar, Banibrata; Giurgiutiu, Victor
2016-04-01
Lamb wave propagation is at the center of attention of researchers for structural health monitoring of thin walled structures. This is due to the fact that Lamb wave modes are natural modes of wave propagation in these structures with long travel distances and without much attenuation. This brings the prospect of monitoring large structure with few sensors/actuators. However the problem of damage detection and identification is an "inverse problem" where we do not have the luxury to know the exact mathematical model of the system. On top of that the problem is more challenging due to the confounding factors of statistical variation of the material and geometric properties. Typically this problem may also be ill posed. Due to all these complexities the direct solution of the problem of damage detection and identification in SHM is impossible. Therefore an indirect method using the solution of the "forward problem" is popular for solving the "inverse problem". This requires a fast forward problem solver. Due to the complexities involved with the forward problem of scattering of Lamb waves from damages researchers rely primarily on numerical techniques such as FEM, BEM, etc. But these methods are slow and practically impossible to be used in structural health monitoring. We have developed a fast and accurate analytical forward problem solver for this purpose. This solver, CMEP (complex modes expansion and vector projection), can simulate scattering of Lamb waves from all types of damages in thin walled structures fast and accurately to assist the inverse problem solver.
Walter, Johannes; Thajudeen, Thaseem; Süss, Sebastian; Segets, Doris; Peukert, Wolfgang
2015-04-21
Analytical centrifugation (AC) is a powerful technique for the characterisation of nanoparticles in colloidal systems. As a direct and absolute technique it requires no calibration or measurements of standards. Moreover, it offers simple experimental design and handling, high sample throughput as well as moderate investment costs. However, the full potential of AC for nanoparticle size analysis requires the development of powerful data analysis techniques. In this study we show how the application of direct boundary models to AC data opens up new possibilities in particle characterisation. An accurate analysis method, successfully applied to sedimentation data obtained by analytical ultracentrifugation (AUC) in the past, was used for the first time in analysing AC data. Unlike traditional data evaluation routines for AC using a designated number of radial positions or scans, direct boundary models consider the complete sedimentation boundary, which results in significantly better statistics. We demonstrate that meniscus fitting, as well as the correction of radius and time invariant noise significantly improves the signal-to-noise ratio and prevents the occurrence of false positives due to optical artefacts. Moreover, hydrodynamic non-ideality can be assessed by the residuals obtained from the analysis. The sedimentation coefficient distributions obtained by AC are in excellent agreement with the results from AUC. Brownian dynamics simulations were used to generate numerical sedimentation data to study the influence of diffusion on the obtained distributions. Our approach is further validated using polystyrene and silica nanoparticles. In particular, we demonstrate the strength of AC for analysing multimodal distributions by means of gold nanoparticles. PMID:25789666
Accurate analytical method for the extraction of solar cell model parameters
NASA Astrophysics Data System (ADS)
Phang, J. C. H.; Chan, D. S. H.; Phillips, J. R.
1984-05-01
Single diode solar cell model parameters are rapidly extracted from experimental data by means of the presently derived analytical expressions. The parameter values obtained have a less than 5 percent error for most solar cells, in light of the extraction of model parameters for two cells of differing quality which were compared with parameters extracted by means of the iterative method.
Accurate analytical modelling of cosmic ray induced failure rates of power semiconductor devices
NASA Astrophysics Data System (ADS)
Bauer, Friedhelm D.
2009-06-01
A new, simple and efficient approach is presented to conduct estimations of the cosmic ray induced failure rate for high voltage silicon power devices early in the design phase. This allows combining common design issues such as device losses and safe operating area with the constraints imposed by the reliability to result in a better and overall more efficient design methodology. Starting from an experimental and theoretical background brought forth a few yeas ago [Kabza H et al. Cosmic radiation as a cause for power device failure and possible countermeasures. In: Proceedings of the sixth international symposium on power semiconductor devices and IC's, Davos, Switzerland; 1994. p. 9-12, Zeller HR. Cosmic ray induced breakdown in high voltage semiconductor devices, microscopic model and possible countermeasures. In: Proceedings of the sixth international symposium on power semiconductor devices and IC's, Davos, Switzerland; 1994. p. 339-40, and Matsuda H et al. Analysis of GTO failure mode during d.c. blocking. In: Proceedings of the sixth international symposium on power semiconductor devices and IC's, Davos, Switzerland; 1994. p. 221-5], an exact solution of the failure rate integral is derived and presented in a form which lends itself to be combined with the results available from commercial semiconductor simulation tools. Hence, failure rate integrals can be obtained with relative ease for realistic two- and even three-dimensional semiconductor geometries. Two case studies relating to IGBT cell design and planar junction termination layout demonstrate the purpose of the method.
NASA Astrophysics Data System (ADS)
RoyChoudhury, Raja; RoyChoudhury, Arundhati
2011-02-01
This paper presents a semi analytical formulation of modal properties of a non linear optical fiber having Kerr non linearity with a three parameter approximation of fundamental modal field. The minimization of core parameter ( U) which involves Kerr nonlinearity through the non-stationary expression of propagation constant, is carried out by Nelder-Mead Simplex method of non linear unconstrained minimization, suitable for problems with non-smooth functions as the method does not require any derivative information. Use of three parameters in modal approximation and implementation of Simplex methods enables our semi analytical description to be an alternative way having less computational burden for calculation of modal parameters than full numerical methods.
Robust Accurate Non-Invasive Analyte Monitor
Robinson, Mark R.
1998-11-03
An improved method and apparatus for determining noninvasively and in vivo one or more unknown values of a known characteristic, particularly the concentration of an analyte in human tissue. The method includes: (1) irradiating the tissue with infrared energy (400 nm-2400 nm) having at least several wavelengths in a given range of wavelengths so that there is differential absorption of at least some of the wavelengths by the tissue as a function of the wavelengths and the known characteristic, the differential absorption causeing intensity variations of the wavelengths incident from the tissue; (2) providing a first path through the tissue; (3) optimizing the first path for a first sub-region of the range of wavelengths to maximize the differential absorption by at least some of the wavelengths in the first sub-region; (4) providing a second path through the tissue; and (5) optimizing the second path for a second sub-region of the range, to maximize the differential absorption by at least some of the wavelengths in the second sub-region. In the preferred embodiment a third path through the tissue is provided for, which path is optimized for a third sub-region of the range. With this arrangement, spectral variations which are the result of tissue differences (e.g., melanin and temperature) can be reduced. At least one of the paths represents a partial transmission path through the tissue. This partial transmission path may pass through the nail of a finger once and, preferably, twice. Also included are apparatus for: (1) reducing the arterial pulsations within the tissue; and (2) maximizing the blood content i the tissue.
Anatomically accurate individual face modeling.
Zhang, Yu; Prakash, Edmond C; Sung, Eric
2003-01-01
This paper presents a new 3D face model of a specific person constructed from the anatomical perspective. By exploiting the laser range data, a 3D facial mesh precisely representing the skin geometry is reconstructed. Based on the geometric facial mesh, we develop a deformable multi-layer skin model. It takes into account the nonlinear stress-strain relationship and dynamically simulates the non-homogenous behavior of the real skin. The face model also incorporates a set of anatomically-motivated facial muscle actuators and underlying skull structure. Lagrangian mechanics governs the facial motion dynamics, dictating the dynamic deformation of facial skin in response to the muscle contraction. PMID:15455936
An accurate analytic representation of the water pair potential.
Cencek, Wojciech; Szalewicz, Krzysztof; Leforestier, Claude; van Harrevelt, Rob; van der Avoird, Ad
2008-08-28
The ab initio water dimer interaction energies obtained from coupled cluster calculations and used in the CC-pol water pair potential (Bukowski et al., Science, 2007, 315, 1249) have been refitted to a site-site form containing eight symmetry-independent sites in each monomer and denoted as CC-pol-8s. Initially, the site-site functions were assumed in a B-spline form, which allowed a precise optimization of the positions of the sites. Next, these functions were assumed in the standard exponential plus inverse powers form. The root mean square error of the CC-pol-8s fit with respect to the 2510 ab initio points is 0.10 kcal mol(-1), compared to 0.42 kcal mol(-1) of the CC-pol fit (0.010 kcal mol(-1) compared to 0.089 kcal mol(-1) for points with negative interaction energies). The energies of the stationary points in the CC-pol-8s potential are considerably more accurate than in the case of CC-pol. The water dimer vibration-rotation-tunneling spectrum predicted by the CC-pol-8s potential agrees substantially and systematically better with experiment than the already very accurate spectrum predicted by CC-pol, while specific features that could not be accurately predicted previously now agree very well with experiment. This shows that the uncertainties of the fit were the largest source of error in the previous predictions and that the present potential sets a new standard of accuracy in investigations of the water dimer. PMID:18688514
Accurate mask model for advanced nodes
NASA Astrophysics Data System (ADS)
Zine El Abidine, Nacer; Sundermann, Frank; Yesilada, Emek; Ndiaye, El Hadji Omar; Mishra, Kushlendra; Paninjath, Sankaranarayanan; Bork, Ingo; Buck, Peter; Toublan, Olivier; Schanen, Isabelle
2014-07-01
Standard OPC models consist of a physical optical model and an empirical resist model. The resist model compensates the optical model imprecision on top of modeling resist development. The optical model imprecision may result from mask topography effects and real mask information including mask ebeam writing and mask process contributions. For advanced technology nodes, significant progress has been made to model mask topography to improve optical model accuracy. However, mask information is difficult to decorrelate from standard OPC model. Our goal is to establish an accurate mask model through a dedicated calibration exercise. In this paper, we present a flow to calibrate an accurate mask enabling its implementation. The study covers the different effects that should be embedded in the mask model as well as the experiment required to model them.
Pre-Modeling Ensures Accurate Solid Models
ERIC Educational Resources Information Center
Gow, George
2010-01-01
Successful solid modeling requires a well-organized design tree. The design tree is a list of all the object's features and the sequential order in which they are modeled. The solid-modeling process is faster and less prone to modeling errors when the design tree is a simple and geometrically logical definition of the modeled object. Few high…
New model accurately predicts reformate composition
Ancheyta-Juarez, J.; Aguilar-Rodriguez, E. )
1994-01-31
Although naphtha reforming is a well-known process, the evolution of catalyst formulation, as well as new trends in gasoline specifications, have led to rapid evolution of the process, including: reactor design, regeneration mode, and operating conditions. Mathematical modeling of the reforming process is an increasingly important tool. It is fundamental to the proper design of new reactors and revamp of existing ones. Modeling can be used to optimize operating conditions, analyze the effects of process variables, and enhance unit performance. Instituto Mexicano del Petroleo has developed a model of the catalytic reforming process that accurately predicts reformate composition at the higher-severity conditions at which new reformers are being designed. The new AA model is more accurate than previous proposals because it takes into account the effects of temperature and pressure on the rate constants of each chemical reaction.
Accurate modeling of parallel scientific computations
NASA Technical Reports Server (NTRS)
Nicol, David M.; Townsend, James C.
1988-01-01
Scientific codes are usually parallelized by partitioning a grid among processors. To achieve top performance it is necessary to partition the grid so as to balance workload and minimize communication/synchronization costs. This problem is particularly acute when the grid is irregular, changes over the course of the computation, and is not known until load time. Critical mapping and remapping decisions rest on the ability to accurately predict performance, given a description of a grid and its partition. This paper discusses one approach to this problem, and illustrates its use on a one-dimensional fluids code. The models constructed are shown to be accurate, and are used to find optimal remapping schedules.
The importance of accurate atmospheric modeling
NASA Astrophysics Data System (ADS)
Payne, Dylan; Schroeder, John; Liang, Pang
2014-11-01
This paper will focus on the effect of atmospheric conditions on EO sensor performance using computer models. We have shown the importance of accurately modeling atmospheric effects for predicting the performance of an EO sensor. A simple example will demonstrated how real conditions for several sites in China will significantly impact on image correction, hyperspectral imaging, and remote sensing. The current state-of-the-art model for computing atmospheric transmission and radiance is, MODTRAN® 5, developed by the US Air Force Research Laboratory and Spectral Science, Inc. Research by the US Air Force, Navy and Army resulted in the public release of LOWTRAN 2 in the early 1970's. Subsequent releases of LOWTRAN and MODTRAN® have continued until the present. Please verify that (1) all pages are present, (2) all figures are correct, (3) all fonts and special characters are correct, and (4) all text and figures fit within the red margin lines shown on this review document. Complete formatting information is available at http://SPIE.org/manuscripts Return to the Manage Active Submissions page at http://spie.org/submissions/tasks.aspx and approve or disapprove this submission. Your manuscript will not be published without this approval. Please contact author_help@spie.org with any questions or concerns. The paper will demonstrate the importance of using validated models and local measured meteorological, atmospheric and aerosol conditions to accurately simulate the atmospheric transmission and radiance. Frequently default conditions are used which can produce errors of as much as 75% in these values. This can have significant impact on remote sensing applications.
University Macro Analytic Simulation Model.
ERIC Educational Resources Information Center
Baron, Robert; Gulko, Warren
The University Macro Analytic Simulation System (UMASS) has been designed as a forecasting tool to help university administrators budgeting decisions. Alternative budgeting strategies can be tested on a computer model and then an operational alternative can be selected on the basis of the most desirable projected outcome. UMASS uses readily…
Analytic Time Depending Galaxy Models
NASA Astrophysics Data System (ADS)
Sala, F.
1990-11-01
RESUMEN. Considerando las hip6tesis de Chandrasekhar para el estudjo de la GalActicaq se han desarrollado varios modelos analiticos integrables con simetria axial y dependientes del . . By considering Chandrasekhar hypotheses +or the study o+ Galactic Dynamics, several integrable analytic axisymmetric time-depending galactic models have been developed. Ke ords; GALAXY-DYNAMICS - GALAXY-STRUCTURE
Highly accurate analytical energy of a two-dimensional exciton in a constant magnetic field
NASA Astrophysics Data System (ADS)
Hoang, Ngoc-Tram D.; Nguyen, Duy-Anh P.; Hoang, Van-Hung; Le, Van-Hoang
2016-08-01
Explicit expressions are given for analytically describing the dependence of the energy of a two-dimensional exciton on magnetic field intensity. These expressions are highly accurate with the precision of up to three decimal places for the whole range of the magnetic field intensity. The results are shown for the ground state and some excited states; moreover, we have all formulae to obtain similar expressions of any excited state. Analysis of numerical results shows that the precision of three decimal places is maintained for the excited states with the principal quantum number of up to n=100.
NASA Astrophysics Data System (ADS)
Colalongo, Luigi; Ghittorelli, Matteo; Torricelli, Fabrizio; Kovács-Vajna, Zsolt Miklos
2015-12-01
Surface-potential-based mathematical models are among the most accurate and physically based compact models of Thin-Film Transistors (TFTs) and, in turn, of Organic Thin-Film Transistors (OTFTs), available today. However, the need for iterative computations of the surface potential limits their computational efficiency and diffusion in CAD applications. The existing closed-form approximations of the surface potential are based on regional approximations and empirical smoothing functions that could result not enough accurate to model OTFTs and, in particular, transconductances and transcapacitances. In this paper we present an accurate and computationally efficient closed-form approximation of the surface potential, based on the Lagrange Reversion Theorem, that can be exploited in advanced surface-potential-based OTFTs and TFTs device models.
Analytical model for ramp compression
NASA Astrophysics Data System (ADS)
Xue, Quanxi; Jiang, Shaoen; Wang, Zhebin; Wang, Feng; Hu, Yun; Ding, Yongkun
2016-08-01
An analytical ramp compression model for condensed matter, which can provide explicit solutions for isentropic compression flow fields, is reported. A ramp compression experiment can be easily designed according to the capability of the loading source using this model. Specifically, important parameters, such as the maximum isentropic region width, material properties, profile of the pressure pulse, and the pressure pulse duration can be reasonably allocated or chosen. To demonstrate and study this model, laser-direct-driven ramp compression experiments and code simulation are performed successively, and the factors influencing the accuracy of the model are studied. The application and simulation show that this model can be used as guidance in the design of a ramp compression experiment. However, it is verified that further optimization work is required for a precise experimental design.
Highly accurate analytic formulae for projectile motion subjected to quadratic drag
NASA Astrophysics Data System (ADS)
Turkyilmazoglu, Mustafa
2016-05-01
The classical phenomenon of motion of a projectile fired (thrown) into the horizon through resistive air charging a quadratic drag onto the object is revisited in this paper. No exact solution is known that describes the full physical event under such an exerted resistance force. Finding elegant analytical approximations for the most interesting engineering features of dynamical behavior of the projectile is the principal target. Within this purpose, some analytical explicit expressions are derived that accurately predict the maximum height, its arrival time as well as the flight range of the projectile at the highest ascent. The most significant property of the proposed formulas is that they are not restricted to the initial speed and firing angle of the object, nor to the drag coefficient of the medium. In combination with the available approximations in the literature, it is possible to gain information about the flight and complete the picture of a trajectory with high precision, without having to numerically simulate the full governing equations of motion.
Accurate astronomical atmospheric dispersion models in ZEMAX
NASA Astrophysics Data System (ADS)
Spanò, P.
2014-07-01
ZEMAX provides a standard built-in atmospheric model to simulate atmospheric refraction and dispersion. This model has been compared with other ones to assess its intrinsic accuracy, critical for very demanding application like ADCs for AO-assisted extremely large telescopes. A revised simple model, based on updated published data of the air refractivity, is proposed by using the "Gradient 5" surface of Zemax. At large zenith angles (65 deg), discrepancies up to 100 mas in the differential refraction are expected near the UV atmospheric transmission cutoff. When high-accuracy modeling is required, the latter model should be preferred.
Analytical modeling of the steady radiative shock
NASA Astrophysics Data System (ADS)
Boireau, L.; Bouquet, S.; Michaut, C.; Clique, C.
2006-06-01
In a paper dated 2000 [1], a fully analytical theory of the radiative shock has been presented. This early model had been used to design [2] radiative shock experiments at the Laboratory for the Use of Intense Lasers (LULI) [3 5]. It became obvious from numerical simulations [6, 7] that this model had to be improved in order to accurately recover experiments. In this communication, we present a new theory in which the ionization rates in the unshocked (bar{Z_1}) and shocked (bar{Z_2} neq bar{Z_1}) material, respectively, are included. Associated changes in excitation energy are also taken into account. We study the influence of these effects on the compression and temperature in the shocked medium.
Modeling Biodegradation and Reactive Transport: Analytical and Numerical Models
Sun, Y; Glascoe, L
2005-06-09
The computational modeling of the biodegradation of contaminated groundwater systems accounting for biochemical reactions coupled to contaminant transport is a valuable tool for both the field engineer/planner with limited computational resources and the expert computational researcher less constrained by time and computer power. There exists several analytical and numerical computer models that have been and are being developed to cover the practical needs put forth by users to fulfill this spectrum of computational demands. Generally, analytical models provide rapid and convenient screening tools running on very limited computational power, while numerical models can provide more detailed information with consequent requirements of greater computational time and effort. While these analytical and numerical computer models can provide accurate and adequate information to produce defensible remediation strategies, decisions based on inadequate modeling output or on over-analysis can have costly and risky consequences. In this chapter we consider both analytical and numerical modeling approaches to biodegradation and reactive transport. Both approaches are discussed and analyzed in terms of achieving bioremediation goals, recognizing that there is always a tradeoff between computational cost and the resolution of simulated systems.
Accurate spectral modeling for infrared radiation
NASA Technical Reports Server (NTRS)
Tiwari, S. N.; Gupta, S. K.
1977-01-01
Direct line-by-line integration and quasi-random band model techniques are employed to calculate the spectral transmittance and total band absorptance of 4.7 micron CO, 4.3 micron CO2, 15 micron CO2, and 5.35 micron NO bands. Results are obtained for different pressures, temperatures, and path lengths. These are compared with available theoretical and experimental investigations. For each gas, extensive tabulations of results are presented for comparative purposes. In almost all cases, line-by-line results are found to be in excellent agreement with the experimental values. The range of validity of other models and correlations are discussed.
Accurate theoretical chemistry with coupled pair models.
Neese, Frank; Hansen, Andreas; Wennmohs, Frank; Grimme, Stefan
2009-05-19
Quantum chemistry has found its way into the everyday work of many experimental chemists. Calculations can predict the outcome of chemical reactions, afford insight into reaction mechanisms, and be used to interpret structure and bonding in molecules. Thus, contemporary theory offers tremendous opportunities in experimental chemical research. However, even with present-day computers and algorithms, we cannot solve the many particle Schrodinger equation exactly; inevitably some error is introduced in approximating the solutions of this equation. Thus, the accuracy of quantum chemical calculations is of critical importance. The affordable accuracy depends on molecular size and particularly on the total number of atoms: for orientation, ethanol has 9 atoms, aspirin 21 atoms, morphine 40 atoms, sildenafil 63 atoms, paclitaxel 113 atoms, insulin nearly 800 atoms, and quaternary hemoglobin almost 12,000 atoms. Currently, molecules with up to approximately 10 atoms can be very accurately studied by coupled cluster (CC) theory, approximately 100 atoms with second-order Møller-Plesset perturbation theory (MP2), approximately 1000 atoms with density functional theory (DFT), and beyond that number with semiempirical quantum chemistry and force-field methods. The overwhelming majority of present-day calculations in the 100-atom range use DFT. Although these methods have been very successful in quantum chemistry, they do not offer a well-defined hierarchy of calculations that allows one to systematically converge to the correct answer. Recently a number of rather spectacular failures of DFT methods have been found-even for seemingly simple systems such as hydrocarbons, fueling renewed interest in wave function-based methods that incorporate the relevant physics of electron correlation in a more systematic way. Thus, it would be highly desirable to fill the gap between 10 and 100 atoms with highly correlated ab initio methods. We have found that one of the earliest (and now
Accurate FDTD modelling for dispersive media using rational function and particle swarm optimisation
NASA Astrophysics Data System (ADS)
Chung, Haejun; Ha, Sang-Gyu; Choi, Jaehoon; Jung, Kyung-Young
2015-07-01
This article presents an accurate finite-difference time domain (FDTD) dispersive modelling suitable for complex dispersive media. A quadratic complex rational function (QCRF) is used to characterise their dispersive relations. To obtain accurate coefficients of QCRF, in this work, we use an analytical approach and a particle swarm optimisation (PSO) simultaneously. In specific, an analytical approach is used to obtain the QCRF matrix-solving equation and PSO is applied to adjust a weighting function of this equation. Numerical examples are used to illustrate the validity of the proposed FDTD dispersion model.
Analytic Model of Reactive Flow
Souers, P C; Vitello, P
2004-11-15
A simple analytic model allows prediction of rate constants and size effect behavior before a hydrocode run if size effect data exists. At infinite radius, it defines not only detonation velocity but also average detonation rate, pressure and energy. This allows the derivation of a generalized radius, which becomes larger as the explosive becomes more non-ideal. The model is applied to near-ideal PBX 9404, in-between ANFO and most non-ideal AN. The power of the pressure declines from 2.3, 1.5 to 0.8 across this set. The power of the burn fraction, F, is 0.8, 0 and 0, so that an F-term is important only for the ideal explosives. The size effect shapes change from concave-down to nearly straight to concave-up. Failure is associated with ideal explosives when the calculated detonation velocity turns in a double-valued way. The effect of the power of the pressure may be simulated by including a pressure cutoff in the detonation rate. The models allows comparison of a wide spectrum of explosives providing that a single detonation rate is feasible.
Analytic Model of Reactive Flow
Souers, P C; Vitello, P
2004-08-02
A simple analytic model allows prediction of rate constants and size effect behavior before a hydrocode run if size effect data exists. At infinite radius, it defines not only detonation velocity but also average detonation rate, pressure and energy. This allows the derivation of a generalized radius, which becomes larger as the explosive becomes more non-ideal. The model is applied to near-ideal PBX 9404, in-between ANFO and most non-ideal AN. The power of the pressure declines from 2.3, 1.5 to 0.8 across this set. The power of the burn fraction, F, is 0.8, 0 and 0, so that an F-term is important only for the ideal explosives. The size effect shapes change from concave-down to nearly straight to concave-up. Failure is associated with ideal explosives when the calculated detonation velocity turns in a double-valued way. The effect of the power of the pressure may be simulated by including a pressure cutoff in the detonation rate. The models allows comparison of a wide spectrum of explosives providing that a single detonation rate is feasible.
Analytic Model of Antenna Sheaths
NASA Astrophysics Data System (ADS)
D'Ippolito, D. A.; Myra, J. R.
2008-11-01
RF sheaths are generated on ICRF antennas whenever the launched fast wave also drives a slow wave, e.g. when the magnetic field is tilted (not perpendicular to the current straps). A new approach to sheath modeling was recently proposed in which the RF waves are computed using a modified boundary condition at the sheath surface to describe the plasma-sheath coupling. Here, we illustrate the use of the sheath BC for antenna sheaths by a model electromagnetic perturbation calculation, treating the B field tilt as a small parameter. Analytic expressions are obtained for the sheath voltage and the rf electric field parallel to B in both sheath and plasma regions, including the Child-Langmuir (self-consistency) constraint. It is shown that the plasma corrections to the sheath voltage (which screen the rf field) can be important. The simple vacuum-field sheath-voltage estimate is obtained as a limiting case. Implications for antenna codes such as TOPICA will be discussed. D.A. D'Ippolito and J.R. Myra, Phys. Plasmas 13, 102508 (2006). V. Lancellotti et al., Nucl. Fusion 46, S476 (2006).
Analytic Modeling of Pressurization and Cryogenic Propellant
NASA Technical Reports Server (NTRS)
Corpening, Jeremy H.
2010-01-01
An analytic model for pressurization and cryogenic propellant conditions during all mission phases of any liquid rocket based vehicle has been developed and validated. The model assumes the propellant tanks to be divided into five nodes and also implements an empirical correlation for liquid stratification if desired. The five nodes include a tank wall node exposed to ullage gas, an ullage gas node, a saturated propellant vapor node at the liquid-vapor interface, a liquid node, and a tank wall node exposed to liquid. The conservation equations of mass and energy are then applied across all the node boundaries and, with the use of perfect gas assumptions, explicit solutions for ullage and liquid conditions are derived. All fluid properties are updated real time using NIST Refprop.1 Further, mass transfer at the liquid-vapor interface is included in the form of evaporation, bulk boiling of liquid propellant, and condensation given the appropriate conditions for each. Model validation has proven highly successful against previous analytic models and various Saturn era test data and reasonably successful against more recent LH2 tank self pressurization ground test data. Finally, this model has been applied to numerous design iterations for the Altair Lunar Lander, Ares V Core Stage, and Ares V Earth Departure Stage in order to characterize Helium and autogenous pressurant requirements, propellant lost to evaporation and thermodynamic venting to maintain propellant conditions, and non-uniform tank draining in configurations utilizing multiple LH2 or LO2 propellant tanks. In conclusion, this model provides an accurate and efficient means of analyzing multiple design configurations for any cryogenic propellant tank in launch, low-acceleration coast, or in-space maneuvering and supplies the user with pressurization requirements, unusable propellants from evaporation and liquid stratification, and general ullage gas, liquid, and tank wall conditions as functions of time.
NASA Technical Reports Server (NTRS)
Schlosser, Herbert; Ferrante, John
1989-01-01
An accurate analytic expression for the nonlinear change of the volume of a solid as a function of applied pressure is of great interest in high-pressure experimentation. It is found that a two-parameter analytic expression, fits the experimental volume-change data to within a few percent over the entire experimentally attainable pressure range. Results are presented for 24 different materials including metals, ceramic semiconductors, polymers, and ionic and rare-gas solids.
Simple analytic model for astrophysical S factors
Yakovlev, D. G.; Beard, M.; Gasques, L. R.; Wiescher, M.
2010-10-15
We propose a physically transparent analytic model of astrophysical S factors as a function of a center-of-mass energy E of colliding nuclei (below and above the Coulomb barrier) for nonresonant fusion reactions. For any given reaction, the S(E) model contains four parameters [two of which approximate the barrier potential, U(r)]. They are easily interpolated along many reactions involving isotopes of the same elements; they give accurate practical expressions for S(E) with only several input parameters for many reactions. The model reproduces the suppression of S(E) at low energies (of astrophysical importance) due to the shape of the low-r wing of U(r). The model can be used to reconstruct U(r) from computed or measured S(E). For illustration, we parametrize our recent calculations of S(E) (using the Sao Paulo potential and the barrier penetration formalism) for 946 reactions involving stable and unstable isotopes of C, O, Ne, and Mg (with nine parameters for all reactions involving many isotopes of the same elements, e.g., C+O). In addition, we analyze astrophysically important {sup 12}C+{sup 12}C reaction, compare theoretical models with experimental data, and discuss the problem of interpolating reliably known S(E) values to low energies (E < or approx. 2-3 MeV).
An analytic performance model of disk arrays and its application
NASA Technical Reports Server (NTRS)
Lee, Edward K.; Katz, Randy H.
1991-01-01
As disk arrays become widely used, tools for understanding and analyzing their performance become increasingly important. In particular, performance models can be invaluable in both configuring and designing disk arrays. Accurate analytic performance models are desirable over other types of models because they can be quickly evaluated, are applicable under a wide range of system and workload parameters, and can be manipulated by a range of mathematical techniques. Unfortunately, analytical performance models of disk arrays are difficult to formulate due to the presence of queuing and fork-join synchronization; a disk array request is broken up into independent disk requests which must all complete to satisfy the original request. We develop, validate, and apply an analytic performance model for disk arrays. We derive simple equations for approximating their utilization, response time, and throughput. We then validate the analytic model via simulation and investigate the accuracy of each approximation used in deriving the analytical model. Finally, we apply the analytical model to derive an equation for the optimal unit of data striping in disk arrays.
ANALYTICAL ELEMENT MODELING OF COASTAL AQUIFERS
Four topics were studied concerning the modeling of groundwater flow in coastal aquifers with analytic elements: (1) practical experience was obtained by constructing a groundwater model of the shallow aquifers below the Delmarva Peninsula USA using the commercial program MVAEM; ...
Analytic barrage attack model. Final report, January 1986-January 1989
St Ledger, J.W.; Naegeli, R.E.; Dowden, N.A.
1989-01-01
An analytic model is developed for a nuclear barrage attack, assuming weapons with no aiming error and a cookie-cutter damage function. The model is then extended with approximations for the effects of aiming error and distance damage sigma. The final result is a fast running model which calculates probability of damage for a barrage attack. The probability of damage is accurate to within seven percent or better, for weapon reliabilities of 50 to 100 percent, distance damage sigmas of 0.5 or less, and zero to very large circular error probabilities. FORTRAN 77 coding is included in the report for the analytic model and for a numerical model used to check the analytic results.
Evaluation of one dimensional analytical models for vegetation canopies
NASA Technical Reports Server (NTRS)
Goel, Narendra S.; Kuusk, Andres
1992-01-01
The SAIL model for one-dimensional homogeneous vegetation canopies has been modified to include the specular reflectance and hot spot effects. This modified model and the Nilson-Kuusk model are evaluated by comparing the reflectances given by them against those given by a radiosity-based computer model, Diana, for a set of canopies, characterized by different leaf area index (LAI) and leaf angle distribution (LAD). It is shown that for homogeneous canopies, the analytical models are generally quite accurate in the visible region, but not in the infrared region. For architecturally realistic heterogeneous canopies of the type found in nature, these models fall short. These shortcomings are quantified.
Leng, Wei; Ju, Lili; Gunzburger, Max; Price, Stephen; Ringler, Todd
2012-01-01
The numerical modeling of glacier and ice sheet evolution is a subject of growing interest, in part because of the potential for models to inform estimates of global sea level change. This paper focuses on the development of a numerical model that determines the velocity and pressure fields within an ice sheet. Our numerical model features a high-fidelity mathematical model involving the nonlinear Stokes system and combinations of no-sliding and sliding basal boundary conditions, high-order accurate finite element discretizations based on variable resolution grids, and highly scalable parallel solution strategies, all of which contribute to a numerical model that can achieve accurate velocity and pressure approximations in a highly efficient manner. We demonstrate the accuracy and efficiency of our model by analytical solution tests, established ice sheet benchmark experiments, and comparisons with other well-established ice sheet models.
Piezoresistive Cantilever Performance-Part I: Analytical Model for Sensitivity.
Park, Sung-Jin; Doll, Joseph C; Pruitt, Beth L
2010-02-01
An accurate analytical model for the change in resistance of a piezoresistor is necessary for the design of silicon piezoresistive transducers. Ion implantation requires a high-temperature oxidation or annealing process to activate the dopant atoms, and this treatment results in a distorted dopant profile due to diffusion. Existing analytical models do not account for the concentration dependence of piezoresistance and are not accurate for nonuniform dopant profiles. We extend previous analytical work by introducing two nondimensional factors, namely, the efficiency and geometry factors. A practical benefit of this efficiency factor is that it separates the process parameters from the design parameters; thus, designers may address requirements for cantilever geometry and fabrication process independently. To facilitate the design process, we provide a lookup table for the efficiency factor over an extensive range of process conditions. The model was validated by comparing simulation results with the experimentally determined sensitivities of piezoresistive cantilevers. We performed 9200 TSUPREM4 simulations and fabricated 50 devices from six unique process flows; we systematically explored the design space relating process parameters and cantilever sensitivity. Our treatment focuses on piezoresistive cantilevers, but the analytical sensitivity model is extensible to other piezoresistive transducers such as membrane pressure sensors. PMID:20336183
Piezoresistive Cantilever Performance—Part I: Analytical Model for Sensitivity
Park, Sung-Jin; Doll, Joseph C.; Pruitt, Beth L.
2010-01-01
An accurate analytical model for the change in resistance of a piezoresistor is necessary for the design of silicon piezoresistive transducers. Ion implantation requires a high-temperature oxidation or annealing process to activate the dopant atoms, and this treatment results in a distorted dopant profile due to diffusion. Existing analytical models do not account for the concentration dependence of piezoresistance and are not accurate for nonuniform dopant profiles. We extend previous analytical work by introducing two nondimensional factors, namely, the efficiency and geometry factors. A practical benefit of this efficiency factor is that it separates the process parameters from the design parameters; thus, designers may address requirements for cantilever geometry and fabrication process independently. To facilitate the design process, we provide a lookup table for the efficiency factor over an extensive range of process conditions. The model was validated by comparing simulation results with the experimentally determined sensitivities of piezoresistive cantilevers. We performed 9200 TSUPREM4 simulations and fabricated 50 devices from six unique process flows; we systematically explored the design space relating process parameters and cantilever sensitivity. Our treatment focuses on piezoresistive cantilevers, but the analytical sensitivity model is extensible to other piezoresistive transducers such as membrane pressure sensors. PMID:20336183
NASA Astrophysics Data System (ADS)
Le Roy, Robert J.; Meshkov, Vladimir V.; Stolyarov, Andrej V.
2009-06-01
We have shown that one and two-parameter analytical mapping functions such as r(y;bar{r}, α)=bar{r}[1 + {1}/{α} tan(π y/2)] and r(y;bar{r})=bar{r} [ {1+ y}/{1-y} ] transform the conventional radial Schrödinger equation into equivalent alternate forms {d^2φ(y)}/{dy^2} = [{π^2}/{4}+({2μ}/ {hbar^2} ) g^2(y) [E - U(r(y))
Fixing a rigorous formalism for the accurate analytic derivation of halo properties
NASA Astrophysics Data System (ADS)
Juan, Enric; Salvador-Solé, Eduard; Domènech, Guillem; Manrique, Alberto
2014-03-01
We establish a one-to-one correspondence between virialized haloes and their seeds, namely peaks with a given density contrast at appropriate Gaussian-filtering radii, in the initial Gaussian random density field. This fixes a rigorous formalism for the analytic derivation of halo properties from the linear power spectrum of density perturbations in any hierarchical cosmology. The typical spherically averaged density profile and mass function of haloes so obtained match those found in numerical simulations.
MPD Thruster Performance Analytic Models
NASA Astrophysics Data System (ADS)
Gilland, James; Johnston, Geoffrey
2003-01-01
Magnetoplasmadynamic (MPD) thrusters are capable of accelerating quasi-neutral plasmas to high exhaust velocities using Megawatts (MW) of electric power. These characteristics make such devices worthy of consideration for demanding, far-term missions such as the human exploration of Mars or beyond. Assessment of MPD thrusters at the system and mission level is often difficult due to their status as ongoing experimental research topics rather than developed thrusters. However, in order to assess MPD thrusters' utility in later missions, some adequate characterization of performance, or more exactly, projected performance, and system level definition are required for use in analyses. The most recent physical models of self-field MPD thrusters have been examined, assessed, and reconfigured for use by systems and mission analysts. The physical models allow for rational projections of thruster performance based on physical parameters that can be measured in the laboratory. The models and their implications for the design of future MPD thrusters are presented.
MPD Thruster Performance Analytic Models
NASA Technical Reports Server (NTRS)
Gilland, James; Johnston, Geoffrey
2007-01-01
Magnetoplasmadynamic (MPD) thrusters are capable of accelerating quasi-neutral plasmas to high exhaust velocities using Megawatts (MW) of electric power. These characteristics make such devices worthy of consideration for demanding, far-term missions such as the human exploration of Mars or beyond. Assessment of MPD thrusters at the system and mission level is often difficult due to their status as ongoing experimental research topics rather than developed thrusters. However, in order to assess MPD thrusters utility in later missions, some adequate characterization of performance, or more exactly, projected performance, and system level definition are required for use in analyses. The most recent physical models of self-field MPD thrusters have been examined, assessed, and reconfigured for use by systems and mission analysts. The physical models allow for rational projections of thruster performance based on physical parameters that can be measured in the laboratory. The models and their implications for the design of future MPD thrusters are presented.
MPD Thruster Performance Analytic Models
NASA Technical Reports Server (NTRS)
Gilland, James; Johnston, Geoffrey
2003-01-01
Magnetoplasmadynamic (MPD) thrusters are capable of accelerating quasi-neutral plasmas to high exhaust velocities using Megawatts (MW) of electric power. These characteristics make such devices worthy of consideration for demanding, far-term missions such as the human exploration of Mars or beyond. Assessment of MPD thrusters at the system and mission level is often difficult due to their status as ongoing experimental research topics rather than developed thrusters. However, in order to assess MPD thrusters utility in later missions, some adequate characterization of performance, or more exactly, projected performance, and system level definition are required for use in analyses. The most recent physical models of self-field MPD thrusters have been examined, assessed, and reconfigured for use by systems and mission analysts. The physical models allow for rational projections of thruster performance based on physical parameters that can be measured in the laboratory. The models and their implications for the design of future MPD thrusters are presented.
Analytic gain in probabilistic decompression sickness models.
Howle, Laurens E
2013-11-01
Decompression sickness (DCS) is a disease known to be related to inert gas bubble formation originating from gases dissolved in body tissues. Probabilistic DCS models, which employ survival and hazard functions, are optimized by fitting model parameters to experimental dive data. In the work reported here, I develop methods to find the survival function gain parameter analytically, thus removing it from the fitting process. I show that the number of iterations required for model optimization is significantly reduced. The analytic gain method substantially improves the condition number of the Hessian matrix which reduces the model confidence intervals by more than an order of magnitude. PMID:24209920
Knowledge Generation Model for Visual Analytics.
Sacha, Dominik; Stoffel, Andreas; Stoffel, Florian; Kwon, Bum Chul; Ellis, Geoffrey; Keim, Daniel A
2014-12-01
Visual analytics enables us to analyze huge information spaces in order to support complex decision making and data exploration. Humans play a central role in generating knowledge from the snippets of evidence emerging from visual data analysis. Although prior research provides frameworks that generalize this process, their scope is often narrowly focused so they do not encompass different perspectives at different levels. This paper proposes a knowledge generation model for visual analytics that ties together these diverse frameworks, yet retains previously developed models (e.g., KDD process) to describe individual segments of the overall visual analytic processes. To test its utility, a real world visual analytics system is compared against the model, demonstrating that the knowledge generation process model provides a useful guideline when developing and evaluating such systems. The model is used to effectively compare different data analysis systems. Furthermore, the model provides a common language and description of visual analytic processes, which can be used for communication between researchers. At the end, our model reflects areas of research that future researchers can embark on. PMID:26356874
Water wave model with accurate dispersion and vertical vorticity
NASA Astrophysics Data System (ADS)
Bokhove, Onno
2010-05-01
Cotter and Bokhove (Journal of Engineering Mathematics 2010) derived a variational water wave model with accurate dispersion and vertical vorticity. In one limit, it leads to Luke's variational principle for potential flow water waves. In the another limit it leads to the depth-averaged shallow water equations including vertical vorticity. Presently, focus will be put on the Hamiltonian formulation of the variational model and its boundary conditions.
Analytic modeling of aerosol size distributions
NASA Technical Reports Server (NTRS)
Deepack, A.; Box, G. P.
1979-01-01
Mathematical functions commonly used for representing aerosol size distributions are studied parametrically. Methods for obtaining best fit estimates of the parameters are described. A catalog of graphical plots depicting the parametric behavior of the functions is presented along with procedures for obtaining analytical representations of size distribution data by visual matching of the data with one of the plots. Examples of fitting the same data with equal accuracy by more than one analytic model are also given.
An Accurate and Dynamic Computer Graphics Muscle Model
NASA Technical Reports Server (NTRS)
Levine, David Asher
1997-01-01
A computer based musculo-skeletal model was developed at the University in the departments of Mechanical and Biomedical Engineering. This model accurately represents human shoulder kinematics. The result of this model is the graphical display of bones moving through an appropriate range of motion based on inputs of EMGs and external forces. The need existed to incorporate a geometric muscle model in the larger musculo-skeletal model. Previous muscle models did not accurately represent muscle geometries, nor did they account for the kinematics of tendons. This thesis covers the creation of a new muscle model for use in the above musculo-skeletal model. This muscle model was based on anatomical data from the Visible Human Project (VHP) cadaver study. Two-dimensional digital images from the VHP were analyzed and reconstructed to recreate the three-dimensional muscle geometries. The recreated geometries were smoothed, reduced, and sliced to form data files defining the surfaces of each muscle. The muscle modeling function opened these files during run-time and recreated the muscle surface. The modeling function applied constant volume limitations to the muscle and constant geometry limitations to the tendons.
Analytical Grid Generation for accurate representation of clearances in CFD for Screw Machines
NASA Astrophysics Data System (ADS)
Rane, S.; Kovačević, A.; Stošić, N.
2015-08-01
One of the major factors affecting the performance prediction of twin screw compressors by use of computational fluid dynamics (CFD) is the accuracy with which the leakage gaps are captured by the discretization methods. The accuracy of mapping leakage flows can be improved by increasing the number of grid points on the profile. However, this method faces limitations when it comes to the complex deforming domains of a twin screw compressor because the computational time increases tremendously. In order to address this problem, an analytical grid distribution procedure is formulated that can independently refine the region of high importance for leakage flows in the interlobe space. This paper describes the procedure of analytical grid generation with the refined mesh in the interlobe area and presents a test case to show the influence of the mesh refinement in that area on the performance prediction. It is shown that by using this method, the flow domains in the vicinity of the interlobe gap and the blowhole area are refined which improves accuracy of leakage flow predictions.
Analytical models of slug tests
Karasaki, K.; Long, J.C.S.; Witherspoon, P.A.
1988-01-01
In the present paper, attempts are made to develop solutions to various models of slug tests that may be applicable in analyzing the results of such tests where existing solutions are inadequate. Various geometries that may be encountered in heterogeneous systems such as fractured rocks are considered. Solutions are presented for linear flow, radial flow with boundaries, two layer, and concentric composite models with different flow geometries between the inner and outer region. Solutions are obtained in Laplace space and inverted back to real space numerically. Type curves are presented for each solution. Analyses of the type curves and derivative response curves reveal that many curves have unique shapes only for certain combination of the flow parameters and the distance. Other sets of type curves are similar in shape, although log-log plots and derivative plots may emphasize some features that may not be apparent in semilog plots. These results show that slug tests suffer problems of nonuniqueness to a greater extent than other well tests.
Analytic models of warm plasma dispersion relations
Seough, J. J.; Yoon, P. H.
2009-09-15
The present paper is concerned with analytic models of warm plasma dispersion relations for electromagnetic waves propagating parallel to the ambient magnetic field. Specifically, effects of finite betas on two slow modes, namely, the left-hand circularly polarized ion-cyclotron mode and the right-hand circularly polarized whistler mode, are investigated. Analytic models of the warm plasma dispersion relations are constructed on the basis of conjecture and upon comparisons with numerically found roots. It is shown that the model solutions are good substitutes for actual roots. The significance of the present work in the context of nonlinear plasma research is discussed.
An analytic model for the Phobos surface
NASA Technical Reports Server (NTRS)
Duxbury, Thomas C.
1991-01-01
Analytic expressions are derived to model the surface topography and the normal to the surface of Phobos. The analytic expressions are comprised of a spherical harmonic expansion for the global figure of Phobos, augmented by addition terms for the large crater Stickney and other craters. Over 300 craters were measured in more than 100 Viking Orbiter images to produce the model. In general, the largest craters were measured since they have a significant effect on topography. The topographic model derived has a global spatial and topographic accuracy ranging from about 100 m in areas having the highest resolution and convergent, stereo coverage, up to 500 m in the poorest areas.
WELLHEAD ANALYTIC ELEMENT MODEL FOR WINDOWS
WhAEM2000 (wellhead analytic element model for Win 98/00/NT/XP) is a public domain, ground-water flow model designed to facilitate capture zone delineation and protection area mapping in support of the State's and Tribe's Wellhead Protection Programs (WHPP) and Source Water Asses...
ESTIMATING UNCERTAINITIES IN FACTOR ANALYTIC MODELS
When interpreting results from factor analytic models as used in receptor modeling, it is important to quantify the uncertainties in those results. For example, if the presence of a species on one of the factors is necessary to interpret the factor as originating from a certain ...
ANALYTIC ELEMENT MODELING OF COASTAL AQUIFERS
Four topics were studied concerning the modeling of ground-water flow in coastal aquifers with analytic elements: (1) practical experience was obtained by constructing a ground-water model of the shallow aquifers below the Delmarva Peninsula USA using the commercial program MVAEM...
Palm: Easing the Burden of Analytical Performance Modeling
Tallent, Nathan R.; Hoisie, Adolfy
2014-06-01
Analytical (predictive) application performance models are critical for diagnosing performance-limiting resources, optimizing systems, and designing machines. Creating models, however, is difficult because they must be both accurate and concise. To ease the burden of performance modeling, we developed Palm, a modeling tool that combines top-down (human-provided) semantic insight with bottom-up static and dynamic analysis. To express insight, Palm defines a source code modeling annotation language. By coordinating models and source code, Palm's models are `first-class' and reproducible. Unlike prior work, Palm formally links models, functions, and measurements. As a result, Palm (a) uses functions to either abstract or express complexity (b) generates hierarchical models (representing an application's static and dynamic structure); and (c) automatically incorporates measurements to focus attention, represent constant behavior, and validate models. We discuss generating models for three different applications.
Analytical models for electrically thin flat lenses and reflectors.
Ruphuy, Miguel; Ramahi, Omar M
2015-04-01
This work presents analytical models for two-dimensional (2D) and three-dimensional (3D) electrically thin lenses and reflectors. The 2D formulation is based on infinite current line sources, whereas the 3D formulation is based on electrically small dipoles. These models emulate the energy convergence of an electrically thin flat lens and reflector when illuminated by a plane wave with specific polarization. The advantages of these models are twofold: first, prediction of the performance of electrically thin flat lenses and reflectors can be made significantly faster than full-wave simulators, and second, providing insight on the performance of these electrically thin devices. The analytic models were validated by comparison with full-wave simulation for several interesting examples. The validation results show that the focal point of the electrically thin flat lenses and reflectors can be accurately predicted through a design that assumes low coupling between different layers of an inhomogeneous media. PMID:26366759
Local Debonding and Fiber Breakage in Composite Materials Modeled Accurately
NASA Technical Reports Server (NTRS)
Bednarcyk, Brett A.; Arnold, Steven M.
2001-01-01
A prerequisite for full utilization of composite materials in aerospace components is accurate design and life prediction tools that enable the assessment of component performance and reliability. Such tools assist both structural analysts, who design and optimize structures composed of composite materials, and materials scientists who design and optimize the composite materials themselves. NASA Glenn Research Center's Micromechanics Analysis Code with Generalized Method of Cells (MAC/GMC) software package (http://www.grc.nasa.gov/WWW/LPB/mac) addresses this need for composite design and life prediction tools by providing a widely applicable and accurate approach to modeling composite materials. Furthermore, MAC/GMC serves as a platform for incorporating new local models and capabilities that are under development at NASA, thus enabling these new capabilities to progress rapidly to a stage in which they can be employed by the code's end users.
A novel fast and accurate pseudo-analytical simulation approach for MOAO
NASA Astrophysics Data System (ADS)
Gendron, É.; Charara, A.; Abdelfattah, A.; Gratadour, D.; Keyes, D.; Ltaief, H.; Morel, C.; Vidal, F.; Sevin, A.; Rousset, G.
2014-08-01
Multi-object adaptive optics (MOAO) is a novel adaptive optics (AO) technique for wide-field multi-object spectrographs (MOS). MOAO aims at applying dedicated wavefront corrections to numerous separated tiny patches spread over a large field of view (FOV), limited only by that of the telescope. The control of each deformable mirror (DM) is done individually using a tomographic reconstruction of the phase based on measurements from a number of wavefront sensors (WFS) pointing at natural and artificial guide stars in the field. We have developed a novel hybrid, pseudo-analytical simulation scheme, somewhere in between the end-to- end and purely analytical approaches, that allows us to simulate in detail the tomographic problem as well as noise and aliasing with a high fidelity, and including fitting and bandwidth errors thanks to a Fourier-based code. Our tomographic approach is based on the computation of the minimum mean square error (MMSE) reconstructor, from which we derive numerically the covariance matrix of the tomographic error, including aliasing and propagated noise. We are then able to simulate the point-spread function (PSF) associated to this covariance matrix of the residuals, like in PSF reconstruction algorithms. The advantage of our approach is that we compute the same tomographic reconstructor that would be computed when operating the real instrument, so that our developments open the way for a future on-sky implementation of the tomographic control, plus the joint PSF and performance estimation. The main challenge resides in the computation of the tomographic reconstructor which involves the inversion of a large matrix (typically 40 000 × 40 000 elements). To perform this computation efficiently, we chose an optimized approach based on the use of GPUs as accelerators and using an optimized linear algebra library: MORSE providing a significant speedup against standard CPU oriented libraries such as Intel MKL. Because the covariance matrix is
An Analytical Model for the Influence of Contact Resistance on Thermoelectric Efficiency
NASA Astrophysics Data System (ADS)
Bjørk, Rasmus
2016-03-01
An analytical model is presented that can account for both electrical and hot and cold thermal contact resistances when calculating the efficiency of a thermoelectric generator. The model is compared to a numerical model of a thermoelectric leg for 16 different thermoelectric materials, as well as to the analytical models of Ebling et al. (J Electron Mater 39:1376, 2010) and Min and Rowe (J Power Sour 38:253, 1992). The model presented here is shown to accurately calculate the efficiency for all systems and all contact resistances considered, with an average difference in efficiency between the numerical model and the analytical model of -0.07 ± 0.35pp. This makes the model more accurate than previously published models. The maximum absolute difference in efficiency between the analytical model and the numerical model is 1.14pp for all materials and all contact resistances considered.
An Accurate Temperature Correction Model for Thermocouple Hygrometers 1
Savage, Michael J.; Cass, Alfred; de Jager, James M.
1982-01-01
Numerous water relation studies have used thermocouple hygrometers routinely. However, the accurate temperature correction of hygrometer calibration curve slopes seems to have been largely neglected in both psychrometric and dewpoint techniques. In the case of thermocouple psychrometers, two temperature correction models are proposed, each based on measurement of the thermojunction radius and calculation of the theoretical voltage sensitivity to changes in water potential. The first model relies on calibration at a single temperature and the second at two temperatures. Both these models were more accurate than the temperature correction models currently in use for four psychrometers calibrated over a range of temperatures (15-38°C). The model based on calibration at two temperatures is superior to that based on only one calibration. The model proposed for dewpoint hygrometers is similar to that for psychrometers. It is based on the theoretical voltage sensitivity to changes in water potential. Comparison with empirical data from three dewpoint hygrometers calibrated at four different temperatures indicates that these instruments need only be calibrated at, e.g. 25°C, if the calibration slopes are corrected for temperature. PMID:16662241
An accurate temperature correction model for thermocouple hygrometers.
Savage, M J; Cass, A; de Jager, J M
1982-02-01
Numerous water relation studies have used thermocouple hygrometers routinely. However, the accurate temperature correction of hygrometer calibration curve slopes seems to have been largely neglected in both psychrometric and dewpoint techniques.In the case of thermocouple psychrometers, two temperature correction models are proposed, each based on measurement of the thermojunction radius and calculation of the theoretical voltage sensitivity to changes in water potential. The first model relies on calibration at a single temperature and the second at two temperatures. Both these models were more accurate than the temperature correction models currently in use for four psychrometers calibrated over a range of temperatures (15-38 degrees C). The model based on calibration at two temperatures is superior to that based on only one calibration.The model proposed for dewpoint hygrometers is similar to that for psychrometers. It is based on the theoretical voltage sensitivity to changes in water potential. Comparison with empirical data from three dewpoint hygrometers calibrated at four different temperatures indicates that these instruments need only be calibrated at, e.g. 25 degrees C, if the calibration slopes are corrected for temperature. PMID:16662241
More-Accurate Model of Flows in Rocket Injectors
NASA Technical Reports Server (NTRS)
Hosangadi, Ashvin; Chenoweth, James; Brinckman, Kevin; Dash, Sanford
2011-01-01
An improved computational model for simulating flows in liquid-propellant injectors in rocket engines has been developed. Models like this one are needed for predicting fluxes of heat in, and performances of, the engines. An important part of predicting performance is predicting fluctuations of temperature, fluctuations of concentrations of chemical species, and effects of turbulence on diffusion of heat and chemical species. Customarily, diffusion effects are represented by parameters known in the art as the Prandtl and Schmidt numbers. Prior formulations include ad hoc assumptions of constant values of these parameters, but these assumptions and, hence, the formulations, are inaccurate for complex flows. In the improved model, these parameters are neither constant nor specified in advance: instead, they are variables obtained as part of the solution. Consequently, this model represents the effects of turbulence on diffusion of heat and chemical species more accurately than prior formulations do, and may enable more-accurate prediction of mixing and flows of heat in rocket-engine combustion chambers. The model has been implemented within CRUNCH CFD, a proprietary computational fluid dynamics (CFD) computer program, and has been tested within that program. The model could also be implemented within other CFD programs.
On the importance of having accurate data for astrophysical modelling
NASA Astrophysics Data System (ADS)
Lique, Francois
2016-06-01
The Herschel telescope and the ALMA and NOEMA interferometers have opened new windows of observation for wavelengths ranging from far infrared to sub-millimeter with spatial and spectral resolutions previously unmatched. To make the most of these observations, an accurate knowledge of the physical and chemical processes occurring in the interstellar and circumstellar media is essential.In this presentation, I will discuss what are the current needs of astrophysics in terms of molecular data and I will show that accurate molecular data are crucial for the proper determination of the physical conditions in molecular clouds.First, I will focus on collisional excitation studies that are needed for molecular lines modelling beyond the Local Thermodynamic Equilibrium (LTE) approach. In particular, I will show how new collisional data for the HCN and HNC isomers, two tracers of star forming conditions, have allowed solving the problem of their respective abundance in cold molecular clouds. I will also present the last collisional data that have been computed in order to analyse new highly resolved observations provided by the ALMA interferometer.Then, I will present the calculation of accurate rate constants for the F+H2 → HF+H and Cl+H2 ↔ HCl+H reactions, which have allowed a more accurate determination of the physical conditions in diffuse molecular clouds. I will also present the recent work on the ortho-para-H2 conversion due to hydrogen exchange that allow more accurate determination of the ortho-to-para-H2 ratio in the universe and that imply a significant revision of the cooling mechanism in astrophysical media.
Accurate method of modeling cluster scaling relations in modified gravity
NASA Astrophysics Data System (ADS)
He, Jian-hua; Li, Baojiu
2016-06-01
We propose a new method to model cluster scaling relations in modified gravity. Using a suite of nonradiative hydrodynamical simulations, we show that the scaling relations of accumulated gas quantities, such as the Sunyaev-Zel'dovich effect (Compton-y parameter) and the x-ray Compton-y parameter, can be accurately predicted using the known results in the Λ CDM model with a precision of ˜3 % . This method provides a reliable way to analyze the gas physics in modified gravity using the less demanding and much more efficient pure cold dark matter simulations. Our results therefore have important theoretical and practical implications in constraining gravity using cluster surveys.
Anisotropic Turbulence Modeling for Accurate Rod Bundle Simulations
Baglietto, Emilio
2006-07-01
An improved anisotropic eddy viscosity model has been developed for accurate predictions of the thermal hydraulic performances of nuclear reactor fuel assemblies. The proposed model adopts a non-linear formulation of the stress-strain relationship in order to include the reproduction of the anisotropic phenomena, and in combination with an optimized low-Reynolds-number formulation based on Direct Numerical Simulation (DNS) to produce correct damping of the turbulent viscosity in the near wall region. This work underlines the importance of accurate anisotropic modeling to faithfully reproduce the scale of the turbulence driven secondary flows inside the bundle subchannels, by comparison with various isothermal and heated experimental cases. The very low scale secondary motion is responsible for the increased turbulence transport which produces a noticeable homogenization of the velocity distribution and consequently of the circumferential cladding temperature distribution, which is of main interest in bundle design. Various fully developed bare bundles test cases are shown for different geometrical and flow conditions, where the proposed model shows clearly improved predictions, in close agreement with experimental findings, for regular as well as distorted geometries. Finally the applicability of the model for practical bundle calculations is evaluated through its application in the high-Reynolds form on coarse grids, with excellent results. (author)
Some analytical models of radiating collapsing spheres
Herrera, L.; Di Prisco, A; Ospino, J.
2006-08-15
We present some analytical solutions to the Einstein equations, describing radiating collapsing spheres in the diffusion approximation. Solutions allow for modeling physical reasonable situations. The temperature is calculated for each solution, using a hyperbolic transport equation, which permits to exhibit the influence of relaxational effects on the dynamics of the system.
Bond indices in solids: extended analytical model.
Ponec, Robert
2011-11-15
The analytical model suggested some time ago for the calculation of bond indices in infinite periodical structures was reconsidered and extended so as to provide not only realistic estimate of the extent of electron sharing localized among individual pairs of the atoms in the lattice but also to detect the eventual presence of multicenter bonding in metallic solids. PMID:21823136
NASA Astrophysics Data System (ADS)
Martínez, M. J.; Marco, F. J.; López, J. A.
2009-02-01
The Hipparcos catalog provides a reference frame at optical wavelengths for the new International Celestial Reference System (ICRS). This new reference system was adopted following the resolution agreed at the 23rd IAU General Assembly held in Kyoto in 1997. Differences in the Hipparcos system of proper motions and the previous materialization of the reference frame, the FK5, are expected to be caused only by the combined effects of the motion of the equinox of the FK5 and the precession of the equator and the ecliptic. Several authors have pointed out an inconsistency between the differences in proper motion of the Hipparcos-FK5 and the correction of the precessional values derived from VLBI and lunar laser ranging (LLR) observations. Most of them have claimed that these discrepancies are due to slightly biased proper motions in the FK5 catalog. The different mathematical models that have been employed to explain these errors have not fully accounted for the discrepancies in the correction of the precessional parameters. Our goal here is to offer an explanation for this fact. We propose the use of independent parametric and nonparametric models. The introduction of a nonparametric model, combined with the inner product in the square integrable functions over the unitary sphere, would give us values which do not depend on the possible interdependencies existing in the data set. The evidence shows that zonal studies are needed. This would lead us to introduce a local nonparametric model. All these models will provide independent corrections to the precessional values, which could then be compared in order to study the reliability in each case. Finally, we obtain values for the precession corrections that are very consistent with those that are currently adopted.
Analytical modeling for microwave and optical metasurfaces
NASA Astrophysics Data System (ADS)
Monti, Alessio; Soric, Jason; Alù, Andrea; Toscano, Alessandro; Bilotti, Filiberto
2016-06-01
A metasurface is an artificial structure composed by an ultrathin surface textured at a subwavelength scale. In the last years, metasurfaces have been revealed to be particularly useful in the design of electromagnetic scattering cancellation devices operating at microwave and optical frequencies. In this contribution we summarize our results about the analytical modelling of microwave and optical metasurfaces composed, respectively, by patterned metallic surfaces and arrays of plasmonic nanoparticles. The analytical results are compared with the numerical ones obtained with a proper set of full-wave simulations showing an excellent agreement.
Analytic Model for the Rototranslational Torsion Pendulum
NASA Astrophysics Data System (ADS)
De Marchi, F.; Bassan, M.; Pucacco, G.; Marconi, L.; Stanga, R.; Visco, M.
2013-01-01
We develop an analytic model to describe the motion of the RotoTranslational Torsion Pendulum PETER in a wide range of frequencies (from 1mHz up to 10-15Hz). We also try to explain some unexpected features we found in the data with only 1 soft degree of freedom and we estimate values for the misalignment angles and other parameters of the model.
Analytical model of internally coupled ears.
Vossen, Christine; Christensen-Dalsgaard, Jakob; van Hemmen, J Leo
2010-08-01
Lizards and many birds possess a specialized hearing mechanism: internally coupled ears where the tympanic membranes connect through a large mouth cavity so that the vibrations of the tympanic membranes influence each other. This coupling enhances the phase differences and creates amplitude differences in the tympanic membrane vibrations. Both cues show strong directionality. The work presented herein sets out the derivation of a three dimensional analytical model of internally coupled ears that allows for calculation of a complete vibration profile of the membranes. The analytical model additionally provides the opportunity to incorporate the effect of the asymmetrically attached columella, which leads to the activation of higher membrane vibration modes. Incorporating this effect, the analytical model can explain measurements taken from the tympanic membrane of a living lizard, for example, data demonstrating an asymmetrical spatial pattern of membrane vibration. As the analytical calculations show, the internally coupled ears increase the directional response, appearing in large directional internal amplitude differences (iAD) and in large internal time differences (iTD). Numerical simulations of the eigenfunctions in an exemplary, realistically reconstructed mouth cavity further estimate the effects of its complex geometry. PMID:20707461
Automated statistical modeling of analytical measurement systems
Jacobson, J J
1992-08-01
The statistical modeling of analytical measurement systems at the Idaho Chemical Processing Plant (ICPP) has been completely automated through computer software. The statistical modeling of analytical measurement systems is one part of a complete quality control program used by the Remote Analytical Laboratory (RAL) at the ICPP. The quality control program is an integration of automated data input, measurement system calibration, database management, and statistical process control. The quality control program and statistical modeling program meet the guidelines set forth by the American Society for Testing Materials and American National Standards Institute. A statistical model is a set of mathematical equations describing any systematic bias inherent in a measurement system and the precision of a measurement system. A statistical model is developed from data generated from the analysis of control standards. Control standards are samples which are made up at precise known levels by an independent laboratory and submitted to the RAL. The RAL analysts who process control standards do not know the values of those control standards. The object behind statistical modeling is to describe real process samples in terms of their bias and precision and, to verify that a measurement system is operating satisfactorily. The processing of control standards gives us this ability.
Fast micromagnetic simulations using an analytic mathematical model
NASA Astrophysics Data System (ADS)
Tsiantos, Vassilios; Miles, Jim
2006-02-01
In this paper an analytic mathematical model is presented for fast micromagnetic simulations. In dynamic micromagnetic simulations the Landau-Lifshitz-Gilbert (LLG) equation is solved for the observation of the reversal magnetisation mechanisms. In stiff micromagnetic simulations the large system of ordinary differential equations has to be solved with an appropriate method, such as the Backward Differentiation Formulas (BDF) method, which leads to the solution of a large linear system. The latter is solved efficiently employing matrix-free techniques, such as Krylov methods with preconditioning. Within the Krylov methods framework a product of a matrix times a vector is involved which is usually approximated with directional differences. This paper provides an analytic mathematical model to calculate efficiently this product, leading to more accurate calculations and consequently faster micromagnetic simulations due to better convergence properties.
Accurate pressure gradient calculations in hydrostatic atmospheric models
NASA Technical Reports Server (NTRS)
Carroll, John J.; Mendez-Nunez, Luis R.; Tanrikulu, Saffet
1987-01-01
A method for the accurate calculation of the horizontal pressure gradient acceleration in hydrostatic atmospheric models is presented which is especially useful in situations where the isothermal surfaces are not parallel to the vertical coordinate surfaces. The present method is shown to be exact if the potential temperature lapse rate is constant between the vertical pressure integration limits. The technique is applied to both the integration of the hydrostatic equation and the computation of the slope correction term in the horizontal pressure gradient. A fixed vertical grid and a dynamic grid defined by the significant levels in the vertical temperature distribution are employed.
Mouse models of human AML accurately predict chemotherapy response.
Zuber, Johannes; Radtke, Ina; Pardee, Timothy S; Zhao, Zhen; Rappaport, Amy R; Luo, Weijun; McCurrach, Mila E; Yang, Miao-Miao; Dolan, M Eileen; Kogan, Scott C; Downing, James R; Lowe, Scott W
2009-04-01
The genetic heterogeneity of cancer influences the trajectory of tumor progression and may underlie clinical variation in therapy response. To model such heterogeneity, we produced genetically and pathologically accurate mouse models of common forms of human acute myeloid leukemia (AML) and developed methods to mimic standard induction chemotherapy and efficiently monitor therapy response. We see that murine AMLs harboring two common human AML genotypes show remarkably diverse responses to conventional therapy that mirror clinical experience. Specifically, murine leukemias expressing the AML1/ETO fusion oncoprotein, associated with a favorable prognosis in patients, show a dramatic response to induction chemotherapy owing to robust activation of the p53 tumor suppressor network. Conversely, murine leukemias expressing MLL fusion proteins, associated with a dismal prognosis in patients, are drug-resistant due to an attenuated p53 response. Our studies highlight the importance of genetic information in guiding the treatment of human AML, functionally establish the p53 network as a central determinant of chemotherapy response in AML, and demonstrate that genetically engineered mouse models of human cancer can accurately predict therapy response in patients. PMID:19339691
An accurate model potential for alkali neon systems.
Zanuttini, D; Jacquet, E; Giglio, E; Douady, J; Gervais, B
2009-12-01
We present a detailed investigation of the ground and lowest excited states of M-Ne dimers, for M=Li, Na, and K. We show that the potential energy curves of these Van der Waals dimers can be obtained accurately by considering the alkali neon systems as one-electron systems. Following previous authors, the model describes the evolution of the alkali valence electron in the combined potentials of the alkali and neon cores by means of core polarization pseudopotentials. The key parameter for an accurate model is the M(+)-Ne potential energy curve, which was obtained by means of ab initio CCSD(T) calculation using a large basis set. For each MNe dimer, a systematic comparison with ab initio computation of the potential energy curve for the X, A, and B states shows the remarkable accuracy of the model. The vibrational analysis and the comparison with existing experimental data strengthens this conclusion and allows for a precise assignment of the vibrational levels. PMID:19968334
Mouse models of human AML accurately predict chemotherapy response
Zuber, Johannes; Radtke, Ina; Pardee, Timothy S.; Zhao, Zhen; Rappaport, Amy R.; Luo, Weijun; McCurrach, Mila E.; Yang, Miao-Miao; Dolan, M. Eileen; Kogan, Scott C.; Downing, James R.; Lowe, Scott W.
2009-01-01
The genetic heterogeneity of cancer influences the trajectory of tumor progression and may underlie clinical variation in therapy response. To model such heterogeneity, we produced genetically and pathologically accurate mouse models of common forms of human acute myeloid leukemia (AML) and developed methods to mimic standard induction chemotherapy and efficiently monitor therapy response. We see that murine AMLs harboring two common human AML genotypes show remarkably diverse responses to conventional therapy that mirror clinical experience. Specifically, murine leukemias expressing the AML1/ETO fusion oncoprotein, associated with a favorable prognosis in patients, show a dramatic response to induction chemotherapy owing to robust activation of the p53 tumor suppressor network. Conversely, murine leukemias expressing MLL fusion proteins, associated with a dismal prognosis in patients, are drug-resistant due to an attenuated p53 response. Our studies highlight the importance of genetic information in guiding the treatment of human AML, functionally establish the p53 network as a central determinant of chemotherapy response in AML, and demonstrate that genetically engineered mouse models of human cancer can accurately predict therapy response in patients. PMID:19339691
Turbulence Models for Accurate Aerothermal Prediction in Hypersonic Flows
NASA Astrophysics Data System (ADS)
Zhang, Xiang-Hong; Wu, Yi-Zao; Wang, Jiang-Feng
Accurate description of the aerodynamic and aerothermal environment is crucial to the integrated design and optimization for high performance hypersonic vehicles. In the simulation of aerothermal environment, the effect of viscosity is crucial. The turbulence modeling remains a major source of uncertainty in the computational prediction of aerodynamic forces and heating. In this paper, three turbulent models were studied: the one-equation eddy viscosity transport model of Spalart-Allmaras, the Wilcox k-ω model and the Menter SST model. For the k-ω model and SST model, the compressibility correction, press dilatation and low Reynolds number correction were considered. The influence of these corrections for flow properties were discussed by comparing with the results without corrections. In this paper the emphasis is on the assessment and evaluation of the turbulence models in prediction of heat transfer as applied to a range of hypersonic flows with comparison to experimental data. This will enable establishing factor of safety for the design of thermal protection systems of hypersonic vehicle.
Analytical steam injection model for layered systems
Abdual-Razzaq; Brigham, W.E.; Castanier, L.M.
1993-08-01
Screening, evaluation and optimization of the steam flooding process in homogeneous reservoirs can be performed by using simple analytical predictive models. In the absence of any analytical model for layered reservoirs, at present, only numerical simulators can be used. And these are expensive. In this study, an analytical model has been developed considering two isolated layers of differing permeabilities. The principle of equal flow potential is applied across the two layers. Gajdica`s (1990) single layer linear steam drive model is extended for the layered system. The formulation accounts for variation of heat loss area in the higher permeability layer, and the development of a hot liquid zone in the lower permeability layer. These calculations also account for effects of viscosity, density, fractional flow curves and pressure drops in the hot liquid zone. Steam injection rate variations in the layers are represented by time weighted average rates. For steam zone calculations, Yortsos and Gavalas`s (1981) upper bound method is used with a correction factor. The results of the model are compared with a numerical simulator. Comparable oil and water flow rates, and breakthrough times were achieved for 100 cp oil. Results with 10 cp and 1000 cp oils indicate the need to improve the formulation to properly handle differing oil viscosities.
NASA Astrophysics Data System (ADS)
Lloyd, N. S.; Bouman, C.; Horstwood, M. S.; Parrish, R. R.; Schwieters, J. B.
2010-12-01
This presentation describes progress in mass spectrometry for analysing very small analyte quantities, illustrated by example applications from nuclear forensics. In this challenging application, precise and accurate (‰) uranium isotope ratios are required from 1 - 2 µm diameter uranium oxide particles, which comprise less than 40 pg of uranium. Traditionally these are analysed using thermal ionisation mass spectrometry (TIMS), and more recently using secondary ionisation mass spectrometry (SIMS). Multicollector inductively-coupled plasma mass spectrometry (MC-ICP-MS) can offer higher productivity compared to these techniques, but is traditionally limited by low efficiency of analyte utilisation (sample through to ion detection). Samples can either be introduced as a solution, or sampled directly from solid using laser ablation. Large multi-isotope ratio datasets can help identify provenance and intended use of anthropogenic uranium and other nuclear materials [1]. The Thermo Scientific NEPTUNE Plus (Bremen, Germany) with ‘Jet Interface’ option offers unparalleled MC-ICP-MS sensitivity. An analyte utilisation of c. 4% has previously been reported for uranium [2]. This high-sensitivity configuration utilises a dry high-capacity (100 m3/h) interface pump, special skimmer and sampler cones and a desolvating nebuliser system. Coupled with new acquisition methodologies, this sensitivity enhancement makes possible the analysis of micro-particles and small sample volumes at higher precision levels than previously achieved. New, high-performance, full-size and compact discrete dynode secondary electron multipliers (SEM) exhibit excellent stability and linearity over a large dynamic range and can be configured to simultaneously measure all of the uranium isotopes. Options for high abundance-sensitivity filters on two ion beams are also available, e.g. for 236U and 234U. Additionally, amplifiers with high ohm (1012 - 1013) feedback resistors have been developed to
A parsimonious analytical model for simulating multispecies plume migration
NASA Astrophysics Data System (ADS)
Chen, J.-S.; Liang, C.-P.; Liu, C.-W.; Li, L. Y.
2015-09-01
A parsimonious analytical model for rapidly predicting the long-term plume behavior of decaying contaminant such as radionuclide and dissolved chlorinated solvent is presented in this study. Generalized analytical solutions in compact format are derived for the two-dimensional advection-dispersion equations coupled with sequential first-order decay reactions involving an arbitrary number of species in groundwater system. The solution techniques involve the sequential applications of the Laplace, finite Fourier cosine, and generalized integral transforms to reduce the coupled partial differential equation system to a set of linear algebraic equations. The system of algebraic equations is next solved for each species in the transformed domain, and the solutions in the original domain are then obtained through consecutive integral transform inversions. Explicit form solutions for a special case are derived using the generalized analytical solutions and are verified against the numerical solutions. The analytical results indicate that the parsimonious analytical solutions are robust and accurate. The solutions are useful for serving as simulation or screening tools for assessing plume behaviors of decaying contaminants including the radionuclides and dissolved chlorinated solvents in groundwater systems.
Some analytical models of anisotropic strange stars
NASA Astrophysics Data System (ADS)
Murad, Mohammad Hassan
2016-01-01
Over the years of the concept of local isotropy has become a too stringent condition in modeling relativistic self-gravitating objects. Taking local anisotropy into consideration, in this work, some analytical models of relativistic anisotropic charged strange stars have been developed. The Einstein-Maxwell gravitational field equations have been solved with a particular form of one of the metric potentials. The radial pressure and the energy density have been assumed to follow the usual linear equation of state of strange quark matter, the MIT bag model.
Combining Modeling and Gaming for Predictive Analytics
Riensche, Roderick M.; Whitney, Paul D.
2012-08-22
Many of our most significant challenges involve people. While human behavior has long been studied, there are recent advances in computational modeling of human behavior. With advances in computational capabilities come increases in the volume and complexity of data that humans must understand in order to make sense of and capitalize on these modeling advances. Ultimately, models represent an encapsulation of human knowledge. One inherent challenge in modeling is efficient and accurate transfer of knowledge from humans to models, and subsequent retrieval. The simulated real-world environment of games presents one avenue for these knowledge transfers. In this paper we describe our approach of combining modeling and gaming disciplines to develop predictive capabilities, using formal models to inform game development, and using games to provide data for modeling.
An analytical model of memristors in plants
Markin, Vladislav S; Volkov, Alexander G; Chua, Leon
2014-01-01
The memristor, a resistor with memory, was postulated by Chua in 1971 and the first solid-state memristor was built in 2008. Recently, we found memristors in vivo in plants. Here we propose a simple analytical model of 2 types of memristors that can be found within plants. The electrostimulation of plants by bipolar periodic waves induces electrical responses in the Aloe vera and Mimosa pudica with fingerprints of memristors. Memristive properties of the Aloe vera and Mimosa pudica are linked to the properties of voltage gated K+ ion channels. The potassium channel blocker TEACl transform plant memristors to conventional resistors. The analytical model of a memristor with a capacitor connected in parallel exhibits different characteristic behavior at low and high frequency of applied voltage, which is the same as experimental data obtained by cyclic voltammetry in vivo. PMID:25482769
Generating Facial Expressions Using an Anatomically Accurate Biomechanical Model.
Wu, Tim; Hung, Alice; Mithraratne, Kumar
2014-11-01
This paper presents a computational framework for modelling the biomechanics of human facial expressions. A detailed high-order (Cubic-Hermite) finite element model of the human head was constructed using anatomical data segmented from magnetic resonance images. The model includes a superficial soft-tissue continuum consisting of skin, the subcutaneous layer and the superficial Musculo-Aponeurotic system. Embedded within this continuum mesh, are 20 pairs of facial muscles which drive facial expressions. These muscles were treated as transversely-isotropic and their anatomical geometries and fibre orientations were accurately depicted. In order to capture the relative composition of muscles and fat, material heterogeneity was also introduced into the model. Complex contact interactions between the lips, eyelids, and between superficial soft tissue continuum and deep rigid skeletal bones were also computed. In addition, this paper investigates the impact of incorporating material heterogeneity and contact interactions, which are often neglected in similar studies. Four facial expressions were simulated using the developed model and the results were compared with surface data obtained from a 3D structured-light scanner. Predicted expressions showed good agreement with the experimental data. PMID:26355331
Analytical Solution of Traffic Cellular Automata Model
NASA Astrophysics Data System (ADS)
Lo, Shih-Ching; Hsu, Chia-Hung
2009-08-01
Complex traffic system seems to be simulated successfully by cellular automaton (CA) models. Various models are developed to understand single-lane traffic, multilane traffic, lane-changing behavior and network traffic situations. However, the result of CA simulation can only be obtained after massive microscopic computation. Although, the mean field theory (MFT) has been studied to be the approximation of CA model, the MFT can only applied to the simple CA rules or small value of parameters. In this study, we simulate traffic flow by the NaSch model under different combination of parameters, which are maximal speed, dawdling probability and density. After that, the position of critical density, the slope of free-flow and congested regime are observed and modeled due to the simulated data. Finally, the coefficients of the model will be calibrated by the simulated data and the analytical solution of traffic CA is obtained.
Accurate verification of the conserved-vector-current and standard-model predictions
Sirlin, A.; Zucchini, R.
1986-10-20
An approximate analytic calculation of O(Z..cap alpha../sup 2/) corrections to Fermi decays is presented. When the analysis of Koslowsky et al. is modified to take into account the new results, it is found that each of the eight accurately studied scrFt values differs from the average by approx. <1sigma, thus significantly improving the comparison of experiments with conserved-vector-current predictions. The new scrFt values are lower than before, which also brings experiments into very good agreement with the three-generation standard model, at the level of its quantum corrections.
An Analytic Function of Lunar Surface Temperature for Exospheric Modeling
NASA Technical Reports Server (NTRS)
Hurley, Dana M.; Sarantos, Menelaos; Grava, Cesare; Williams, Jean-Pierre; Retherford, Kurt D.; Siegler, Matthew; Greenhagen, Benjamin; Paige, David
2014-01-01
We present an analytic expression to represent the lunar surface temperature as a function of Sun-state latitude and local time. The approximation represents neither topographical features nor compositional effects and therefore does not change as a function of selenographic latitude and longitude. The function reproduces the surface temperature measured by Diviner to within +/-10 K at 72% of grid points for dayside solar zenith angles of less than 80, and at 98% of grid points for nightside solar zenith angles greater than 100. The analytic function is least accurate at the terminator, where there is a strong gradient in the temperature, and the polar regions. Topographic features have a larger effect on the actual temperature near the terminator than at other solar zenith angles. For exospheric modeling the effects of topography on the thermal model can be approximated by using an effective longitude for determining the temperature. This effective longitude is randomly redistributed with 1 sigma of 4.5deg. The resulting ''roughened'' analytical model well represents the statistical dispersion in the Diviner data and is expected to be generally useful for future models of lunar surface temperature, especially those implemented within exospheric simulations that address questions of volatile transport.
An Improved Analytic Model for Microdosimeter Response
NASA Technical Reports Server (NTRS)
Shinn, Judy L.; Wilson, John W.; Xapsos, Michael A.
2001-01-01
An analytic model used to predict energy deposition fluctuations in a microvolume by ions through direct events is improved to include indirect delta ray events. The new model can now account for the increase in flux at low lineal energy when the ions are of very high energy. Good agreement is obtained between the calculated results and available data for laboratory ion beams. Comparison of GCR (galactic cosmic ray) flux between Shuttle TEPC (tissue equivalent proportional counter) flight data and current calculations draws a different assessment of developmental work required for the GCR transport code (HZETRN) than previously concluded.
Inverter Modeling For Accurate Energy Predictions Of Tracking HCPV Installations
NASA Astrophysics Data System (ADS)
Bowman, J.; Jensen, S.; McDonald, Mark
2010-10-01
High efficiency high concentration photovoltaic (HCPV) solar plants of megawatt scale are now operational, and opportunities for expanded adoption are plentiful. However, effective bidding for sites requires reliable prediction of energy production. HCPV module nameplate power is rated for specific test conditions; however, instantaneous HCPV power varies due to site specific irradiance and operating temperature, and is degraded by soiling, protective stowing, shading, and electrical connectivity. These factors interact with the selection of equipment typically supplied by third parties, e.g., wire gauge and inverters. We describe a time sequence model accurately accounting for these effects that predicts annual energy production, with specific reference to the impact of the inverter on energy output and interactions between system-level design decisions and the inverter. We will also show two examples, based on an actual field design, of inverter efficiency calculations and the interaction between string arrangements and inverter selection.
NASA Astrophysics Data System (ADS)
Roy Choudhury, Raja; Roy Choudhury, Arundhati; Kanti Ghose, Mrinal
2013-01-01
A semi-analytical model with three optimizing parameters and a novel non-Gaussian function as the fundamental modal field solution has been proposed to arrive at an accurate solution to predict various propagation parameters of graded-index fibers with less computational burden than numerical methods. In our semi analytical formulation the optimization of core parameter U which is usually uncertain, noisy or even discontinuous, is being calculated by Nelder-Mead method of nonlinear unconstrained minimizations as it is an efficient and compact direct search method and does not need any derivative information. Three optimizing parameters are included in the formulation of fundamental modal field of an optical fiber to make it more flexible and accurate than other available approximations. Employing variational technique, Petermann I and II spot sizes have been evaluated for triangular and trapezoidal-index fibers with the proposed fundamental modal field. It has been demonstrated that, the results of the proposed solution identically match with the numerical results over a wide range of normalized frequencies. This approximation can also be used in the study of doped and nonlinear fiber amplifier.
Analytical modeling of worldwide medical radiation use
Mettler, F.A. Jr.; Davis, M.; Kelsey, C.A.; Rosenberg, R.; Williams, A.
1987-02-01
An analytical model was developed to estimate the availability and frequency of medical radiation use on a worldwide basis. This model includes medical and dental x-ray, nuclear medicine, and radiation therapy. The development of an analytical model is necessary as the first step in estimating the radiation dose to the world's population from this source. Since there is no data about the frequency of medical radiation use in more than half the countries in the world and only fragmentary data in an additional one-fourth of the world's countries, such a model can be used to predict the uses of medical radiation in these countries. The model indicates that there are approximately 400,000 medical x-ray machines worldwide and that approximately 1.2 billion diagnostic medical x-ray examinations are performed annually. Dental x-ray examinations are estimated at 315 million annually and approximately 22 million in-vivo diagnostic nuclear medicine examinations. Approximately 4 million radiation therapy procedures or courses of treatment are undertaken annually.
Accurate, low-cost 3D-models of gullies
NASA Astrophysics Data System (ADS)
Onnen, Nils; Gronz, Oliver; Ries, Johannes B.; Brings, Christine
2015-04-01
Soil erosion is a widespread problem in arid and semi-arid areas. The most severe form is the gully erosion. They often cut into agricultural farmland and can make a certain area completely unproductive. To understand the development and processes inside and around gullies, we calculated detailed 3D-models of gullies in the Souss Valley in South Morocco. Near Taroudant, we had four study areas with five gullies different in size, volume and activity. By using a Canon HF G30 Camcorder, we made varying series of Full HD videos with 25fps. Afterwards, we used the method Structure from Motion (SfM) to create the models. To generate accurate models maintaining feasible runtimes, it is necessary to select around 1500-1700 images from the video, while the overlap of neighboring images should be at least 80%. In addition, it is very important to avoid selecting photos that are blurry or out of focus. Nearby pixels of a blurry image tend to have similar color values. That is why we used a MATLAB script to compare the derivatives of the images. The higher the sum of the derivative, the sharper an image of similar objects. MATLAB subdivides the video into image intervals. From each interval, the image with the highest sum is selected. E.g.: 20min. video at 25fps equals 30.000 single images. The program now inspects the first 20 images, saves the sharpest and moves on to the next 20 images etc. Using this algorithm, we selected 1500 images for our modeling. With VisualSFM, we calculated features and the matches between all images and produced a point cloud. Then, MeshLab has been used to build a surface out of it using the Poisson surface reconstruction approach. Afterwards we are able to calculate the size and the volume of the gullies. It is also possible to determine soil erosion rates, if we compare the data with old recordings. The final step would be the combination of the terrestrial data with the data from our aerial photography. So far, the method works well and we
Towards Accurate Molecular Modeling of Plastic Bonded Explosives
NASA Astrophysics Data System (ADS)
Chantawansri, T. L.; Andzelm, J.; Taylor, D.; Byrd, E.; Rice, B.
2010-03-01
There is substantial interest in identifying the controlling factors that influence the susceptibility of polymer bonded explosives (PBXs) to accidental initiation. Numerous Molecular Dynamics (MD) simulations of PBXs using the COMPASS force field have been reported in recent years, where the validity of the force field in modeling the solid EM fill has been judged solely on its ability to reproduce lattice parameters, which is an insufficient metric. Performance of the COMPASS force field in modeling EMs and the polymeric binder has been assessed by calculating structural, thermal, and mechanical properties, where only fair agreement with experimental data is obtained. We performed MD simulations using the COMPASS force field for the polymer binder hydroxyl-terminated polybutadiene and five EMs: cyclotrimethylenetrinitramine, 1,3,5,7-tetranitro-1,3,5,7-tetra-azacyclo-octane, 2,4,6,8,10,12-hexantirohexaazazisowurzitane, 2,4,6-trinitro-1,3,5-benzenetriamine, and pentaerythritol tetranitate. Predicted EM crystallographic and molecular structural parameters, as well as calculated properties for the binder will be compared with experimental results for different simulation conditions. We also present novel simulation protocols, which improve agreement between experimental and computation results thus leading to the accurate modeling of PBXs.
Analytic Models of Plausible Gravitational Lens Potentials
Baltz, Edward A.; Marshall, Phil; Oguri, Masamune
2007-05-04
Gravitational lenses on galaxy scales are plausibly modeled as having ellipsoidal symmetry and a universal dark matter density profile, with a Sersic profile to describe the distribution of baryonic matter. Predicting all lensing effects requires knowledge of the total lens potential: in this work we give analytic forms for that of the above hybrid model. Emphasizing that complex lens potentials can be constructed from simpler components in linear combination, we provide a recipe for attaining elliptical symmetry in either projected mass or lens potential.We also provide analytic formulae for the lens potentials of Sersic profiles for integer and half-integer index. We then present formulae describing the gravitational lensing effects due to smoothly-truncated universal density profiles in cold dark matter model. For our isolated haloes the density profile falls off as radius to the minus fifth or seventh power beyond the tidal radius, functional forms that allow all orders of lens potential derivatives to be calculated analytically, while ensuring a non-divergent total mass. We show how the observables predicted by this profile differ from that of the original infinite-mass NFW profile. Expressions for the gravitational flexion are highlighted. We show how decreasing the tidal radius allows stripped haloes to be modeled, providing a framework for a fuller investigation of dark matter substructure in galaxies and clusters. Finally we remark on the need for finite mass halo profiles when doing cosmological ray-tracing simulations, and the need for readily-calculable higher order derivatives of the lens potential when studying catastrophes in strong lenses.
Analytic models of plausible gravitational lens potentials
Baltz, Edward A.; Marshall, Phil; Oguri, Masamune E-mail: pjm@physics.ucsb.edu
2009-01-15
Gravitational lenses on galaxy scales are plausibly modelled as having ellipsoidal symmetry and a universal dark matter density profile, with a Sersic profile to describe the distribution of baryonic matter. Predicting all lensing effects requires knowledge of the total lens potential: in this work we give analytic forms for that of the above hybrid model. Emphasising that complex lens potentials can be constructed from simpler components in linear combination, we provide a recipe for attaining elliptical symmetry in either projected mass or lens potential. We also provide analytic formulae for the lens potentials of Sersic profiles for integer and half-integer index. We then present formulae describing the gravitational lensing effects due to smoothly-truncated universal density profiles in cold dark matter model. For our isolated haloes the density profile falls off as radius to the minus fifth or seventh power beyond the tidal radius, functional forms that allow all orders of lens potential derivatives to be calculated analytically, while ensuring a non-divergent total mass. We show how the observables predicted by this profile differ from that of the original infinite-mass NFW profile. Expressions for the gravitational flexion are highlighted. We show how decreasing the tidal radius allows stripped haloes to be modelled, providing a framework for a fuller investigation of dark matter substructure in galaxies and clusters. Finally we remark on the need for finite mass halo profiles when doing cosmological ray-tracing simulations, and the need for readily-calculable higher order derivatives of the lens potential when studying catastrophes in strong lenses.
Accurate modelling of flow induced stresses in rigid colloidal aggregates
NASA Astrophysics Data System (ADS)
Vanni, Marco
2015-07-01
A method has been developed to estimate the motion and the internal stresses induced by a fluid flow on a rigid aggregate. The approach couples Stokesian dynamics and structural mechanics in order to take into account accurately the effect of the complex geometry of the aggregates on hydrodynamic forces and the internal redistribution of stresses. The intrinsic error of the method, due to the low-order truncation of the multipole expansion of the Stokes solution, has been assessed by comparison with the analytical solution for the case of a doublet in a shear flow. In addition, it has been shown that the error becomes smaller as the number of primary particles in the aggregate increases and hence it is expected to be negligible for realistic reproductions of large aggregates. The evaluation of internal forces is performed by an adaptation of the matrix methods of structural mechanics to the geometric features of the aggregates and to the particular stress-strain relationship that occurs at intermonomer contacts. A preliminary investigation on the stress distribution in rigid aggregates and their mode of breakup has been performed by studying the response to an elongational flow of both realistic reproductions of colloidal aggregates (made of several hundreds monomers) and highly simplified structures. A very different behaviour has been evidenced between low-density aggregates with isostatic or weakly hyperstatic structures and compact aggregates with highly hyperstatic configuration. In low-density clusters breakup is caused directly by the failure of the most stressed intermonomer contact, which is typically located in the inner region of the aggregate and hence originates the birth of fragments of similar size. On the contrary, breakup of compact and highly cross-linked clusters is seldom caused by the failure of a single bond. When this happens, it proceeds through the removal of a tiny fragment from the external part of the structure. More commonly, however
An accurate and simple quantum model for liquid water.
Paesani, Francesco; Zhang, Wei; Case, David A; Cheatham, Thomas E; Voth, Gregory A
2006-11-14
The path-integral molecular dynamics and centroid molecular dynamics methods have been applied to investigate the behavior of liquid water at ambient conditions starting from a recently developed simple point charge/flexible (SPC/Fw) model. Several quantum structural, thermodynamic, and dynamical properties have been computed and compared to the corresponding classical values, as well as to the available experimental data. The path-integral molecular dynamics simulations show that the inclusion of quantum effects results in a less structured liquid with a reduced amount of hydrogen bonding in comparison to its classical analog. The nuclear quantization also leads to a smaller dielectric constant and a larger diffusion coefficient relative to the corresponding classical values. Collective and single molecule time correlation functions show a faster decay than their classical counterparts. Good agreement with the experimental measurements in the low-frequency region is obtained for the quantum infrared spectrum, which also shows a higher intensity and a redshift relative to its classical analog. A modification of the original parametrization of the SPC/Fw model is suggested and tested in order to construct an accurate quantum model, called q-SPC/Fw, for liquid water. The quantum results for several thermodynamic and dynamical properties computed with the new model are shown to be in a significantly better agreement with the experimental data. Finally, a force-matching approach was applied to the q-SPC/Fw model to derive an effective quantum force field for liquid water in which the effects due to the nuclear quantization are explicitly distinguished from those due to the underlying molecular interactions. Thermodynamic and dynamical properties computed using standard classical simulations with this effective quantum potential are found in excellent agreement with those obtained from significantly more computationally demanding full centroid molecular dynamics
Analytical Performance Models for Geologic Repositories
Chambre, P.L.; Pigford, T.H.; Fujita, A.; Kanki, T.; Kobayashi,A.; Lung, H.; Ting, D.; Sato, Y.; Savoshy, S.J.
1982-10-01
This report presents analytical solutions of the dissolution and hydrogeologic transport of radionuclides in geologic repositories. Numerical examples are presented to demonstrate the equations resulting from these analyses. The subjects treated in the present report are: (a) Solubility-limited transport with transverse dispersion (Chapter 2); (b) Transport of a radionuclide chain with nonequilibrium chemical reactions (Chapter 3); (c) Advective transport in a two-dimensional flow field (Chapter 4); (d) Radionuclide.transport in fractured media (Chapter 5); (e) A mathematical model for EPA's analysis of generic repositories (Chapter 6); and (f) Dissolution of radionuclides from solid waste (Chapter 7).
NASA Astrophysics Data System (ADS)
Ukaegbu, Ikechi Augustine; Sangirov, Jamshid; Cho, Mu Hee; Lee, Tae-Woo; Park, Hyo-Hoon
2011-07-01
In this paper, a crosstalk expression and equivalent circuit model have been proposed based on RLC line model and interconnect parameters for wire-bonded and flip-chip bonded multichannel optoelectronic modules. The analytical expression and model are accurate for computing crosstalk of interconnects used in chip packaging. In addition, full-wave simulation and experimental results from total crosstalk measurement are discussed.
Analytical Ion Thruster Discharge Performance Model
NASA Technical Reports Server (NTRS)
Goebel, Dan M.; Wirz, Richard E.; Katz, Ira
2006-01-01
A particle and energy balance model of the plasma discharge in magnetic ring-cusp ion thrusters has been developed. The model follows the original work of Brophy in the development of global 0-D discharge models that utilize conservation of particles into and out of the thruster and conservation of energy into the discharge and out of the plasma in the form of charged particles to the walls and beam and plasma radiation. The present model is significantly expanded over Brophy's original work by including self-consistent calculations of the internal neutral pressure, electron temperature, primary electron density, electrostatic ion confinement (due to the ring-cusp fields), plasma potential, discharge stability, and time dependent behavior during recycling. The model only requires information on the thruster geometry, ion optics performance and electrical inputs such as discharge voltage and currents, etc. to produce accurate performance curves of discharge loss versus mass utilization efficiency. The model has been benchmarked against the NEXIS Laboratory Model (LM) and Development Model (DM) thrusters, and successfully predicts the thruster discharge loss as a function of mass utilization efficiency for a variety of thrusters. The discharge performance model will be presented and results showing ion thruster performance and stability given.
Multilayer analytic element modeling of radial collector wells.
Bakker, Mark; Kelson, Victor A; Luther, Kenneth H
2005-01-01
A new multilayer approach is presented for the modeling of ground water flow to radial collector wells. The approach allows for the inclusion of all aspects of the unique boundary condition along the lateral arms of a collector well, including skin effect and internal friction losses due to flow in the arms. The hydraulic conductivity may differ between horizontal layers within the aquifer, and vertical anisotropy can be taken into account. The approach is based on the multilayer analytic element method, such that regional flow and local three-dimensional detail may be simulated simultaneously and accurately within one regional model. Horizontal flow inside a layer is computed analytically, while vertical flow is approximated with a standard finite-difference scheme. Results obtained with the proposed approach compare well to results obtained with three-dimensional analytic element solutions for flow in unconfined aquifers. The presented approach may be applied to predict the yield of a collector well in a regional setting and to compute the origin and residence time, and thus the quality, of water pumped by the collector well. As an example, the addition of three lateral arms to a collector well that already has three laterals is investigated. The new arms are added at an elevation of 2 m above the existing laterals. The yield increase of the collector well is computed as a function of the lengths of the three new arms. PMID:16324013
Analytical model of batch magnetophoretic separation
NASA Astrophysics Data System (ADS)
Kashevsky, S. B.; Kashevsky, B. E.
2013-06-01
Magnetophoresis (the motion of magnetic particles driven by the nonuniform magnetic field), that for a long time has been used for extracting magnetically susceptible objects in diverse industries, now attracts interest to develop more sophisticated microfluidic and batch techniques for separation and manipulation of biological particles, and magnetically assisted absorption and catalysis in organic chemistry, biochemistry, and petrochemistry. A deficiency of magnetic separation science is the lack of simple analytical models imitating real processes of magnetic separation. We have studied the motion of superparamagnetic (generally, soft magnetic) particles in liquid in the three-dimensional field of the diametrically polarized permanent cylindrical magnet; this geometry is basically representative of the batch separation mode. In the limit of the infinite-length magnet, we found the particle magnetophoresis proceeds independently of the magnet polarization direction, following the simple analytical relation incorporating all the relevant physical and geometrical parameters of the particle-magnet system. In experiments with a finite-length magnet we have shown applicability of the developed theory as to analyze the performance of the real batch separation systems in the noncooperative mode, and finally, we have presented an example of such analysis for the case of immunomagnetic cell separation and developed a criterion of the model limitation imposed by the magnetic aggregation of particles.
Analytical models for complex swirling flows
NASA Astrophysics Data System (ADS)
Borissov, A.; Hussain, V.
1996-11-01
We develops a new class of analytical solutions of the Navier-Stokes equations for swirling flows, and suggests ways to predict and control such flows occurring in various technological applications. We view momentum accumulation on the axis as a key feature of swirling flows and consider vortex-sink flows on curved axisymmetric surfaces with an axial flow. We show that these solutions model swirling flows in a cylindrical can, whirlpools, tornadoes, and cosmic swirling jets. The singularity of these solutions on the flow axis is removed by matching them with near-axis Schlichting and Long's swirling jets. The matched solutions model flows with very complex patterns, consisting of up to seven separation regions with recirculatory 'bubbles' and vortex rings. We apply the matched solutions for computing flows in the Ranque-Hilsch tube, in the meniscus of electrosprays, in vortex breakdown, and in an industrial vortex burner. The simple analytical solutions allow a clear understanding of how different control parameters affect the flow and guide selection of optimal parameter values for desired flow features. These solutions permit extension to other problems (such as heat transfer and chemical reaction) and have the potential of being significantly useful for further detailed investigation by direct or large-eddy numerical simulations as well as laboratory experimentation.
An analytical model of a longitudinal-torsional ultrasonic transducer
NASA Astrophysics Data System (ADS)
Al-Budairi, Hassan; Lucas, Margaret
2012-08-01
The combination of longitudinal and torsional (LT) vibrations at high frequencies finds many applications such as ultrasonic drilling, ultrasonic welding, and ultrasonic motors. The LT mode can be obtained by modifications to the design of a standard bolted Langevin ultrasonic transducer driven by an axially poled piezoceramic stack, by a technique that degenerates the longitudinal mode to an LT motion by a geometrical alteration of the wave path. The transducer design is developed and optimised through numerical modelling which can represent the geometry and mechanical properties of the transducer and its vibration response to an electrical input applied across the piezoceramic stack. However, although these models can allow accurate descriptions of the mechanical behaviour, they do not generally provide adequate insights into the electrical characteristics of the transducer. In this work, an analytical model is developed to present the LT transducer based on the equivalent circuit method. This model can represent both the mechanical and electrical aspects and is used to extract many of the design parameters, such as resonance and anti-resonance frequencies, the impedance spectra and the coupling coefficient of the transducer. The validity of the analytical model is demonstrated by close agreement with experimental results.
Coarse-grained red blood cell model with accurate mechanical properties, rheology and dynamics.
Fedosov, Dmitry A; Caswell, Bruce; Karniadakis, George E
2009-01-01
We present a coarse-grained red blood cell (RBC) model with accurate and realistic mechanical properties, rheology and dynamics. The modeled membrane is represented by a triangular mesh which incorporates shear inplane energy, bending energy, and area and volume conservation constraints. The macroscopic membrane elastic properties are imposed through semi-analytic theory, and are matched with those obtained in optical tweezers stretching experiments. Rheological measurements characterized by time-dependent complex modulus are extracted from the membrane thermal fluctuations, and compared with those obtained from the optical magnetic twisting cytometry results. The results allow us to define a meaningful characteristic time of the membrane. The dynamics of RBCs observed in shear flow suggests that a purely elastic model for the RBC membrane is not appropriate, and therefore a viscoelastic model is required. The set of proposed analyses and numerical tests can be used as a complete model testbed in order to calibrate the modeled viscoelastic membranes to accurately represent RBCs in health and disease. PMID:19965026
An analytically enriched finite element method for cohesive crack modeling.
Cox, James V.
2010-04-01
Meaningful computational investigations of many solid mechanics problems require accurate characterization of material behavior through failure. A recent approach to fracture modeling has combined the partition of unity finite element method (PUFEM) with cohesive zone models. Extension of the PUFEM to address crack propagation is often referred to as the extended finite element method (XFEM). In the PUFEM, the displacement field is enriched to improve the local approximation. Most XFEM studies have used simplified enrichment functions (e.g., generalized Heaviside functions) to represent the strong discontinuity but have lacked an analytical basis to represent the displacement gradients in the vicinity of the cohesive crack. As such, the mesh had to be sufficiently fine for the FEM basis functions to capture these gradients.In this study enrichment functions based upon two analytical investigations of the cohesive crack problem are examined. These functions have the potential of representing displacement gradients in the vicinity of the cohesive crack with a relatively coarse mesh and allow the crack to incrementally advance across each element. Key aspects of the corresponding numerical formulation are summarized. Analysis results for simple model problems are presented to evaluate if quasi-static crack propagation can be accurately followed with the proposed formulation. A standard finite element solution with interface elements is used to provide the accurate reference solution, so the model problems are limited to a straight, mode I crack in plane stress. Except for the cohesive zone, the material model for the problems is homogenous, isotropic linear elasticity. The effects of mesh refinement, mesh orientation, and enrichment schemes that enrich a larger region around the cohesive crack are considered in the study. Propagation of the cohesive zone tip and crack tip, time variation of the cohesive zone length, and crack profiles are presented. The analysis
Comparison between analytical and numerical solution of mathematical drying model
NASA Astrophysics Data System (ADS)
Shahari, N.; Rasmani, K.; Jamil, N.
2016-02-01
Drying is often related to the food industry as a process of shifting heat and mass inside food, which helps in preserving food. Previous research using a mass transfer equation showed that the results were mostly concerned with the comparison between the simulation model and the experimental data. In this paper, the finite difference method was used to solve a mass equation during drying using different kinds of boundary condition, which are equilibrium and convective boundary conditions. The results of these two models provide a comparison between the analytical and the numerical solution. The result shows a close match between the two solution curves. It is concluded that the two proposed models produce an accurate solution to describe the moisture distribution content during the drying process. This analysis indicates that we have confidence in the behaviour of moisture in the numerical simulation. This result demonstrated that a combined analytical and numerical approach prove that the system is behaving physically. Based on this assumption, the model of mass transfer was extended to include the temperature transfer, and the result shows a similar trend to those presented in the simpler case.
The Immediate Exchange model: an analytical investigation
NASA Astrophysics Data System (ADS)
Katriel, Guy
2015-01-01
We study the Immediate Exchange model, recently introduced by Heinsalu and Patriarca [Eur. Phys. J. B 87, 170 (2014)], who showed by simulations that the wealth distribution in this model converges to a Gamma distribution with shape parameter 2. Here we justify this conclusion analytically, in the infinite-population limit. An infinite-population version of the model is derived, describing the evolution of the wealth distribution in terms of iterations of a nonlinear operator on the space of probability densities. It is proved that the Gamma distributions with shape parameter 2 are fixed points of this operator, and that, starting with an arbitrary wealth distribution, the process converges to one of these fixed points. We also discuss the mixed model introduced in the same paper, in which exchanges are either bidirectional or unidirectional with fixed probability. We prove that, although, as found by Heinsalu and Patriarca, the equilibrium distribution can be closely fit by Gamma distributions, the equilibrium distribution for this model is not a Gamma distribution.
ANALYTIC MODELING OF THE MORETON WAVE KINEMATICS
Temmer, M.; Veronig, A. M.
2009-09-10
The issue whether Moreton waves are flare-ignited or coronal mass ejection (CME)-driven, or a combination of both, is still a matter of debate. We develop an analytical model describing the evolution of a large-amplitude coronal wave emitted by the expansion of a circular source surface in order to mimic the evolution of a Moreton wave. The model results are confronted with observations of a strong Moreton wave observed in association with the X3.8/3B flare/CME event from 2005 January 17. Using different input parameters for the expansion of the source region, either derived from the real CME observations (assuming that the upward moving CME drives the wave), or synthetically generated scenarios (expanding flare region, lateral expansion of the CME flanks), we calculate the kinematics of the associated Moreton wave signature. Those model input parameters are determined which fit the observed Moreton wave kinematics best. Using the measured kinematics of the upward moving CME as the model input, we are not able to reproduce the observed Moreton wave kinematics. The observations of the Moreton wave can be reproduced only by applying a strong and impulsive acceleration for the source region expansion acting in a piston mechanism scenario. Based on these results we propose that the expansion of the flaring region or the lateral expansion of the CME flanks is more likely the driver of the Moreton wave than the upward moving CME front.
Analytic Modeling of Collector Current and Delay Time in Hbts
NASA Astrophysics Data System (ADS)
Jung, Hee-Bum
1992-01-01
Collector current in abrupt Al_ {0.48}In_{0.52} As/In_{0.53}Ga _{0.47}As HBTs is investigated. Because tunneling plays an important role for abrupt heterojunctions, thermionic field emission (TF) mechanism is included, as a part of the model, in addition to thermionic emission (TE) theory. To model the modulation of the effective barrier height correctly, non-ideal doping profile across the heterojunction is considered. Calculations showed that under nominal operating conditions, TF is dominant over TE in determining the collector current. Furthermore, modulation of the effective barrier height manifests itself in the collector ideality factor that is greater than unity. It is shown that, by calculating the above mentioned transport mechanisms and including the barrier height modulation, the collector current and its temperature dependence in abrupt AlInAs/InGaAs HBTs can be predicted correctly. The detailed calculation is reduced to an analytical closed -form model by assuming a Gaussian energy spectrum for TF current. The model is determined to be accurate over a wide range of bias and temperatures. A simple TE/TF Ebers -Moll model for abrupt HBTs is derived. The classical expression for collector small signal delay time is inadequate for vertically scaled transistors where transient velocity effects can no longer be ignored. Analytical expressions for collector transit time and small signal delay time are proposed for circuit simulation. These models use a general non-uniform velocity profile described entirely in terms of five physical parameters: momentum and energy relaxation times, and initial, peak, and saturated velocities. A C_infty-continuous function approximation for the transit time is used to obtain analytical closed-form expressions for collector small signal delay time in terms of physically meaningful transport parameters. An accurate empirical two-piece model is also proposed. As the collector thickness is scaled down, the ratio of small signal
Analytical model of Europa's O2 exosphere
NASA Astrophysics Data System (ADS)
Milillo, Anna; Plainaki, Christina; Orsini, Stefano; De Angelis, Elisabetta; Mangano, Valeria; Massetti, Stefano; Mura, Alessandro; Rispoli, Rosanna; Colasanti, Luca
2015-04-01
The origin of the exosphere of Europa is its water ice surface. The existing exosphere models, assuming either a collisionless environment (simple Monte Carlo techniques) or a kinetic approach (Direct Monte Carlo Method) both predict that the major constituent of the exosphere is molecular oxygen. Specifically, O2 is generated at the surface through radiolysis and chemical interactions of the water dissociation products. The non-escaping O2 molecules circulate around the moon impacting the surface several times, due to their long lifetime and due to their non- sticking, suffering thermalization to the surface temperature after each impact. In fact, the HST observations of the O emission lines have manifested the presence of an asymmetric atomic Oxygen envelope, evidencing the possible existence of a thin asymmetric molecular Oxygen atmosphere. The existing Monte Carlo models are not easily applicable as input of simulations devoted to the study of the plasma interactions with the moon. On the other hand, the simple exponential density profiles cannot well depict the higher temperature/higher altitudes component originating by radiolysis. On the contrary, it would be important to have a suitable and user-friendly model to use as a tool. This study presents an analytical 3D model that is able to describe the molecular Oxygen exosphere by reproducing the asymmetries due to two configurations among Europa, Jupiter and the Sun, that is illumination at leading and at trailing side. This model is obtained by a non-linear fit procedure of the EGEON Monte Carlo model to a Chamberlain density profile. Different parameters of the model are able to describe various exosphere properties thus allowing a detailed investigation of the exospheric characteristics.
Improved analytical model for residual stress prediction in orthogonal cutting
NASA Astrophysics Data System (ADS)
Qi, Zhaoxu; Li, Bin; Xiong, Liangshan
2014-09-01
The analytical model of residual stress in orthogonal cutting proposed by Jiann is an important tool for residual stress prediction in orthogonal cutting. In application of the model, a problem of low precision of the surface residual stress prediction is found. By theoretical analysis, several shortages of Jiann's model are picked out, including: inappropriate boundary conditions, unreasonable calculation method of thermal stress, ignorance of stress constraint and cyclic loading algorithm. These shortages may directly lead to the low precision of the surface residual stress prediction. To eliminate these shortages and make the prediction more accurate, an improved model is proposed. In this model, a new contact boundary condition between tool and workpiece is used to make it in accord with the real cutting process; an improved calculation method of thermal stress is adopted; a stress constraint is added according to the volumeconstancy of plastic deformation; and the accumulative effect of the stresses during cyclic loading is considered. At last, an experiment for measuring residual stress in cutting AISI 1045 steel is conducted. Also, Jiann's model and the improved model are simulated under the same conditions with cutting experiment. The comparisons show that the surface residual stresses predicted by the improved model is closer to the experimental results than the results predicted by Jiann's model.
Improved analytical model for residual stress prediction in orthogonal cutting
NASA Astrophysics Data System (ADS)
Qi, Zhaoxu; Li, Bin; Xiong, Liangshan
2014-09-01
The analytical model of residual stress in orthogonal cutting proposed by Jiann is an important tool for residual stress prediction in orthogonal cutting. In application of the model, a problem of low precision of the surface residual stress prediction is found. By theoretical analysis, several shortages of Jiann's model are picked out, including: inappropriate boundary conditions, unreasonable calculation method of thermal stress, ignorance of stress constraint and cyclic loading algorithm. These shortages may directly lead to the low precision of the surface residual stress prediction. To eliminate these shortages and make the prediction more accurate, an improved model is proposed. In this model, a new contact boundary condition between tool and workpiece is used to make it in accord with the real cutting process; an improved calculation method of thermal stress is adopted; a stress constraint is added according to the volume-constancy of plastic deformation; and the accumulative effect of the stresses during cyclic loading is considered. At last, an experiment for measuring residual stress in cutting AISI 1045 steel is conducted. Also, Jiann's model and the improved model are simulated under the same conditions with cutting experiment. The comparisons show that the surface residual stresses predicted by the improved model is closer to the experimental results than the results predicted by Jiann's model.
WHAEM: PROGRAM DOCUMENTATION FOR THE WELLHEAD ANALYTIC ELEMENT MODEL
The Wellhead Analytic Element Model (WhAEM) demonstrates a new technique for the definition of time-of-travel capture zones in relatively simple geohydrologic settings. he WhAEM package includes an analytic element model that uses superposition of (many) analytic solutions to gen...
Lumb, Matthew P.; Steiner, Myles A.; Geisz, John F.; Walters, Robert J.
2014-11-21
The analytical drift-diffusion formalism is able to accurately simulate a wide range of solar cell architectures and was recently extended to include those with back surface reflectors. However, as solar cells approach the limits of material quality, photon recycling effects become increasingly important in predicting the behavior of these cells. In particular, the minority carrier diffusion length is significantly affected by the photon recycling, with consequences for the solar cell performance. In this paper, we outline an approach to account for photon recycling in the analytical Hovel model and compare analytical model predictions to GaAs-based experimental devices operating close to the fundamental efficiency limit.
Applying an accurate spherical model to gamma-ray burst afterglow observations
NASA Astrophysics Data System (ADS)
Leventis, K.; van der Horst, A. J.; van Eerten, H. J.; Wijers, R. A. M. J.
2013-05-01
We present results of model fits to afterglow data sets of GRB 970508, GRB 980703 and GRB 070125, characterized by long and broad-band coverage. The model assumes synchrotron radiation (including self-absorption) from a spherical adiabatic blast wave and consists of analytic flux prescriptions based on numerical results. For the first time it combines the accuracy of hydrodynamic simulations through different stages of the outflow dynamics with the flexibility of simple heuristic formulas. The prescriptions are especially geared towards accurate description of the dynamical transition of the outflow from relativistic to Newtonian velocities in an arbitrary power-law density environment. We show that the spherical model can accurately describe the data only in the case of GRB 970508, for which we find a circumburst medium density n ∝ r-2. We investigate in detail the implied spectra and physical parameters of that burst. For the microphysics we show evidence for equipartition between the fraction of energy density carried by relativistic electrons and magnetic field. We also find that for the blast wave to be adiabatic, the fraction of electrons accelerated at the shock has to be smaller than 1. We present best-fitting parameters for the afterglows of all three bursts, including uncertainties in the parameters of GRB 970508, and compare the inferred values to those obtained by different authors.
NASA Astrophysics Data System (ADS)
Smith, R.; Flynn, C.; Candlish, G. N.; Fellhauer, M.; Gibson, B. K.
2015-04-01
We present accurate models of the gravitational potential produced by a radially exponential disc mass distribution. The models are produced by combining three separate Miyamoto-Nagai discs. Such models have been used previously to model the disc of the Milky Way, but here we extend this framework to allow its application to discs of any mass, scalelength, and a wide range of thickness from infinitely thin to near spherical (ellipticities from 0 to 0.9). The models have the advantage of simplicity of implementation, and we expect faster run speeds over a double exponential disc treatment. The potentials are fully analytical, and differentiable at all points. The mass distribution of our models deviates from the radial mass distribution of a pure exponential disc by <0.4 per cent out to 4 disc scalelengths, and <1.9 per cent out to 10 disc scalelengths. We tabulate fitting parameters which facilitate construction of exponential discs for any scalelength, and a wide range of disc thickness (a user-friendly, web-based interface is also available). Our recipe is well suited for numerical modelling of the tidal effects of a giant disc galaxy on star clusters or dwarf galaxies. We consider three worked examples; the Milky Way thin and thick disc, and a discy dwarf galaxy.
Analytic Ballistic Performance Model of Whipple Shields
NASA Technical Reports Server (NTRS)
Miller, J. E.; Bjorkman, M. D.; Christiansen, E. L.; Ryan, S. J.
2015-01-01
The dual-wall, Whipple shield is the shield of choice for lightweight, long-duration flight. The shield uses an initial sacrificial wall to initiate fragmentation and melt an impacting threat that expands over a void before hitting a subsequent shield wall of a critical component. The key parameters to this type of shield are the rear wall and its mass which stops the debris, as well as the minimum shock wave strength generated by the threat particle impact of the sacrificial wall and the amount of room that is available for expansion. Ensuring the shock wave strength is sufficiently high to achieve large scale fragmentation/melt of the threat particle enables the expansion of the threat and reduces the momentum flux of the debris on the rear wall. Three key factors in the shock wave strength achieved are the thickness of the sacrificial wall relative to the characteristic dimension of the impacting particle, the density and material cohesion contrast of the sacrificial wall relative to the threat particle and the impact speed. The mass of the rear wall and the sacrificial wall are desirable to minimize for launch costs making it important to have an understanding of the effects of density contrast and impact speed. An analytic model is developed here, to describe the influence of these three key factors. In addition this paper develops a description of a fourth key parameter related to fragmentation and its role in establishing the onset of projectile expansion.
Analytic wave model of Stark deceleration dynamics
Gubbels, Koos; Meijer, Gerard; Friedrich, Bretislav
2006-06-15
Stark deceleration relies on time-dependent inhomogeneous electric fields which repetitively exert a decelerating force on polar molecules. Fourier analysis reveals that such fields, generated by an array of field stages, consist of a superposition of partial waves with well-defined phase velocities. Molecules whose velocities come close to the phase velocity of a given wave get a ride from that wave. For a square-wave temporal dependence of the Stark field, the phase velocities of the waves are found to be odd-fraction multiples of a fundamental phase velocity {lambda}/{tau}, with {lambda} and {tau} the spatial and temporal periods of the field. Here we study explicitly the dynamics due to any of the waves as well as due to their mutual perturbations. We first solve the equations of motion for the case of single-wave interactions and exploit their isomorphism with those for the biased pendulum. Next we analyze the perturbations of the single-wave dynamics by other waves and find that these have no net effect on the phase stability of the acceleration or deceleration process. Finally, we find that a packet of molecules can also ride a wave which results from an interference of adjacent waves. In this case, small phase stability areas form around phase velocities that are even-fraction multiples of the fundamental velocity. A detailed comparison with classical trajectory simulations and with experiment demonstrates that the analytic 'wave model' encompasses all the longitudinal physics encountered in a Stark decelerator.
Digitalized accurate modeling of SPCB with multi-spiral surface based on CPC algorithm
NASA Astrophysics Data System (ADS)
Huang, Yanhua; Gu, Lizhi
2015-09-01
The main methods of the existing multi-spiral surface geometry modeling include spatial analytic geometry algorithms, graphical method, interpolation and approximation algorithms. However, there are some shortcomings in these modeling methods, such as large amount of calculation, complex process, visible errors, and so on. The above methods have, to some extent, restricted the design and manufacture of the premium and high-precision products with spiral surface considerably. This paper introduces the concepts of the spatially parallel coupling with multi-spiral surface and spatially parallel coupling body. The typical geometry and topological features of each spiral surface forming the multi-spiral surface body are determined, by using the extraction principle of datum point cluster, the algorithm of coupling point cluster by removing singular point, and the "spatially parallel coupling" principle based on the non-uniform B-spline for each spiral surface. The orientation and quantitative relationships of datum point cluster and coupling point cluster in Euclidean space are determined accurately and in digital description and expression, coupling coalescence of the surfaces with multi-coupling point clusters under the Pro/E environment. The digitally accurate modeling of spatially parallel coupling body with multi-spiral surface is realized. The smooth and fairing processing is done to the three-blade end-milling cutter's end section area by applying the principle of spatially parallel coupling with multi-spiral surface, and the alternative entity model is processed in the four axis machining center after the end mill is disposed. And the algorithm is verified and then applied effectively to the transition area among the multi-spiral surface. The proposed model and algorithms may be used in design and manufacture of the multi-spiral surface body products, as well as in solving essentially the problems of considerable modeling errors in computer graphics and
Analytical model of diffuse reflectance spectrum of skin tissue
NASA Astrophysics Data System (ADS)
Lisenko, S. A.; Kugeiko, M. M.; Firago, V. A.; Sobchuk, A. N.
2014-01-01
We have derived simple analytical expressions that enable highly accurate calculation of diffusely reflected light signals of skin in the spectral range from 450 to 800 nm at a distance from the region of delivery of exciting radiation. The expressions, taking into account the dependence of the detected signals on the refractive index, transport scattering coefficient, absorption coefficient and anisotropy factor of the medium, have been obtained in the approximation of a two-layer medium model (epidermis and dermis) for the same parameters of light scattering but different absorption coefficients of layers. Numerical experiments on the retrieval of the skin biophysical parameters from the diffuse reflectance spectra simulated by the Monte Carlo method show that commercially available fibre-optic spectrophotometers with a fixed distance between the radiation source and detector can reliably determine the concentration of bilirubin, oxy- and deoxyhaemoglobin in the dermis tissues and the tissue structure parameter characterising the size of its effective scatterers. We present the examples of quantitative analysis of the experimental data, confirming the correctness of estimates of biophysical parameters of skin using the obtained analytical expressions.
Analytical model of diffuse reflectance spectrum of skin tissue
Lisenko, S A; Kugeiko, M M; Firago, V A; Sobchuk, A N
2014-01-31
We have derived simple analytical expressions that enable highly accurate calculation of diffusely reflected light signals of skin in the spectral range from 450 to 800 nm at a distance from the region of delivery of exciting radiation. The expressions, taking into account the dependence of the detected signals on the refractive index, transport scattering coefficient, absorption coefficient and anisotropy factor of the medium, have been obtained in the approximation of a two-layer medium model (epidermis and dermis) for the same parameters of light scattering but different absorption coefficients of layers. Numerical experiments on the retrieval of the skin biophysical parameters from the diffuse reflectance spectra simulated by the Monte Carlo method show that commercially available fibre-optic spectrophotometers with a fixed distance between the radiation source and detector can reliably determine the concentration of bilirubin, oxy- and deoxyhaemoglobin in the dermis tissues and the tissue structure parameter characterising the size of its effective scatterers. We present the examples of quantitative analysis of the experimental data, confirming the correctness of estimates of biophysical parameters of skin using the obtained analytical expressions. (biophotonics)
Analytical modeling of Cosmic Winds and Jets
NASA Astrophysics Data System (ADS)
Vlahakis, Nektarios
1998-11-01
stellar wind and the Blandford and Payne (1982) model of a disk-wind; it also contains nonpolytropic models, such as those of winds/jets in Sauty and Tsinganos (1994), Lima et al (1996) and Trussoni et al (1997). Besides the unification of all known cases under a common scheme, several new classes emerge and some are briefly analyzed; they could be explored for a further understanding of the physical properties of MHD outflows from various magnetized astrophysical rotators. We also propose a new class of exact and self-consistent MHD solutions which describe steady and axisymmetric hydromagnetic outflows from the magnetized atmosphere of a rotating gravitating central object with possibly an orbiting accretion disk. The plasma is driven by a thermal pressure gradient, as well as by magnetic rotator and radiative forces. At the Alfvenic and fast critical points the appropriate criticality conditions are applied. The outflows start almost radially but after the Alfven transition and before the fast critical surface is encountered the magnetic pinching force bends the poloidal streamlines into a cylindrical jet-type shape. The terminal speed, Alfven number, cross-sectional area of the jet, as well as its final pressure and density obtain uniform values at large distances from the source. The goal of the study is to give an analytical discussion of the two-dimensional interplay of the thermal pressure gradient, gravitational, Lorentz and inertial forces in accelerating and collimating an MHD flow. A parametric study of the model is given, as well as a brief sketch of its applicability to a self-consistent modeling of collimated outflows from various astrophysical objects. For example, the obtained characteristics of the collimated outflow in agreement with those in jets associated with YSO's. General theoretical arguments and various analytic self-similar solutions have recently shown that magnetized and rotating astrophysical outflows may become asymptotically cylindrical
New process model proves accurate in tests on catalytic reformer
Aguilar-Rodriguez, E.; Ancheyta-Juarez, J. )
1994-07-25
A mathematical model has been devised to represent the process that takes place in a fixed-bed, tubular, adiabatic catalytic reforming reactor. Since its development, the model has been applied to the simulation of a commercial semiregenerative reformer. The development of mass and energy balances for this reformer led to a model that predicts both concentration and temperature profiles along the reactor. A comparison of the model's results with experimental data illustrates its accuracy at predicting product profiles. Simple steps show how the model can be applied to simulate any fixed-bed catalytic reformer.
Numerical-analytical modeling of the Earth's pole oscillations
NASA Astrophysics Data System (ADS)
Markov, Y.; Filippova, A.
2015-08-01
For the purpose of more accurate forecasting the oscillatory process of the Earth pole in time periods with significant anomalies (irregular deviations) a numerical-analytical approach is presented for the combined modeling of the interdependent dynamical processes - the oscillatory-rotational motion of the Earth and the time dependent coefficients of the geopotential. The oscillations of the inertia tensor components of the Earth depend on various factors such as mechanical and physical parameters of the planet, the motions of the tide-generating bodies and observed large scale natural events. Time variations of these and some other factors affect the Earth orientation parameters. The generalization of the previously researched mathematical model of Chandler and annual oscillations of the Earth pole is being held with the use of celestial mechanics methods and the mathematical description of the Earth gravitational field's temporal variations. The latter makes possible to improve the forecast precision of the Earth pole trajectory. Also the more precise model is to have small number of parameters and to agree with the previously developed one (to have the same structural features and to have a correspondence between the averaged dynamical parameters and the parameters of the basic model).
Analytical model for CO(2) laser ablation of fused quartz.
Nowak, Krzysztof M; Baker, Howard J; Hall, Denis R
2015-10-10
This paper reports the development of an analytical model, with supporting experimental data, which quite accurately describes the key features of CO_{2} laser ablation of fused silica glass. The quantitative model of nonexplosive, evaporative material removal is shown to match the experimental data very well, to the extent that it can be used as a tool for ablative measurements of absorption coefficient and vaporization energy. The experimental results indicated that a minimum of 12 MJ kg^{-1} is required to fully vaporize fused quartz initially held at room temperature, which is in good agreement with the prediction of the model supplied with input data available in the literature. An optimal window for the machining of fused quartz was revealed in terms of pulse duration 20-80 μs and CO_{2} laser wavelength optimized for maximum absorption coefficient. Material removal rates of 0.33 μm per J cm^{-2} allow for a high-precision depth control with modest laser stability. The model may also be used as a parameter selection guide for CO_{2} laser ablation of fused silica or other materials of similar thermophysical properties. PMID:26479800
Coupling Efforts to the Accurate and Efficient Tsunami Modelling System
NASA Astrophysics Data System (ADS)
Son, S.
2015-12-01
In the present study, we couple two different types of tsunami models, i.e., nondispersive shallow water model of characteristic form(MOST ver.4) and dispersive Boussinesq model of non-characteristic form(Son et al. (2011)) in an attempt to improve modelling accuracy and efficiency. Since each model deals with different type of primary variables, additional care on matching boundary condition is required. Using an absorbing-generating boundary condition developed by Van Dongeren and Svendsen(1997), model coupling and integration is achieved. Characteristic variables(i.e., Riemann invariants) in MOST are converted to non-characteristic variables for Boussinesq solver without any loss of physical consistency. Established modelling system has been validated through typical test problems to realistic tsunami events. Simulated results reveal good performance of developed modelling system. Since coupled modelling system provides advantageous flexibility feature during implementation, great efficiencies and accuracies are expected to be gained through spot-focusing application of Boussinesq model inside the entire domain of tsunami propagation.
Accurate modelling of anisotropic effects in austenitic stainless steel welds
NASA Astrophysics Data System (ADS)
Nowers, O. D.; Duxbury, D. J.; Drinkwater, B. W.
2014-02-01
The ultrasonic inspection of austenitic steel welds is challenging due to the formation of highly anisotropic and heterogeneous structures post-welding. This is due to the intrinsic crystallographic structure of austenitic steel, driving the formation of dendritic grain structures on cooling. The anisotropy is manifested as both a `steering' of the ultrasonic beam and the back-scatter of energy due to the macroscopic granular structure of the weld. However, the quantitative effects and relative impacts of these phenomena are not well-understood. A semi-analytical simulation framework has been developed to allow the study of anisotropic effects in austenitic stainless steel welds. Frequency-dependent scatterers are allocated to a weld-region to approximate the coarse grain-structures observed within austenitic welds and imaged using a simulated array. The simulated A-scans are compared against an equivalent experimental setup demonstrating excellent agreement of the Signal to Noise (S/N) ratio. Comparison of images of the simulated and experimental data generated using the Total Focusing Method (TFM) indicate a prominent layered effect in the simulated data. A superior grain allocation routine is required to improve upon this.
Accurate modelling of anisotropic effects in austenitic stainless steel welds
Nowers, O. D.; Duxbury, D. J.; Drinkwater, B. W.
2014-02-18
The ultrasonic inspection of austenitic steel welds is challenging due to the formation of highly anisotropic and heterogeneous structures post-welding. This is due to the intrinsic crystallographic structure of austenitic steel, driving the formation of dendritic grain structures on cooling. The anisotropy is manifested as both a ‘steering’ of the ultrasonic beam and the back-scatter of energy due to the macroscopic granular structure of the weld. However, the quantitative effects and relative impacts of these phenomena are not well-understood. A semi-analytical simulation framework has been developed to allow the study of anisotropic effects in austenitic stainless steel welds. Frequency-dependent scatterers are allocated to a weld-region to approximate the coarse grain-structures observed within austenitic welds and imaged using a simulated array. The simulated A-scans are compared against an equivalent experimental setup demonstrating excellent agreement of the Signal to Noise (S/N) ratio. Comparison of images of the simulated and experimental data generated using the Total Focusing Method (TFM) indicate a prominent layered effect in the simulated data. A superior grain allocation routine is required to improve upon this.
Analytical Model of Nano-Electromechanical (NEM) Nonvolatile Memory Cells
NASA Astrophysics Data System (ADS)
Han, Boram; Choi, Woo Young
The fringe field effects of nano-electromechanical (NEM) nonvolatile memory cells have been investigated analytically for the accurate evaluation of NEM memory cells. As the beam width is scaled down, fringe field effect becomes more severe. It has been observed that pull-in, release and hysteresis voltage decrease more than our prediction. Also, the fringe field on cell characteristics has been discussed.
Linaro, Daniele; Storace, Marco; Giugliano, Michele
2011-01-01
Stochastic channel gating is the major source of intrinsic neuronal noise whose functional consequences at the microcircuit- and network-levels have been only partly explored. A systematic study of this channel noise in large ensembles of biophysically detailed model neurons calls for the availability of fast numerical methods. In fact, exact techniques employ the microscopic simulation of the random opening and closing of individual ion channels, usually based on Markov models, whose computational loads are prohibitive for next generation massive computer models of the brain. In this work, we operatively define a procedure for translating any Markov model describing voltage- or ligand-gated membrane ion-conductances into an effective stochastic version, whose computer simulation is efficient, without compromising accuracy. Our approximation is based on an improved Langevin-like approach, which employs stochastic differential equations and no Montecarlo methods. As opposed to an earlier proposal recently debated in the literature, our approximation reproduces accurately the statistical properties of the exact microscopic simulations, under a variety of conditions, from spontaneous to evoked response features. In addition, our method is not restricted to the Hodgkin-Huxley sodium and potassium currents and is general for a variety of voltage- and ligand-gated ion currents. As a by-product, the analysis of the properties emerging in exact Markov schemes by standard probability calculus enables us for the first time to analytically identify the sources of inaccuracy of the previous proposal, while providing solid ground for its modification and improvement we present here. PMID:21423712
Health Informatics for Neonatal Intensive Care Units: An Analytical Modeling Perspective.
Khazaei, Hamzeh; Mench-Bressan, Nadja; McGregor, Carolyn; Pugh, James Edward
2015-01-01
The effective use of data within intensive care units (ICUs) has great potential to create new cloud-based health analytics solutions for disease prevention or earlier condition onset detection. The Artemis project aims to achieve the above goals in the area of neonatal ICUs (NICU). In this paper, we proposed an analytical model for the Artemis cloud project which will be deployed at McMaster Children's Hospital in Hamilton. We collect not only physiological data but also the infusion pumps data that are attached to NICU beds. Using the proposed analytical model, we predict the amount of storage, memory, and computation power required for the system. Capacity planning and tradeoff analysis would be more accurate and systematic by applying the proposed analytical model in this paper. Numerical results are obtained using real inputs acquired from McMaster Children's Hospital and a pilot deployment of the system at The Hospital for Sick Children (SickKids) in Toronto. PMID:27170907
Analytic Models of High-Temperature Hohlraums
Stygar, W.A.; Olson, R.E.; Spielman, R.B.; Leeper, R.J.
2000-11-29
A unified set of high-temperature-hohlraum models has been developed. For a simple hohlraum, P{sub s} = [A{sub s}+(1{minus}{alpha}{sub W})A{sub W}+A{sub H}]{sigma}T{sub R}{sup 4} + (4V{sigma}/c)(dT{sub R}{sup r}/dt) where P{sub S} is the total power radiated by the source, A{sub s} is the source area, A{sub W} is the area of the cavity wall excluding the source and holes in the wall, A{sub H} is the area of the holes, {sigma} is the Stefan-Boltzmann constant, T{sub R} is the radiation brightness temperature, V is the hohlraum volume, and c is the speed of light. The wall albedo {alpha}{sub W} {triple_bond} (T{sub W}/T{sub R}){sup 4} where T{sub W} is the brightness temperature of area A{sub W}. The net power radiated by the source P{sub N} = P{sub S}-A{sub S}{sigma}T{sub R}{sup 4}, which suggests that for laser-driven hohlraums the conversion efficiency {eta}{sub CE} be defined as P{sub N}/P{sub LASER}. The characteristic time required to change T{sub R}{sup 4} in response to a change in P{sub N} is 4V/C[(l{minus}{alpha}{sub W})A{sub W}+A{sub H}]. Using this model, T{sub R}, {alpha}{sub W}, and {eta}{sub CE} can be expressed in terms of quantities directly measurable in a hohlraum experiment. For a steady-state hohlraum that encloses a convex capsule, P{sub N} = {l_brace}(1{minus}{alpha}{sub W})A{sub W}+A{sub H}+[(1{minus}{alpha}{sub C})(A{sub S}+A{sub W}{alpha}{sub W})A{sub C}/A{sub T}]{r_brace}{sigma}T{sub RC}{sup 4} where {alpha}{sub C} is the capsule albedo, A{sub C} is the capsule area, A{sub T} {triple_bond} (A{sub S}+A{sub W}+A{sub H}), and T{sub RC} is the brightness temperature of the radiation that drives the capsule. According to this relation, the capsule-coupling efficiency of the baseline National-Ignition-Facility (NIF) hohlraum is 15% higher than predicted by previous analytic expressions. A model of a hohlraum that encloses a z pinch is also presented.
Magnetic field models of nine CP stars from "accurate" measurements
NASA Astrophysics Data System (ADS)
Glagolevskij, Yu. V.
2013-01-01
The dipole models of magnetic fields in nine CP stars are constructed based on the measurements of metal lines taken from the literature, and performed by the LSD method with an accuracy of 10-80 G. The model parameters are compared with the parameters obtained for the same stars from the hydrogen line measurements. For six out of nine stars the same type of structure was obtained. Some parameters, such as the field strength at the poles B p and the average surface magnetic field B s differ considerably in some stars due to differences in the amplitudes of phase dependences B e (Φ) and B s (Φ), obtained by different authors. It is noted that a significant increase in the measurement accuracy has little effect on the modelling of the large-scale structures of the field. By contrast, it is more important to construct the shape of the phase dependence based on a fairly large number of field measurements, evenly distributed by the rotation period phases. It is concluded that the Zeeman component measurement methods have a strong effect on the shape of the phase dependence, and that the measurements of the magnetic field based on the lines of hydrogen are more preferable for modelling the large-scale structures of the field.
Project Summary. ANALYTICAL ELEMENT MODELING OF COASTAL AQUIFERS
Four topics were studied concerning the modeling of groundwater flow in coastal aquifers with analytic elements: (1) practical experience was obtained by constructing a groundwater model of the shallow aquifers below the Delmarva Peninsula USA using the commercial program MVAEM; ...
An Accurate In Vitro Model of the E. coli Envelope
Clifton, Luke A; Holt, Stephen A; Hughes, Arwel V; Daulton, Emma L; Arunmanee, Wanatchaporn; Heinrich, Frank; Khalid, Syma; Jefferies, Damien; Charlton, Timothy R; Webster, John R P; Kinane, Christian J; Lakey, Jeremy H
2015-01-01
Gram-negative bacteria are an increasingly serious source of antibiotic-resistant infections, partly owing to their characteristic protective envelope. This complex, 20 nm thick barrier includes a highly impermeable, asymmetric bilayer outer membrane (OM), which plays a pivotal role in resisting antibacterial chemotherapy. Nevertheless, the OM molecular structure and its dynamics are poorly understood because the structure is difficult to recreate or study in vitro. The successful formation and characterization of a fully asymmetric model envelope using Langmuir–Blodgett and Langmuir–Schaefer methods is now reported. Neutron reflectivity and isotopic labeling confirmed the expected structure and asymmetry and showed that experiments with antibacterial proteins reproduced published in vivo behavior. By closely recreating natural OM behavior, this model provides a much needed robust system for antibiotic development. PMID:26331292
Leidenfrost effect: accurate drop shape modeling and new scaling laws
NASA Astrophysics Data System (ADS)
Sobac, Benjamin; Rednikov, Alexey; Dorbolo, Stéphane; Colinet, Pierre
2014-11-01
In this study, we theoretically investigate the shape of a drop in a Leidenfrost state, focusing on the geometry of the vapor layer. The drop geometry is modeled by numerically matching the solution of the hydrostatic shape of a superhydrophobic drop (for the upper part) with the solution of the lubrication equation of the vapor flow underlying the drop (for the bottom part). The results highlight that the vapor layer, fed by evaporation, forms a concave depression in the drop interface that becomes increasingly marked with the drop size. The vapor layer then consists of a gas pocket in the center and a thin annular neck surrounding it. The film thickness increases with the size of the drop, and the thickness at the neck appears to be of the order of 10--100 μm in the case of water. The model is compared to recent experimental results [Burton et al., Phys. Rev. Lett., 074301 (2012)] and shows an excellent agreement, without any fitting parameter. New scaling laws also emerge from this model. The geometry of the vapor pocket is only weakly dependent on the superheat (and thus on the evaporation rate), this weak dependence being more pronounced in the neck region. In turn, the vapor layer characteristics strongly depend on the drop size.
Automated dynamic analytical model improvement for damped structures
NASA Technical Reports Server (NTRS)
Fuh, J. S.; Berman, A.
1985-01-01
A method is described to improve a linear nonproportionally damped analytical model of a structure. The procedure finds the smallest changes in the analytical model such that the improved model matches the measured modal parameters. Features of the method are: (1) ability to properly treat complex valued modal parameters of a damped system; (2) applicability to realistically large structural models; and (3) computationally efficiency without involving eigensolutions and inversion of a large matrix.
Analytic model for the bispectrum of galaxies in redshift space
Smith, Robert E.; Sheth, Ravi K.; Scoccimarro, Roman
2008-07-15
We develop an analytic theory for the redshift space bispectrum of dark matter, haloes, and galaxies. This is done within the context of the halo model of structure formation, as this allows for the self-consistent inclusion of linear and nonlinear redshift-space distortions and also for the nonlinearity of the halo bias. The model is applicable over a wide range of scales: on the largest scales the predictions reduce to those of the standard perturbation theory (PT); on smaller scales they are determined primarily by the nonlinear virial velocities of galaxies within haloes, and this gives rise to the U-shaped anisotropy in the reduced bispectrum--a finger print of the Finger-Of-God distortions. We then confront the predictions with measurements of the redshift-space bispectrum of dark matter from an ensemble of numerical simulations. On very large scales, k=0.05h Mpc{sup -1}, we find reasonably good agreement between our halo model, PT and the data, to within the errors. On smaller scales, k=0.1h Mpc{sup -1}, the measured bispectra differ from the PT at the level of {approx}10%-20%, especially for colinear triangle configurations. The halo model predictions improve over PT, but are accurate to no better than 10%. On smaller scales k=0.5-1.0h Mpc{sup -1}, our model provides a significant improvement over PT, which breaks down. This implies that studies which use the lowest order PT to extract galaxy bias information are not robust on scales k > or approx. 0.1h Mpc{sup -1}. The analytic and simulation results also indicate that there is no observable scale for which the configuration dependence of the reduced bispectrum is constant--hierarchical models for the higher-order correlation functions in redshift space are unlikely to be useful. It is hoped that our model will facilitate extraction of information from large-scale structure surveys of the Universe, because different galaxy populations are naturally included into our description.
An analytically linearized helicopter model with improved modeling accuracy
NASA Technical Reports Server (NTRS)
Jensen, Patrick T.; Curtiss, H. C., Jr.; Mckillip, Robert M., Jr.
1991-01-01
An analytically linearized model for helicopter flight response including rotor blade dynamics and dynamic inflow, that was recently developed, was studied with the objective of increasing the understanding, the ease of use, and the accuracy of the model. The mathematical model is described along with a description of the UH-60A Black Hawk helicopter and flight test used to validate the model. To aid in utilization of the model for sensitivity analysis, a new, faster, and more efficient implementation of the model was developed. It is shown that several errors in the mathematical modeling of the system caused a reduction in accuracy. These errors in rotor force resolution, trim force and moment calculation, and rotor inertia terms were corrected along with improvements to the programming style and documentation. Use of a trim input file to drive the model is examined. Trim file errors in blade twist, control input phase angle, coning and lag angles, main and tail rotor pitch, and uniform induced velocity, were corrected. Finally, through direct comparison of the original and corrected model responses to flight test data, the effect of the corrections on overall model output is shown.
NASA Astrophysics Data System (ADS)
Mead, A. J.; Peacock, J. A.; Heymans, C.; Joudaki, S.; Heavens, A. F.
2015-12-01
We present an optimized variant of the halo model, designed to produce accurate matter power spectra well into the non-linear regime for a wide range of cosmological models. To do this, we introduce physically motivated free parameters into the halo-model formalism and fit these to data from high-resolution N-body simulations. For a variety of Λ cold dark matter (ΛCDM) and wCDM models, the halo-model power is accurate to ≃ 5 per cent for k ≤ 10h Mpc-1 and z ≤ 2. An advantage of our new halo model is that it can be adapted to account for the effects of baryonic feedback on the power spectrum. We demonstrate this by fitting the halo model to power spectra from the OWLS (OverWhelmingly Large Simulations) hydrodynamical simulation suite via parameters that govern halo internal structure. We are able to fit all feedback models investigated at the 5 per cent level using only two free parameters, and we place limits on the range of these halo parameters for feedback models investigated by the OWLS simulations. Accurate predictions to high k are vital for weak-lensing surveys, and these halo parameters could be considered nuisance parameters to marginalize over in future analyses to mitigate uncertainty regarding the details of feedback. Finally, we investigate how lensing observables predicted by our model compare to those from simulations and from HALOFIT for a range of k-cuts and feedback models and quantify the angular scales at which these effects become important. Code to calculate power spectra from the model presented in this paper can be found at https://github.com/alexander-mead/hmcode.
Fast Analytical Methods for Macroscopic Electrostatic Models in Biomolecular Simulations*
Xu, Zhenli; Cai, Wei
2013-01-01
We review recent developments of fast analytical methods for macroscopic electrostatic calculations in biological applications, including the Poisson–Boltzmann (PB) and the generalized Born models for electrostatic solvation energy. The focus is on analytical approaches for hybrid solvation models, especially the image charge method for a spherical cavity, and also the generalized Born theory as an approximation to the PB model. This review places much emphasis on the mathematical details behind these methods. PMID:23745011
Analytic Model For Estimation Of Cold Bulk Metal Forming Simulations
Skunca, Marko; Keran, Zdenka; Math, Miljenko
2007-05-17
Numerical simulation of bulk metal forming plays an important role in predicting a key parameters in cold forging. Comparison of numerical and experimental data is of great importance, but there is always a need of more universal analytical tools. Therefore, many papers besides experiment and simulation of a particular bulk metal forming technology, include an analytic model. In this paper an analytical model for evaluation of commercially available simulation program packages is proposed. Based on elementary theory of plasticity, being only geometry dependent, model represents a good analytical reference to estimate given modeling preferences like; element types, solver, remeshing influence and many others. Obtained, geometry dependent, stress fields compared with numerical data give a clear picture of numerical possibilities and limitations of particular modeling program package.
Feedbacks Between Numerical and Analytical Models in Hydrogeology
NASA Astrophysics Data System (ADS)
Zlotnik, V. A.; Cardenas, M. B.; Toundykov, D.; Cohn, S.
2012-12-01
Hydrogeology is a relatively young discipline which combines elements of Earth science and engineering. Mature fundamental disciplines (e.g., physics, chemistry, fluid mechanics) have centuries-long history of mathematical modeling even prior to discovery of Darcy's law. Thus, in hydrogeology, relatively few classic analytical models (such those by Theis, Polubarinova-Kochina, Philip, Toth, Henry, Dagan, Neuman) were developed by the early 1970's. The advent of computers and practical demands refocused mathematical models towards numerical techniques. With more diverse but less mathematically-oriented training, most hydrogeologists shifted from analytical methods to use of standardized computational software. Spatial variability in internal properties and external boundary conditions and geometry, and the added complexity of chemical and biological processes will remain major challenges for analytical modeling. Possibly, analytical techniques will play a subordinate role to numerical approaches in many applications. On the other hand, the rise of analytical element modeling of groundwater flow is a strong alternative to numerical models when data demand and computational efficiency is considered. The hallmark of analytical models - transparency and accuracy - will remain indispensable for scientific exploration of complex phenomena and for benchmarking numerical models. Therefore, there will always be feedbacks and complementarities between numerical and analytical techniques, as well as a certain ideological schism among various views to modeling. We illustrate the idea of feedbacks by reviewing evolution of Joszef Toth's analytical model of gravity driven flow systems. Toth's (1963) approach was to reduce the flow domain to a rectangle which allowed for closed-form solution of the governing equations. Succeeding numerical finite-element models by Freeze and Witherspoon (1966-1968) explored the effects of geometry and heterogeneity on regional groundwater flow
Towards an analytical model of water: The octupolar model
NASA Astrophysics Data System (ADS)
Blum, L.; Vericat, F.; Bratko, D.
1995-01-01
A simple potential for the water intermolecular potential, consisting of hard spheres with a point dipole and a potential well with the symmetry of a tetrahedral octupole was proposed by Bratko, Blum, and Luzar some time ago. This structural model was formally solved by Blum, Cummings, and Bratko, and explicit solutions have been recently obtained. We show here that this very simple model agrees surprisingly well with the experimental pair correlation functions of Soper and Phillips. The agreement with the gOH(r) and gHH(r) functions is quite good. For the gOO(r) the agreement is not as good, but this has to do with the hard core nature of the potential. Analytical solutions for soft spherical cores exist.
Computational and analytical modeling of eye refractive surgery
NASA Astrophysics Data System (ADS)
Cabrera, Delia
As the number of corneal refractive procedures increases annually, concerns about their long-term stability and predictability have become the center of attention in the ophthalmic community. This thesis focuses on developing quantitative biomechanical models of the cornea that will overcome shortcomings of previous models and incorporate new observations of corneal elastic properties. Our intent is to provide a more accurate model of the corneal structure to guide current and future developments. The second chapter shows that neural networks could rapidly prototype practical solutions to obtain a better estimate of the average corneal power using the contrast and image size parameters provided by the topographic systems. After establishing improved measurements of the corneal shape the thesis focuses on the development of various corneal models. The analytical model proposed shows that geometric optics, corneal structural properties and surgical nomograms could be used to gain a better understanding of corneal response to surgical interventions. The predictions of this model are closer to the values provided by the published nomograms and clinical data than that obtained by the traditional geometric model. Three surgical procedures (Ultrafast Laser-Automated Lamellar Keratomileusis, Corneal Transplant and Intrastromal Refractive Keratectomy) were simulated using the finite element method. A new formulation was developed that simulates the changes on corneal curvature after refractive surgery when the stiffness inhomogeneities across the corneal thickness are considered. It has been shown that the predictability of the surgical outcome is improved when the stiffness inhomogeneities and nonlinearities of the deformations are included in the finite element simulations. Moreover, a finite element formulation has been developed first time to characterize the intrastromal refractive keratectomy procedure. An inhomogeneous (small displacements) model was identified as an
Combined experimental/analytical modeling using component mode synthesis
Martinez, D.R.; Carne, T.G.; Gregory, D.L.; Miller, A.K.
1984-01-01
This study evaluates the accuracy of computed model frequencies and mode shapes obtained from a combined experimental/analytical model for a simple beam structure. The structure was divided into two subsystems, and one subsystem was tested to obtain its free-free modes. Using a Component Mode Synthesis (CMS) technique, the experimental model data base for one subsystem was directly coupled with a finite element model of the other subsystem to create an experimental/analytical model of the total structure. Both the translational and rotational elements of the residual flexibilities and mode shapes at the interface of the experimental subsystem were measured and used in the coupling. The modal frequencies and mode shapes obtained for the combined experimental/analytical model are compared to those for a reference finite element model of the entire structure. The sensitivity of the CMS model predictions to errors in the modal parameters and residual flexibilities, which are required to define a subsystem, is also examined.
Analytical modelling of Thirty Meter Telescope optics polarization
NASA Astrophysics Data System (ADS)
Anche, Ramya M.; Anupama, G. C.; Reddy, Krishna; Sen, Asoke; Sankarasubramanian, K.; Ramaprakash, A. N.; Sengupta, Sujan; Skidmore, Warren; Atwood, Jenny; Tirupathi, Sivarani; Pandey, Shashi Bhushan
2015-06-01
The polarization introduced due to Thirty Meter Telescope (TMT) optics is calculated using an analytical model. Mueller matrices are also generated for each optical element using Zemax, based on which the instrumental polarization due to the entire system at the focal plane is estimated and compared with the analytical model. This study is significant in the estimation of the telescope sensitivity and also has great implications for future instruments.
Analytical and numerical modeling of non-collinear shear wave mixing at an imperfect interface.
Zhang, Ziyin; Nagy, Peter B; Hassan, Waled
2016-02-01
Non-collinear shear wave mixing at an imperfect interface between two solids can be exploited for nonlinear ultrasonic assessment of bond quality. In this study we developed two analytical models for nonlinear imperfect interfaces. The first model uses a finite nonlinear interfacial stiffness representation of an imperfect interface of vanishing thickness, while the second model relies on a thin nonlinear interphase layer to represent an imperfect interface region. The second model is actually a derivative of the first model obtained by calculating the equivalent interfacial stiffness of a thin isotropic nonlinear interphase layer in the quasi-static approximation. The predictions of both analytical models were numerically verified by comparison to COMSOL finite element simulations. These models can accurately predict the additional nonlinearity caused by interface imperfections based on the strength of the reflected and transmitted mixed longitudinal waves produced by them under non-collinear shear wave interrogation. PMID:26482394
Analytical and numerical modeling of non-collinear shear wave mixing at an imperfect interface
NASA Astrophysics Data System (ADS)
Zhang, Ziyin; Nagy, Peter B.; Hassan, Waled
2016-02-01
Non-collinear shear wave mixing at an imperfect interface between two solids can be exploited for nonlinear ultrasonic assessment of bond quality. In this study we developed two analytical models for nonlinear imperfect interfaces. The first model uses a finite nonlinear interfacial stiffness representation of an imperfect interface of vanishing thickness, while the second model relies on a thin nonlinear interphase layer to represent an imperfect interface region. The second model is actually a derivative of the first model obtained by calculating the equivalent interfacial stiffness of a thin isotropic nonlinear interphase layer in the quasi-static approximation. The predictions of both analytical models were numerically verified by comparison to COMSOL finite element simulations. These models can accurately predict the excess nonlinearity caused by interface imperfections based on the strength of the reflected and transmitted mixed longitudinal waves produced by them under non-collinear shear wave interrogation.
Analytical model for orbital debris environmental management
NASA Technical Reports Server (NTRS)
Talent, David L.
1990-01-01
A differential equation, also referred to as the PIB (particle-in-a-box) model, expressing the time rate of change of the number of objects in orbit, is developed, and its applicability is illustrated. The model can be used as a tool for the assessment of LEO environment stability, and as a starting point for the development of numerical evolutionary models. Within the context of the model, evolutionary scenarios are examined, and found to be sensitive to the growth rate. It is determined that the present environment is slightly unstable to catastrophic growth, and that the number of particles on orbit will continue to increase until approximately 2250-2350 AD, with a maximum of 2,000,000. The model is expandable to the more realistic (complex) case of multiple species in a multiple-tier system.
Experimental and analytical generic space station dynamic models
NASA Technical Reports Server (NTRS)
Belvin, W. K.; Edighoffer, H. H.
1986-01-01
A dynamic model used for verification of analytical and experimental methods is documented. The model consists of five substructures to simulate the multibody, low frequency nature of large space structures. Design considerations which led to a fundamental vibration frequency of less than one Hz are described. Finite element analysis used to predict the vibration modes and frequencies of the experimental model is presented. In addition, modeling of cable suspension effects using prestressed vibration analysis is described. Details of the expermental and analytical models are included to permit replication of the study. Results of the modal vibration tests and analysis are presented in a separate document.
Combined experimental/analytical modeling using component mode synthesis
Martinez, D.R.; Carne, T.G.; Gregory, D.L.; Miller, A.K.
1984-04-01
This study evaluates the accuracy of computed modal frequencies and mode shapes obtained from a combined experimental/analytical model for a simple beam structure. The structure was divided into two subsystems and one subsystem was tested to obtain its free-free modes. Using a Component Mode Synthesis (CMS) technique, the experimental modal data base for one subsystem was directly coupled with a finite element model of the other subsystem to create an experimental/analytical model of the total structure. Both the translational and rotational elements of the residual flexibilities and mode shapes at the interface of the experimental subsystem were measured and used in the coupling. The modal frequencies and mode shapes obtained for the combined experimental/analytical model are compared to those for a reference finite element model of the entire structure. The sensitivity of the CMS model predictions to errors in the modal parameters and residual flexibilities, which are required to define a subsystem, is also examined.
Johnson, Timothy C.; Wellman, Dawn M.
2015-06-26
Electrical resistivity tomography (ERT) has been widely used in environmental applications to study processes associated with subsurface contaminants and contaminant remediation. Anthropogenic alterations in subsurface electrical conductivity associated with contamination often originate from highly industrialized areas with significant amounts of buried metallic infrastructure. The deleterious influence of such infrastructure on imaging results generally limits the utility of ERT where it might otherwise prove useful for subsurface investigation and monitoring. In this manuscript we present a method of accurately modeling the effects of buried conductive infrastructure within the forward modeling algorithm, thereby removing them from the inversion results. The method is implemented in parallel using immersed interface boundary conditions, whereby the global solution is reconstructed from a series of well-conditioned partial solutions. Forward modeling accuracy is demonstrated by comparison with analytic solutions. Synthetic imaging examples are used to investigate imaging capabilities within a subsurface containing electrically conductive buried tanks, transfer piping, and well casing, using both well casings and vertical electrode arrays as current sources and potential measurement electrodes. Results show that, although accurate infrastructure modeling removes the dominating influence of buried metallic features, the presence of metallic infrastructure degrades imaging resolution compared to standard ERT imaging. However, accurate imaging results may be obtained if electrodes are appropriately located.
An Analytical Model of Tribocharging in Regolith
NASA Astrophysics Data System (ADS)
Carter, D. P.; Hartzell, C. M.
2015-12-01
Nongravitational forces, including electrostatic forces and cohesion, can drive the behavior of regolith in low gravity environments such as the Moon and asteroids. Regolith is the 'skin' of solid planetary bodies: it is the outer coating that is observed by orbiters and the first material contacted by landers. Triboelectric charging, the phenomenon by which electrical charge accumulates during the collision or rubbing of two surfaces, has been found to occur in initially electrically neutral granular mixtures. Although charge transfer is often attributed to chemical differences between the different materials, charge separation has also been found to occur in mixtures containing grains of a single material, but with a variety of grain sizes. In such cases, the charge always separates according to grain size; typically the smaller grains acquire a more negative charge than the larger grains. Triboelectric charging may occur in a variety of planetary phenomena (including mass wasting and dust storms) as well as during spacecraft-surface interactions (including sample collection and wheel motion). Interactions between charged grains or with the solar wind plasma could produce regolith motion. However, a validated, predictive model of triboelectric charging between dielectric grains has not yet been developed. A model for such size-dependent charge separation will be presented, demonstrating how random collisions between initially electrically neutral grains lead to net migration of electrons toward the smaller grains. The model is applicable to a wide range of single-material granular mixtures, including those with unusual or wildly varying size distributions, and suggests a possible mechanism for the reversal of the usual size-dependent charge polarity described above. This is a significant improvement over existing charge exchange models, which are restricted to two discrete grains sizes and provide severely limited estimates for charge magnitude. We will also
Analytical modelling of soil effects on electromagnetic induction sensor for humanitarian demining
NASA Astrophysics Data System (ADS)
Vasić, D.; Ambruš, D.; Bilas, V.
2013-06-01
Accurate compensation of the soil effect is essential for a new generation of sensitive classification-based electromagnetic induction landmine detectors. We present an analytical model for evaluation of the soil effect suitable for straightforward numerical implementation. The modelled soil consists of arbitrary number of conductive and magnetic layers. The solution region is truncated leading to the solution in form of a series rather than infinite integrals. Frequency-dependent permeability is inherent to the model, and time domain analysis can be made using DFT. In order to illustrate the model usage, we evaluate performances of three metal detector designs.
Hybrid experimental/analytical models of structural dynamics - Creation and use for predictions
NASA Technical Reports Server (NTRS)
Balmes, Etienne
1993-01-01
An original complete methodology for the construction of predictive models of damped structural vibrations is introduced. A consistent definition of normal and complex modes is given which leads to an original method to accurately identify non-proportionally damped normal mode models. A new method to create predictive hybrid experimental/analytical models of damped structures is introduced, and the ability of hybrid models to predict the response to system configuration changes is discussed. Finally a critical review of the overall methodology is made by application to the case of the MIT/SERC interferometer testbed.
NASA Astrophysics Data System (ADS)
Lachaume, Regis; Rabus, Markus; Jordan, Andres
2015-08-01
In stellar interferometry, the assumption that the observables can be seen as Gaussian, independent variables is the norm. In particular, neither the optical interferometry FITS (OIFITS) format nor the most popular fitting software in the field, LITpro, offer means to specify a covariance matrix or non-Gaussian uncertainties. Interferometric observables are correlated by construct, though. Also, the calibration by an instrumental transfer function ensures that the resulting observables are not Gaussian, even if uncalibrated ones happened to be so.While analytic frameworks have been published in the past, they are cumbersome and there is no generic implementation available. We propose here a relatively simple way of dealing with correlated errors without the need to extend the OIFITS specification or making some Gaussian assumptions. By repeatedly picking at random which interferograms, which calibrator stars, and which are the errors on their diameters, and performing the data processing on the bootstrapped data, we derive a sampling of p(O), the multivariate probability density function (PDF) of the observables O. The results can be stored in a normal OIFITS file. Then, given a model m with parameters P predicting observables O = m(P), we can estimate the PDF of the model parameters f(P) = p(m(P)) by using a density estimation of the observables' PDF p.With observations repeated over different baselines, on nights several days apart, and with a significant set of calibrators systematic errors are de facto taken into account. We apply the technique to a precise and accurate assessment of stellar diameters obtained at the Very Large Telescope Interferometer with PIONIER.
Analytical and numerical modeling for flexible pipes
NASA Astrophysics Data System (ADS)
Wang, Wei; Chen, Geng
2011-12-01
The unbonded flexible pipe of eight layers, in which all the layers except the carcass layer are assumed to have isotropic properties, has been analyzed. Specifically, the carcass layer shows the orthotropic characteristics. The effective elastic moduli of the carcass layer have been developed in terms of the influence of deformation to stiffness. With consideration of the effective elastic moduli, the structure can be properly analyzed. Also the relative movements of tendons and relative displacements of wires in helical armour layer have been investigated. A three-dimensional nonlinear finite element model has been presented to predict the response of flexible pipes under axial force and torque. Further, the friction and contact of interlayer have been considered. Comparison between the finite element model and experimental results obtained in literature has been given and discussed, which might provide practical and technical support for the application of unbonded flexible pipes.
Analytical model for screening potential CO2 repositories
Okwen, R.T.; Stewart, M.T.; Cunningham, J.A.
2011-01-01
Assessing potential repositories for geologic sequestration of carbon dioxide using numerical models can be complicated, costly, and time-consuming, especially when faced with the challenge of selecting a repository from a multitude of potential repositories. This paper presents a set of simple analytical equations (model), based on the work of previous researchers, that could be used to evaluate the suitability of candidate repositories for subsurface sequestration of carbon dioxide. We considered the injection of carbon dioxide at a constant rate into a confined saline aquifer via a fully perforated vertical injection well. The validity of the analytical model was assessed via comparison with the TOUGH2 numerical model. The metrics used in comparing the two models include (1) spatial variations in formation pressure and (2) vertically integrated brine saturation profile. The analytical model and TOUGH2 show excellent agreement in their results when similar input conditions and assumptions are applied in both. The analytical model neglects capillary pressure and the pressure dependence of fluid properties. However, simulations in TOUGH2 indicate that little error is introduced by these simplifications. Sensitivity studies indicate that the agreement between the analytical model and TOUGH2 depends strongly on (1) the residual brine saturation, (2) the difference in density between carbon dioxide and resident brine (buoyancy), and (3) the relationship between relative permeability and brine saturation. The results achieved suggest that the analytical model is valid when the relationship between relative permeability and brine saturation is linear or quasi-linear and when the irreducible saturation of brine is zero or very small. ?? 2011 Springer Science+Business Media B.V.
Analytical model for fast-shock ignition
Ghasemi, S. A. Farahbod, A. H.; Sobhanian, S.
2014-07-15
A model and its improvements are introduced for a recently proposed approach to inertial confinement fusion, called fast-shock ignition (FSI). The analysis is based upon the gain models of fast ignition, shock ignition and considerations for the fast electrons penetration into the pre-compressed fuel to examine the formation of an effective central hot spot. Calculations of fast electrons penetration into the dense fuel show that if the initial electron kinetic energy is of the order ∼4.5 MeV, the electrons effectively reach the central part of the fuel. To evaluate more realistically the performance of FSI approach, we have used a quasi-two temperature electron energy distribution function of Strozzi (2012) and fast ignitor energy formula of Bellei (2013) that are consistent with 3D PIC simulations for different values of fast ignitor laser wavelength and coupling efficiency. The general advantages of fast-shock ignition in comparison with the shock ignition can be estimated to be better than 1.3 and it is seen that the best results can be obtained for the fuel mass around 1.5 mg, fast ignitor laser wavelength ∼0.3 micron and the shock ignitor energy weight factor about 0.25.
NASA Astrophysics Data System (ADS)
Shi, De-Heng; Liu, Yu-Fang; Sun, Jin-Feng; Zhu, Zun-Lue; Yang, Xiang-Dong
2006-12-01
The reasonable dissociation limit of the second excited singlet state B1Π of 7LiH molecule is obtained. The accurate dissociation energy and equilibrium geometry of the B1Π state are calculated using a symmetry-adapted-cluster configuration-interaction method in full active space. The whole potential energy curve for the B1Π state is obtained over the internuclear distance ranging from about 0.10 nm to 0.54 nm, and has a least-square fit to the analytic Murrell-Sorbie function form. The vertical excitation energy is calculated from the ground state to the B1Π state and compared with previous theoretical results. The equilibrium internuclear distance obtained by geometry optimization is found to be quite different from that obtained by single-point energy scanning under the same calculation condition. Based on the analytic potential energy function, the harmonic frequency value of the B1Π state is estimated. A comparison of the theoretical calculations of dissociation energies, equilibrium interatomic distances and the analytic potential energy function with those obtained by previous theoretical results clearly shows that the present work is more comprehensive and in better agreement with experiments than previous theories, thus it is an improvement on previous theories.
A new analytic model for fracture-dominated reservoirs
Karasaki, K.; Long, J.C.S.; Witherspoon, P.A.
1988-03-01
A new analytic model for analyzing well test data from fracture-dominated reservoirs is presented. It is a concentric composite model with a finite-radius well located in the center. In the inner region the flow is assumed to be linear, and in the outer region the flow is assumed to be radial. Solutions are obtained analytically and type curves for ranges of dimensionless parameters are presented. The model can be used to find the extent and the flow parameters of the fractures near the well and the average values for the entire system provided that wellbore-storage effects do not mask the early-time data.
Analytical modelling of no-vent fill process
NASA Technical Reports Server (NTRS)
Vaughan, David A.; Schmidt, George R.
1990-01-01
An analytical model called FILL is presented which represents the first step in attaining the capability for no-vent fill of cryogens in space. The model's analytical structure is described, including the equations used to calculate transient thermodynamic behavior in different regions of the tank. The code predictions are compared with data from recent no-vent fill ground tests using Freon-114. The results are used to validate the FILL model to evaluate the viability of universal submerged jet theory in predicting system-level condensation effects.
MONA: An accurate two-phase well flow model based on phase slippage
Asheim, H.
1984-10-01
In two phase flow, holdup and pressure loss are related to interfacial slippage. A model based on the slippage concept has been developed and tested using production well data from Forties, the Ekofisk area, and flowline data from Prudhoe Bay. The model developed turned out considerably more accurate than the standard models used for comparison.
ANALYTIC ELEMENT MODELING OF COASTAL AQUIFERS - PROJECT SUMMARY
Four topics were studied concerning the modeling of groundwater flow in coastal aquifers with analytic elements: (1) practical experience was obtained by constructing a groundwater model
of the shallow aquifers below the Delmarva Peninsula USA using the commercial program MVAE...
Analytical development and optimization of a graphene–solution interface capacitance model
Karimi, Hediyeh; Rahmani, Rasoul; Mashayekhi, Reza; Ranjbari, Leyla; Shirdel, Amir H; Haghighian, Niloofar; Movahedi, Parisa; Hadiyan, Moein
2014-01-01
Summary Graphene, which as a new carbon material shows great potential for a range of applications because of its exceptional electronic and mechanical properties, becomes a matter of attention in these years. The use of graphene in nanoscale devices plays an important role in achieving more accurate and faster devices. Although there are lots of experimental studies in this area, there is a lack of analytical models. Quantum capacitance as one of the important properties of field effect transistors (FETs) is in our focus. The quantum capacitance of electrolyte-gated transistors (EGFETs) along with a relevant equivalent circuit is suggested in terms of Fermi velocity, carrier density, and fundamental physical quantities. The analytical model is compared with the experimental data and the mean absolute percentage error (MAPE) is calculated to be 11.82. In order to decrease the error, a new function of E composed of α and β parameters is suggested. In another attempt, the ant colony optimization (ACO) algorithm is implemented for optimization and development of an analytical model to obtain a more accurate capacitance model. To further confirm this viewpoint, based on the given results, the accuracy of the optimized model is more than 97% which is in an acceptable range of accuracy. PMID:24991496
Analytical development and optimization of a graphene-solution interface capacitance model.
Karimi, Hediyeh; Rahmani, Rasoul; Mashayekhi, Reza; Ranjbari, Leyla; Shirdel, Amir H; Haghighian, Niloofar; Movahedi, Parisa; Hadiyan, Moein; Ismail, Razali
2014-01-01
Graphene, which as a new carbon material shows great potential for a range of applications because of its exceptional electronic and mechanical properties, becomes a matter of attention in these years. The use of graphene in nanoscale devices plays an important role in achieving more accurate and faster devices. Although there are lots of experimental studies in this area, there is a lack of analytical models. Quantum capacitance as one of the important properties of field effect transistors (FETs) is in our focus. The quantum capacitance of electrolyte-gated transistors (EGFETs) along with a relevant equivalent circuit is suggested in terms of Fermi velocity, carrier density, and fundamental physical quantities. The analytical model is compared with the experimental data and the mean absolute percentage error (MAPE) is calculated to be 11.82. In order to decrease the error, a new function of E composed of α and β parameters is suggested. In another attempt, the ant colony optimization (ACO) algorithm is implemented for optimization and development of an analytical model to obtain a more accurate capacitance model. To further confirm this viewpoint, based on the given results, the accuracy of the optimized model is more than 97% which is in an acceptable range of accuracy. PMID:24991496
Experimental verification of the SP-100 TEM pump analytical models
Salamah, S.A.; Miller, D.D.; Sinha, U.N.; Narkiewicz, R.S. )
1993-01-15
Validation of the TEM pump analytical model is conducted via experimental verification of the model prediction. Two key tests that have provided essential information toward this objective are the Magnetic Bench Test (MBT) and the ElectroMagnetic Integration Test (EMIT). The tests are briefly described and experimental results are compared with predictions of simulation models that form part of overall TEM pump performance model.
Creation of Anatomically Accurate Computer-Aided Design (CAD) Solid Models from Medical Images
NASA Technical Reports Server (NTRS)
Stewart, John E.; Graham, R. Scott; Samareh, Jamshid A.; Oberlander, Eric J.; Broaddus, William C.
1999-01-01
Most surgical instrumentation and implants used in the world today are designed with sophisticated Computer-Aided Design (CAD)/Computer-Aided Manufacturing (CAM) software. This software automates the mechanical development of a product from its conceptual design through manufacturing. CAD software also provides a means of manipulating solid models prior to Finite Element Modeling (FEM). Few surgical products are designed in conjunction with accurate CAD models of human anatomy because of the difficulty with which these models are created. We have developed a novel technique that creates anatomically accurate, patient specific CAD solids from medical images in a matter of minutes.
NASA Astrophysics Data System (ADS)
Grasso, Robert J.; Russo, Leonard P.; Barrett, John L.; Odhner, Jefferson E.; Egbert, Paul I.
2007-09-01
BAE Systems presents the results of a program to model the performance of Raman LIDAR systems for the remote detection of atmospheric gases, air polluting hydrocarbons, chemical and biological weapons, and other molecular species of interest. Our model, which integrates remote Raman spectroscopy, 2D and 3D LADAR, and USAF atmospheric propagation codes permits accurate determination of the performance of a Raman LIDAR system. The very high predictive performance accuracy of our model is due to the very accurate calculation of the differential scattering cross section for the specie of interest at user selected wavelengths. We show excellent correlation of our calculated cross section data, used in our model, with experimental data obtained from both laboratory measurements and the published literature. In addition, the use of standard USAF atmospheric models provides very accurate determination of the atmospheric extinction at both the excitation and Raman shifted wavelengths.
Semi-analytical modelling of piezoelectric excitation of guided waves
NASA Astrophysics Data System (ADS)
Kalkowski, Michał K.; Rustighi, Emiliano; Waters, Timothy P.
2015-03-01
Piezoelectric elements are a key component of modern non-destructive testing (NDT) and structural health monitoring (SHM) systems and play a significant role in many other areas involving dynamic interaction with the structure such as energy harvesting, active control, power ultrasonics or removal of surface accretions using structural waves. In this paper we present a wave-based technique for modelling waveguides equipped with piezoelectric actuators in which there is no need for common simplifications regarding their dynamic behaviour or mutual interaction with the structure. The proposed approach is based on the semi-analytical finite element (SAFE) method. We developed a new piezoelectric semi-analytical element and employed the analytical wave approach to model the distributed electric excitation and scattering of the waves at discontinuities. The model is successfully validated against an experiment on a beam-like waveguide with emulated anechoic terminations.
Combined experimental/analytical modeling of shell/payload structures
Martinez, D.R.; Miller, A.K.; Carne, T.G.
1985-12-01
This study evaluates the accuracy of computed modal frequencies obtained from a combined experimental/analytical model of a shell/payload structure. A component mode synthesis technique was used which incorporated free modes and residual effects. The total structure is physically divided into the two subsystems which are connected through stiff joints. The payload was tested to obtain its free-free modes, while a finite element model of the shell was analyzed to obtain its modal description. Both the translational and rotational components of the experimental mode shapes at the payload interface were used in the coupling. Sensitivity studies were also performed to determine the effect of neglecting the residual terms of the payload. Results from a previous study of a combined experimental/analytical model for a beam structure are also given. The beam structure was used to examine the basic procedures and difficulties in experimentally measuring, and analytically accounting for the rotational and residual quantities.
Analytical modeling of the radial pn junction nanowire solar cells
NASA Astrophysics Data System (ADS)
Ali, Nouran M.; Allam, Nageh K.; Abdel Haleem, Ashraf M.; Rafat, Nadia H.
2014-07-01
In photovoltaic solar cells, radial p-n junctions have been considered a very promising structure to improve the carrier collection efficiency and accordingly the conversion efficiency. In the present study, the semiconductor equations, namely Poisson's and continuity equations for a cylindrical p-n junction solar cell, have been solved analytically. The analytical model is based on Green's function theory to calculate the current density, open circuit voltage, fill factor, and conversion efficiency. The model has been used to simulate p-n and p-i-n silicon radial solar cells. The validity and accuracy of the present simulator were confirmed through a comparison with previously published experimental and numerical reports.
An analytical model of the HINT performance metric
Snell, Q.O.; Gustafson, J.L.
1996-10-01
The HINT benchmark was developed to provide a broad-spectrum metric for computers and to measure performance over the full range of memory sizes and time scales. We have extended our understanding of why HINT performance curves look the way they do and can now predict the curves using an analytical model based on simple hardware specifications as input parameters. Conversely, by fitting the experimental curves with the analytical model, hardware specifications such as memory performance can be inferred to provide insight into the nature of a given computer system.
Torrengo, E; Cigliutti, R; Bosco, G; Carena, A; Curri, V; Poggiolini, P; Nespola, A; Zeolla, D; Forghieri, F
2011-12-12
Link design for optical communication systems requires accurate modeling of nonlinear propagation in fibers. This topic has been widely analyzed in last decades with partial successes in special conditions, but without a comprehensive solution. Since the introduction of coherent detection with electronic signal processing the scenario completely changed because this category of systems shows better performances in links without in-line dispersion management. This change to uncompensated transmission allowed to modify the approach in the study of nonlinear fiber propagation and in recent years a series of promising analytical models have been proposed. In this paper, we present an experimental validation over different fiber types of an analytical model for nonlinear propagation over uncompensated optical transmission links. Considering an ultra-dense WDM system, we transmitted ten 120-Gb/s PM-QPSK signals over a multi-span system probing different fiber types: SSMF, PSCF and NZDSF. A good matching was found in all cases showing the potential of the analytical model for accurate performance estimation that could lead to powerful tools for link design. PMID:22274104
Helicopter derivative identification from analytic models and flight test data
NASA Technical Reports Server (NTRS)
Molusis, J. H.; Briczinski, S.
1974-01-01
Recent results of stability derivative identification from helicopter analytic models and flight test data are presented. Six and nine degree-of-freedom (DOF) linear models are identified from an analytic nonlinear helicopter simulation using a least square technique. The identified models are compared with the convectional partial differentiation method for obtaining derivatives to form the basis for interpretation of derivatives identified from flight data. Six degree-of-freedom models are identified from CH-53A and CH-54B flight data, using an extended Kalman filter modified to process several maneuvers simultaneously. The a priori derivative estimate is obtained by optimal filtering of the data and then using a least square method. The results demonstrate that a six DOF identified model is sufficient to determine the low frequency modes of motion, but a nine DOF rotor/body model is necessary for proper representation of short-term response.
Analytic models of the chemical evolution of galaxies
NASA Technical Reports Server (NTRS)
Clayton, Donald D.
1986-01-01
Techniques are described for constructing analytic models of the chemical evolution of galaxies subject to infall of metal-poor material onto a maturing disk. A class of linear models is discussed which takes the star-formation rate within a defined region to be proportional to the mass of interstellar gas within that region, and the instantaneous recycling approximation is adopted. The solutions are obtained by approximately matching the infall rate to parametrized familiies of functions for which the equations are exactly soluble. The masses, the primary and secondary metallicities, and the gas concentrations of radioactive chronometers can all then be analytically expressed. Surveys of galactic abundances in location and in time can be compared to the parameter spaces of the analytic representations.
Analytical modeling for gamma radiation damage on silicon photodiodes
NASA Astrophysics Data System (ADS)
Jafari, H.; Feghhi, S. A. H.
2016-04-01
Radiation-induced damage in PIN silicon photodiode induces degradation of the photodiode parameters. In this work, by presenting an analytical model, the effect of gamma dose on the dark current in a PIN photodiode array was investigated. Geant4 was used to obtain the damage constant as a result of primary incident particle fluence and NIEL distribution calculations. Experimental measurements as well as numerical simulation of the semiconductor with ATLAS were carried out to verify and parameterize the analytical model calculations. A reasonable agreement has been found between analytical results and experimental data for BPX65 silicon photodiodes irradiated by a Co-60 gamma source at total doses up to 500 krad under different reverse voltages. Moreover, the results showed that the dark current of each photodiode array pixel has considerably increased by gamma dose irradiation.
An analytical model of anisotropic low-field electron mobility in wurtzite indium nitride
NASA Astrophysics Data System (ADS)
Wang, Shulong; Liu, Hongxia; Song, Xin; Guo, Yulong; Yang, Zhaonian
2014-03-01
This paper presents a theoretical analysis of anisotropic transport properties and develops an anisotropic low-field electron analytical mobility model for wurtzite indium nitride (InN). For the different effective masses in the Γ-A and Γ-M directions of the lowest valley, both the transient and steady state transport behaviors of wurtzite InN show different transport characteristics in the two directions. From the relationship between velocity and electric field, the difference is more obvious when the electric field is low in the two directions. To make an accurate description of the anisotropic transport properties under low field, for the first time, we present an analytical model of anisotropic low-field electron mobility in wurtzite InN. The effects of different ionized impurity scattering models on the low-field mobility calculated by Monte Carlo method (Conwell-Weisskopf and Brooks-Herring method) are also considered.
Large-signal numerical and analytical HBT models
Teeter, D.A.; East, J.R.; Mains, R.K.; Haddad, G.I. )
1993-05-01
Several large-signal HBT models are investigated in this paper to determine their usefulness at millimeter-wave frequencies. The most detailed model involves numerically solving moments of the Boltzmann Transport Equation. A description of the numerical model is given along with several simulated results. The numerical model is then used to evaluate two analytical HBT models, the conventional Gummel-Poon model and a modified Ebers-Moll model. It is found that the commonly used Gummel-Poon model exhibits poor agreement with numerical and experimental data at millimeter-wave frequencies due to neglect of transit-time delays. Improved agreement between measured and modeled data result by including transit-time effects in an Ebers-Moll model. This simple model has direct application to millimeter-wave power amplifier and oscillator design. Several measured results are presented to help verify the simple model.
A new analytical model for wind farm power prediction
NASA Astrophysics Data System (ADS)
Niayifar, Amin; Porte-Agel, Fernando
2015-04-01
In this study, a new analytical approach is presented and validated to predict wind farm power production. The new model assumes a Gaussian distribution for the velocity deficit in the wake which has been recently proposed by Bastankhah and Porté-Agel (2014). To estimate the velocity deficit in the wake, this model needs the local wake growth rate parameter which is calculated based on the local turbulence intensity in the wind farm. The interaction of the wakes is modeled by use of the velocity deficit superposition principle. Finally, the power curve is used to estimate the power production from the wind turbines. The wind farm model is compared to large-eddy simulation (LES) data of Horns Rev wind farm for a wide range of wind directions. Reasonable agreement between the proposed analytical model and LES data is obtained. This prediction is substantially better than the one obtained with common wind farm softwares such as WAsP.
Modeling Analyte Transport and Capture in Porous Bead Sensors
Chou, Jie; Lennart, Alexis; Wong, Jorge; Ali, Mehnaaz F.; Floriano, Pierre N.; Christodoulides, Nicolaos; Camp, James; McDevitt, John T.
2013-01-01
Porous agarose microbeads, with high surface to volume ratios and high binding densities, are attracting attention as highly sensitive, affordable sensor elements for a variety of high performance bioassays. While such polymer microspheres have been extensively studied and reported on previously and are now moving into real-world clinical practice, very little work has been completed to date to model the convection, diffusion, and binding kinetics of soluble reagents captured within such fibrous networks. Here, we report the development of a three-dimensional computational model and provide the initial evidence for its agreement with experimental outcomes derived from the capture and detection of representative protein and genetic biomolecules in 290μm porous beads. We compare this model to antibody-mediated capture of C-reactive protein and bovine serum albumin, along with hybridization of oligonucleotide sequences to DNA probes. These results suggest that due to the porous interior of the agarose bead, internal analyte transport is both diffusion- and convection-based, and regardless of the nature of analyte, the bead interiors reveal an interesting trickle of convection-driven internal flow. Based on this model, the internal to external flow rate ratio is found to be in the range of 1:3100 to 1:170 for beads with agarose concentration ranging from 0.5% to 8% for the sensor ensembles here studied. Further, both model and experimental evidence suggest that binding kinetics strongly affect analyte distribution of captured reagents within the beads. These findings reveal that high association constants create a steep moving boundary in which unbound analytes are held back at the periphery of the bead sensor. Low association constants create a more shallow moving boundary in which unbound analytes diffuse further into the bead before binding. These models agree with experimental evidence and thus serve as a new tool set for the study of bio-agent transport processes
An analytic spin chain model with fractional revival
NASA Astrophysics Data System (ADS)
Lemay, Jean-Michel; Vinet, Luc; Zhedanov, Alexei
2016-08-01
New analytic spin chains with fractional revival are introduced. Their nearest-neighbor couplings and local magnetic fields correspond to the recurrence coefficients of para-Racah polynomials which are orthogonal on quadratic bi-lattices. These models generalize the spin chain associated to the dual-Hahn polynomials. Instances where perfect state transfer also occurs are identified.
FACTOR ANALYTIC MODELS OF CLUSTERED MULTIVARIATE DATA WITH INFORMATIVE CENSORING
This paper describes a general class of factor analytic models for the analysis of clustered multivariate data in the presence of informative missingness. We assume that there are distinct sets of cluster-level latent variables related to the primary outcomes and to the censorin...
Analytical Models of Legislative Texts for Muslim Scholars
ERIC Educational Resources Information Center
Alwan, Ammar Abdullah Naseh; Yusoff, Mohd Yakubzulkifli Bin Mohd; Al-Hami, Mohammad Said M.
2011-01-01
The significance of the analytical models in traditional Islamic studies is that they contribute in sharpening the intellectual capacity of the students of Islamic studies. Research literature in Islamic studies has descriptive side predominantly; the information is gathered and compiled and rarely analyzed properly. This weakness is because of…
Analytical properties of the anisotropic cubic Ising model
Hansel, D.; Maillard, J.M.; Oitmaa, J.; Velgakis, M.J.
1987-07-01
The authors combine an exact functional relation, the inversion relation, with conventional high-temperature expansions to explore the analytic properties of the anisotropic Ising model on both the square and simple cubic lattice. In particular, they investigate the nature of the singularities that occur in partially resummed expansions of the partition function and of the susceptibility.
Improving light propagation Monte Carlo simulations with accurate 3D modeling of skin tissue
Paquit, Vincent C; Price, Jeffery R; Meriaudeau, Fabrice; Tobin Jr, Kenneth William
2008-01-01
In this paper, we present a 3D light propagation model to simulate multispectral reflectance images of large skin surface areas. In particular, we aim to simulate more accurately the effects of various physiological properties of the skin in the case of subcutaneous vein imaging compared to existing models. Our method combines a Monte Carlo light propagation model, a realistic three-dimensional model of the skin using parametric surfaces and a vision system for data acquisition. We describe our model in detail, present results from the Monte Carlo modeling and compare our results with those obtained with a well established Monte Carlo model and with real skin reflectance images.
Yi, S. A.; Khudik, V.; Siemon, C.; Shvets, G.
2012-12-21
Self-injection of background electrons in plasma wakefield accelerators in the highly nonlinear bubble regime is analyzed using particle-in-cell and semi-analytic modeling. It is shown that the return current in the bubble sheath layer is crucial for accurate determination of the trapped particle trajectories.
An analytic model of toroidal half-wave oscillations: Implication on plasma density estimates
NASA Astrophysics Data System (ADS)
Bulusu, Jayashree; Sinha, A. K.; Vichare, Geeta
2015-06-01
The developed analytic model for toroidal oscillations under infinitely conducting ionosphere ("Rigid-end") has been extended to "Free-end" case when the conjugate ionospheres are infinitely resistive. The present direct analytic model (DAM) is the only analytic model that provides the field line structures of electric and magnetic field oscillations associated with the "Free-end" toroidal wave for generalized plasma distribution characterized by the power law ρ = ρo(ro/r)m, where m is the density index and r is the geocentric distance to the position of interest on the field line. This is important because different regions in the magnetosphere are characterized by different m. Significant improvement over standard WKB solution and an excellent agreement with the numerical exact solution (NES) affirms validity and advancement of DAM. In addition, we estimate the equatorial ion number density (assuming H+ atom as the only species) using DAM, NES, and standard WKB for Rigid-end as well as Free-end case and illustrate their respective implications in computing ion number density. It is seen that WKB method overestimates the equatorial ion density under Rigid-end condition and underestimates the same under Free-end condition. The density estimates through DAM are far more accurate than those computed through WKB. The earlier analytic estimates of ion number density were restricted to m = 6, whereas DAM can account for generalized m while reproducing the density for m = 6 as envisaged by earlier models.
Terahertz microstructured optical fibers: An analytical field model
NASA Astrophysics Data System (ADS)
Sharma, Dinesh Kumar; Sharma, Anurag; Varshney, R. K.; Pal, B. P.
2014-10-01
Microstructured optical fibers (MOFs) have wavelength scale periodic microstructure running along their length. Their core and two-dimensional microstructured cladding might be based on varied geometries and materials, enabling light guidance due to different propagation mechanisms over an extremely large wavelength range, extending to the terahertz (THz) frequency region. As a result, these fibers have revolutionized the optical fiber technology by means of creating new degrees of freedom in the fiber design, fabrication and applicability. We analytically study the modal properties of terahertz microstructured optical fiber (THz MOF), by using our analytical field model, developed for optical waveguides.
Improvements to the analytical linescan model for SEM metrology
NASA Astrophysics Data System (ADS)
Mack, Chris A.; Bunday, Benjamin D.
2016-03-01
Critical dimension scanning electron microscope (CD-SEM) metrology has long used empirical approaches to determine edge locations. While such solutions are very flexible, physics-based models offer the potential for improved accuracy and precision for specific applications. Here, Monte Carlo simulation is used to generate theoretical linescans from single step and line/space targets in order to build a physics-based analytical model, including the presence of bottom footing and top corner rounding. The resulting analytical linescan model fits the Monte Carlo simulation results for different feature heights, widths, pitches, sidewall angles, bottom footing, and top corner rounding. This model has also been successfully applied to asymetric features such as sidewall spacers encountered in self-aligned double patterning.
On the Development of Parameterized Linear Analytical Longitudinal Airship Models
NASA Technical Reports Server (NTRS)
Kulczycki, Eric A.; Johnson, Joseph R.; Bayard, David S.; Elfes, Alberto; Quadrelli, Marco B.
2008-01-01
In order to explore Titan, a moon of Saturn, airships must be able to traverse the atmosphere autonomously. To achieve this, an accurate model and accurate control of the vehicle must be developed so that it is understood how the airship will react to specific sets of control inputs. This paper explains how longitudinal aircraft stability derivatives can be used with airship parameters to create a linear model of the airship solely by combining geometric and aerodynamic airship data. This method does not require system identification of the vehicle. All of the required data can be derived from computational fluid dynamics and wind tunnel testing. This alternate method of developing dynamic airship models will reduce time and cost. Results are compared to other stable airship dynamic models to validate the methods. Future work will address a lateral airship model using the same methods.
Analytic model of an IR radiation heat pipe
NASA Technical Reports Server (NTRS)
Hoffman, Pamela J.
1990-01-01
An IR radiation heat pipe made from multilayer insulation blankets and proposed to be used aboard spacecraft to transfer waste heat was modeled analytically. A circular cross section pipe 9-in. in diameter, 10-ft long, with a specular reflectivity of 0.94 was found to have an efficiency of 58.6 percent. Several key parameters were varied for the circular model to understand their significance. In addition, square and triangular cross section pipes were investigated.
Analytical modeling of thermoluminescent albedo detectors for neutron dosimetry.
Glickstein, S S
1983-02-01
In order to gain an in-depth understanding of the neutron physics of a 6LiF TLD when used as an albedo neutron dosimeter, an analytical model was developed to simulate the response of a 6LiF chip. The analytical model was used to examine the sensitivity of the albedo TLD response to incident monoenergetic neutrons and to evaluate a multiple chip TLD neutron dosimeter. Contrary to initial experimental studies, which were hampered by statistical uncertainties, the analytical evaluation revealed that a three-energy-group detector could not reliably measure the dose equivalent to personnel exposed to multiple neutron spectra. The analysis clearly illustrates that there may be order of magnitude errors in the measured neutron dose if the dosimeter has not been calibrated for the same flux spectrum to which it is exposed. As a result of this analysis, it was concluded that, for personnel neutron monitoring, a present TLD badge must be calibrated for the neutron spectrum into which the badge is to be introduced. The analytical model used in this study can readily be adopted for evaluating other possible detectors and shield material that might be proposed in the future as suitable for use in neutron dosimetry applications. PMID:6826377
Analytic model for assessing thermal performance of SCUBA divers
NASA Technical Reports Server (NTRS)
Montgomery, L. D.
1975-01-01
To assist design of adequate protective clothing, mathematical model of man's thermoregulatory system has been developed so that body thermal responses under immersed conditions can be predicted accurately. Experimental data encompassed wide range of water temperatures, protective clothing, breathing-gas mixtures, and durations of immersion.
Yang, Li; Tunega, Daniel; Xu, Lai; Govind, Niranjan; Sun, Rui; Taylor, Ramona; Lischka, Hans; De Jong, Wibe A.; Hase, William L.
2013-08-29
In a previous study (J. Phys. Chem. C 2011, 115, 12403) cluster models for the TiO2 rutile (110) surface and MP2 calculations were used to develop an analytic potential energy function for dimethyl methylphosphonate (DMMP) interacting with this surface. In the work presented here, this analytic potential and MP2 cluster models are compared with DFT "slab" calculations for DMMP interacting with the TiO2 (110) surface and with DFT cluster models for the TiO2 (110) surface. The DFT slab calculations were performed with the PW91 and PBE functionals. The analytic potential gives DMMP/ TiO2 (110) potential energy curves in excellent agreement with those obtained from the slab calculations. The cluster models for the TiO2 (110) surface, used for the MP2 calculations, were extended to DFT calculations with the B3LYP, PW91, and PBE functional. These DFT calculations do not give DMMP/TiO2 (110) interaction energies which agree with those from the DFT slab calculations. Analyses of the wave functions for these cluster models show that they do not accurately represent the HOMO and LUMO for the surface, which should be 2p and 3d orbitals, respectively, and the models also do not give an accurate band gap. The MP2 cluster models do not accurately represent the LUMO and that they give accurate DMMP/TiO2 (110) interaction energies is apparently fortuitous, arising from their highly inaccurate band gaps. Accurate cluster models, consisting of 7, 10, and 15 Ti-atoms and which have the correct HOMO and LUMO properties, are proposed. The work presented here illustrates the care that must be taken in "constructing" cluster models which accurately model surfaces.
An analytical model for permeability of isotropic porous media
NASA Astrophysics Data System (ADS)
Yang, Xiaohu; Lu, Tian Jian; Kim, Tongbeum
2014-06-01
We demonstrate that permeability of isotropic porous media e.g., open-cell foams can be analytically presented as a function of two morphological parameters: porosity and pore size. Adopting a cubic unit cell model, an existing tortuosity model from the branching algorithm method is incorporated into a generalized permeability model. The present model shows that dimensionless permeability significantly increases as the porosity of isotropic porous media and unifies the previously reported data in a wide range of porosity (ɛ=0.55-0.98) and pore size (Dp=0.254 mm-5.08 mm).
NASA Astrophysics Data System (ADS)
O'Brien, Edward P.; Morrison, Greg; Brooks, Bernard R.; Thirumalai, D.
2009-03-01
Single molecule Förster resonance energy transfer (FRET) experiments are used to infer the properties of the denatured state ensemble (DSE) of proteins. From the measured average FRET efficiency, ⟨E⟩, the distance distribution P(R ) is inferred by assuming that the DSE can be described as a polymer. The single parameter in the appropriate polymer model (Gaussian chain, wormlike chain, or self-avoiding walk) for P(R ) is determined by equating the calculated and measured ⟨E⟩. In order to assess the accuracy of this "standard procedure," we consider the generalized Rouse model (GRM), whose properties [⟨E⟩ and P(R )] can be analytically computed, and the Molecular Transfer Model for protein L for which accurate simulations can be carried out as a function of guanadinium hydrochloride (GdmCl) concentration. Using the precisely computed ⟨E⟩ for the GRM and protein L, we infer P(R ) using the standard procedure. We find that the mean end-to-end distance can be accurately inferred (less than 10% relative error) using ⟨E⟩ and polymer models for P(R ). However, the value extracted for the radius of gyration (Rg) and the persistence length (lp) are less accurate. For protein L, the errors in the inferred properties increase as the GdmCl concentration increases for all polymer models. The relative error in the inferred Rg and lp, with respect to the exact values, can be as large as 25% at the highest GdmCl concentration. We propose a self-consistency test, requiring measurements of ⟨E⟩ by attaching dyes to different residues in the protein, to assess the validity of describing DSE using the Gaussian model. Application of the self-consistency test to the GRM shows that even for this simple model, which exhibits an order→disorder transition, the Gaussian P(R ) is inadequate. Analysis of experimental data of FRET efficiencies with dyes at several locations for the cold shock protein, and simulations results for protein L, for which accurate FRET
An analytical model of thermal mechanical stress induced by through silicon via
NASA Astrophysics Data System (ADS)
Dong, Gang; Shi, Tao; Zhao, Ying-Bo; Yang, Yin-Tang
2015-05-01
We present an accurate through silicon via (TSV) thermal mechanical stress analytical model which is verified by using finite element method (FEM). The results show only a very small error. By using the proposed analytical model, we also study the impacts of the TSV radius size, the thickness, the material of Cu diffusion barrier, and liner on the stress. It is found that the liner can absorb the stress effectively induced by coefficient of thermal expansion mismatch. The stress decreases with the increase of liner thickness. Benzocyclobutene (BCB) as a liner material is better than SiO2. However, the Cu diffusion barrier has little effect on the stress. The stress with a smaller TSV has a smaller value. Based on the analytical model, we explore and validate the linear superposition principle of stress tensors and demonstrate the accuracy of this method against detailed FEM simulations. The analytic solutions of stress of two TSVs and three TSVs have high precision against the finite element result. Project supported by the National Natural Science Foundation of China (Grant No. 61334003) and the Kunshan Innovation Institute of Xidian University.
Health Informatics for Neonatal Intensive Care Units: An Analytical Modeling Perspective
Mench-Bressan, Nadja; McGregor, Carolyn; Pugh, James Edward
2015-01-01
The effective use of data within intensive care units (ICUs) has great potential to create new cloud-based health analytics solutions for disease prevention or earlier condition onset detection. The Artemis project aims to achieve the above goals in the area of neonatal ICUs (NICU). In this paper, we proposed an analytical model for the Artemis cloud project which will be deployed at McMaster Children’s Hospital in Hamilton. We collect not only physiological data but also the infusion pumps data that are attached to NICU beds. Using the proposed analytical model, we predict the amount of storage, memory, and computation power required for the system. Capacity planning and tradeoff analysis would be more accurate and systematic by applying the proposed analytical model in this paper. Numerical results are obtained using real inputs acquired from McMaster Children’s Hospital and a pilot deployment of the system at The Hospital for Sick Children (SickKids) in Toronto. PMID:27170907
Comprehensive analytical model to characterize randomness in optical waveguides.
Zhou, Junhe; Gallion, Philippe
2016-04-01
In this paper, the coupled mode theory (CMT) is used to derive the corresponding stochastic differential equations (SDEs) for the modal amplitude evolution inside optical waveguides with random refractive index variations. Based on the SDEs, the ordinary differential equations (ODEs) are derived to analyze the statistics of the modal amplitudes, such as the optical power and power variations as well as the power correlation coefficients between the different modal powers. These ODEs can be solved analytically and therefore, it greatly simplifies the analysis. It is demonstrated that the ODEs for the power evolution of the modes are in excellent agreement with the Marcuse' coupled power model. The higher order statistics, such as the power variations and power correlation coefficients, which are not exactly analyzed in the Marcuse' model, are discussed afterwards. Monte-Carlo simulations are performed to demonstrate the validity of the analytical model. PMID:27136981
Analytical Modeling for the Grating Eddy Current Displacement Sensors
NASA Astrophysics Data System (ADS)
Lv, Chunfeng; Tao, Wei; Lei, Huaming; Jiang, Yingying; Zhao, Hui
2015-02-01
As a new type of displacement sensor, grating eddy current displacement sensor (GECDS) combines traditional eddy current sensors and grating structure in one. The GECDS performs a wide range displacement measurement without precision reduction. This paper proposes an analytical modeling approach for the GECDS. The solution model is established in the Cartesian coordinate system, and the solving domain is limited to finite extents by using the truncated region eigenfunction expansion method. Based on the second order vector potential, expressions for the electromagnetic field as well as coil impedance related to the displacement can be expressed in closed-form. Theoretical results are then confirmed by experiments, which prove the suitability and effectiveness of the analytical modeling approach.
Analytical Model for Thermal Elastoplastic Stresses of Functionally Graded Materials
Zhai, P. C.; Chen, G.; Liu, L. S.; Fang, C.; Zhang, Q. J.
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 geometries 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.
Galactic chemical evolution and nucleocosmochronology - Analytic quadratic models
NASA Technical Reports Server (NTRS)
Clayton, D. D.
1985-01-01
Quadratic models of the chemical evolution of the Galaxy for a star formation rate proportional to the square of the gas mass are studied. The search for analytic solutions to the gas mass and star mass for time-dependent rates of gaseous infall onto the disk is examined. The quadratic models are compared to models having linear star formation rates. The mass, metallicity, number of stars, and U-235/U-238 isotopic ratio for the models which are subjected to the same infall rate, the same initial disk mass, and the same final gas fraction are compared. The results of the comparison indicate that: (1) the average dwarf age is greater in the quadratic model, (2) the metallicity grows initially faster in the quadratic model, (3) the quadratic model has a smaller percentage of low-Z dwarfs, and (4) the U-235/U-238 isotopic ratio indicates a younger quadratic model.
Identification of accurate nonlinear rainfall-runoff models with unique parameters
NASA Astrophysics Data System (ADS)
Schoups, G.; Vrugt, J. A.; Fenicia, F.; van de Giesen, N.
2009-04-01
We propose a strategy to identify models with unique parameters that yield accurate streamflow predictions, given a time-series of rainfall inputs. The procedure consists of five general steps. First, an a priori range of model structures is specified based on prior general and site-specific hydrologic knowledge. To this end, we rely on a flexible model code that allows a specification of a wide range of model structures, from simple to complex. Second, using global optimization each model structure is calibrated to a record of rainfall-runoff data, yielding optimal parameter values for each model structure. Third, accuracy of each model structure is determined by estimating model prediction errors using independent validation and statistical theory. Fourth, parameter identifiability of each calibrated model structure is estimated by means of Monte Carlo Markov Chain simulation. Finally, an assessment is made about each model structure in terms of its accuracy of mimicking rainfall-runoff processes (step 3), and the uniqueness of its parameters (step 4). The procedure results in the identification of the most complex and accurate model supported by the data, without causing parameter equifinality. As such, it provides insight into the information content of the data for identifying nonlinear rainfall-runoff models. We illustrate the method using rainfall-runoff data records from several MOPEX basins in the US.
Secondary metallicity in analytic models of chemical evolution of galaxies
NASA Technical Reports Server (NTRS)
Clayton, D. D.; Pantelaki, I.
1986-01-01
Analytic models of the chemical evolution of galactic regions that grow in mass owing to the continuous infall of matter are characterized, emphasizing the solutions for secondary nuclei (defined as those nuclei whose stellar yields are proportional to the abundance of a primary seed nucleus) in the families of models described by Clayton (1984 and 1985). Wide variations in time dependence of both primary and secondary nuclei as well as in the ratio of secondary to primary are displayed by these model families, confirming again the usefulness of these families as interpretive guides if galaxies do in fact evolve with substantial infall. Additionally, analytic solutions are presented for two other possible interesting systems: the evolution of abundances if the primary metallicity in the infall is increasing in time, and the evolution of abundances if the primary yield changes linearly with time owing to continuous changes in the stellar mass function, the opacity, or other astrophysical agents. Finally, test evaluations of the instantaneous recycling approximation on which these analytic models rely are presented.
Human performance modeling for system of systems analytics.
Dixon, Kevin R.; Lawton, Craig R.; Basilico, Justin Derrick; Longsine, Dennis E.; Forsythe, James Chris; Gauthier, John Henry; Le, Hai D.
2008-10-01
A Laboratory-Directed Research and Development project was initiated in 2005 to investigate Human Performance Modeling in a System of Systems analytic environment. SAND2006-6569 and SAND2006-7911 document interim results from this effort; this report documents the final results. The problem is difficult because of the number of humans involved in a System of Systems environment and the generally poorly defined nature of the tasks that each human must perform. A two-pronged strategy was followed: one prong was to develop human models using a probability-based method similar to that first developed for relatively well-understood probability based performance modeling; another prong was to investigate more state-of-art human cognition models. The probability-based modeling resulted in a comprehensive addition of human-modeling capability to the existing SoSAT computer program. The cognitive modeling resulted in an increased understanding of what is necessary to incorporate cognition-based models to a System of Systems analytic environment.
Material Models for Accurate Simulation of Sheet Metal Forming and Springback
NASA Astrophysics Data System (ADS)
Yoshida, Fusahito
2010-06-01
For anisotropic sheet metals, modeling of anisotropy and the Bauschinger effect is discussed in the framework of Yoshida-Uemori kinematic hardening model incorporating with anisotropic yield functions. The performances of the models in predicting yield loci, cyclic stress-strain responses on several types of steel and aluminum sheets are demonstrated by comparing the numerical simulation results with the corresponding experimental observations. From some examples of FE simulation of sheet metal forming and springback, it is concluded that modeling of both the anisotropy and the Bauschinger effect is essential for the accurate numerical simulation.
Progress on Analytical Modeling of Coherent Electron Cooling
Wang, G.; Blaskiewicz, M.; Litvinenko, V.; Webb, S.
2010-05-23
We report recent progresses on analytical studies of Coherent Electron Cooling. The phase space electron beam distribution obtained from the 1D FEL amplifier is applied to an infinite electron plasma model and the electron density evolution inside the kicker is derived. We also investigate the velocity modulation in the modulator and obtain a closed form solution for the current density evolution for infinite homogeneous electron plasma.
Development of modified cable models to simulate accurate neuronal active behaviors
2014-01-01
In large network and single three-dimensional (3-D) neuron simulations, high computing speed dictates using reduced cable models to simulate neuronal firing behaviors. However, these models are unwarranted under active conditions and lack accurate representation of dendritic active conductances that greatly shape neuronal firing. Here, realistic 3-D (R3D) models (which contain full anatomical details of dendrites) of spinal motoneurons were systematically compared with their reduced single unbranched cable (SUC, which reduces the dendrites to a single electrically equivalent cable) counterpart under passive and active conditions. The SUC models matched the R3D model's passive properties but failed to match key active properties, especially active behaviors originating from dendrites. For instance, persistent inward currents (PIC) hysteresis, frequency-current (FI) relationship secondary range slope, firing hysteresis, plateau potential partial deactivation, staircase currents, synaptic current transfer ratio, and regional FI relationships were not accurately reproduced by the SUC models. The dendritic morphology oversimplification and lack of dendritic active conductances spatial segregation in the SUC models caused significant underestimation of those behaviors. Next, SUC models were modified by adding key branching features in an attempt to restore their active behaviors. The addition of primary dendritic branching only partially restored some active behaviors, whereas the addition of secondary dendritic branching restored most behaviors. Importantly, the proposed modified models successfully replicated the active properties without sacrificing model simplicity, making them attractive candidates for running R3D single neuron and network simulations with accurate firing behaviors. The present results indicate that using reduced models to examine PIC behaviors in spinal motoneurons is unwarranted. PMID:25277743
NASA Astrophysics Data System (ADS)
Rumple, Christopher; Krane, Michael; Richter, Joseph; Craven, Brent
2013-11-01
The mammalian nose is a multi-purpose organ that houses a convoluted airway labyrinth responsible for respiratory air conditioning, filtering of environmental contaminants, and chemical sensing. Because of the complexity of the nasal cavity, the anatomy and function of these upper airways remain poorly understood in most mammals. However, recent advances in high-resolution medical imaging, computational modeling, and experimental flow measurement techniques are now permitting the study of respiratory airflow and olfactory transport phenomena in anatomically-accurate reconstructions of the nasal cavity. Here, we focus on efforts to manufacture an anatomically-accurate transparent model for stereoscopic particle image velocimetry (SPIV) measurements. Challenges in the design and manufacture of an index-matched anatomical model are addressed. PIV measurements are presented, which are used to validate concurrent computational fluid dynamics (CFD) simulations of mammalian nasal airflow. Supported by the National Science Foundation.
Analytical models quantify the military benefit of collaborative search
NASA Astrophysics Data System (ADS)
Friedman, Melvin H.; Du Bosq, Todd W.; Flug, Eric A.
2010-04-01
Analytical Model 1 describes how long it takes the first observer to find a target when multiple observers search a field of regard using imagery provided by a single sensor. This model, developed using probability concepts, suggests considerable benefits accrue from collaborative search: when P is near one and with ten observers the mean detection time (in reduced time) is reduced by almost an order of magnitude when compared to that of a single observer. To get the instant of detection in clock time we add the delay time td to the reduced time. Empirical fits for td and are also given in the paper. Model 1 was verified/validated by computer simulation and perception experiments. Here ten observers searched sixty computer generated fields of regard (each one was 60 x 20 degrees) for a single military vehicle. Analytical Model 2 describes how the probability of target acquisition increases with the number of observers. The results of Model 2 suggest that probability of target acquisition increases considerably when multiple observers independently search a field of regard. Model 2 was verified by simulation but not by perception experiment. Models 1 and 2 are pertinent to development of search strategies with multiple observers and are expected to find use in wargaming for evaluating the efficacy of networked imaging sensors.
Optimization of Analytical Potentials for Coarse-Grained Biopolymer Models.
Mereghetti, Paolo; Maccari, Giuseppe; Spampinato, Giulia Lia Beatrice; Tozzini, Valentina
2016-08-25
The increasing trend in the recent literature on coarse grained (CG) models testifies their impact in the study of complex systems. However, the CG model landscape is variegated: even considering a given resolution level, the force fields are very heterogeneous and optimized with very different parametrization procedures. Along the road for standardization of CG models for biopolymers, here we describe a strategy to aid building and optimization of statistics based analytical force fields and its implementation in the software package AsParaGS (Assisted Parameterization platform for coarse Grained modelS). Our method is based on the use and optimization of analytical potentials, optimized by targeting internal variables statistical distributions by means of the combination of different algorithms (i.e., relative entropy driven stochastic exploration of the parameter space and iterative Boltzmann inversion). This allows designing a custom model that endows the force field terms with a physically sound meaning. Furthermore, the level of transferability and accuracy can be tuned through the choice of statistical data set composition. The method-illustrated by means of applications to helical polypeptides-also involves the analysis of two and three variable distributions, and allows handling issues related to the FF term correlations. AsParaGS is interfaced with general-purpose molecular dynamics codes and currently implements the "minimalist" subclass of CG models (i.e., one bead per amino acid, Cα based). Extensions to nucleic acids and different levels of coarse graining are in the course. PMID:27150459
AN ANALYTIC RADIATIVE-CONVECTIVE MODEL FOR PLANETARY ATMOSPHERES
Robinson, Tyler D.; Catling, David C.
2012-09-20
We present an analytic one-dimensional radiative-convective model of the thermal structure of planetary atmospheres. Our model assumes that thermal radiative transfer is gray and can be represented by the two-stream approximation. Model atmospheres are assumed to be in hydrostatic equilibrium, with a power-law scaling between the atmospheric pressure and the gray thermal optical depth. The convective portions of our models are taken to follow adiabats that account for condensation of volatiles through a scaling parameter to the dry adiabat. By combining these assumptions, we produce simple, analytic expressions that allow calculations of the atmospheric-pressure-temperature profile, as well as expressions for the profiles of thermal radiative flux and convective flux. We explore the general behaviors of our model. These investigations encompass (1) worlds where atmospheric attenuation of sunlight is weak, which we show tend to have relatively high radiative-convective boundaries; (2) worlds with some attenuation of sunlight throughout the atmosphere, which we show can produce either shallow or deep radiative-convective boundaries, depending on the strength of sunlight attenuation; and (3) strongly irradiated giant planets (including hot Jupiters), where we explore the conditions under which these worlds acquire detached convective regions in their mid-tropospheres. Finally, we validate our model and demonstrate its utility through comparisons to the average observed thermal structure of Venus, Jupiter, and Titan, and by comparing computed flux profiles to more complex models.
NASA Astrophysics Data System (ADS)
Hurlbatt, A.; O’Connell, D.; Gans, T.
2016-08-01
Analytical and numerical models allow investigation of complicated discharge phenomena and the interplay that makes plasmas such a complex environment. Global models are quick to implement and can have almost negligible computation cost, but provide only bulk or spatially averaged values. Full fluid models take longer to develop, and can take days to solve, but provide accurate spatio-temporal profiles of the whole plasma. The work presented here details a different type of model, analytically similar to fluid models, but computationally closer to a global model, and able to give spatially resolved solutions for the challenging environment of electronegative plasmas. Included are non-isothermal electrons, gas heating, and coupled neutral dynamics. Solutions are reached in seconds to minutes, and spatial profiles are given for densities, fluxes, and temperatures. This allows the semi-analytical model to fill the gap that exists between global and full fluid models, extending the tools available to researchers. The semi-analytical model can perform broad parameter sweeps that are not practical with more computationally expensive models, as well as exposing non-trivial trends that global models cannot capture. Examples are given for a low pressure oxygen CCP. Excellent agreement is shown with a full fluid model, and comparisons are drawn with the corresponding global model.
Cimpoesu, Dorin Stoleriu, Laurentiu; Stancu, Alexandru
2013-12-14
We propose a generalized Stoner-Wohlfarth (SW) type model to describe various experimentally observed angular dependencies of the switching field in non-single-domain magnetic particles. Because the nonuniform magnetic states are generally characterized by complicated spin configurations with no simple analytical description, we maintain the macrospin hypothesis and we phenomenologically include the effects of nonuniformities only in the anisotropy energy, preserving as much as possible the elegance of SW model, the concept of critical curve and its geometric interpretation. We compare the results obtained with our model with full micromagnetic simulations in order to evaluate the performance and limits of our approach.
A new analytical model for wind farm power prediction
NASA Astrophysics Data System (ADS)
Niayifar, Amin; Porté-Agel, Fernando
2015-06-01
In this study, a new analytical approach is presented and validated to predict wind farm power production. The new model is an extension of the recently proposed by Bastankhah and Porté-Agel for a single wake. It assumes a self-similar Gaussian shape of the velocity deficit and satisfies conservation of mass and momentum. To estimate the velocity deficit in the wake, this model needs the local wake growth rate parameter which is calculated based on the local turbulence intensity in the wind farm. The interaction of the wakes is modeled by use of the velocity deficit superposition principle. Finally, the power curve is used to estimate the power production from the wind turbines. The wind farm model is compared to large-eddy simulation (LES) data and measurments of Horns Rev wind farm for a wide range of wind directions. Reasonable agreement between the proposed analytical model, LES data and measurments is obtained. This prediction is also found to be substantially better than the one obtained with a commonly used wind farm wake model.
Methodology to set up accurate OPC model using optical CD metrology and atomic force microscopy
NASA Astrophysics Data System (ADS)
Shim, Yeon-Ah; Kang, Jaehyun; Lee, Sang-Uk; Kim, Jeahee; Kim, Keeho
2007-03-01
For the 90nm node and beyond, smaller Critical Dimension(CD) control budget is required and the ways to control good CD uniformity are needed. Moreover Optical Proximity Correction(OPC) for the sub-90nm node demands more accurate wafer CD data in order to improve accuracy of OPC model. Scanning Electron Microscope (SEM) is the typical method for measuring CD until ArF process. However SEM can give serious attack such as shrinkage of Photo Resist(PR) by burning of weak chemical structure of ArF PR due to high energy electron beam. In fact about 5nm CD narrowing occur when we measure CD by using CD-SEM in ArF photo process. Optical CD Metrology(OCD) and Atomic Force Microscopy(AFM) has been considered to the method for measuring CD without attack of organic materials. Also the OCD and AFM measurement system have the merits of speed, easiness and accurate data. For model-based OPC, the model is generated using CD data of test patterns transferred onto the wafer. In this study we discuss to generate accurate OPC model using OCD and AFM measurement system.
Model and Analytic Processes for Export License Assessments
Thompson, Sandra E.; Whitney, Paul D.; Weimar, Mark R.; Wood, Thomas W.; Daly, Don S.; Brothers, Alan J.; Sanfilippo, Antonio P.; Cook, Diane; Holder, Larry
2011-09-29
This paper represents the Department of Energy Office of Nonproliferation Research and Development (NA-22) Simulations, Algorithms and Modeling (SAM) Program's first effort to identify and frame analytical methods and tools to aid export control professionals in effectively predicting proliferation intent; a complex, multi-step and multi-agency process. The report focuses on analytical modeling methodologies that alone, or combined, may improve the proliferation export control license approval process. It is a follow-up to an earlier paper describing information sources and environments related to international nuclear technology transfer. This report describes the decision criteria used to evaluate modeling techniques and tools to determine which approaches will be investigated during the final 2 years of the project. The report also details the motivation for why new modeling techniques and tools are needed. The analytical modeling methodologies will enable analysts to evaluate the information environment for relevance to detecting proliferation intent, with specific focus on assessing risks associated with transferring dual-use technologies. Dual-use technologies can be used in both weapons and commercial enterprises. A decision-framework was developed to evaluate which of the different analytical modeling methodologies would be most appropriate conditional on the uniqueness of the approach, data availability, laboratory capabilities, relevance to NA-22 and Office of Arms Control and Nonproliferation (NA-24) research needs and the impact if successful. Modeling methodologies were divided into whether they could help micro-level assessments (e.g., help improve individual license assessments) or macro-level assessment. Macro-level assessment focuses on suppliers, technology, consumers, economies, and proliferation context. Macro-level assessment technologies scored higher in the area of uniqueness because less work has been done at the macro level. An approach to
Aspirating Seal Development: Analytical Modeling and Seal Test Rig
NASA Technical Reports Server (NTRS)
Bagepalli, Bharat
1996-01-01
This effort is to develop large diameter (22 - 36 inch) Aspirating Seals for application in aircraft engines. Stein Seal Co. will be fabricating the 36-inch seal(s) for testing. GE's task is to establish a thorough understanding of the operation of Aspirating Seals through analytical modeling and full-scale testing. The two primary objectives of this project are to develop the analytical models of the aspirating seal system, to upgrade using GE's funds, GE's 50-inch seal test rig for testing the Aspirating Seal (back-to-back with a corresponding brush seal), test the aspirating seal(s) for seal closure, tracking and maneuver transients (tilt) at operating pressures and temperatures, and validate the analytical model. The objective of the analytical model development is to evaluate the transient and steady-state dynamic performance characteristics of the seal designed by Stein. The transient dynamic model uses a multi-body system approach: the Stator, Seal face and the rotor are treated as individual bodies with relative degrees of freedom. Initially, the thirty-six springs are represented as a single one trying to keep open the aspirating face. Stops (Contact elements) are provided between the stator and the seal (to compensate the preload in the fully-open position) and between the rotor face and Seal face (to detect rub). The secondary seal is considered as part of the stator. The film's load, damping and stiffness characteristics as functions of pressure and clearance are evaluated using a separate (NASA) code GFACE. Initially, a laminar flow theory is used. Special two-dimensional interpolation routines are written to establish exact film load and damping values at each integration time step. Additionally, other user-routines are written to read-in actual pressure, rpm, stator-growth and rotor growth data and, later, to transfer these as appropriate loads/motions in the system-dynamic model. The transient dynamic model evaluates the various motions, clearances
Peat pyrolysis and the analytical semi-empirical model
Feng, J.; Green, A.E.S.
2007-07-01
Pyrolysis of peat could convert this material into useful fuels and valuable hydrocarbons. A study of peat pyrolysis can also serve as a useful bridge between studies of coal pyrolysis and biomass pyrolysis. Using an analytical model of pyrolysis that has previously been applied to biomass and to coal, we present here the results of applications of this model to a representative peat. The analysis suggests means of organizing and processing rate and yield data that should be useful in applications of pyrolysis for the production of fuels and chemicals.
Analytical properties of a three-compartmental dynamical demographic model
NASA Astrophysics Data System (ADS)
Postnikov, E. B.
2015-07-01
The three-compartmental demographic model by Korotaeyv-Malkov-Khaltourina, connecting population size, economic surplus, and education level, is considered from the point of view of dynamical systems theory. It is shown that there exist two integrals of motion, which enables the system to be reduced to one nonlinear ordinary differential equation. The study of its structure provides analytical criteria for the dominance ranges of the dynamics of Malthus and Kremer. Additionally, the particular ranges of parameters enable the derived general ordinary differential equations to be reduced to the models of Gompertz and Thoularis-Wallace.
An Analytical Thermal Model for Autonomous Soaring Research
NASA Technical Reports Server (NTRS)
Allen, Michael
2006-01-01
A viewgraph presentation describing an analytical thermal model used to enable research on autonomous soaring for a small UAV aircraft is given. The topics include: 1) Purpose; 2) Approach; 3) SURFRAD Data; 4) Convective Layer Thickness; 5) Surface Heat Budget; 6) Surface Virtual Potential Temperature Flux; 7) Convective Scaling Velocity; 8) Other Calculations; 9) Yearly trends; 10) Scale Factors; 11) Scale Factor Test Matrix; 12) Statistical Model; 13) Updraft Strength Calculation; 14) Updraft Diameter; 15) Updraft Shape; 16) Smoothed Updraft Shape; 17) Updraft Spacing; 18) Environment Sink; 19) Updraft Lifespan; 20) Autonomous Soaring Research; 21) Planned Flight Test; and 22) Mixing Ratio.
An Analytic Model for Buoyancy Resonances in Protoplanetary Disks
NASA Astrophysics Data System (ADS)
Lubow, Stephen H.; Zhu, Zhaohuan
2014-04-01
Zhu et al. found in three-dimensional shearing box simulations a new form of planet-disk interaction that they attributed to a vertical buoyancy resonance in the disk. We describe an analytic linear model for this interaction. We adopt a simplified model involving azimuthal forcing that produces the resonance and permits an analytic description of its structure. We derive an analytic expression for the buoyancy torque and show that the vertical torque distribution agrees well with the results of the Athena simulations and a Fourier method for linear numerical calculations carried out with the same forcing. The buoyancy resonance differs from the classic Lindblad and corotation resonances in that the resonance lies along tilted planes. Its width depends on damping effects and is independent of the gas sound speed. The resonance does not excite propagating waves. At a given large azimuthal wavenumber ky > h -1 (for disk thickness h), the buoyancy resonance exerts a torque over a region that lies radially closer to the corotation radius than the Lindblad resonance. Because the torque is localized to the region of excitation, it is potentially subject to the effects of nonlinear saturation. In addition, the torque can be reduced by the effects of radiative heat transfer between the resonant region and its surroundings. For each azimuthal wavenumber, the resonance establishes a large scale density wave pattern in a plane within the disk.
Analytical dynamics models for space missions around minor bodies
NASA Astrophysics Data System (ADS)
Cardoso dos Santos, Josué; dos Santos Carvalho, Jean Paulo; Vilhena de Moraes, Rodolpho; Bertachini de Almeida Prado, Antônio Fernando
2015-08-01
In recent years, the dynamics of orbits around minor bodies and icy moons in our solar system has become important in planning future missions that intend to visit dwarf planets, planetary moons, asteroids and comets. Due to their special characteristics, Europa, Ganymede, Callisto, Enceladus, Titan and Triton are among the group of objects with greater potential to receive missions in a near future. In order to provide a semi-analytical theory for tuture space exploration of these celestial bodies, this work aims to present two analytical models to describe and evaluate gravitational disturbances over a spacecrat's orbit around a minor body. A search for these less perturbed orbits is performed. An analytical model for the third-body perturbation is presented and consideres it in an eccentric-inclined orbit. Some harmonic terms due to the non-uniform distribuition of mass are considered according they are available in the literature. The dynamic of these orbits is explored by numerical simulations. The results are in accordance with the requirements for missions present in the literature.
An analytic model for buoyancy resonances in protoplanetary disks
Lubow, Stephen H.; Zhu, Zhaohuan E-mail: zhzhu@astro.princeton.edu
2014-04-10
Zhu et al. found in three-dimensional shearing box simulations a new form of planet-disk interaction that they attributed to a vertical buoyancy resonance in the disk. We describe an analytic linear model for this interaction. We adopt a simplified model involving azimuthal forcing that produces the resonance and permits an analytic description of its structure. We derive an analytic expression for the buoyancy torque and show that the vertical torque distribution agrees well with the results of the Athena simulations and a Fourier method for linear numerical calculations carried out with the same forcing. The buoyancy resonance differs from the classic Lindblad and corotation resonances in that the resonance lies along tilted planes. Its width depends on damping effects and is independent of the gas sound speed. The resonance does not excite propagating waves. At a given large azimuthal wavenumber k{sub y} > h {sup –1} (for disk thickness h), the buoyancy resonance exerts a torque over a region that lies radially closer to the corotation radius than the Lindblad resonance. Because the torque is localized to the region of excitation, it is potentially subject to the effects of nonlinear saturation. In addition, the torque can be reduced by the effects of radiative heat transfer between the resonant region and its surroundings. For each azimuthal wavenumber, the resonance establishes a large scale density wave pattern in a plane within the disk.
Can phenological models predict tree phenology accurately under climate change conditions?
NASA Astrophysics Data System (ADS)
Chuine, Isabelle; Bonhomme, Marc; Legave, Jean Michel; García de Cortázar-Atauri, Inaki; Charrier, Guillaume; Lacointe, André; Améglio, Thierry
2014-05-01
The onset of the growing season of trees has been globally earlier by 2.3 days/decade during the last 50 years because of global warming and this trend is predicted to continue according to climate forecast. The effect of temperature on plant phenology is however not linear because temperature has a dual effect on bud development. On one hand, low temperatures are necessary to break bud dormancy, and on the other hand higher temperatures are necessary to promote bud cells growth afterwards. Increasing phenological changes in temperate woody species have strong impacts on forest trees distribution and productivity, as well as crops cultivation areas. Accurate predictions of trees phenology are therefore a prerequisite to understand and foresee the impacts of climate change on forests and agrosystems. Different process-based models have been developed in the last two decades to predict the date of budburst or flowering of woody species. They are two main families: (1) one-phase models which consider only the ecodormancy phase and make the assumption that endodormancy is always broken before adequate climatic conditions for cell growth occur; and (2) two-phase models which consider both the endodormancy and ecodormancy phases and predict a date of dormancy break which varies from year to year. So far, one-phase models have been able to predict accurately tree bud break and flowering under historical climate. However, because they do not consider what happens prior to ecodormancy, and especially the possible negative effect of winter temperature warming on dormancy break, it seems unlikely that they can provide accurate predictions in future climate conditions. It is indeed well known that a lack of low temperature results in abnormal pattern of bud break and development in temperate fruit trees. An accurate modelling of the dormancy break date has thus become a major issue in phenology modelling. Two-phases phenological models predict that global warming should delay
NASA Astrophysics Data System (ADS)
Pop, Eric; Dutton, Robert W.; Goodson, Kenneth E.
2004-11-01
We describe the implementation of a Monte Carlo model for electron transport in silicon. The model uses analytic, nonparabolic electron energy bands, which are computationally efficient and sufficiently accurate for future low-voltage (<1V) nanoscale device applications. The electron-lattice scattering is incorporated using an isotropic, analytic phonon-dispersion model, which distinguishes between the optical/acoustic and the longitudinal/transverse phonon branches. We show that this approach avoids introducing unphysical thresholds in the electron distribution function, and that it has further applications in computing detailed phonon generation spectra from Joule heating. A set of deformation potentials for electron-phonon scattering is introduced and shown to yield accurate transport simulations in bulk silicon across a wide range of electric fields and temperatures. The shear deformation potential is empirically determined at Ξu=6.8eV, and consequently, the isotropically averaged scattering potentials with longitudinal and transverse acoustic phonons are DLA=6.39eV and DTA=3.01eV, respectively, in reasonable agreement with previous studies. The room-temperature electron mobility in strained silicon is also computed and shown to be in better agreement with the most recent phonon-limited data available. As a result, we find that electron coupling with g-type phonons is about 40% lower, and the coupling with f-type phonons is almost twice as strong as previously reported.
Universal analytic model for tunnel FET circuit simulation
NASA Astrophysics Data System (ADS)
Lu, Hao; Esseni, David; Seabaugh, Alan
2015-06-01
A simple analytic model based on the Kane-Sze formula is used to describe the current-voltage characteristics of tunnel field-effect transistors (TFETs). This model captures the unique features of the TFET including the decrease in subthreshold swing with drain current and the superlinear onset of the output characteristic. The model also captures the ambipolar current characteristic at negative gate-source bias and the negative differential resistance for negative drain-source biases. A simple empirical capacitance model is also included to enable circuit simulation. The model has fairly general validity and is not specific to a particular TFET geometry. Good agreement is shown with published atomistic simulations of an InAs double-gate TFET with gate perpendicular to the tunnel junction and with numerical simulations of a broken-gap AlGaSb/InAs TFET with gate in parallel with the tunnel junction.
Seth A Veitzer
2008-10-21
Effects of stray electrons are a main factor limiting performance of many accelerators. Because heavy-ion fusion (HIF) accelerators will operate in regimes of higher current and with walls much closer to the beam than accelerators operating today, stray electrons might have a large, detrimental effect on the performance of an HIF accelerator. A primary source of stray electrons is electrons generated when halo ions strike the beam pipe walls. There is some research on these types of secondary electrons for the HIF community to draw upon, but this work is missing one crucial ingredient: the effect of grazing incidence. The overall goal of this project was to develop the numerical tools necessary to accurately model the effect of grazing incidence on the behavior of halo ions in a HIF accelerator, and further, to provide accurate models of heavy ion stopping powers with applications to ICF, WDM, and HEDP experiments.
Building an accurate 3D model of a circular feature for robot vision
NASA Astrophysics Data System (ADS)
Li, L.
2012-06-01
In this paper, an accurate 3D model analysis of a circular feature is built with error compensation for robot vision. We propose an efficient method of fitting ellipses to data points by minimizing the algebraic distance subject to the constraint that a conic should be an ellipse and solving the ellipse parameters through a direct ellipse fitting method by analysing the 3D geometrical representation in a perspective projection scheme, the 3D position of a circular feature with known radius can be obtained. A set of identical circles, machined on a calibration board whose centres were known, was calibrated with a camera and did the model analysis that our method developed. Experimental results show that our method is more accurate than other methods.
Analytical model of reactive transport processes with spatially variable coefficients.
Simpson, Matthew J; Morrow, Liam C
2015-05-01
Analytical solutions of partial differential equation (PDE) models describing reactive transport phenomena in saturated porous media are often used as screening tools to provide insight into contaminant fate and transport processes. While many practical modelling scenarios involve spatially variable coefficients, such as spatially variable flow velocity, v(x), or spatially variable decay rate, k(x), most analytical models deal with constant coefficients. Here we present a framework for constructing exact solutions of PDE models of reactive transport. Our approach is relevant for advection-dominant problems, and is based on a regular perturbation technique. We present a description of the solution technique for a range of one-dimensional scenarios involving constant and variable coefficients, and we show that the solutions compare well with numerical approximations. Our general approach applies to a range of initial conditions and various forms of v(x) and k(x). Instead of simply documenting specific solutions for particular cases, we present a symbolic worksheet, as supplementary material, which enables the solution to be evaluated for different choices of the initial condition, v(x) and k(x). We also discuss how the technique generalizes to apply to models of coupled multispecies reactive transport as well as higher dimensional problems. PMID:26064648
Analytical modeling of glucose biosensors based on carbon nanotubes
2014-01-01
In recent years, carbon nanotubes have received widespread attention as promising carbon-based nanoelectronic devices. Due to their exceptional physical, chemical, and electrical properties, namely a high surface-to-volume ratio, their enhanced electron transfer properties, and their high thermal conductivity, carbon nanotubes can be used effectively as electrochemical sensors. The integration of carbon nanotubes with a functional group provides a good and solid support for the immobilization of enzymes. The determination of glucose levels using biosensors, particularly in the medical diagnostics and food industries, is gaining mass appeal. Glucose biosensors detect the glucose molecule by catalyzing glucose to gluconic acid and hydrogen peroxide in the presence of oxygen. This action provides high accuracy and a quick detection rate. In this paper, a single-wall carbon nanotube field-effect transistor biosensor for glucose detection is analytically modeled. In the proposed model, the glucose concentration is presented as a function of gate voltage. Subsequently, the proposed model is compared with existing experimental data. A good consensus between the model and the experimental data is reported. The simulated data demonstrate that the analytical model can be employed with an electrochemical glucose sensor to predict the behavior of the sensing mechanism in biosensors. PMID:24428818
Analytical model of reactive transport processes with spatially variable coefficients
Simpson, Matthew J.; Morrow, Liam C.
2015-01-01
Analytical solutions of partial differential equation (PDE) models describing reactive transport phenomena in saturated porous media are often used as screening tools to provide insight into contaminant fate and transport processes. While many practical modelling scenarios involve spatially variable coefficients, such as spatially variable flow velocity, v(x), or spatially variable decay rate, k(x), most analytical models deal with constant coefficients. Here we present a framework for constructing exact solutions of PDE models of reactive transport. Our approach is relevant for advection-dominant problems, and is based on a regular perturbation technique. We present a description of the solution technique for a range of one-dimensional scenarios involving constant and variable coefficients, and we show that the solutions compare well with numerical approximations. Our general approach applies to a range of initial conditions and various forms of v(x) and k(x). Instead of simply documenting specific solutions for particular cases, we present a symbolic worksheet, as supplementary material, which enables the solution to be evaluated for different choices of the initial condition, v(x) and k(x). We also discuss how the technique generalizes to apply to models of coupled multispecies reactive transport as well as higher dimensional problems. PMID:26064648
Długosz, Maciej; Antosiewicz, Jan M
2015-07-01
Proper treatment of hydrodynamic interactions is of importance in evaluation of rigid-body mobility tensors of biomolecules in Stokes flow and in simulations of their folding and solution conformation, as well as in simulations of the translational and rotational dynamics of either flexible or rigid molecules in biological systems at low Reynolds numbers. With macromolecules conveniently modeled in calculations or in dynamic simulations as ensembles of spherical frictional elements, various approximations to hydrodynamic interactions, such as the two-body, far-field Rotne-Prager approach, are commonly used, either without concern or as a compromise between the accuracy and the numerical complexity. Strikingly, even though the analytical Rotne-Prager approach fails to describe (both in the qualitative and quantitative sense) mobilities in the simplest system consisting of two spheres, when the distance between their surfaces is of the order of their size, it is commonly applied to model hydrodynamic effects in macromolecular systems. Here, we closely investigate hydrodynamic effects in two and three-body systems, consisting of bead-shell molecular models, using either the analytical Rotne-Prager approach, or an accurate numerical scheme that correctly accounts for the many-body character of hydrodynamic interactions and their short-range behavior. We analyze mobilities, and translational and rotational velocities of bodies resulting from direct forces acting on them. We show, that with the sufficient number of frictional elements in hydrodynamic models of interacting bodies, the far-field approximation is able to provide a description of hydrodynamic effects that is in a reasonable qualitative as well as quantitative agreement with the description resulting from the application of the virtually exact numerical scheme, even for small separations between bodies. PMID:26068580
Research on bathymetry estimation by Worldview-2 based with the semi-analytical model
NASA Astrophysics Data System (ADS)
Sheng, L.; Bai, J.; Zhou, G.-W.; Zhao, Y.; Li, Y.-C.
2015-04-01
South Sea Islands of China are far away from the mainland, the reefs takes more than 95% of south sea, and most reefs scatter over interested dispute sensitive area. Thus, the methods of obtaining the reefs bathymetry accurately are urgent to be developed. Common used method, including sonar, airborne laser and remote sensing estimation, are limited by the long distance, large area and sensitive location. Remote sensing data provides an effective way for bathymetry estimation without touching over large area, by the relationship between spectrum information and bathymetry. Aimed at the water quality of the south sea of China, our paper develops a bathymetry estimation method without measured water depth. Firstly the semi-analytical optimization model of the theoretical interpretation models has been studied based on the genetic algorithm to optimize the model. Meanwhile, OpenMP parallel computing algorithm has been introduced to greatly increase the speed of the semi-analytical optimization model. One island of south sea in China is selected as our study area, the measured water depth are used to evaluate the accuracy of bathymetry estimation from Worldview-2 multispectral images. The results show that: the semi-analytical optimization model based on genetic algorithm has good results in our study area;the accuracy of estimated bathymetry in the 0-20 meters shallow water area is accepted.Semi-analytical optimization model based on genetic algorithm solves the problem of the bathymetry estimation without water depth measurement. Generally, our paper provides a new bathymetry estimation method for the sensitive reefs far away from mainland.
Compact analytical modeling of SOI partially depleted MOSFETs with LETISOI
NASA Astrophysics Data System (ADS)
Faynot, O.; Poiroux, T.; Pelloie, J. L.
2001-04-01
As SOI technology becomes very attractive for ULSI CMOS, a dedicated and accurate SOI model has to be developed in order to take into account all specific electrical effects related to the SOI structure. The LETISOI model has been developed for partially depleted SOI devices; current and charge equations are built on a physical basis. In addition to the classical MOS conduction, FB effects, self-heating and bipolar transistor action have to be accurately modeled. The transient behavior, very different from the static one, has to be analyzed and kept in mind by designers to take full benefit of SOI devices in circuits, while avoiding any design issue due to SOI. Both FB and body-contacted devices can be described with this model. This paper describes how the model is built with an emphasis on the specific SOI needs and the related strategy for parameter extraction. Typical simulation results are also presented, outlining the capability of the model to simulate SOI specific dynamic behavior like bipolar activation and history effects.
Oksel, Ceyda; Winkler, David A; Ma, Cai Y; Wilkins, Terry; Wang, Xue Z
2016-09-01
The number of engineered nanomaterials (ENMs) being exploited commercially is growing rapidly, due to the novel properties they exhibit. Clearly, it is important to understand and minimize any risks to health or the environment posed by the presence of ENMs. Data-driven models that decode the relationships between the biological activities of ENMs and their physicochemical characteristics provide an attractive means of maximizing the value of scarce and expensive experimental data. Although such structure-activity relationship (SAR) methods have become very useful tools for modelling nanotoxicity endpoints (nanoSAR), they have limited robustness and predictivity and, most importantly, interpretation of the models they generate is often very difficult. New computational modelling tools or new ways of using existing tools are required to model the relatively sparse and sometimes lower quality data on the biological effects of ENMs. The most commonly used SAR modelling methods work best with large datasets, are not particularly good at feature selection, can be relatively opaque to interpretation, and may not account for nonlinearity in the structure-property relationships. To overcome these limitations, we describe the application of a novel algorithm, a genetic programming-based decision tree construction tool (GPTree) to nanoSAR modelling. We demonstrate the use of GPTree in the construction of accurate and interpretable nanoSAR models by applying it to four diverse literature datasets. We describe the algorithm and compare model results across the four studies. We show that GPTree generates models with accuracies equivalent to or superior to those of prior modelling studies on the same datasets. GPTree is a robust, automatic method for generation of accurate nanoSAR models with important advantages that it works with small datasets, automatically selects descriptors, and provides significantly improved interpretability of models. PMID:26956430
Accurate and efficient halo-based galaxy clustering modelling with simulations
NASA Astrophysics Data System (ADS)
Zheng, Zheng; Guo, Hong
2016-06-01
Small- and intermediate-scale galaxy clustering can be used to establish the galaxy-halo connection to study galaxy formation and evolution and to tighten constraints on cosmological parameters. With the increasing precision of galaxy clustering measurements from ongoing and forthcoming large galaxy surveys, accurate models are required to interpret the data and extract relevant information. We introduce a method based on high-resolution N-body simulations to accurately and efficiently model the galaxy two-point correlation functions (2PCFs) in projected and redshift spaces. The basic idea is to tabulate all information of haloes in the simulations necessary for computing the galaxy 2PCFs within the framework of halo occupation distribution or conditional luminosity function. It is equivalent to populating galaxies to dark matter haloes and using the mock 2PCF measurements as the model predictions. Besides the accurate 2PCF calculations, the method is also fast and therefore enables an efficient exploration of the parameter space. As an example of the method, we decompose the redshift-space galaxy 2PCF into different components based on the type of galaxy pairs and show the redshift-space distortion effect in each component. The generalizations and limitations of the method are discussed.
NASA Astrophysics Data System (ADS)
Blackman, Jonathan; Field, Scott E.; Galley, Chad R.; Szilágyi, Béla; Scheel, Mark A.; Tiglio, Manuel; Hemberger, Daniel A.
2015-09-01
Simulating a binary black hole coalescence by solving Einstein's equations is computationally expensive, requiring days to months of supercomputing time. Using reduced order modeling techniques, we construct an accurate surrogate model, which is evaluated in a millisecond to a second, for numerical relativity (NR) waveforms from nonspinning binary black hole coalescences with mass ratios in [1, 10] and durations corresponding to about 15 orbits before merger. We assess the model's uncertainty and show that our modeling strategy predicts NR waveforms not used for the surrogate's training with errors nearly as small as the numerical error of the NR code. Our model includes all spherical-harmonic -2Yℓm waveform modes resolved by the NR code up to ℓ=8 . We compare our surrogate model to effective one body waveforms from 50 M⊙ to 300 M⊙ for advanced LIGO detectors and find that the surrogate is always more faithful (by at least an order of magnitude in most cases).
Chuine, Isabelle; Bonhomme, Marc; Legave, Jean-Michel; García de Cortázar-Atauri, Iñaki; Charrier, Guillaume; Lacointe, André; Améglio, Thierry
2016-10-01
The onset of the growing season of trees has been earlier by 2.3 days per decade during the last 40 years in temperate Europe because of global warming. The effect of temperature on plant phenology is, however, not linear because temperature has a dual effect on bud development. On one hand, low temperatures are necessary to break bud endodormancy, and, on the other hand, higher temperatures are necessary to promote bud cell growth afterward. Different process-based models have been developed in the last decades to predict the date of budbreak of woody species. They predict that global warming should delay or compromise endodormancy break at the species equatorward range limits leading to a delay or even impossibility to flower or set new leaves. These models are classically parameterized with flowering or budbreak dates only, with no information on the endodormancy break date because this information is very scarce. Here, we evaluated the efficiency of a set of phenological models to accurately predict the endodormancy break dates of three fruit trees. Our results show that models calibrated solely with budbreak dates usually do not accurately predict the endodormancy break date. Providing endodormancy break date for the model parameterization results in much more accurate prediction of this latter, with, however, a higher error than that on budbreak dates. Most importantly, we show that models not calibrated with endodormancy break dates can generate large discrepancies in forecasted budbreak dates when using climate scenarios as compared to models calibrated with endodormancy break dates. This discrepancy increases with mean annual temperature and is therefore the strongest after 2050 in the southernmost regions. Our results claim for the urgent need of massive measurements of endodormancy break dates in forest and fruit trees to yield more robust projections of phenological changes in a near future. PMID:27272707