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.
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.
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.
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.
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 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.
Analytic model of the radiation-dominated decay of a compact toroid
Auerbach, S.P.
1981-11-11
The coaxial-gun, compact-torus experiments at LLNL and LASNL are believed to be radiation-dominated, in the sense that most or all of the input energy is lost by impurity radiation. This paper presents a simple analytic model of the radiation-dominated decay of a compact torus, and demonstrates that several striking features of the experiment (finite lifetime, linear current decay, insensitivity of the lifetime to density or stored magnetic energy) may also be explained by the hypothesis that impurity radiation dominates the energy loss. The model incorporates the essential features of the more elaborate 1 1/2-D simulations of Shumaker et al., yet is simple enough to be solved exactly. Based on the analytic results, a simple criterion is given for the maximum tolerable impurity density.
Analytical-HZETRN model for rapid assessment of active magnetic radiation shielding
NASA Astrophysics Data System (ADS)
Washburn, S. A.; Blattnig, S. R.; Singleterry, R. C.; Westover, S. C.
2014-01-01
The use of active radiation shielding designs has the potential to reduce the radiation exposure received by astronauts on deep-space missions at a significantly lower mass penalty than designs utilizing only passive shielding. Unfortunately, the determination of the radiation exposure inside these shielded environments often involves lengthy and computationally intensive Monte Carlo analysis. In order to evaluate the large trade space of design parameters associated with a magnetic radiation shield design, an analytical model was developed for the determination of flux inside a solenoid magnetic field due to the Galactic Cosmic Radiation (GCR) radiation environment. This analytical model was then coupled with NASA's radiation transport code, HZETRN, to account for the effects of passive/structural shielding mass. The resulting model can rapidly obtain results for a given configuration and can therefore be used to analyze an entire trade space of potential variables in less time than is required for even a single Monte Carlo run. Analyzing this trade space for a solenoid magnetic shield design indicates that active shield bending powers greater than ˜15 Tm and passive/structural shielding thicknesses greater than 40 g/cm2 have a limited impact on reducing dose equivalent values. Also, it is shown that higher magnetic field strengths are more effective than thicker magnetic fields at reducing dose equivalent.
Analytical-HZETRN Model for Rapid Assessment of Active Magnetic Radiation Shielding
NASA Technical Reports Server (NTRS)
Washburn, S. A.; Blattnig, S. R.; Singleterry, R. C.; Westover, S. C.
2014-01-01
The use of active radiation shielding designs has the potential to reduce the radiation exposure received by astronauts on deep-space missions at a significantly lower mass penalty than designs utilizing only passive shielding. Unfortunately, the determination of the radiation exposure inside these shielded environments often involves lengthy and computationally intensive Monte Carlo analysis. In order to evaluate the large trade space of design parameters associated with a magnetic radiation shield design, an analytical model was developed for the determination of flux inside a solenoid magnetic field due to the Galactic Cosmic Radiation (GCR) radiation environment. This analytical model was then coupled with NASA's radiation transport code, HZETRN, to account for the effects of passive/structural shielding mass. The resulting model can rapidly obtain results for a given configuration and can therefore be used to analyze an entire trade space of potential variables in less time than is required for even a single Monte Carlo run. Analyzing this trade space for a solenoid magnetic shield design indicates that active shield bending powers greater than 15 Tm and passive/structural shielding thicknesses greater than 40 g/cm2 have a limited impact on reducing dose equivalent values. Also, it is shown that higher magnetic field strengths are more effective than thicker magnetic fields at reducing dose equivalent.
Analytical model for radiative transfer including the effects of a rough material interface.
Giddings, Thomas E; Kellems, Anthony R
2016-08-20
The reflected and transmitted radiance due to a source located above a water surface is computed based on models for radiative transfer in continuous optical media separated by a discontinuous air-water interface with random surface roughness. The air-water interface is described as the superposition of random, unresolved roughness on a deterministic realization of a stochastic wave surface at resolved scales. Under the geometric optics assumption, the bidirectional reflection and transmission functions for the air-water interface are approximated by applying regular perturbation methods to Snell's law and including the effects of a random surface roughness component. Formal analytical solutions to the radiative transfer problem under the small-angle scattering approximation account for the effects of scattering and absorption as light propagates through the atmosphere and water and also capture the diffusive effects due to the interaction of light with the rough material interface that separates the two optical media. Results of the analytical models are validated against Monte Carlo simulations, and the approximation to the bidirectional reflection function is also compared to another well-known analytical model. PMID:27556978
Coefficients of an analytical aerosol forcing equation determined with a Monte-Carlo radiation model
NASA Astrophysics Data System (ADS)
Hassan, Taufiq; Moosmüller, H.; Chung, Chul E.
2015-10-01
Simple analytical equations for global-average direct aerosol radiative forcing are useful to quickly estimate aerosol forcing changes as function of key atmosphere, surface and aerosol parameters. The surface and atmosphere parameters in these analytical equations are the globally uniform atmospheric transmittance and surface albedo, and have so far been estimated from simplified observations under untested assumptions. In the present study, we take the state-of-the-art analytical equation and write the aerosol forcing as a linear function of the single scattering albedo (SSA) and replace the average upscatter fraction with the asymmetry parameter (ASY). Then we determine the surface and atmosphere parameter values of this equation using the output from the global MACR (Monte-Carlo Aerosol Cloud Radiation) model, as well as testing the validity of the equation. The MACR model incorporated spatio-temporally varying observations for surface albedo, cloud optical depth, water vapor, stratosphere column ozone, etc., instead of assuming as in the analytical equation that the atmosphere and surface parameters are globally uniform, and should thus be viewed as providing realistic radiation simulations. The modified analytical equation needs globally uniform aerosol parameters that consist of AOD (Aerosol Optical Depth), SSA, and ASY. The MACR model is run here with the same globally uniform aerosol parameters. The MACR model is also run without cloud to test the cloud effect. In both cloudy and cloud-free runs, the equation fits in the model output well whether SSA or ASY varies. This means the equation is an excellent approximation for the atmospheric radiation. On the other hand, the determined parameter values are somewhat realistic for the cloud-free runs but unrealistic for the cloudy runs. The global atmospheric transmittance, one of the determined parameters, is found to be around 0.74 in case of the cloud-free conditions and around 1.03 with cloud. The surface
Fast radiative transfer of dust reprocessing in semi-analytic models with artificial neural networks
NASA Astrophysics Data System (ADS)
Silva, Laura; Fontanot, Fabio; Granato, Gian Luigi
2012-06-01
A serious concern for semi-analytical galaxy formation models, aiming to simulate multiwavelength surveys and to thoroughly explore the model parameter space, is the extremely time-consuming numerical solution of the radiative transfer of stellar radiation through dusty media. To overcome this problem, we have implemented an artificial neural network (ANN) algorithm in the radiative transfer code GRASIL, in order to significantly speed up the computation of the infrared (IR) spectral energy distribution (SED). The ANN we have implemented is of general use, in that its input neurons are defined as those quantities effectively determining the shape of the IR SED. Therefore, the training of the ANN can be performed with any model and then applied to other models. We made a blind test to check the algorithm, by applying a net trained with a standard chemical evolution model (i.e. CHE_EVO) to a mock catalogue extracted from the semi-analytic model MORGANA, and compared galaxy counts and evolution of the luminosity functions in several near-IR to sub-millimetre (sub-mm) bands, and also the spectral differences for a large subset of randomly extracted models. The ANN is able to excellently approximate the full computation, but with a gain in CPU time by ˜2 orders of magnitude. It is only advisable that the training covers reasonably well the range of values of the input neurons in the application. Indeed in the sub-mm at high redshift, a tiny fraction of models with some sensible input neurons out of the range of the trained net cause wrong answer by the ANN. These are extreme starbursting models with high optical depths, favourably selected by sub-mm observations, and are difficult to predict a priori.
NASA Astrophysics Data System (ADS)
Heng, Kevin; Mendonça, João M.; Lee, Jae-Min
2014-11-01
We present a comprehensive analytical study of radiative transfer using the method of moments and include the effects of non-isotropic scattering in the coherent limit. Within this unified formalism, we derive the governing equations and solutions describing two-stream radiative transfer (which approximates the passage of radiation as a pair of outgoing and incoming fluxes), flux-limited diffusion (which describes radiative transfer in the deep interior), and solutions for the temperature-pressure profiles. Generally, the problem is mathematically underdetermined unless a set of closures (Eddington coefficients) is specified. We demonstrate that the hemispheric (or hemi-isotropic) closure naturally derives from the radiative transfer equation if energy conservation is obeyed, while the Eddington closure produces spurious enhancements of both reflected light and thermal emission. We concoct recipes for implementing two-stream radiative transfer in stand-alone numerical calculations and general circulation models. We use our two-stream solutions to construct toy models of the runaway greenhouse effect. We present a new solution for temperature-pressure profiles with a non-constant optical opacity and elucidate the effects of non-isotropic scattering in the optical and infrared. We derive generalized expressions for the spherical and Bond albedos and the photon deposition depth. We demonstrate that the value of the optical depth corresponding to the photosphere is not always 2/3 (Milne's solution) and depends on a combination of stellar irradiation, internal heat, and the properties of scattering in both the optical and infrared. Finally, we derive generalized expressions for the total, net, outgoing, and incoming fluxes in the convective regime.
Heng, Kevin; Mendonça, João M.; Lee, Jae-Min E-mail: joao.mendonca@csh.unibe.ch
2014-11-01
We present a comprehensive analytical study of radiative transfer using the method of moments and include the effects of non-isotropic scattering in the coherent limit. Within this unified formalism, we derive the governing equations and solutions describing two-stream radiative transfer (which approximates the passage of radiation as a pair of outgoing and incoming fluxes), flux-limited diffusion (which describes radiative transfer in the deep interior), and solutions for the temperature-pressure profiles. Generally, the problem is mathematically underdetermined unless a set of closures (Eddington coefficients) is specified. We demonstrate that the hemispheric (or hemi-isotropic) closure naturally derives from the radiative transfer equation if energy conservation is obeyed, while the Eddington closure produces spurious enhancements of both reflected light and thermal emission. We concoct recipes for implementing two-stream radiative transfer in stand-alone numerical calculations and general circulation models. We use our two-stream solutions to construct toy models of the runaway greenhouse effect. We present a new solution for temperature-pressure profiles with a non-constant optical opacity and elucidate the effects of non-isotropic scattering in the optical and infrared. We derive generalized expressions for the spherical and Bond albedos and the photon deposition depth. We demonstrate that the value of the optical depth corresponding to the photosphere is not always 2/3 (Milne's solution) and depends on a combination of stellar irradiation, internal heat, and the properties of scattering in both the optical and infrared. Finally, we derive generalized expressions for the total, net, outgoing, and incoming fluxes in the convective regime.
NASA Astrophysics Data System (ADS)
Bostater, Charles R., Jr.; Huddleston, Lisa H.
2000-12-01
A remote sensing reflectance model, which describes the transfer of irradiant light within a homogeneous water column has previously been used to simulate the nadir viewing reflectance just above or below the water surface by Bostater, et al. Wavelength dependent features in the water surface reflectance depend upon the nature of the down welling irradiance, bottom reflectance and the water absorption and backscatter coefficients. The latter are very important coefficients, and depend upon the constituents in water and both vary as a function of the water depth and wavelength in actual water bodies. This paper describes a preliminary approach for the analytical solution of the radiative transfer equations in a two-stream representation of the irradiance field with variable coefficients due to the depth dependent water concentrations of substances such as chlorophyl pigments, dissolved organic matter and suspended particulate matter. The analytical model formulation makes use of analytically based solutions to the 2-flow equations. However, in this paper we describe the use of the unique Cauchy boundary conditions previously used, along with a matrix solution to allow for the prediction of the synthetic water surface reflectance signatures within a nonhomogeneous medium. Observed reflectance signatures as well as model derived 'synthetic signatures' are processed using efficient algorithms which demonstrate the error induced using the layered matrix approach is much less than 1 percent when compared to the analytical homogeneous water column solution. The influence of vertical gradients of water constituents may be extremely important in remote sensing of coastal water constituents as well as in remote sensing of submerged targets and different bottom types such as corals, sea grasses and sand.
Sazykina, Tatiana G; Kryshev, Alexander I
2016-01-01
A dynamic mathematical model is formulated, predicting the development of radiation effects in a generic animal population, inhabiting an elemental ecosystem 'population-limiting resource'. Differential equations of the model describe the dynamic responses to radiation damage of the following population characteristics: gross biomass; intrinsic fractions of healthy and reversibly damaged tissues in biomass; intrinsic concentrations of the self-repairing pool and the growth factor; and amount of the limiting resource available in the environment. Analytical formulae are found for the steady states of model variables as non-linear functions of the dose rate of chronic radiation exposure. Analytical solutions make it possible to predict the expected severity of radiation effects in a model ecosystem, including such endpoints as morbidity, mortality, life shortening, biosynthesis, and population biomass. Model parameters are selected from species data on lifespan, physiological growth and mortality rates, and individual radiosensitivity. Thresholds for population extinction can be analytically calculated for different animal species, examples are provided for generic mice and wolf populations. The ecosystem model demonstrates a compensatory effect of the environment on the development of radiation effects in wildlife. The model can be employed to construct a preliminary scale 'radiation exposure-population effects' for different animal species; species can be identified, which are vulnerable at a population level to chronic radiation exposure. PMID:26408836
A new analytical solar radiation pressure model for current BeiDou satellites: IGGBSPM.
Tan, Bingfeng; Yuan, Yunbin; Zhang, Baocheng; Hsu, Hou Ze; Ou, Jikun
2016-01-01
An analytical solar radiation pressure (SRP) model, IGGBSPM (an abbreviation for Institute of Geodesy and Geophysics BeiDou Solar Pressure Model), has been developed for three BeiDou satellite types, namely, geostationary orbit (GEO), inclined geosynchronous orbit (IGSO) and medium earth orbit (MEO), based on a ray-tracing method. The performance of IGGBSPM was assessed based on numerical integration, SLR residuals and analyses of empirical SRP parameters (except overlap computations). The numerical results show that the integrated orbit resulting from IGGBSPM differs from the precise ephemerides by approximately 5 m and 2 m for GEO and non-GEO satellites, respectively. Moreover, when IGGBSPM is used as an a priori model to enhance the ECOM (5-parameter) model with stochastic pulses, named ECOM + APR, for precise orbit determination, the SLR RMS residual improves by approximately 20-25 percent over the ECOM-only solution during the yaw-steering period and by approximately 40 percent during the yaw-fixed period. For the BeiDou GEO01 satellite, improvements of 18 and 32 percent can be achieved during the out-of-eclipse season and during the eclipse season, respectively. An investigation of the estimated ECOM D0 parameters indicated that the β-angle dependence that is evident in the ECOM-only solution is no longer present in the ECOM + APR solution. PMID:27595795
A new analytical solar radiation pressure model for current BeiDou satellites: IGGBSPM
Tan, Bingfeng; Yuan, Yunbin; Zhang, Baocheng; Hsu, Hou Ze; Ou, Jikun
2016-01-01
An analytical solar radiation pressure (SRP) model, IGGBSPM (an abbreviation for Institute of Geodesy and Geophysics BeiDou Solar Pressure Model), has been developed for three BeiDou satellite types, namely, geostationary orbit (GEO), inclined geosynchronous orbit (IGSO) and medium earth orbit (MEO), based on a ray-tracing method. The performance of IGGBSPM was assessed based on numerical integration, SLR residuals and analyses of empirical SRP parameters (except overlap computations). The numerical results show that the integrated orbit resulting from IGGBSPM differs from the precise ephemerides by approximately 5 m and 2 m for GEO and non-GEO satellites, respectively. Moreover, when IGGBSPM is used as an a priori model to enhance the ECOM (5-parameter) model with stochastic pulses, named ECOM + APR, for precise orbit determination, the SLR RMS residual improves by approximately 20–25 percent over the ECOM-only solution during the yaw-steering period and by approximately 40 percent during the yaw-fixed period. For the BeiDou GEO01 satellite, improvements of 18 and 32 percent can be achieved during the out-of-eclipse season and during the eclipse season, respectively. An investigation of the estimated ECOM D0 parameters indicated that the β-angle dependence that is evident in the ECOM-only solution is no longer present in the ECOM + APR solution. PMID:27595795
Development of PARMA: PHITS-based analytical radiation model in the atmosphere.
Sato, Tatsuhiko; Yasuda, Hiroshi; Niita, Koji; Endo, Akira; Sihver, Lembit
2008-08-01
Estimation of cosmic-ray spectra in the atmosphere has been essential for the evaluation of aviation doses. We therefore calculated these spectra by performing Monte Carlo simulation of cosmic-ray propagation in the atmosphere using the PHITS code. The accuracy of the simulation was well verified by experimental data taken under various conditions, even near sea level. Based on a comprehensive analysis of the simulation results, we proposed an analytical model for estimating the cosmic-ray spectra of neutrons, protons, helium ions, muons, electrons, positrons and photons applicable to any location in the atmosphere at altitudes below 20 km. Our model, named PARMA, enables us to calculate the cosmic radiation doses rapidly with a precision equivalent to that of the Monte Carlo simulation, which requires much more computational time. With these properties, PARMA is capable of improving the accuracy and efficiency of the cosmic-ray exposure dose estimations not only for aircrews but also for the public on the ground. PMID:18666812
Verley, Jason C.; Axness, Carl L.; Hembree, Charles Edward; Keiter, Eric Richard; Kerr, Bert
2012-04-01
Photocurrent generated by ionizing radiation represents a threat to microelectronics in radiation environments. Circuit simulation tools such as SPICE [1] can be used to analyze these threats, and typically rely on compact models for individual electrical components such as transistors and diodes. Compact models consist of a handful of differential and/or algebraic equations, and are derived by making simplifying assumptions to any of the many semiconductor transport equations. Historically, many photocurrent compact models have suffered from accuracy issues due to the use of qualitative approximation, rather than mathematically correct solutions to the ambipolar diffusion equation. A practical consequence of this inaccuracy is that a given model calibration is trustworthy over only a narrow range of operating conditions. This report describes work to produce improved compact models for photocurrent. Specifically, an analytic model is developed for epitaxial diode structures that have a highly doped subcollector. The analytic model is compared with both numerical TCAD calculations, as well as the compact model described in reference [2]. The new analytic model compares well against TCAD over a wide range of operating conditions, and is shown to be superior to the compact model from reference [2].
NASA Astrophysics Data System (ADS)
Gaetano, A.; Roncolato, J.; Montorfano, D.; Barbato, M. C.; Ambrosetti, G.; Pedretti, A.
2016-05-01
The employment of new gaseous heat transfer fluids as air or CO2, which are cheaper and environmentally friendly, is drawing more and more attention within the field of Concentrated Solar Power applications. However, despite the advantages, their use requires receivers with a larger heat transfer area and flow cross section with a consequent greater volume of thermal insulation. Solid thermal insulations currently used present high thermal inertia which is energetically penalizing during the daily transient phases faced by the main plant components (e.g. receivers). With the aim of overcoming this drawback a thermal insulation based on radiative shields is presented in this study. Starting from an initial layout comprising a solid thermal insulation layer, the geometry was optimized avoiding the use of the solid insulation keeping performance and fulfilling the geometrical constraints. An analytical Matlab model was implemented to assess the system thermal behavior in terms of heat loss taking into account conductive, convective and radiative contributions. Accurate 2D Computational Fluid Dynamics (CFD) simulations were run to validate the Matlab model which was then used to select the most promising among three new different designs.
NASA Astrophysics Data System (ADS)
Upcraft, L. M.
2010-09-01
There is an ongoing need for numerically efficient algorithms that are capable of calculating the radiative and collisional rates of arbitrarily complex ions that are present in hot plasmas to a level of accuracy that surpasses that available in many existing approximations. Hydrogen-like solutions for determining these rates in more general ions by use of an effective (and generally non-integer) atomic number frequently give poor results and are of limited validity. This paper illustrates that results accurate to of order 20% can be obtained for matrix elements of both rates for arbitrarily complex ions by use of hydrogenic wavefunctions that use different effective atomic numbers for the initial and final sub-shells. Not only does this allow for the realistic modelling of inner shell transitions, it naturally allows for the physical effect of orbital relaxation. It is shown that the integral of the generalised oscillator strength used by the Plane-wave Born approximation has an analytic solution that can be reduced to a form suitable for efficient numerical integration over an arbitrary electron distribution. Extensive use of the computer algebra package Mathematica ® has generated a unique formula for each transition and the results have been transformed to efficient fortran 90 code for all transitions between non-relativistic sub-shells with principal quantum numbers n ≤ 10. In the case of the collisional matrix elements these are typically two to three orders of magnitude faster to calculate than by direct numerical evaluation. The fortran code is available upon request from the author.
ERIC Educational Resources Information Center
James, W. G. G.
1970-01-01
Discusses the historical development of both the wave and the corpuscular photon model of light. Suggests that students should be informed that the two models are complementary and that each model successfully describes a wide range of radiation phenomena. Cites 19 references which might be of interest to physics teachers and students. (LC)
Analytic approximate radiation effects due to Bremsstrahlung
Ben-Zvi I.
2012-02-01
The purpose of this note is to provide analytic approximate expressions that can provide quick estimates of the various effects of the Bremsstrahlung radiation produced relatively low energy electrons, such as the dumping of the beam into the beam stop at the ERL or field emission in superconducting cavities. The purpose of this work is not to replace a dependable calculation or, better yet, a measurement under real conditions, but to provide a quick but approximate estimate for guidance purposes only. These effects include dose to personnel, ozone generation in the air volume exposed to the radiation, hydrogen generation in the beam dump water cooling system and radiation damage to near-by magnets. These expressions can be used for other purposes, but one should note that the electron beam energy range is limited. In these calculations the good range is from about 0.5 MeV to 10 MeV. To help in the application of this note, calculations are presented as a worked out example for the beam dump of the R&D Energy Recovery Linac.
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
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.
Boyd, Kathleen A; Jones, Rob J; Paul, Jim; Birrell, Fiona; Briggs, Andrew H; Leung, Hing Y
2015-01-01
Objective To determine the cost-effectiveness of salvage cryotherapy (SC) in men with radiation recurrent prostate cancer (RRPC). Design Cost-utility analysis using decision analytic modelling by a Markov model. Setting and methods Compared SC and androgen deprivation therapy (ADT) in a cohort of patients with RRPC (biopsy proven local recurrence, no evidence of metastatic disease). A literature review captured published data to inform the decision model, and resource use data were from the Scottish Prostate Cryotherapy Service. The model was run in monthly cycles for RRPC men, mean age of 70 years. The model was run over the patient lifetime, to assess changes in patient health states and the associated quality of life, survival and cost impacts. Results are reported in terms of the discounted incremental costs and discounted incremental quality-adjusted life years (QALYs) gained between the 2 alternative interventions. Probabilistic sensitivity analysis used a 10 000 iteration Monte Carlo simulation. Results SC has a high upfront treatment cost, but delays the ongoing monthly cost of ADT. SC is the dominant strategy over the patient lifetime; it is more effective with an incremental 0.56 QALY gain (95% CI 0.28 to 0.87), and less costly with a reduced lifetime cost of £29 719 (€37 619) (95% CI −51 985 to −9243). For a ceiling ratio of £30 000, SC has a 100% probability to be cost-effective. The cost neutral point was at 3.5 years, when the upfront cost of SC (plus any subsequent cumulative cost of side effects and ADT) equates the cumulative cost in the ADT arm. Limitations of our model may arise from its insensitivity to parameter or structural uncertainty. Conclusions The platform for SC versus ADT cost-effective analysis can be employed to evaluate other treatment modalities or strategies in RRPC. SC is the dominant strategy, costing less over a patient's lifetime with improvements in QALYs. Trial registration number This economic analysis
Analytic expressions for ULF wave radiation belt radial diffusion coefficients
Ozeke, Louis G; Mann, Ian R; Murphy, Kyle R; Jonathan Rae, I; Milling, David K
2014-01-01
We present analytic expressions for ULF wave-derived radiation belt radial diffusion coefficients, as a function of L and Kp, which can easily be incorporated into global radiation belt transport models. The diffusion coefficients are derived from statistical representations of ULF wave power, electric field power mapped from ground magnetometer data, and compressional magnetic field power from in situ measurements. We show that the overall electric and magnetic diffusion coefficients are to a good approximation both independent of energy. We present example 1-D radial diffusion results from simulations driven by CRRES-observed time-dependent energy spectra at the outer boundary, under the action of radial diffusion driven by the new ULF wave radial diffusion coefficients and with empirical chorus wave loss terms (as a function of energy, Kp and L). There is excellent agreement between the differential flux produced by the 1-D, Kp-driven, radial diffusion model and CRRES observations of differential electron flux at 0.976 MeV—even though the model does not include the effects of local internal acceleration sources. Our results highlight not only the importance of correct specification of radial diffusion coefficients for developing accurate models but also show significant promise for belt specification based on relatively simple models driven by solar wind parameters such as solar wind speed or geomagnetic indices such as Kp. Key Points Analytic expressions for the radial diffusion coefficients are presented The coefficients do not dependent on energy or wave m value The electric field diffusion coefficient dominates over the magnetic PMID:26167440
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).
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; ...
NASA Astrophysics Data System (ADS)
Stupakov, Gennady; Zhou, Demin
2016-04-01
We develop a general model of coherent synchrotron radiation (CSR) impedance with shielding provided by two parallel conducting plates. This model allows us to easily reproduce all previously known analytical CSR wakes and to expand the analysis to situations not explored before. It reduces calculations of the impedance to taking integrals along the trajectory of the beam. New analytical results are derived for the radiation impedance with shielding for the following orbits: a kink, a bending magnet, a wiggler of finite length, and an infinitely long wiggler. All our formulas are benchmarked against numerical simulations with the CSRZ computer code.
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
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.
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.
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.
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...
Analytical Characterization of CFRP Laser Treated by Excimer Laser Radiation
NASA Astrophysics Data System (ADS)
Kreling, S.; Fischer, F.; Delmdahl, R.; Gäbler, F.; Dilger, K.
Due to the increasing interest in lightweight structures, carbon-fiber reinforced plastics are increasingly applied, especially in the transportation industry. An interesting technology for joining these materials is adhesive bonding due to numerous advantages compared to conventional techniques like riveting. However, to achieve a strong and durable bond, surface pre-treatment is necessary to remove residues of release agents that are transferred to the surface during manufacturing. This paper describes analytical experiments, namely SEM and XPS, performed on CFRP surfaces pre-treated with 308 nm excimer laser radiation.
Analytical properties of the radiance in atmospheric radiative transfer theory
NASA Astrophysics Data System (ADS)
Otto, Sebastian
2014-01-01
It is demonstrated mathematically strictly that state density functions, as the radiance (specific intensity), exist to describe certain state properties of transported photons on microscopic and the state of the radiation field on macroscopic scale, which have independent physical meanings. Analytical properties as boundedness, continuity, differentiability and integrability of these functions to describe the photon transport are discussed. It is shown that the density functions may be derived based on the assumption of photons as real particles of non-zero and finite size, independently of usual electrodynamics, and certain historically postulated functional relationships between them were proved, that is, these functions can be derived mathematically strictly and consistently within the framework of the theory of the phenomenological radiative transfer if one takes the theory seriously by really assuming photons as particles. In this sense these functions may be treated as fundamental physical quantities within the scope of this theory, if one considers the possibility of the existence of photons.
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
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.
Saturn Radiation (SATRAD) Model
NASA Technical Reports Server (NTRS)
Garrett, H. B.; Ratliff, J. M.; Evans, R. W.
2005-01-01
The Saturnian radiation belts have not received as much attention as the Jovian radiation belts because they are not nearly as intense-the famous Saturnian particle rings tend to deplete the belts near where their peak would occur. As a result, there has not been a systematic development of engineering models of the Saturnian radiation environment for mission design. A primary exception is that of Divine (1990). That study used published data from several charged particle experiments aboard the Pioneer 1 1, Voyager 1, and Voyager 2 spacecraft during their flybys at Saturn to generate numerical models for the electron and proton radiation belts between 2.3 and 13 Saturn radii. The Divine Saturn radiation model described the electron distributions at energies between 0.04 and 10 MeV and the proton distributions at energies between 0.14 and 80 MeV. The model was intended to predict particle intensity, flux, and fluence for the Cassini orbiter. Divine carried out hand calculations using the model but never formally developed a computer program that could be used for general mission analyses. This report seeks to fill that void by formally developing a FORTRAN version of the model that can be used as a computer design tool for missions to Saturn that require estimates of the radiation environment around the planet. The results of that effort and the program listings are presented here along with comparisons with the original estimates carried out by Divine. In addition, Pioneer and Voyager data were scanned in from the original references and compared with the FORTRAN model s predictions. The results were statistically analyzed in a manner consistent with Divine s approach to provide estimates of the ability of the model to reproduce the original data. Results of a formal review of the model by a panel of experts are also presented. Their recommendations for further tests, analyses, and extensions to the model are discussed.
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.
Shumway, R.W.
1987-10-01
The ATHENA computer program has many features that make it desirable to use as a space reactor evaluation tool. One of the missing features was a surface-to-surface thermal radiation model. A model was developed that allows any of the regular ATHENA heat slabs to radiate to any other heat slab. The view factors and surface emissivities must be specified by the user. To verify that the model was properly accounting for radiant energy transfer, two different types of test calculations were performed. Both calculations have excellent results. The updates have been used on both the INEL CDC-176 and the Livermore Cray. 7 refs., 2 figs., 6 tabs.
Radiation risk estimation models
Hoel, D.G.
1987-11-01
Cancer risk models and their relationship to ionizing radiation are discussed. There are many model assumptions and risk factors that have a large quantitative impact on the cancer risk estimates. Other health end points such as mental retardation may be an even more serious risk than cancer for those with in utero exposures. 8 references.
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
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.
NASA Astrophysics Data System (ADS)
Smirnova, Olga
Biologically motivated mathematical models, which describe the dynamics of the major hematopoietic lineages (the thrombocytopoietic, lymphocytopoietic, granulocytopoietic, and erythropoietic systems) in acutely/chronically irradiated humans are developed. These models are implemented as systems of nonlinear differential equations, which variables and constant parameters have clear biological meaning. It is shown that the developed models are capable of reproducing clinical data on the dynamics of these systems in humans exposed to acute radiation in the result of incidents and accidents, as well as in humans exposed to low-level chronic radiation. Moreover, the averaged value of the "lethal" dose rates of chronic irradiation evaluated within models of these four major hematopoietic lineages coincides with the real minimal dose rate of lethal chronic irradiation. The demonstrated ability of the models of the human thrombocytopoietic, lymphocytopoietic, granulocytopoietic, and erythropoietic systems to predict the dynamical response of these systems to acute/chronic irradiation in wide ranges of doses and dose rates implies that these mathematical models form an universal tool for the investigation and prediction of the dynamics of the major human hematopoietic lineages for a vast pattern of irradiation scenarios. In particular, these models could be applied for the radiation risk assessment for health of astronauts exposed to space radiation during long-term space missions, such as voyages to Mars or Lunar colonies, as well as for health of people exposed to acute/chronic irradiation due to environmental radiological events.
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
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.
Radiation torque exerted on a spheroid: Analytical solution
NASA Astrophysics Data System (ADS)
Xu, Feng; Lock, James A.; Gouesbet, Gérard; Tropea, Cameron
2008-07-01
As a companion work to our previous study [F. Xu, , Phys. Rev. E. 75, 026613 (2007)] on radiation force prediction for a homogeneous spheroid, we provide in this paper the analytical solution to the radiation torque exerted by an arbitrarily shaped beam on a spheroid, which can be prolate or oblate, transparent or absorbing. Calculations based upon this theoretical development are performed for both linearly and circularly polarized incident beams, and the results are compared to those of a sphere. Stable orientations of spheroids inside a linearly and a circularly polarized Gaussian beam are predicted. We analyze two physical mechanisms, the polarization torque and the reaction force torque, which do not exist or have no contribution to the torque on a sphere but cause rotation of a spheroid. As verification, the dipole method is also developed for the torque calculation for spheroids of size much less than the wavelength, and geometrical optics is developed to qualitatively analyze the torque exerted on spheroids of large size.
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.
NASA Astrophysics Data System (ADS)
Déau, Estelle
2015-06-01
The opposition effect in Saturn's main rings is characterized by a surge in ring brightness, when the phase angle approaches zero degree. This effect can be used to derive: physical properties of the ring particles and the ring layer, via the shadow hiding mechanism; and physical properties of the regolith grains that cover the ring particles, via the coherent backscattering mechanism. Since the exact origin of this effect is still a matter of debate, we try different combinations of the physical mechanisms cited above to derive constraints on the nature, the texture, and the disposition of the ring particles. In particular, we derive regolith grain sizes, particle sizes, differential power law indices, filling factors, and vertical thicknesses; and we compare them with independent works to validate or invalidate the assumptions of the opposition effect models used. Our coherent backscattering model provides grain sizes similar to the sizes estimated from water ice band depth modeling in the near infrared. Our shadow hiding model assuming a power law size distribution provides vertical thickness consistent with previous estimates from density waves measurements and N-body simulations. We show that the assumption of an homogeneous medium is a key parameter in the shadow hiding modeling. In the case of the B ring, we demonstrate that all previous photometric models assuming an homogeneous ring layer (i.e. uniform particle size distribution, random spacing of the particles and small filling factor) have led to a set of unconfirmed solutions. This result reinforces the idea that the Saturn's main rings should be modeled as an heterogeneous medium.
Status of LDEF radiation modeling
NASA Technical Reports Server (NTRS)
Watts, John W.; Armstrong, T. W.; Colborn, B. L.
1995-01-01
The current status of model prediction and comparison with LDEF radiation dosimetry measurements is summarized with emphasis on major results obtained in evaluating the uncertainties of present radiation environment model. The consistency of results and conclusions obtained from model comparison with different sets of LDEF radiation data (dose, activation, fluence, LET spectra) is discussed. Examples where LDEF radiation data and modeling results can be utilized to provide improved radiation assessments for planned LEO missions (e.g., Space Station) are given.
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.
A semi-analytic model of magnetized liner inertial fusion
McBride, Ryan D.; Slutz, Stephen A.
2015-05-15
Presented is a semi-analytic model of magnetized liner inertial fusion (MagLIF). This model accounts for several key aspects of MagLIF, including: (1) preheat of the fuel (optionally via laser absorption); (2) pulsed-power-driven liner implosion; (3) liner compressibility with an analytic equation of state, artificial viscosity, internal magnetic pressure, and ohmic heating; (4) adiabatic compression and heating of the fuel; (5) radiative losses and fuel opacity; (6) magnetic flux compression with Nernst thermoelectric losses; (7) magnetized electron and ion thermal conduction losses; (8) end losses; (9) enhanced losses due to prescribed dopant concentrations and contaminant mix; (10) deuterium-deuterium and deuterium-tritium primary fusion reactions for arbitrary deuterium to tritium fuel ratios; and (11) magnetized α-particle fuel heating. We show that this simplified model, with its transparent and accessible physics, can be used to reproduce the general 1D behavior presented throughout the original MagLIF paper [S. A. Slutz et al., Phys. Plasmas 17, 056303 (2010)]. We also discuss some important physics insights gained as a result of developing this model, such as the dependence of radiative loss rates on the radial fraction of the fuel that is preheated.
Radiative transfer model: matrix operator method.
Liu, Q; Ruprecht, E
1996-07-20
A radiative transfer model, the matrix operator method, is discussed here. The matrix operator method is applied to a plane-parallel atmosphere within three spectral ranges: the visible, the infrared, and the microwave. For a homogeneous layer with spherical scattering, the radiative transfer equation can be solved analytically. The vertically inhomogeneous atmosphere can be subdivided into a set of homogeneous layers. The solution of the radiative transfer equation for the vertically inhomogeneous atmosphere is obtained recurrently from the analytical solutions for the subdivided layers. As an example for the application of the matrix operator method, the effects of the cirrus and the stratocumulus clouds on the net radiation at the surface and at the top of the atmosphere are investigated. The relationship between the polarization in the microwave range and the rain rates is also studied. Copies of the FORTRAN program and the documentation of the FORTRAN program on a diskette are available. PMID:21102832
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)
Beth, A.; Garnier, P.; Toublanc, D.; Dandouras, I. S.; Mazelle, C. X.
2013-12-01
a) The planetary exospheres are poorly known in their outer parts, since the neutral densities are low compared with the instruments detection capabilities. The exospheric models are thus often the main source of information at such high altitudes. We revisit here the importance of a specific exospheric population, i.e. the satellite particles, which is usually neglected in the models (see also Beth et al., Icarus, accepted). These particles are indeed produced through rare collisions in the exospheres, and may either be negligible or dominate the exospheres of all planets with dense atmospheres in our solar system, depending on the balance between their sources and losses. Richter et al. (1979) were the first to propose, beyond the Chamberlain's (Chamberlain (1963)) rough approximation, a rigorous approach for these particles by using the Boltzmann equation in the Earth exosphere below 3000 km altitude. We here further investigate this approach and determine the contribution of satellite populations to the densities of light elements at Titan (H2 species) and Mars (H species). The results confirm that the Chamberlain approximation overestimates the satellite particles densities at high altitudes, but that there may be enough collisions to produce a significant amount of satellite particles in some conditions, up to more than 50% of the contributions due to ballistic and escaping populations (i.e. those considered in the collisionless models) in the case of Mars or even 36% at Titan. This suggests that considering collisionless exospheric profiles for light species can lead to an underestimation of the total densities at high altitudes. b) We also present a new analytical approach to understand the structure of the exospheres submitted to the radiation pressure and determine the densities of others populations, i.e. ballistic and escaping particles, using the Liouville theorem and Hamiltonian mechanics. Our modeling work will in particular enable to better
NASA Technical Reports Server (NTRS)
Horwitz, James L.
1992-01-01
The purpose of this work was to assist with the development of analytical techniques for the interpretation of infrared observations. We have done the following: (1) helped to develop models for continuum absorption calculations for water vapor in the far infrared spectral region; (2) worked on models for pressure-induced absorption for O2 and N2 and their comparison with available observations; and (3) developed preliminary studies of non-local thermal equilibrium effects in the upper stratosphere and mesosphere for infrared gases. These new techniques were employed for analysis of balloon-borne far infrared data by a group at the Harvard-Smithsonian Center for Astrophysics. The empirical continuum absorption model for water vapor in the far infrared spectral region and the pressure-induced N2 absorption model were found to give satisfactory results in the retrieval of the mixing ratios of a number of stratospheric trace constituents from balloon-borne far infrared observations.
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.
A rapid radiative transfer model for reflection of solar radiation
NASA Technical Reports Server (NTRS)
Xiang, X.; Smith, E. A.; Justus, C. G.
1994-01-01
A rapid analytical radiative transfer model for reflection of solar radiation in plane-parallel atmospheres is developed based on the Sobolev approach and the delta function transformation technique. A distinct advantage of this model over alternative two-stream solutions is that in addition to yielding the irradiance components, which turn out to be mathematically equivalent to the delta-Eddington approximation, the radiance field can also be expanded in a mathematically consistent fashion. Tests with the model against a more precise multistream discrete ordinate model over a wide range of input parameters demonstrate that the new approximate method typically produces average radiance differences of less than 5%, with worst average differences of approximately 10%-15%. By the same token, the computational speed of the new model is some tens to thousands times faster than that of the more precise model when its stream resolution is set to generate precise calculations.
r-mode instability: Analytical solution with gravitational radiation reaction
Dias, Oscar J.C.; Sa, Paulo M.
2005-07-15
Analytical r-mode solutions are investigated within the linearized theory in the case of a slowly rotating, Newtonian, barotropic, nonmagnetized, perfect-fluid star in which the gravitational radiation (GR) reaction force is present. For the GR reaction term we use the 3.5 post-Newtonian order expansion of the GR force, in order to include the contribution of the current quadrupole moment. We find the explicit expression for the r-mode velocity perturbations and we conclude that they are sinusoidal with the same frequency as the well-known GR force-free linear r-mode solution, and that the GR force drives the r-modes unstable with a growth time scale that agrees with the expression first found by Lindblom, Owen, and Morsink. We also show that the amplitude of these velocity perturbations is corrected, relatively to the GR force-free case, by a term of order {omega}{sup 6}, where {omega} is the angular velocity of the star.
NASA Technical Reports Server (NTRS)
Jin, Zhonghai; Charlock, Thomas P.; Rutledge, Ken; Knut Stamnes; Wang, Yingjian
2006-01-01
Using the efficient discrete-ordinate method, we present an analytical solution for radiative transfer in the coupled atmosphere-ocean system with rough air-water interface. The theoretical formulations of the radiative transfer equation and solution are described. The effects of surface roughness on radiation field in the atmosphere and ocean are studied and compared with measurements. The results show that ocean surface roughness has significant effects on the upwelling radiation in the atmosphere and the downwelling radiation in the ocean. As wind speed increases, the angular domain of sunglint broadens, the surface albedo decreases, and the transmission to ocean increases. The downward radiance field in the upper ocean is highly anisotropic, but this anisotropy decreases rapidly as surface wind increases and as depth in ocean increases. The effects of surface roughness on radiation also depend greatly on both wavelength and angle of incidence (i.e., solar elevation); these effects are significantly smaller throughout the spectrum at high sun. The model-observation discrepancies may indicate that the Cox-Munk surface roughness model is not sufficient for high wind conditions.
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.
Analytical investigation of adaptive control of radiated inlet noise from turbofan engines
NASA Technical Reports Server (NTRS)
Risi, John D.; Burdisso, Ricardo A.
1994-01-01
An analytical model has been developed to predict the resulting far field radiation from a turbofan engine inlet. A feedforward control algorithm was simulated to predict the controlled far field radiation from the destructive combination of fan noise and secondary control sources. Numerical results were developed for two system configurations, with the resulting controlled far field radiation patterns showing varying degrees of attenuation and spillover. With one axial station of twelve control sources and error sensors with equal relative angular positions, nearly global attenuation is achieved. Shifting the angular position of one error sensor resulted in an increase of spillover to the extreme sidelines. The complex control inputs for each configuration was investigated to identify the structure of the wave pattern created by the control sources, giving an indication of performance of the system configuration. It is deduced that the locations of the error sensors and the control source configuration are equally critical to the operation of the active noise control system.
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.
A physically based analytical spatial air temperature and humidity model
NASA Astrophysics Data System (ADS)
Yang, Yang; Endreny, Theodore A.; Nowak, David J.
2013-09-01
Spatial variation of urban surface air temperature and humidity influences human thermal comfort, the settling rate of atmospheric pollutants, and plant physiology and growth. Given the lack of observations, we developed a Physically based Analytical Spatial Air Temperature and Humidity (PASATH) model. The PASATH model calculates spatial solar radiation and heat storage based on semiempirical functions and generates spatially distributed estimates based on inputs of topography, land cover, and the weather data measured at a reference site. The model assumes that for all grids under the same mesoscale climate, grid air temperature and humidity are modified by local variation in absorbed solar radiation and the partitioning of sensible and latent heat. The model uses a reference grid site for time series meteorological data and the air temperature and humidity of any other grid can be obtained by solving the heat flux network equations. PASATH was coupled with the USDA iTree-Hydro water balance model to obtain evapotranspiration terms and run from 20 to 29 August 2010 at a 360 m by 360 m grid scale and hourly time step across a 285 km2 watershed including the urban area of Syracuse, NY. PASATH predictions were tested at nine urban weather stations representing variability in urban topography and land cover. The PASATH model predictive efficiency R2 ranged from 0.81 to 0.99 for air temperature and 0.77 to 0.97 for dew point temperature. PASATH is expected to have broad applications on environmental and ecological models.
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).
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.
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.
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.
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.
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.
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.
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...
Chandra Radiation Environment Modeling
NASA Technical Reports Server (NTRS)
Minow, Joseph I.; Blackwell, W. C.
2003-01-01
CRMFLX (Chandra Radiation Model of ion FluX) is a radiation environment risk mitigation tool for use as a decision aid in planning the operations times for Chandra's Advanced CCD Imaging Spectrometer (ACIS) detector. The accurate prediction of the proton flux environment with energies of 100 - 200 keV is needed in order to protect the ACIS detector against proton degradation. Unfortunately, protons of this energy are abundant in the region of space Chandra must operate, and on-board particle detectors do not measure proton flux levels of the required energy range. This presentation will describe the plasma environment data analysis and modeling basis of the CRMFLX engineering environment model developed to predict the proton flux in the solar wind, magnetosheath, and magnetosphere phenomenological regions of geospace. The recently released CRMFLX Version 2 implementation includes an algorithm that propagates flux from an observation location to other regions of the magnetosphere based on convective ExB and VB-curvature particle drift motions. This technique has the advantage of more completely filling out the database and makes maximum use of limited data obtained during high Kp periods or in areas of the magnetosphere with poor satellite flux measurement coverage.
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)
NASA Astrophysics Data System (ADS)
McMahon, Jay W.
Solar radiation pressure is the largest non-gravitational perturbation for most satellites in the solar system, and can therefore have a significant influence on their orbital dynamics. This work presents a new method for representing the solar radiation pressure force acting on a satellite, and applies this theory to natural and artificial satellites. The solar radiation pressure acceleration is modeled as a Fourier series which depends on the Sun's location in a body-fixed frame; a new set of Fourier coefficients are derived for every latitude of the Sun in this frame, and the series is expanded in terms of the longitude of the Sun. The secular effects due to the solar radiation pressure perturbations are given analytically through the application of averaging theory when the satellite is in a synchronous orbit. This theory is then applied to binary asteroid systems to explain the Binary YORP effect. Long term predictions of the evolution of the near-Earth asteroid 1999 KW4 are discussed under the influence of solar radiation pressure, J2, and 3rd body gravitational effects from the Sun. Secular effects are shown to remain when the secondary asteroid becomes non-synchronous due to a librational motion. The theory is also applied to Earth orbiting spacecraft, and is shown to be a valuable tool for improved orbit determination. The Fourier series solar radiation pressure model derived here is shown to give comparable results for orbit determination of the GPS IIR-M satellites as JPL's solar radiation pressure model. The theory is also extended to incorporate the effects of the Earth's shadow analytically. This theory is briefly applied to the evolution of orbital debris to explain the assumptions that are necessary in order to use the cannonball model for debris orbit evolution, as is common in the literature. Finally, the averaging theory methodology is applied to a class of Earth orbiting solar sail spacecraft to show the orbital effects when the sails are made
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.
Quasi-analytical treatment of spatially averaged radiation transfer in complex terrain
NASA Astrophysics Data System (ADS)
LöWe, H.; Helbig, N.
2012-10-01
We provide a new quasi-analytical method to compute the subgrid topographic influences on the shortwave radiation fluxes and the effective albedo in complex terrain as required for large-scale meteorological, land surface, or climate models. We investigate radiative transfer in complex terrain via the radiosity equation on isotropic Gaussian random fields. Under controlled approximations we derive expressions for domain-averaged fluxes of direct, diffuse, and terrain radiation and the sky view factor. Domain-averaged quantities can be related to a type of level-crossing probability of the random field, which is approximated by long-standing results developed for acoustic scattering at ocean boundaries. This allows us to express all nonlocal horizon effects in terms of a local terrain parameter, namely, the mean-square slope. Emerging integrals are computed numerically, and fit formulas are given for practical purposes. As an implication of our approach, we provide an expression for the effective albedo of complex terrain in terms of the Sun elevation angle, mean-square slope, the area-averaged surface albedo, and the ratio of atmospheric direct beam to diffuse radiation. For demonstration we compute the decrease of the effective albedo relative to the area-averaged albedo in Switzerland for idealized snow-covered and clear-sky conditions at noon in winter. We find an average decrease of 5.8% and spatial patterns which originate from characteristics of the underlying relief. Limitations and possible generalizations of the method are discussed.
Analytic Shielding Optimization to Reduce Crew Exposure to Ionizing Radiation Inside Space Vehicles
NASA Technical Reports Server (NTRS)
Gaza, Razvan; Cooper, Tim P.; Hanzo, Arthur; Hussein, Hesham; Jarvis, Kandy S.; Kimble, Ryan; Lee, Kerry T.; Patel, Chirag; Reddell, Brandon D.; Stoffle, Nicholas; Zapp, E. Neal; Shelfer, Tad D.
2009-01-01
A sustainable lunar architecture provides capabilities for leveraging out-of-service components for alternate uses. Discarded architecture elements may be used to provide ionizing radiation shielding to the crew habitat in case of a Solar Particle Event. The specific location relative to the vehicle where the additional shielding mass is placed, as corroborated with particularities of the vehicle design, has a large influence on protection gain. This effect is caused by the exponential- like decrease of radiation exposure with shielding mass thickness, which in turn determines that the most benefit from a given amount of shielding mass is obtained by placing it so that it preferentially augments protection in under-shielded areas of the vehicle exposed to the radiation environment. A novel analytic technique to derive an optimal shielding configuration was developed by Lockheed Martin during Design Analysis Cycle 3 (DAC-3) of the Orion Crew Exploration Vehicle (CEV). [1] Based on a detailed Computer Aided Design (CAD) model of the vehicle including a specific crew positioning scenario, a set of under-shielded vehicle regions can be identified as candidates for placement of additional shielding. Analytic tools are available to allow capturing an idealized supplemental shielding distribution in the CAD environment, which in turn is used as a reference for deriving a realistic shielding configuration from available vehicle components. While the analysis referenced in this communication applies particularly to the Orion vehicle, the general method can be applied to a large range of space exploration vehicles, including but not limited to lunar and Mars architecture components. In addition, the method can be immediately applied for optimization of radiation shielding provided to sensitive electronic components.
NASA Astrophysics Data System (ADS)
Beth, A.; Garnier, P.; Toublanc, D.; Dandouras, I. S.; Mazelle, C. X.
2014-12-01
Because of rare collisions, the motion of light species (H, H2) in the planetary exospheres is essentially determined by the external forces: the gravitation from the planet and the radiation pressure, ... Currently, the only analytical model used to model exospheric neutral density profiles is the well-known Chamberlain model which takes into account only the gravity. In this work and in the same way as Chamberlain, we solve rigorously and analytically, based on the Hamiltonian mechanics and Liouville theorem, the additional effect of the radiation pressure in particular for hydrogen (the model works for any species sensitive to the radiation pressure) on the structure of the exosphere and on the density profiles of ballistic particles. This approach was initially developed by Bishop and Chamberlain (1989) only in the Sun-planet direction. We extend it here to the whole exosphere with a 2D model. Also, we determine analytically the escape flux on the dayside at SZA=0, which can be compared with the Jeans' escape flux. We thus show that the radiation pressure induces : strong density asymmetries at high altitudes in the planetary exospheres, leading to the phenomenon of geotail at Earth for example the natural existence of an external limit (or exopause) for the exosphere, whose location is analytically determined an increase of the exospheric densities compared with Chamberlain profiles without radiation pressure (e.g. up to +150% at 5 Martian radius) a significant increase of the thermal escape flux (up to 30/35% for Earth/Mars today), until a «blow-off » regime with a constant escape flux for an extreme radiation pressure. The influence of the radiation pressure on the escape flux may thus bring conditions on the size of primary atmospheres, because of a strong radiation pressure in the Sun's young years. Finally, we show that this model may be applied to exoplanets, in particular to the hot Jupiters that are also subject to additional effects: centrifugal
Analytical Modeling of Variable Density Multilayer Insulation for Cryogenic Storage
NASA Technical Reports Server (NTRS)
Hedayat, A.; Hastings, L. J.; Brown, T.; Cruit, Wendy (Technical Monitor)
2001-01-01
A unique foam/Multilayer Insulation (MLI) combination concept for orbital cryogenic storage was experimentally evaluated at NASA Marshall Space Flight Center (MSFC) using the Multipurpose Hydrogen Test Bed (MHTB). The MLI was designed for an on-orbit storage period of 45 days and included several unique features such as: a variable layer density and larger but fewer perforations for venting during ascent to orbit. Test results with liquid hydrogen indicated that the MLI weight or heat leak is reduced by about half in comparison with standard MLI. The focus of this paper is on analytical modeling of the Variable Density MLI (VD-MLI) on-orbit performance (i.e. vacuum/low pressure environment). The foam/VD-MLI combination model is considered to have five segments. The first segment represents the optional foam layer. The second, third, and fourth segments represent three MLI segments with different layer densities. The last segment is considered to be a shroud that surrounds the last MLI layer. Two approaches are considered. In the first approach, the variable density MLI is modeled layer by layer while in the second approach, a semi-empirical model is applied. Both models account for thermal radiation between shields, gas conduction, and solid conduction through the layer separator materials.
Modifiers of radiation action on DNA screened by analytical ultracentrifugation
Cobreros, G.; Lopez Zumel, M.C.; Usobiaga, P.
1982-11-01
The effect of daunomycin, chromomycin A/sub 3/, anthramycin, cis-dichlorodiammineplatinum(II) (cisplatin), mitomycin C, and chloroquine on the frequency of radioninduced strand breaks in X-irradiated aqueous solutions of DNA was studied primarily by analytical ultracentrifugation using the sedimentation velocity method. The results show a potent radiosensitizing effect for daunomycin and chromomycin, a protective action for chloroquine and mitomycin, and a nonmodifying effect for anthramycin. Cisplatin forms a highly aggregated complex with DNA which prevents this kind of study.
Early test facilities and analytic methods for radiation shielding: Proceedings
Ingersoll, D.T. ); Ingersoll, J.K. )
1992-11-01
This report represents a compilation of eight papers presented at the 1992 American Nuclear Society/European Nuclear Society International Meeting. The meeting is of special significance since it commemorates the fiftieth anniversary of the first controlled nuclear chain reaction. The papers contained in this report were presented in a special session organized by the Radiation Protection and Shielding Division in keeping with the historical theme of the meeting. The paper titles are good indicators of their content and are: (1) The origin of radiation shielding research: The Oak Ridge experience, (2) Shielding research at the hanford site, (3) Aircraft shielding experiments at General Dynamics Fort Worth, 1950-1962, (4) Where have the neutrons gone , a history of the tower shielding facility, (5) History and evolution of buildup factors, (6) Early shielding research at Bettis atomic power laboratory, (7) UK reactor shielding: then and now, (8) A very personal view of the development of radiation shielding theory.
Analytical Investigation of Pumped Fluid Loop Radiators for Orion Spacecraft
NASA Technical Reports Server (NTRS)
Reavis, Gretchen
2007-01-01
This viewgraph presentation reviews the history of pumped fluid loop radiators used in Apollo spacecraft, and the problems and challenges for using them in the Orion Spacecraft. Included in this presentation are the issues of Flow stagnation, flow stability, for single panels and multi-panels.
Analytical models for total dose ionization effects in MOS devices.
Campbell, Phillip Montgomery; Bogdan, Carolyn W.
2008-08-01
MOS devices are susceptible to damage by ionizing radiation due to charge buildup in gate, field and SOI buried oxides. Under positive bias holes created in the gate oxide will transport to the Si / SiO{sub 2} interface creating oxide-trapped charge. As a result of hole transport and trapping, hydrogen is liberated in the oxide which can create interface-trapped charge. The trapped charge will affect the threshold voltage and degrade the channel mobility. Neutralization of oxidetrapped charge by electron tunneling from the silicon and by thermal emission can take place over long periods of time. Neutralization of interface-trapped charge is not observed at room temperature. Analytical models are developed that account for the principal effects of total dose in MOS devices under different gate bias. The intent is to obtain closed-form solutions that can be used in circuit simulation. Expressions are derived for the aging effects of very low dose rate radiation over long time periods.
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; ...
Characterization of uniform scanning proton beams with analytical models
NASA Astrophysics Data System (ADS)
Demez, Nebi
Tissue equivalent phantoms have an important place in radiation therapy planning and delivery. They have been manufactured for use in conventional radiotherapy. Their tissue equivalency for proton beams is currently in active investigation. The Bragg-Kleeman rule was used to calculate water equivalent thickness (WET) for available tissue equivalent phantoms from CIRS (Norfolk, VA, USA). WET's of those phantoms were also measured using proton beams at Hampton University Proton Therapy Institute (HUPTI). WET measurements and calculations are in good agreement within ˜1% accuracy except for high Z phantoms. Proton beams were also characterized with an analytical proton dose calculation model, Proton Loss Model (PLM) [26], to investigate protons interactions in water and those phantoms. Depth-dose and lateral dose profiles of protons in water and in those phantoms were calculated, measured, and compared. Water Equivalent Spreadness (WES) was also investigated for those phantoms using the formula for scattering power ratio. Because WES is independent of incident energy of protons, it is possible to estimate spreadness of protons in different media by just knowing WES. Measurements are usually taken for configuration of the treatment planning system (TPS). This study attempted to achieve commissioning data for uniform scanning proton planning with analytical methods, PLM, which have been verified with published measurements and Monte Carlo calculations. Depth doses and lateral profiles calculated by PLM were compared with measurements via the gamma analysis method. While gamma analysis shows that depth doses are in >90% agreement with measured depth doses, the agreement falls to <80% for some lateral profiles. PLM data were imported into the TPS (PLM-TPS). PLM-TPS was tested with different patient cases. The PLM-TPS treatment plans for 5 prostate cases show acceptable agreement. The Planning Treatment Volume (PTV) coverage was 100 % with PLM-TPS except for one case in
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.
NASA Astrophysics Data System (ADS)
Bars, Itzhak; Chen, Shih-Hung; Steinhardt, Paul J.; Turok, Neil
2012-10-01
We study a model of a scalar field minimally coupled to gravity, with a specific potential energy for the scalar field, and include curvature and radiation as two additional parameters. Our goal is to obtain analytically the complete set of configurations of a homogeneous and isotropic universe as a function of time. This leads to a geodesically complete description of the Universe, including the passage through the cosmological singularities, at the classical level. We give all the solutions analytically without any restrictions on the parameter space of the model or initial values of the fields. We find that for generic solutions the Universe goes through a singular (zero-size) bounce by entering a period of antigravity at each big crunch and exiting from it at the following big bang. This happens cyclically again and again without violating the null-energy condition. There is a special subset of geodesically complete nongeneric solutions which perform zero-size bounces without ever entering the antigravity regime in all cycles. For these, initial values of the fields are synchronized and quantized but the parameters of the model are not restricted. There is also a subset of spatial curvature-induced solutions that have finite-size bounces in the gravity regime and never enter the antigravity phase. These exist only within a small continuous domain of parameter space without fine-tuning the initial conditions. To obtain these results, we identified 25 regions of a 6-parameter space in which the complete set of analytic solutions are explicitly obtained.
Quasi-analytical treatment of spatially averaged radiation transfer in complex terrain
NASA Astrophysics Data System (ADS)
Löwe, H.; Helbig, N.
2012-04-01
We provide a new quasi-analytical method to compute the topographic influence on the effective albedo of complex topography as required for meteorological, land-surface or climate models. We investigate radiative transfer in complex terrain via the radiosity equation on isotropic Gaussian random fields. Under controlled approximations we derive expressions for domain averages of direct, diffuse and terrain radiation and the sky view factor. Domain averaged quantities are related to a type of level-crossing probability of the random field which is approximated by longstanding results developed for acoustic scattering at ocean boundaries. This allows us to express all non-local horizon effects in terms of a local terrain parameter, namely the mean squared slope. Emerging integrals are computed numerically and fit formulas are given for practical purposes. As an implication of our approach we provide an expression for the effective albedo of complex terrain in terms of the sun elevation angle, mean squared slope, the area averaged surface albedo, and the direct-to-diffuse ratio of solar radiation. As an application, we compute the effective albedo for the Swiss Alps and discuss possible generalizations of the method.
Analytical Radiation Transport Benchmarks for The Next Century
B.D. Ganapol
2005-01-19
Verification of large-scale computational algorithms used in nuclear engineering and radiological applications is an essential element of reliable code performance. For this reason, the development of a suite of multidimensional semi-analytical benchmarks has been undertaken to provide independent verification of proper operation of codes dealing with the transport of neutral particles. The benchmarks considered cover several one-dimensional, multidimensional, monoenergetic and multigroup, fixed source and critical transport scenarios. The first approach, called the Green's Function. In slab geometry, the Green's function is incorporated into a set of integral equations for the boundary fluxes. Through a numerical Fourier transform inversion and subsequent matrix inversion for the boundary fluxes, a semi-analytical benchmark emerges. Multidimensional solutions in a variety of infinite media are also based on the slab Green's function. In a second approach, a new converged SN method is developed. In this method, the SN solution is ''minded'' to bring out hidden high quality solutions. For this case multigroup fixed source and criticality transport problems are considered. Remarkably accurate solutions can be obtained with this new method called the Multigroup Converged SN (MGCSN) method as will be demonstrated.
Putting synchrotron radiation to work for technology: Analytic methods
Not Available
1992-02-01
This report contains viewgraphs on: Advanced Light Source; Ultra-ESCA: Advanced Capabilities of XPS with High-Brightness Synchrotron Radiation; High-Resolution (20 nm) XPS and XANES with the ALS; Photoelectron Spectroscopy in Industry: Current Capabilities, Needs, and Possible Roles for the ALS; Materials Analysis by Photoemission: Is This Practical at ALS?; Applications of Long-Wavelength X-Ray Fluorescence Spectrometry and X-Ray Powder Diffractometry.
Putting synchrotron radiation to work for technology: Analytic methods
Not Available
1992-02-01
This report contains viewgraphs on: Advanced Light Source; Ultra-ESCA: Advanced Capabilities of XPS with High-Brightness Synchrotron Radiation; High-Resolution (20 nm) XPS and XANES with the ALS; Photoelectron Spectroscopy in Industry: Current Capabilities, Needs, and Possible Roles for the ALS; Materials Analysis by Photoemission: Is This Practical at ALS ; Applications of Long-Wavelength X-Ray Fluorescence Spectrometry and X-Ray Powder Diffractometry.
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.
Econometric model for age- and population-dependent radiation exposures
Sandquist, G.M.; Slaughter, D.M. ); Rogers, V.C.
1991-01-01
The economic impact associated with ionizing radiation exposures in a given human population depends on numerous factors including the individual's mean economic status as a function age, the age distribution of the population, the future life expectancy at each age, and the latency period for the occurrence of radiation-induced health effects. A simple mathematical model has been developed that provides an analytical methodology for estimating the societal econometrics associated with radiation effects are to be assessed and compared for economic evaluation.
Analytical thermal model validation for Cassini radioisotope thermoelectric generator
Lin, E.I.
1997-12-31
The Saturn-bound Cassini spacecraft is designed to rely, without precedent, on the waste heat from its three radioisotope thermoelectric generators (RTGs) to warm the propulsion module subsystem, and the RTG end dome temperature is a key determining factor of the amount of waste heat delivered. A previously validated SINDA thermal model of the RTG was the sole guide to understanding its complex thermal behavior, but displayed large discrepancies against some initial thermal development test data. A careful revalidation effort led to significant modifications and adjustments of the model, which result in a doubling of the radiative heat transfer from the heat source support assemblies to the end domes and bring up the end dome and flange temperature predictions to within 2 C of the pertinent test data. The increased inboard end dome temperature has a considerable impact on thermal control of the spacecraft central body. The validation process offers an example of physically-driven analytical model calibration with test data from not only an electrical simulator but also a nuclear-fueled flight unit, and has established the end dome temperatures of a flight RTG where no in-flight or ground-test data existed before.
Analytical discrete ordinate method for radiative transfer in dense vegetation canopies
NASA Astrophysics Data System (ADS)
Picca, Paolo; Furfaro, Roberto
2013-03-01
The radiative transfer (RT) in dense vegetation canopies can be approximated via linear Boltzmann equation. However, the directionality of the basic scattering element (i.e. the canopy leaf) makes the medium inherently anisotropic and introduces special features in the definition of both scattering kernel and total cross section. In this paper, a classical methodology for the solution of transport problem, namely the analytical discrete ordinate (ADO) method, is extended to account for the peculiarities of photon transport into dense vegetation canopies. It is demonstrated that the special symmetries arising from modeling the leaf as a bi-Lambertian scatterer, enable the derivation of the ADO equations for canopy transport. Several numerical tests have been performed to evaluate the accuracy of ADO against numerical benchmarks available in the literature. The results show that the proposed methodology is highly accurate, computationally efficient and may set future standards for numerical transport in dense vegetation canopies.
Modeling the Space Radiation Environment
NASA Technical Reports Server (NTRS)
Xapsos, Michael A.
2006-01-01
There has been a renaissance of interest in space radiation environment modeling. This has been fueled by the growing need to replace long time standard AP-9 and AE-8 trapped particle models, the interplanetary exploration initiative, the modern satellite instrumentation that has led to unprecedented measurement accuracy, and the pervasive use of Commercial off the Shelf (COTS) microelectronics that require more accurate predictive capabilities. The objective of this viewgraph presentation was to provide basic understanding of the components of the space radiation environment and their variations, review traditional radiation effects application models, and present recent developments.
Early Test Facilities and Analytic Methods for Radiation Shielding
Ingersoll, D.T.
1992-01-01
This report represents a compilation of eight papers presented at the 1992 American Nuclear Society/European Nuclear Society International Meeting held in Chicago, Illinois on November 15 20,1992. The meeting is of special significance since it commemorates the 50th anniversary of the first controlled nuclear chain reaction, which occurred, not coincidentally, in Chicago. The papers contained in this report were presented in a special session organized by the Radiation Protection and Shielding Division in keeping with the historical theme of the meeting.
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
Rabacus: A Python package for analytic cosmological radiative transfer calculations
NASA Astrophysics Data System (ADS)
Altay, G.; Wise, J. H.
2015-04-01
We describe RABACUS, a Python package for calculating the transfer of hydrogen ionizing radiation in simplified geometries relevant to astronomy and cosmology. We present example solutions for three specific cases: (1) a semi-infinite slab gas distribution in a homogeneous isotropic background, (2) a spherically symmetric gas distribution with a point source at the center, and (3) a spherically symmetric gas distribution in a homogeneous isotropic background. All problems can accommodate arbitrary spectra and density profiles as input. The solutions include a treatment of both hydrogen and helium, a self-consistent calculation of equilibrium temperatures, and the transfer of recombination radiation. The core routines are written in Fortran 90 and then wrapped in Python leading to execution speeds thousands of times faster than equivalent routines written in pure Python. In addition, all variables have associated units for ease of analysis. The software is part of the Python Package Index and the source code is available on Bitbucket at https://bitbucket.org/galtay/rabacus. In addition, installation instructions and a detailed users guide are available at http://pythonhosted.org//rabacus.
Analytical solutions for radiation-driven winds in massive stars. I. The fast regime
Araya, I.; Curé, M.; Cidale, L. S.
2014-11-01
Accurate mass-loss rate estimates are crucial keys in the study of wind properties of massive stars and for testing different evolutionary scenarios. From a theoretical point of view, this implies solving a complex set of differential equations in which the radiation field and the hydrodynamics are strongly coupled. The use of an analytical expression to represent the radiation force and the solution of the equation of motion has many advantages over numerical integrations. Therefore, in this work, we present an analytical expression as a solution of the equation of motion for radiation-driven winds in terms of the force multiplier parameters. This analytical expression is obtained by employing the line acceleration expression given by Villata and the methodology proposed by Müller and Vink. On the other hand, we find useful relationships to determine the parameters for the line acceleration given by Müller and Vink in terms of the force multiplier parameters.
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.
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.
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.
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.
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.
Linear radiation model for phase of thermal emission spectroscopy
NASA Astrophysics Data System (ADS)
Bennett, Ted D.; Yu, Fengling
2005-11-01
A linear radiation model is developed that overcomes the analytical complexity in phase of thermal emission spectroscopy. It is shown that the linear radiation model can result in a simple algebraic relation between the phase of thermal emission and four coating properties, enabling these properties to be determined by nonlinear regression analysis of experimental measurements. Suitability of the linear radiation model to various measurement conditions is explored, and the model is applied to the phase of thermal emission measurements performed on a thermal barrier coating.
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.
A multigroup radiation diffusion test problem: Comparison of code results with analytic solution
Shestakov, A I; Harte, J A; Bolstad, J H; Offner, S R
2006-12-21
We consider a 1D, slab-symmetric test problem for the multigroup radiation diffusion and matter energy balance equations. The test simulates diffusion of energy from a hot central region. Opacities vary with the cube of the frequency and radiation emission is given by a Wien spectrum. We compare results from two LLNL codes, Raptor and Lasnex, with tabular data that define the analytic solution.
Acoustic black holes: massless scalar field analytic solutions and analogue Hawking radiation
NASA Astrophysics Data System (ADS)
Vieira, H. S.; Bezerra, V. B.
2016-07-01
We obtain the analytic solutions of the radial part of the massless Klein-Gordon equation in the spacetime of both three dimensional rotating and four dimensional canonical acoustic black holes, which are given in terms of the confluent Heun functions. From these solutions, we obtain the scalar waves near the acoustic horizon. We discuss the analogue Hawking radiation of massless scalar particles and the features of the spectrum associated with the radiation emitted by these acoustic black holes.
Comparison of a semi-analytic and a CFD model uranium combustion to experimental data.
Clarksean, R.
1998-04-01
Two numerical models were developed and compared for the analysis of uranium combustion and ignition in a furnace. Both a semi-analytical solution and a computational fluid dynamics (CFD) numerical solution were obtained. Prediction of uranium oxidation rates is important for fuel storage applications, fuel processing, and the development of spent fuel metal waste forms. The semi-analytical model was based on heat transfer correlations, a semi-analytical model of flow over a flat surface, and simple radiative heat transfer from the material surface. The CFD model numerically determined the flowfield over the object of interest, calculated the heat and mass transfer to the material of interest, and calculated the radiative heat exchange of the material with the furnace. The semi-analytical model is much less detailed than the CFD model, but yields reasonable results and assists in understanding the physical process. Short computation times allowed the analyst to study numerous scenarios. The CFD model had significantly longer run times, was found to have some physical limitations that were not easily modified, but was better able to yield details of the heat and mass transfer and flow field once code limitations were overcome.
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.
NASA Astrophysics Data System (ADS)
Boss, Alan P.
2009-03-01
The disk instability mechanism for giant planet formation is based on the formation of clumps in a marginally gravitationally unstable protoplanetary disk, which must lose thermal energy through a combination of convection and radiative cooling if they are to survive and contract to become giant protoplanets. While there is good observational support for forming at least some giant planets by disk instability, the mechanism has become theoretically contentious, with different three-dimensional radiative hydrodynamics codes often yielding different results. Rigorous code testing is required to make further progress. Here we present two new analytical solutions for radiative transfer in spherical coordinates, suitable for testing the code employed in all of the Boss disk instability calculations. The testing shows that the Boss code radiative transfer routines do an excellent job of relaxing to and maintaining the analytical results for the radial temperature and radiative flux profiles for a spherical cloud with high or moderate optical depths, including the transition from optically thick to optically thin regions. These radial test results are independent of whether the Eddington approximation, diffusion approximation, or flux-limited diffusion approximation routines are employed. The Boss code does an equally excellent job of relaxing to and maintaining the analytical results for the vertical (θ) temperature and radiative flux profiles for a disk with a height proportional to the radial distance. These tests strongly support the disk instability mechanism for forming giant planets.
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.
Several numerical and analytical solutions of the radiative transfer equation (RTE) for plane albedo were compared for solar light reflection by sea water. The study incorporated the simplest case, that being a semi-infinite one-dimensional plane-parallel absorbing and scattering...
Pérez-Andújar, Angélica; Zhang, Rui; Newhauser, Wayne
2013-12-15
Purpose: Stray neutron radiation is of concern after radiation therapy, especially in children, because of the high risk it might carry for secondary cancers. Several previous studies predicted the stray neutron exposure from proton therapy, mostly using Monte Carlo simulations. Promising attempts to develop analytical models have also been reported, but these were limited to only a few proton beam energies. The purpose of this study was to develop an analytical model to predict leakage neutron equivalent dose from passively scattered proton beams in the 100-250-MeV interval.Methods: To develop and validate the analytical model, the authors used values of equivalent dose per therapeutic absorbed dose (H/D) predicted with Monte Carlo simulations. The authors also characterized the behavior of the mean neutron radiation-weighting factor, w{sub R}, as a function of depth in a water phantom and distance from the beam central axis.Results: The simulated and analytical predictions agreed well. On average, the percentage difference between the analytical model and the Monte Carlo simulations was 10% for the energies and positions studied. The authors found that w{sub R} was highest at the shallowest depth and decreased with depth until around 10 cm, where it started to increase slowly with depth. This was consistent among all energies.Conclusion: Simple analytical methods are promising alternatives to complex and slow Monte Carlo simulations to predict H/D values. The authors' results also provide improved understanding of the behavior of w{sub R} which strongly depends on depth, but is nearly independent of lateral distance from the beam central axis.
Simulation and analytic analysis of radiation driven islands at the density limit
NASA Astrophysics Data System (ADS)
Brennan, D. P.; Liu, C.; Gates, D. A.; Delgado-Aparicio, L.; White, R.
2014-10-01
The effect of radiative cooling on the onset and evolution of magnetic islands is investigated with nonlinear resistive MHD simulations and reduced theoretical analysis. The configuration is a cylindrical tokamak with a m/n = 2/1 island and includes three dimensional resistivity and anisotropic heat conduction in the simulations. The radiative cooling is implemented as a temperature perturbation inside the island, which modifies the island structure and drives the island more unstable. Analytic reduction of the saturated island size and structure supports the simulation results. The results offer intuitive understanding of experimental observations of radiation driven magnetic islands, which may explain density limit disruptions.
Diffusion model for lightning radiative transfer
NASA Technical Reports Server (NTRS)
Koshak, William J.; Solakiewicz, Richard J.; Phanord, Dieudonne D.; Blakeslee, Richard J.
1994-01-01
A one-speed Boltzmann transport theory, with diffusion approximations, is applied to study the radiative transfer properties of lightning in optically thick thunderclouds. Near-infrared (lambda = 0.7774 micrometers) photons associated with a prominent oxygen emission triplet in the lightning spectrum are considered. Transient and spatially complex lightning radiation sources are placed inside a rectangular parallelepiped thundercloud geometry and the effects of multiple scattering are studied. The cloud is assumed to be composed of a homogeneous collection of identical spherical water droplets, each droplet a nearly conservative, anisotropic scatterer. Conceptually, we treat the thundercloud like a nuclear reactor, with photons replaced by neutrons, and utilize standard one-speed neutron diffusion techniques common in nuclear reactor analyses. Valid analytic results for the intensity distribution (expanded in spherical harmonics) are obtained for regions sufficiently far from sources. Model estimates of the arrival-time delay and pulse width broadening of lightning signals radiated from within the cloud are determined and the results are in good agreement with both experimental data and previous Monte Carlo estimates. Additional model studies of this kind will be used to study the general information content of cloud top lightning radiation signatures.
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
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.
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.
A toy model linking atmospheric thermal radiation and sea ice growth
NASA Technical Reports Server (NTRS)
Thorndike, A. S.
1992-01-01
A simplified analytical model of sea ice growth is presented where the atmosphere is in thermal radiative equilibrium with the ice. This makes the downwelling longwave radiation reaching the ice surface an internal variable rather than a specified forcing. Analytical results demonstrate how the ice state depends on properties of the ice and on the externally specified climate.
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.
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.
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...
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.
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.
Advantages of Analytical Transformations in Monte Carlo Methods for Radiation Transport
McKinley, M S; Brooks III, E D; Daffin, F
2004-12-13
Monte Carlo methods for radiation transport typically attempt to solve an integral by directly sampling analog or weighted particles, which are treated as physical entities. Improvements to the methods involve better sampling, probability games or physical intuition about the problem. We show that significant improvements can be achieved by recasting the equations with an analytical transform to solve for new, non-physical entities or fields. This paper looks at one such transform, the difference formulation for thermal photon transport, showing a significant advantage for Monte Carlo solution of the equations for time dependent transport. Other related areas are discussed that may also realize significant benefits from similar analytical transformations.
Delgado-Aparicio, L.; Tritz, K.; Kramer, T.; Stutman, D.; Finkentha, M.; Hill, K.; Bitter, M.
2010-08-26
A new set of analytic formulae describes the transmission of soft X-ray (SXR) continuum radiation through a metallic foil for its application to fast electron temperature measurements in fusion plasmas. This novel approach shows good agreement with numerical calculations over a wide range of plasma temperatures in contrast with the solutions obtained when using a transmission approximated by a single-Heaviside function [S. von Goeler, Rev. Sci. Instrum., 20, 599, (1999)]. The new analytic formulae can improve the interpretation of the experimental results and thus contribute in obtaining fast teperature measurements in between intermittent Thomson Scattering data.
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.
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.
An urban radiation obstruction model
NASA Astrophysics Data System (ADS)
Frank, Randall S.; Gerding, R. Bruce; O'Rourke, Patricia A.; Terjung, Werner H.
1981-03-01
An urban street canyon radiation obstruction model has been developed. The model can describe community structure in terms of the type and dimensions of every building, block, road, park, etc. The need for massive data acquisition in regard to obstruction modeling calls for computerized algorithms, relieving the researcher of the needless tedium of hand calculations and the accompanying high degree of error and labor costs. The model program OBSTRUCT was written in FORTRAN IV for use on the IBM 3033. To facilitate changes or modifications, OBSTRUCT was written in modular form.
NASA Astrophysics Data System (ADS)
Täger, Olaf; Dannemann, Martin; Hufenbach, Werner A.
2015-04-01
Lightweight structures for high-technology applications are designed to meet the increasing demands on low structural weight and maximum stiffness. These classical lightweight properties result in lower inertial forces that consequently lead to higher vibration amplitudes thereby increasing sound radiation. Here, special anisotropic multilayered composites offer a high vibro-acoustic lightweight potential. The authors developed analytical vibro-acoustic simulation models, which allow a material-adapted structural-dynamic and sound radiation analysis of anisotropic multilayered composite plates. Compared to numerical methods FEM/BEM these analytical models allow a quick and physically based analysis of the vibro-acoustic properties of anisotropic composite plates. This advantage can be seen by the presented extensive parameter studies, which have been performed in order to analyse the influence of composite-specific design variables on the resulting vibro-acoustic behaviour. Here, it was found that the vibro-acoustic parameters like eigenfrequency and modal damping show direction-dependent properties. Furthermore, the investigations reveal that laminated composites show a so-called damping-dominated sound radiation behaviour. Based on these studies, a vibro-acoustic design procedure is proposed and design guidelines are derived.
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.
Fraidenraich, N
1995-08-01
Knowledge of optical and radiative properties is often essential for the design and evaluation of V-trough solar energy collectors. Using the concept of reflection modes, we derived a set of functions associated with each mode; this allowed us to calculate the optical and radiative properties for rejected light radiation. These expressions, together with those for accepted light radiation published previously, were used to calculate the optical efficiency for beam radiation and the exchange factors (diffuse radiation) between aperture and absorber (accepted light) and between aperture and aperture (rejected light). Numerical results of these factors were obtained for various combinations of concentration ratio and vertex angle. Results are compared between a case in which the reflectivity is constant and one in which the reflectivity varies with incidence angle; the difference does not exceed 1% for a reflectivity of 0.8. Considering the reflectivity as a constant allows us to obtain analytic solutions for the exchange factors, expressed as a sum of trigonometric functions. PMID:21052319
Analytical solution for boundary heat fluxes from a radiating rectangular medium
NASA Technical Reports Server (NTRS)
Siegel, R.
1991-01-01
Reference is made to the work of Shah (1979) which demonstrated the possibility of partially integrating the radiative equations analytically to obtain an 'exact' solution. Shah's solution was given as a double integration of the modified Bessel function of order zero. Here, it is shown that the 'exact' solution for a rectangular region radiating to cold black walls can be conveniently derived, and expressed in simple form, by using an integral function, Sn, analogous to the exponential integral function appearing in plane-layer solutions.
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.
Analytical solution of beam spread function for ocean light radiative transfer.
Xu, Zao; Yue, Dick K P
2015-07-13
We develop a new method to analytically obtain the beam spread function (BSF) for light radiative transfer in oceanic environments. The BSF, which is defined as the lateral distribution of the (scalar) irradiance with increasing depth in response to a uni-directional beam emanating from a point source in an infinite ocean, must in general be obtained by solving the three-dimensional (3D) radiative transfer equation (RTE). By taking advantage of the highly forward-peaked scattering property of the ocean particles, we assume, for a narrow beam source, the dependence of radiance on polar angle and azimuthal angle is deliberately separated; only single scattering takes place in the azimuthal direction while multiple scattering still occurs in the polar direction. This assumption enables us to reduce the five-variable 3D RTE to a three-variable two-dimensional (2D) RTE. With this simplification, we apply Fourier spectral method to both spatial and angular variables so that we are able to analytically solve the 2D RTE and obtain the 2D BSF accordingly. Using the relations between 2D and 3D solutions acquired during the process of simplification, we are able to obtain the 3D BSF in explicit form. The 2D and 3D analytical solutions are validated by comparing with Monte Carlo radiative transfer simulations. The 2D analytical BSF agrees excellently with the Monte Carlo result. Despite assumptions of axial symmetry and spike-like azimuthal profile of the radiance in deriving the 3D BSF, the comparisons to numerical simulations are very satisfactory especially for limited optical depths (< O(5)) for single scattering albedo values typical in the ocean. The explicit form of the analytical BSF and the significant gain in computational efficiency (several orders higher) relative to RTE simulations make many forward and inverse problems in ocean optics practical for routine applications. PMID:26191856
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...
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
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…
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
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.
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.
On Improving Analytical Models of Cosmic Reionization for Matching Numerical Simulations
Kaurov, Alexander A.
2016-01-01
The methods for studying the epoch of cosmic reionization vary from full radiative transfer simulations to purely analytical models. While numerical approaches are computationally expensive and are not suitable for generating many mock catalogs, analytical methods are based on assumptions and approximations. We explore the interconnection between both methods. First, we ask how the analytical framework of excursion set formalism can be used for statistical analysis of numerical simulations and visual representation of the morphology of ionization fronts. Second, we explore the methods of training the analytical model on a given numerical simulation. We present a new code which emerged from this study. Its main application is to match the analytical model with a numerical simulation. Then, it allows one to generate mock reionization catalogs with volumes exceeding the original simulation quickly and computationally inexpensively, meanwhile reproducing large scale statistical properties. These mock catalogs are particularly useful for CMB polarization and 21cm experiments, where large volumes are required to simulate the observed signal.
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.
Models for infrared atmospheric radiation
NASA Technical Reports Server (NTRS)
Tiwari, S. N.
1976-01-01
Line and band models for infrared spectral absorption are discussed. Radiative transmittance and integrated absorptance of Lorentz, Doppler, and voigt line profiles were compared for a range of parameters. It was found that, for the intermediate path lengths, the combined Lorentz-Doppler (Voigt) profile is essential in calculating the atmospheric transmittance. Narrow band model relations for absorptance were used to develop exact formulations for total absorption by four wide band models. Several continuous correlations for the absorption of a wide band model were compared with the numerical solutions of the wide band models. By employing the line-by-line and quasi-random band model formulations, computational procedures were developed for evaluating transmittance and upwelling atmospheric radiance. Homogeneous path transmittances were calculated for selected bands of CO, CO2, and N2O and compared with experimental measurements. The upwelling radiance and signal change in the wave number interval of the CO fundamental band were also calculated.
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.
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
Slot Region Radiation Environment Models
NASA Astrophysics Data System (ADS)
Sandberg, Ingmar; Daglis, Ioannis; Heynderickx, Daniel; Evans, Hugh; Nieminen, Petteri
2013-04-01
Herein we present the main characteristics and first results of the Slot Region Radiation Environment Models (SRREMs) project. The statistical models developed in SRREMs aim to address the variability of trapped electron and proton fluxes in the region between the inner and the outer electron radiation belt. The energetic charged particle fluxes in the slot region are highly dynamic and are known to vary by several orders of magnitude on both short and long timescales. During quiet times, the particle fluxes are much lower than those found at the peak of the inner and outer belts and the region is considered benign. During geospace magnetic storms, though, this region can fill with energetic particles as the peak of the outer belt is pushed Earthwards and the fluxes can increase drastically. There has been a renewed interest in the potential operation of commercial satellites in orbits that are at least partially contained within the Slot Region. Hence, there is a need to improve the current radiation belt models, most of which do not model the extreme variability of the slot region and instead provide long-term averages between the better-known low and medium Earth orbits (LEO and MEO). The statistical models developed in the SRREMs project are based on the analysis of a large volume of available data and on the construction of a virtual database of slot region particle fluxes. The analysis that we have followed retains the long-term temporal, spatial and spectral variations in electron and proton fluxes as well as the short-term enhancement events at altitudes and inclinations relevant for satellites in the slot region. A large number of datasets have been used for the construction, evaluation and inter-calibration of the SRREMs virtual dataset. Special emphasis has been given on the use and analysis of ESA Standard Radiation Environment Monitor (SREM) data from the units on-board PROBA-1, INTEGRAL, and GIOVE-B due to the sufficient spatial and long temporal
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
Radiation dosimetry and biophysical models of space radiation effects
NASA Technical Reports Server (NTRS)
Cucinotta, Francis A.; Wu, Honglu; Shavers, Mark R.; George, Kerry
2003-01-01
Estimating the biological risks from space radiation remains a difficult problem because of the many radiation types including protons, heavy ions, and secondary neutrons, and the absence of epidemiology data for these radiation types. Developing useful biophysical parameters or models that relate energy deposition by space particles to the probabilities of biological outcomes is a complex problem. Physical measurements of space radiation include the absorbed dose, dose equivalent, and linear energy transfer (LET) spectra. In contrast to conventional dosimetric methods, models of radiation track structure provide descriptions of energy deposition events in biomolecules, cells, or tissues, which can be used to develop biophysical models of radiation risks. In this paper, we address the biophysical description of heavy particle tracks in the context of the interpretation of both space radiation dosimetry and radiobiology data, which may provide insights into new approaches to these problems.
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.
Canopy radiation transmission for an energy balance snowmelt model
NASA Astrophysics Data System (ADS)
Mahat, Vinod; Tarboton, David G.
2012-01-01
To better estimate the radiation energy within and beneath the forest canopy for energy balance snowmelt models, a two stream radiation transfer model that explicitly accounts for canopy scattering, absorption and reflection was developed. Upward and downward radiation streams represented by two differential equations using a single path assumption were solved analytically to approximate the radiation transmitted through or reflected by the canopy with multiple scattering. This approximation results in an exponential decrease of radiation intensity with canopy depth, similar to Beer's law for a deep canopy. The solution for a finite canopy is obtained by applying recursive superposition of this two stream single path deep canopy solution. This solution enhances capability for modeling energy balance processes of the snowpack in forested environments, which is important when quantifying the sensitivity of hydrologic response to input changes using physically based modeling. The radiation model was included in a distributed energy balance snowmelt model and results compared with observations made in three different vegetation classes (open, coniferous forest, deciduous forest) at a forest study area in the Rocky Mountains in Utah, USA. The model was able to capture the sensitivity of beneath canopy net radiation and snowmelt to vegetation class consistent with observations and achieve satisfactory predictions of snowmelt from forested areas from parsimonious practically available information. The model is simple enough to be applied in a spatially distributed way, but still relatively rigorously and explicitly represent variability in canopy properties in the simulation of snowmelt over a watershed.
Comparing analytical and Monte Carlo optical diffusion models in phosphor-based X-ray detectors
NASA Astrophysics Data System (ADS)
Kalyvas, N.; Liaparinos, P.
2014-03-01
Luminescent materials are employed as radiation to light converters in detectors of medical imaging systems, often referred to as phosphor screens. Several processes affect the light transfer properties of phosphors. Amongst the most important is the interaction of light. Light attenuation (absorption and scattering) can be described either through "diffusion" theory in theoretical models or "quantum" theory in Monte Carlo methods. Although analytical methods, based on photon diffusion equations, have been preferentially employed to investigate optical diffusion in the past, Monte Carlo simulation models can overcome several of the analytical modelling assumptions. The present study aimed to compare both methodologies and investigate the dependence of the analytical model optical parameters as a function of particle size. It was found that the optical photon attenuation coefficients calculated by analytical modeling are decreased with respect to the particle size (in the region 1- 12 μm). In addition, for particles sizes smaller than 6μm there is no simultaneous agreement between the theoretical modulation transfer function and light escape values with respect to the Monte Carlo data.
Analytical model of non-line-of-sight single-scatter propagation.
Yin, Hongwei; Chang, Shengli; Wang, Xiaofeng; Yang, Jiankun; Yang, Juncai; Tan, Jichun
2010-07-01
An analytical model of non-line-of-sight (NLOS) single-scatter propagation is presented that has no integral form and is intended for performance analysis and system design of NLOS UV communication. Based on isotropic scattering and a continuous wave transmitter, the analytical model is verified by the current NLOS single-scatter propagation model, with consistent results. Several rules concerning NLOS UV communication are put forward on the basis of this analytical model, which are shown as follows: on condition that the minimum single-scatter optical depth is less than 0.1, the path loss factor should be 1; to maintain the NLOS UV communication link, the transmitter needs to radiate neither a continuous wave nor a huge pulse but a low-power wave whose duration is approximately the duration of impulse response; the "best" extinction coefficient is approximately the inverse ratio of the efficient single-scatter range; on condition that the radiation intensity of the transmitter is fixed, the half field of views (FOVs) are positive factors, while the elevation angles are negative factors; on condition that the power of the transmitter is fixed, the conclusions mentioned above remain valid with the exception that the half FOV of the transmitter is a negative factor. These rules also apply to anisotropic scattering. PMID:20596134
Cavity radiation model for solar central receivers
Lipps, F.W.
1981-01-01
The Energy Laboratory of the University of Houston has developed a computer simulation program called CREAM (i.e., Cavity Radiations Exchange Analysis Model) for application to the solar central receiver system. The zone generating capability of CREAM has been used in several solar re-powering studies. CREAM contains a geometric configuration factor generator based on Nusselt's method. A formulation of Nusselt's method provides support for the FORTRAN subroutine NUSSELT. Numerical results from NUSSELT are compared to analytic values and values from Sparrow's method. Sparrow's method is based on a double contour integral and its reduction to a single integral which is approximated by Guassian methods. Nusselt's method is adequate for the intended engineering applications, but Sparrow's method is found to be an order of magnitude more efficient in many situations.
Analytical procedures to identify foods that have been treated with ionizing radiation
Morehouse, K.M.
1994-12-31
Foods can be treated with ionizing radiation to reduce microbial infection and insect infestations, inhibit sprouting, and delay maturation, thereby extending the shelf life of foods. The treatment of different types of foods with ionizing radiation for specific purposes is accepted in several countries, although it is prohibited in others. The U.S. Food and Drug Administration (FDA) has established regulations to allow the treatment of several different foods with ionizing radiation and has received petitions for the approval of radiation treatment of additional foods. When carried out according to established good manufacturing practices, food irradiation yields safe, wholesome foods. Often, the irradiated product may be chemically and/or microbiologically {open_quotes}safer{close_quotes} than the nonirradiated product. An area of great interest in the last several years has been the development of analytical techniques to monitor foods that have been treated with ionizing radiation. A method for the identification of irradiated foods will help to foster compliance with labeling regulations, strengthen national and international regulations for the irradiation of specific foods, and enhance consumer confidence in the safety of food commodities that have been treated with ionizing radiation.
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).
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.
Space shuttle main engine plume radiation model
NASA Technical Reports Server (NTRS)
Reardon, J. E.; Lee, Y. C.
1978-01-01
The methods are described which are used in predicting the thermal radiation received by space shuttles, from the plumes of the main engines. Radiation to representative surface locations were predicted using the NASA gaseous plume radiation GASRAD program. The plume model is used with the radiative view factor (RAVFAC) program to predict sea level radiation at specified body points. The GASRAD program is described along with the predictions. The RAVFAC model is also discussed.
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.
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.
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.
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.
Lee, H. C.; Jiang, T. F.
2010-11-15
We analytically solve the relativistic equation of motion for an electron in ion plasma channels and calculate the corresponding trajectory as well as the synchrotron radiation. The relativistic effect on a trajectory is strong, i.e., many high-order harmonic terms in the trajectory, when the ratio of the initial transverse velocity (v{sub x0}) to the longitudinal velocity (v{sub z0}) of the electron injected to ion plasma channels is high. Interestingly, these high-order harmonic terms result in a quite broad and intense radiation spectrum, especially at an oblique angle, in contrast to an earlier understanding. As the initial velocity ratio (v{sub x0}:v{sub z0}) decreases, the relativistic effect becomes weak; only the first and second harmonic terms remain in the transverse and longitudinal trajectories, respectively, which coincides with the result of Esarey et al. [Phys. Rev. E 65, 056505 (2002)]. Our formalism also allows the description of electron's trajectory in the presence of an applied magnetic field. Critical magnetic fields for cyclotron motions are figured out and compared with semiclassical results. The cyclotron motion leads to more high-order harmonic terms than the trajectory without magnetic fields and causes an immensely broad spectrum with vastly large radiation amplitude for high initial velocity ratios (v{sub x0}:v{sub z0}). The radiation from hard x-ray to gamma-ray regions can be generated with a broad radiation angle, thus available for applications.
Interaction of airborne and structureborne noise radiated by plates. Volume 1: Analytical study
NASA Technical Reports Server (NTRS)
Mcgary, M. C.
1986-01-01
The interaction of airborne and structureborne noise radiated by aircraft materials was examined. The theory and results of several computer simulations of the noise radiated by thin, isotropic, rectangular aluminum plates due to fully coherent combined acoustic and vibrational inputs is presented. The most significant finding was the extremely large influence that the relative phase between inputs has on the combined noise radiation of the plates. Phase dependent effects manifest themselves as cross terms in both the dynamic and acoustic portions of the analysis. Computer simulations show that these cross terms can radically alter the combined sound power radiated by plates constructed of aircraft-type materials. The results suggest that airborne-structureborne interactive effects could be responsible for a significant portion of the overall noise radiated by aircraft-type structures in the low frequency regime. This implies that previous analytical and experimental studies may have neglected an important physical phenomenon in the analayses of the interior noise of propeller dirven aircraft.
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.
Zlatoustov, N N; Bizin, I A; Bondarenko, I M; Kaganov, V M; Grenkova, G K
2001-07-01
The developed information and analytic system (IAS) of decision making support in liquidation medical-and-sanitary consequences of the large-scale radiation accidents permits on the base of computer models to evaluate the size and degree of territorial radioactive pollution, the radiation doses, the number and structure of radiation affection, the forces and means of medical service required for rendering the casualties the first medical aid as well as to compare the available and necessary medical forces and means. IAS is intended for equipping the toxicologist-radiologist's working place. To solve the private problems the IAS computer methods can be used both in complex and autonomously. PMID:11561428
A stochastic model of radiation-induced bone marrow damage
Cotlet, G.; Blue, T.E.
2000-03-01
A stochastic model, based on consensus principles from radiation biology, is used to estimate bone-marrow stem cell pool survival (CFU-S and stroma cells) after irradiation. The dose response model consists of three coupled first order linear differential equations which quantitatively describe time dependent cellular damage, repair, and killing of red bone marrow cells. This system of differential equations is solved analytically through the use of a matrix approach for continuous and fractionated irradiations. The analytic solutions are confirmed through the dynamical solution of the model equations using SIMULINK. Rate coefficients describing the cellular processes of radiation damage and repair, extrapolated to humans from animal data sets and adjusted for neutron-gamma mixed fields, are employed in a SIMULINK analysis of criticality accidents. The results show that, for the time structures which may occur in criticality accidents, cell survival is established mainly by the average dose and dose rate.
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
Pérez-Andújar, Angélica; Zhang, Rui; Newhauser, Wayne
2013-01-01
Purpose: Stray neutron radiation is of concern after radiation therapy, especially in children, because of the high risk it might carry for secondary cancers. Several previous studies predicted the stray neutron exposure from proton therapy, mostly using Monte Carlo simulations. Promising attempts to develop analytical models have also been reported, but these were limited to only a few proton beam energies. The purpose of this study was to develop an analytical model to predict leakage neutron equivalent dose from passively scattered proton beams in the 100-250-MeV interval. Methods: To develop and validate the analytical model, the authors used values of equivalent dose per therapeutic absorbed dose (H/D) predicted with Monte Carlo simulations. The authors also characterized the behavior of the mean neutron radiation-weighting factor, \\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{upgreek} \\usepackage{mathrsfs} \\setlength{\\oddsidemargin}{-69pt} \\begin{document} \\overline {w_R }\\end{document}wR¯, as a function of depth in a water phantom and distance from the beam central axis. Results: The simulated and analytical predictions agreed well. On average, the percentage difference between the analytical model and the Monte Carlo simulations was 10% for the energies and positions studied. The authors found that \\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{upgreek} \\usepackage{mathrsfs} \\setlength{\\oddsidemargin}{-69pt} \\begin{document} \\overline {w_R }\\end{document}wR¯ was highest at the shallowest depth and decreased with depth until around 10 cm, where it started to increase slowly with depth. This was consistent among all energies. Conclusion: Simple analytical methods are promising alternatives to complex and slow Monte Carlo simulations
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.
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.
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 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.
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.
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
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.
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
Analytical model for contaminant mass removal by air sparging
Rabideau, A.J.; Blayden, J.M.
1998-12-31
An analytical model was developed to predict the removal of volatile organic compounds (VOCs) from ground water by air sparging (AS). The model treats the air sparging zone as a completely mixed reactor subject to the removal of dissolved contaminants by volatilization, advection, and first-order decay. Nonequilibrium desorption is approximated as a first-order mass transfer process. The model reproduces the tailing and rebound behavior often observed at AS sites, and would normally require the estimation of three site-specific parameters. Dimensional analysis demonstrates that predicting tailing can be interpreted in terms of kinetic desorption or diffusion of aqueous phase contaminants into discrete air channels. Related work is ongoing to test the model against field data.
Laser satellite constellations for strategic defense - an analytic model
Parmentola, J.A.; Milton, A.F.
1987-10-01
Using mainly geometric reasoning, an analytic model is constructed that predicts the required characteristics of an orbiting constellation of laser battle stations, each of which is designed to destroy ballistic missiles during their boost phase. The geometry of the constellation configuration and some general aspects of the coverage problem are discussed. The determination of the absentee ratio falls into two main categories that depend upon whether the Soviet ICBM threat is concentrated at a single location or whether it is distributed as it is now. A point-threat model and a distributive threat model are considered, the determination of the respective kill rates for these models is discussed, and the scaling properties of the laser constellation with respect to a change in the quantitative nature of the two types of ICBM threats are considered.
A two-dimensional analytical model of petroleum vapor intrusion
NASA Astrophysics Data System (ADS)
Yao, Yijun; Verginelli, Iason; Suuberg, Eric M.
2016-02-01
In this study we present an analytical solution of a two-dimensional petroleum vapor intrusion model, which incorporates a steady-state diffusion-dominated vapor transport in a homogeneous soil and piecewise first-order aerobic biodegradation limited by oxygen availability. This new model can help practitioners to easily generate two-dimensional soil gas concentration profiles for both hydrocarbons and oxygen and estimate hydrocarbon indoor air concentrations as a function of site-specific conditions such as source strength and depth, reaction rate constant, soil characteristics and building features. The soil gas concentration profiles generated by this new model are shown in good agreement with three-dimensional numerical simulations and two-dimensional measured soil gas data from a field study. This implies that for cases involving diffusion dominated soil gas transport, steady state conditions and homogenous source and soil, this analytical model can be used as a fast and easy-to-use risk screening tool by replicating the results of 3-D numerical simulations but with much less computational effort.
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.
Generalized Analytical Model for the Radio-Frequency Sheath
NASA Astrophysics Data System (ADS)
Czarnetzki, Uwe
2014-10-01
An analytical model for the planar radio frequency (RF) sheath in capacitive discharges is developed based on the applied RF voltage as the boundary condition. The model applies to all kind of waveforms for the applied RF voltage, includes both sheaths in a discharge of arbitrary symmetry, and allows for an arbitrary degree of ion collisionallity in the sheaths (charge-exchange collisions). Further, effects of the finite floating potential during sheath collapse are included. The model can even be extended to electronegative plasmas with low bulk conductivity. The individual sheath voltages, the self-bias, and the RF floating potentials are explicitly calculated by a voltage balance equation using a cubic-charge voltage relation for the sheaths. In particular, the RF-phase as a function of the sheath voltage is determined. This is an input for a single second order non-linear integro-differential equation which is governing the ion flow velocity in the sheath. Fast numerical integration is straight forward and in many cases approximate analytical solutions can be obtained. Based on the solution for the ion flow velocity, densities, electric fields, currents, and charge-voltage relations are calculated. Further, the Child-Langmuir laws for the collisionless as well as the highly collisional case are derived. Very good agreement between model and experiments is obtained.
Estimating solar radiation for plant simulation models
NASA Technical Reports Server (NTRS)
Hodges, T.; French, V.; Leduc, S.
1985-01-01
Five algorithms producing daily solar radiation surrogates using daily temperatures and rainfall were evaluated using measured solar radiation data for seven U.S. locations. The algorithms were compared both in terms of accuracy of daily solar radiation estimates and terms of response when used in a plant growth simulation model (CERES-wheat). Requirements for accuracy of solar radiation for plant growth simulation models are discussed. One algorithm is recommended as being best suited for use in these models when neither measured nor satellite estimated solar radiation values are available.
Analytical model for flow duration curves in seasonally dry climates
NASA Astrophysics Data System (ADS)
Müller, Marc F.; Dralle, David N.; Thompson, Sally E.
2014-07-01
Flow duration curves (FDC) display streamflow values against their relative exceedance time. They provide critical information for watershed management by representing the variation in the availability and reliability of surface water to supply ecosystem services and satisfy anthropogenic needs. FDCs are particularly revealing in seasonally dry climates, where surface water supplies are highly variable. While useful, the empirical computation of FDCs is data intensive and challenging in sparsely gauged regions, meaning that there is a need for robust, predictive models to evaluate FDCs with simple parameterization. Here, we derive a process-based analytical expression for FDCs in seasonally dry climates. During the wet season, streamflow is modeled as a stochastic variable driven by rainfall, following the stochastic analytical model of Botter et al. (2007a). During the dry season, streamflow is modeled as a deterministic recession with a stochastic initial condition that accounts for the carryover of catchment storage across seasons. The resulting FDC model is applied to 38 catchments in Nepal, coastal California, and Western Australia, where FDCs are successfully modeled using five physically meaningful parameters with minimal calibration. A Monte Carlo analysis revealed that the model is robust to deviations from its assumptions of Poissonian rainfall, exponentially distributed response times and constant seasonal timing. The approach successfully models period-of-record FDCs and allows interannual and intra-annual sources of variations in dry season streamflow to be separated. The resulting median annual FDCs and confidence intervals allow the simulation of the consequences of interannual flow variations for infrastructure projects. We present an example using run-of-river hydropower in Nepal as a case study.
Analytical model of solute transport by unsteady unsaturated gravitational infiltration.
Lessoff, S C; Indelman, P
2004-08-01
Penetration of reactive solute into a soil during a cycle of water infiltration and redistribution is investigated by deriving analytical closed form solutions for fluid flux, moisture content and contaminant concentration. The solution is developed for gravitational flow and advective transport and is applied to two scenarios of solute applications encountered in the applications: a finite pulse of solute dissolved in irrigation water and an instantaneous pulse broadcasted onto the soil surface. Through comparison to simulations of Richards' flow, capillary suction is shown to have contrasting effects on the upper and lower boundaries of the fluid pulse, speeding penetration of the wetting front and reducing the rate of drying. This leads to agreement between the analytical and numerical solutions for typical field and experimental conditions. The analytical solution is further incorporated into a stochastic column model of flow and transport to compute mean solute concentration in a heterogeneous field. An unusual phenomenon of plume contraction is observed at long times of solute propagation during the drying stage. The mean concentration profiles match those of the Monte-Carlo simulations for capillary length scales typical of sandy soils. PMID:15240168
Analytic Thermoelectric Couple Modeling: Variable Material Properties and Transient Operation
NASA Technical Reports Server (NTRS)
Mackey, Jonathan A.; Sehirlioglu, Alp; Dynys, Fred
2015-01-01
To gain a deeper understanding of the operation of a thermoelectric couple a set of analytic solutions have been derived for a variable material property couple and a transient couple. Using an analytic approach, as opposed to commonly used numerical techniques, results in a set of useful design guidelines. These guidelines can serve as useful starting conditions for further numerical studies, or can serve as design rules for lab built couples. The analytic modeling considers two cases and accounts for 1) material properties which vary with temperature and 2) transient operation of a couple. The variable material property case was handled by means of an asymptotic expansion, which allows for insight into the influence of temperature dependence on different material properties. The variable property work demonstrated the important fact that materials with identical average Figure of Merits can lead to different conversion efficiencies due to temperature dependence of the properties. The transient couple was investigated through a Greens function approach; several transient boundary conditions were investigated. The transient work introduces several new design considerations which are not captured by the classic steady state analysis. The work helps to assist in designing couples for optimal performance, and also helps assist in material selection.
A hybrid finite-difference and analytic element groundwater model.
Haitjema, H M; Feinstein, D T; Hunt, R J; Gusyev, M A
2010-01-01
Regional finite-difference models tend to have large cell sizes, often on the order of 1-2 km on a side. Although the regional flow patterns in deeper formations may be adequately represented by such a model, the intricate surface water and groundwater interactions in the shallower layers are not. Several stream reaches and nearby wells may occur in a single cell, precluding any meaningful modeling of the surface water and groundwater interactions between the individual features. We propose to replace the upper MODFLOW layer or layers, in which the surface water and groundwater interactions occur, by an analytic element model (GFLOW) that does not employ a model grid; instead, it represents wells and surface waters directly by the use of point-sinks and line-sinks. For many practical cases it suffices to provide GFLOW with the vertical leakage rates calculated in the original coarse MODFLOW model in order to obtain a good representation of surface water and groundwater interactions. However, when the combined transmissivities in the deeper (MODFLOW) layers dominate, the accuracy of the GFLOW solution diminishes. For those cases, an iterative coupling procedure, whereby the leakages between the GFLOW and MODFLOW model are updated, appreciably improves the overall solution, albeit at considerable computational cost. The coupled GFLOW-MODFLOW model is applicable to relatively large areas, in many cases to the entire model domain, thus forming an attractive alternative to local grid refinement or inset models. PMID:20132324
Digital forensics: an analytical crime scene procedure model (ACSPM).
Bulbul, Halil Ibrahim; Yavuzcan, H Guclu; Ozel, Mesut
2013-12-10
In order to ensure that digital evidence is collected, preserved, examined, or transferred in a manner safeguarding the accuracy and reliability of the evidence, law enforcement and digital forensic units must establish and maintain an effective quality assurance system. The very first part of this system is standard operating procedures (SOP's) and/or models, conforming chain of custody requirements, those rely on digital forensics "process-phase-procedure-task-subtask" sequence. An acceptable and thorough Digital Forensics (DF) process depends on the sequential DF phases, and each phase depends on sequential DF procedures, respectively each procedure depends on tasks and subtasks. There are numerous amounts of DF Process Models that define DF phases in the literature, but no DF model that defines the phase-based sequential procedures for crime scene identified. An analytical crime scene procedure model (ACSPM) that we suggest in this paper is supposed to fill in this gap. The proposed analytical procedure model for digital investigations at a crime scene is developed and defined for crime scene practitioners; with main focus on crime scene digital forensic procedures, other than that of whole digital investigation process and phases that ends up in a court. When reviewing the relevant literature and interrogating with the law enforcement agencies, only device based charts specific to a particular device and/or more general perspective approaches to digital evidence management models from crime scene to courts are found. After analyzing the needs of law enforcement organizations and realizing the absence of crime scene digital investigation procedure model for crime scene activities we decided to inspect the relevant literature in an analytical way. The outcome of this inspection is our suggested model explained here, which is supposed to provide guidance for thorough and secure implementation of digital forensic procedures at a crime scene. In digital forensic
NASA Astrophysics Data System (ADS)
Flynn, Edward M.; Choi, Paul S.
1990-06-01
Measurements and an analytical investigation of the thermal contact conductance of spotwelded joints are presented. Pairs of .081-cm thick sheets of Ti-6Al-4V alloy were resistance welded to form test coupons representing skins and stiffeners. The spotweld number density was varied. Measurements of contact conductance of these coupons were made at temperatures from 120 to 232 C. In addition, a thermal model of the spotweld interface was constructed which yielded good predictions of spotweld conduction, gaseous conduction, and gap radiation.
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.
Modeling of Laser-generated Radiative Blast Waves
Keilty, K. A.; Liang, E. P.; Ditmire, T.; Remington, B. A.; Shigemori, K.; Rubenchik, A. M.
2000-08-01
We simulate experiments performed with the Falcon laser at Lawrence Livermore National Laboratory to generate strong, cylindrically diverging blast waves of relevance to astrophysics. In particular, we are interested in producing and modeling radiative shocks. We compare numerical simulations with the data and with an analytic approximation to blast-wave propagation with a radiative-loss term included. Our goal is to develop a laboratory setting for studying radiative shocks of relevance to supernova remnants, gamma-ray burst afterglows, and other high-energy astrophysics phenomena. We will show that a good degree of agreement exists between the experimental data and the numerical simulations, demonstrating that it is indeed possible to generate radiative shocks in the laboratory using tabletop femtosecond lasers. In addition, we show how we can determine the energy-loss rate from the blast-wave evolution. This analytic method is independent of the exact mechanism of radiative cooling and is scalable to both the laboratory and astrophysical radiative blast waves. (c) 2000 The American Astronomical Society.
HTS axial flux induction motor with analytic and FEA modeling
NASA Astrophysics Data System (ADS)
Li, S.; Fan, Y.; Fang, J.; Qin, W.; Lv, G.; Li, J. H.
2013-11-01
This paper presents a high-temperature superconductor (HTS) axial-flux induction motor, which can output levitation force and torque simultaneously. In order to analyze the character of the force, analytic method and finite element method are adopted to model the motor. To make sure the HTS can carry sufficiently large current and work well, the magnetic field distribution in HTS coil is calculated. An effective method to improve the critical current of HTS coil is presented. Then, AC losses in HTS windings in the motor are estimated and tested.
"Violent Intent Modeling: Incorporating Cultural Knowledge into the Analytical Process
Sanfilippo, Antonio P.; Nibbs, Faith G.
2007-08-24
While culture has a significant effect on the appropriate interpretation of textual data, the incorporation of cultural considerations into data transformations has not been systematic. Recognizing that the successful prevention of terrorist activities could hinge on the knowledge of the subcultures, Anthropologist and DHS intern Faith Nibbs has been addressing the need to incorporate cultural knowledge into the analytical process. In this Brown Bag she will present how cultural ideology is being used to understand how the rhetoric of group leaders influences the likelihood of their constituents to engage in violent or radicalized behavior, and how violent intent modeling can benefit from understanding that process.
An analytical model of accretion onto white dwarfs
NASA Astrophysics Data System (ADS)
Ospina, N.; Hernanz, M.
2013-05-01
The analytical model of Frank et al. (2002) has been used to investigate the structure of the accretion stream onto white dwarfs (WD). In particular, the post-shock region (temperature, density and gas velocity distributions) and X-ray spectrum emitted by this region. We have obtained the temperature, density and gas velocity distributions of the emission region for different masses of white dwarfs and at different positions in the shock coordinate. Also, we calculated the emitted spectrum for different WD masses and at different positions of the shock with the principal objective of study the accretion at different points of the emission region.
Analytical and numerical modeling of surface morphologies in thin films
Genin, F.Y.
1995-05-01
Experimental studies have show that strains due to thermal expansion mismatch between a film and its substrate can produce very large stresses in the film that can lead to the formation of holes and hillocks. Based on a phenomenological description of the evolution of a solid surface under both capillary and stress driving forces and for surface and grain boundary self-diffusion, this article provides analytical and numerical solutions for surface profiles of model geometries in polycrystalline thin films. Results can explain a variety of surface morphologies commonly observed experimentally and are discussed to give some practical insights on how to control the growth of holes and hillocks in thin films.
Analytical model for the radio-frequency sheath
NASA Astrophysics Data System (ADS)
Czarnetzki, Uwe
2013-12-01
A simple analytical model for the planar radio-frequency (rf) sheath in capacitive discharges is developed that is based on the assumptions of a step profile for the electron front, charge exchange collisions with constant cross sections, negligible ionization within the sheath, and negligible ion dynamics. The continuity, momentum conservation, and Poisson equations are combined in a single integro-differential equation for the square of the ion drift velocity, the so called sheath equation. Starting from the kinetic Boltzmann equation, special attention is paid to the derivation and the validity of the approximate fluid equation for momentum balance. The integrals in the sheath equation appear in the screening function which considers the relative contribution of the temporal mean of the electron density to the space charge in the sheath. It is shown that the screening function is quite insensitive to variations of the effective sheath parameters. The two parameters defining the solution are the ratios of the maximum sheath extension to the ion mean free path and the Debye length, respectively. A simple general analytic expression for the screening function is introduced. By means of this expression approximate analytical solutions are obtained for the collisionless as well as the highly collisional case that compare well with the exact numerical solution. A simple transition formula allows application to all degrees of collisionality. In addition, the solutions are used to calculate all static and dynamic quantities of the sheath, e.g., the ion density, fields, and currents. Further, the rf Child-Langmuir laws for the collisionless as well as the collisional case are derived. An essential part of the model is the a priori knowledge of the wave form of the sheath voltage. This wave form is derived on the basis of a cubic charge-voltage relation for individual sheaths, considering both sheaths and the self-consistent self-bias in a discharge with arbitrary
Shear mechanical properties of the spleen: experiment and analytical modelling.
Nicolle, S; Noguer, L; Palierne, J-F
2012-05-01
This paper aims at providing the first shear mechanical properties of spleen tissue. Rheometric tests on porcine splenic tissues were performed in the linear and nonlinear regime, revealing a weak frequency dependence of the dynamic moduli in linear regime and a distinct strain-hardening effect in nonlinear regime. These behaviours are typical of soft tissues such as kidney and liver, with however a less pronounced strain-hardening for the spleen. An analytical model based on power laws is then proposed to describe the general shear viscoelastic behaviour of the spleen. PMID:22498291
Analytical model of ionospheric convection at subauroral latitudes
NASA Astrophysics Data System (ADS)
Deminov, M. G.; Kim, V. P.; Shubin, V. N.; Khegai, V. V.
An analytical model of plasma convection in the subauroral ionosphere is developed, assuming that the electric shielding of the inner magnetosphere is controlled by polarization of the hot ion zone situated behind the inner boundary of the plasma sheet. It is shown that, at subauroral latitudes at night, the plasma drifts eastward, while during the day it shifts westward. Thus, in the predmidnight sector, the direction of convection in the subauroral ionosphere is opposite to that of auroral convection. In general, the electric field, with a strength of 10 mV/m, has a meridional direction.
Zhou, Dong; Zhang, Hui; Ye, Peiqing
2015-01-01
Background Penumbra characteristics play a significant role in dose delivery accuracy for radiation therapy. For treatment planning, penumbra width and radiation field offset strongly influence target dose conformity and organ at risk sparing. Methods In this study, we present an analytical and numerical approach for evaluation of the rounded leaf end effect on penumbra characteristics. Based on the rule of half-value layer, algorithms for leaf position calculation and radiation field offset correction were developed, which were advantageous particularly in dealing with large radius leaf end. Computer simulation was performed based on the Monte Carlo codes of EGSnrc/BEAMnrc, with groups of leaf end radii and source sizes. Data processing technique of curve fitting was employed for deriving penumbra width and radiation field offset. Results Results showed that penumbra width increased with source size. Penumbra width curves for large radius leaf end were U-shaped. This observation was probably related to the fact that radiation beams penetrated through the proximal and distal leaf sides. In contrast, source size had negligible impact on radiation field offset. Radiation field offsets were found to be constant both for analytical method and numerical simulation. However, the overall resulting values of radiation field offset obtained by analytical method were slightly smaller compared with Monte Carlo simulation. Conclusions The method we proposed could provide insight into the investigation of rounded leaf end effects on penumbra characteristics. Penumbra width and radiation field offset calibration should be carefully performed to commission multileaf collimator for intensity modulated radiotherapy. PMID:26401137
Experimental, numerical and analytical models of mantle starting plumes
NASA Astrophysics Data System (ADS)
Coulliette, D. L.; Loper, D. E.
1995-12-01
The results of a combined experimental, numerical and analytical investigation of starting thermal plumes are described, to obtain a better perspective on plumes within the Earth's mantle. Thermal plumes were generated experimentally in a tank of corn syrup by means of an electrical heater. Viscosity ratios of 400, 30 000, and 10 8 were generated by varying the temperature of the tank. Plumes for the smaller ratios had the traditional 'balloon-on-astring' shape, but that at the highest ratio had a novel morphology. The plume heads in the first two cases were observed to rise at roughly a constant speed, in contrast to most previous studies which found the plume heads to accelerate. Loss of buoyancy from the plume head owing to heat loss is believed to be responsible for this difference. Starting plumes were simulated numerically using an axisymmetric, finite-element code to solve the Boussinesq equations at finite Prandtl numbers. The constant rise speed observed experimentally was confirmed by the numerical simulation for the viscosity ratios of 400 and 30 000, but numerical instability prevented simulation of the case with a viscosity ratio of 10 8. There was very good agreement between the experimental and numerical rise speeds. An analytical model was developed which reduces to previous models in limiting cases. This parameterization gives better agreement with the experimental and numerical results than does any previous model.
An Analytical Model for the Radio-Frequency Sheath
NASA Astrophysics Data System (ADS)
Czarnetzki, Uwe
2013-09-01
An analytical model for the planar radio frequency (RF) sheath in capacitive discharges is developed based on the applied RF voltage as the boundary condition. In a first step, the individual sheath voltages and the self-bias are calculated using a cubic-charge voltage relation. In the second step, a single integro-differential equation is derived to describe the ion flow velocity in the sheath under all conditions of collisionality. Central to the model is the screening function that describes the screening of the ion density by the mean electron density in the sheath. Numerical integration of the sheath equation is straight forward. However, for the collisionless as well as the collisional case explicit, simple, and precise analytical approximations can be found. Drift velocities, densities, fields, currents, and charge-voltage relations are calculated. Further, the Child-Langmuir laws for both cases of collisonality are derived. These solutions are in very good agreement with experimental data from the literature based on laser electric field measurements, the Brinkmann sheath model, and PIC simulations. The technique works well also for other waveforms, e.g. the electrical asymmetry effect or tailored pulse waveforms.
Analytical modeling of turbine wakes in yawed conditions
NASA Astrophysics Data System (ADS)
Bastankhah, Majid; Porté-Agel, Fernando
2016-04-01
Increasing wind energy production has become a unanimous plan for virtually all the developed countries. In addition to constructing new wind farms, this goal can be achieved by making wind farms more efficient. Control strategies in wind farms, such as manipulating the yaw angle of the turbines, have the potential to make wind farms more efficient. Costly numerical simulations or measurements cannot be, however, employed to assess the viability of this strategy in the numerous different scenarios happening in real wind farms. In this study, we aim to develop an inexpensive and simple analytical model that is able for the first time to predict the whole wake of a yawed turbine with an acceptable accuracy. The proposed analytical model is built upon the simplified version of the Reynolds-averaged Navier-Stokes equations. Apart from the ability of the model to predict wake flows in yawed conditions, it can provide a better understanding of turbine wakes in this complex situation. For example, it can give valuable insights on how the wake deflection varies by changing turbine and incoming flow characteristics, such as the thrust coefficient of the turbine or the ambient turbulence.
Urban stormwater management planning with analytical probabilistic models
Adams, B.J.
2000-07-01
Understanding how to properly manage urban stormwater is a critical concern to civil and environmental engineers the world over. Mismanagement of stormwater and urban runoff results in flooding, erosion, and water quality problems. In an effort to develop better management techniques, engineers have come to rely on computer simulation and advanced mathematical modeling techniques to help plan and predict water system performance. This important book outlines a new method that uses probability tools to model how stormwater behaves and interacts in a combined- or single-system municipal water system. Complete with sample problems and case studies illustrating how concepts really work, the book presents a cost-effective, easy-to-master approach to analytical modeling of stormwater management systems.
Analytic modeling of antibody versus nanocell delivery of photosensitizer
NASA Astrophysics Data System (ADS)
Pogue, Brian W.; Srinivasan, Subhadra; Samkoe, Kimberley; Zheng, Lei Zak; Rai, Prakash; Mai, Zhiming; Verma, Sarika; Hasan, Tayyaba
2010-02-01
Delivery of therapeutic agents to solid tumors is challenging, and the issues that govern this can be distilled down into parameters which allow computational modeling. In this paper, the basic rate equations and diffusion kernel for the time and space modeling of delivery are developed, along with an analytical solution to this equation. The model is then used to compare delivery of Avastin antibody relative to delivery encapsulated in a nanocell delivery vehicle. The key factors are the plasma clearance or excretion rates, and the binding, or not, as it transports into the tumor tissue. A reduction in the plasma clearance rate inherently increases available delivery over time, and additionally the reduction in binding from antibody to nanocell allows higher penetration into the tumor at the longer circulation times.
Semi-analytical model of cosmic ray electron transport
NASA Astrophysics Data System (ADS)
Ivascenko, A.; Spanier, F.
2011-07-01
We present a numerical extension to the analytical propagation model introduced in Hein and Spanier (2008) to describe the leptonic population in the galactic disc. The model is used to derive a possible identification of the components that contribute to the leptonic cosmic ray spectrum, as measured by PAMELA, Fermi and HESS, with an emphasis on secondary e+-e- production in collisions of cosmic ray particles with ambient interstellar medium (ISM). We find that besides secondaries, an additional source symmetric in e+ and e- production is needed to explain both the PAMELA anomaly and the Fermi bump, assuming a power-law primary electron spectrum. Our model also allows us to derive constraints for some properties of the ISM.
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