jacobian-free newton krylov: Topics by WorldWideScience.org

Linear Multifrequency-Grey Acceleration Recast for Preconditioned Krylov Iterations

Fully coupled, Newton-Krylov algorithms are investigated for solving strongly coupled, nonlinear systems of partial differential equations arising in the field of computational fluid dynamics. Primitive variable forms of the steady incompressible and compressible Navier-Stokes and energy equations that describe the flow of a laminar Newtonian fluid in two-dimensions are specifically considered. Numerical solutions are obtained by first integrating over discrete finite volumes that compose the computational mesh. The resulting system of nonlinear algebraic equations are linearized using Newton`s method. Preconditioned Krylov subspace based iterative algorithms then solve these linear systems on each Newton iteration. Selected Krylov algorithms include the Arnoldi-based Generalized Minimal RESidual (GMRES) algorithm, and the Lanczos-based Conjugate Gradients Squared (CGS), Bi-CGSTAB, ...

1995-10-01

High Order Finite Volume Nonlinear Schemes for the Boltzmann Transport Equation

Abstract not provided

2007-03-19

Newton's first law and the existence of free tachyons

The authors apply the nonlinear WENO (Weighted Essentially Nonoscillatory) scheme to the spatial discretization of the Boltzmann Transport Equation modeling linear particle transport. The method is a finite volume scheme which ensures not only conservation, but also provides for a more natural handling of boundary conditions, material properties and source terms, as well as an easier parallel implementation and post processing. It is nonlinear in the sense that the stencil depends on the solution at each time step or iteration level. By biasing the gradient calculation towards the stencil with smaller derivatives, the scheme eliminates the Gibb's phenomenon with oscillations of size O(1) and reduces them to O(h{sup r}), where h is the mesh size and r is the order of accuracy. The current implementation is three-dimensional, generalized for unequally spaced meshes, fully parallelized, and up to fifth order accurate (WENO5) in space. For unsteady problems, the resulting ...

2005-03-29

Control area performance improvement by extended security constrained economic dispatch

NEWTON's first law prohibits the coupling of free tachyons to ordinary matter, if it is interpreted as selection rule interdicting emission without change of the rest-mass of the emitting system.

On multitemporal generalization of Newton's gravitational law

An algorithm for solving the extended security constrained economic dispatch (ESCED) problem with real-time economic dispatch grade speed and reliability is presented. The ESCED problem is formulated by adding regulating margin and ramp rate constraints to the network security constrained economic dispatch problem previously solved by the CEDC algorithm. Starting with Newton`s method to optimize the Lagrangian, the ESCED is developed by superimposing on Newton`s method eight major components called Tracking Start Initialization, Hessian Pre-Elimination, Implicit Dual Variable Calculations, Regulating Margin Sensitivity Coefficient Calculations, Traumatic Event Evaluation, Constraint Relaxation, Implicit Ramp Rate Constraint Implementation, and Relaxed Incremental Cost Calculations. Test results are also presented.

1997-02-01

Investigating the complex X-ray spectrum of a broad-line 2MASS red quasar: XMM-Newton observation of FTM 0830+3759

A n-time generalization of Newton's law (of universal gravitation) formula in N =n + d + 1-dimensional space-time is conjectured. This formula implies a relation for effective N-dimensional gravitational constant G_{eff} = G cos^2 \\theta, where \\theta is the angle between the direction of motion of two particles in n-dimensional time manifold R^n.

2009-01-01

Analysis and design of high speed thrust ball bearings. (Reannouncement with new availability information)

We report results from a 50 ks XMM-Newton observation of the dust-reddened broad-line quasar FTM 0830+3759 (z=0.413) selected from the FIRST/2MASS Red Quasar survey. For this AGN, a very short 9 ks Chandra exposure had suggested a feature-rich X-ray spectrum and HST images revealed a very disturbed host galaxy morphology. Contrary to classical, optically-selected quasars, the X-ray properties of red (i.e. with J-Ks> 1.7 and R-Ks> 4) broad line quasars are still quite unexplored, although there is a growing consensus that, due to moderate obscuration, these objects can offer a unique view of spectral components typically swamped by the AGN light in normal, blue quasars. The XMM-Newton observation discussed here has definitely confirmed the complexity of the X-ray spectrum revealing the presence of a cold (or mildly-ionized) absorber with Nh ~10^{22} cm^-2 along the line of sight to the nucleus and a Compton reflection component accompanied ...

2010-01-01

An inexact Newton method for fully-coupled solution of the Navier-Stokes equations with heat and mass transport

In this paper, a complete solution of the steady-state Newton-Euler equations describing the motion of the balls in high speed thrust bill bearings is presented. Kinematic equations are derived and used as ancillary equations to obtain the numerical solution of the Newton-Euler equations. This differs from previously reported work where extraneous assumptions (instead of kinematic analysis) are used to allow solutions of these equations. It is shown that the no-slip condition for the ball motion reported in the literature is based on an unnecessary assumption. The correct form of the condition is developed and used in the analysis of the ball motion. The analysis shows the importance of contact angle on both ball motion and contact force. An approximate theory is then developed for determining the optimum contact angle that would result in minimum contact force in high speed ball bearings. The results are compared with experimental data, and ...

1991-12-31

WITNESSING THE KEY EARLY PHASE OF QUASAR EVOLUTION: AN OBSCURED ACTIVE GALACTIC NUCLEUS PAIR IN THE INTERACTING GALAXY IRAS 20210+1121

The solution of the governing steady transport equations for momentum, heat and mass transfer in flowing fluids can be very difficult. These difficulties arise from the nonlinear, coupled, nonsymmetric nature of the system of algebraic equations that results from spatial discretization of the PDEs. In this manuscript the authors focus on evaluating a proposed nonlinear solution method based on an inexact Newton method with backtracking. In this context they use a particular spatial discretization based on a pressure stabilized Petrov-Galerkin finite element formulation of the low Mach number Navier-Stokes equations with heat and mass transport. The discussion considers computational efficiency, robustness and some implementation issues related to the proposed nonlinear solution scheme. Computational results are presented for several challenging CFD benchmark problems as well as two large scale 3D flow simulations.

1997-02-01

The intriguing nature of the high energy gamma ray source XSSJ12270-4859

We report the discovery of an active galactic nucleus (AGN) pair in the interacting galaxy system IRAS 20210+1121 at z = 0.056. An XMM-Newton observation reveals the presence of an obscured (N _H #approx# 5 x 10"2"3 cm"-"2), Seyfert-like (L _2_-_1_0_k_e_V = 4.7 x 10"4"2 erg s"-"1) nucleus in the northern galaxy, which lacks unambiguous optical AGN signatures. Our spectral analysis also provides strong evidence that the IR-luminous southern galaxy hosts a Type 2 quasar embedded in a bright starburst emission. In particular, the X-ray primary continuum from the nucleus appears totally depressed in the XMM-Newton band as expected in the case of a Compton-thick absorber, and only the emission produced by Compton scattering ('reflection') of the continuum from circumnuclear matter is seen. As such, IRAS 20210+1121 seems to provide an excellent opportunity to witness a key, early phase in the quasar evolution predicted by the theoretical models of ...

2010-10-20

Lubrication of an 85-mm ball bearing with RP-1

The nature of the hard X-ray source XSSJ12270-4859 is still unclear though it was claimed to be a magnetic Cataclysmic Variable. We here present a broad-band X-ray and gamma ray study based on a recent XMM-Newton observation and archival INTEGRAL and RXTE data. From the Fermi/LAT 1-year point source catalogue, we tentatively associate XSSJ12270-4859 with 1FGLJ1227.9-4852, a source of high energy gamma rays with emission up to 10GeV. We complement the study with UV photometry from XMM-Newton and ground-based optical and near-IR photometry. The X-ray emission is highly variable showing flares and intensity dips. The X-ray flares consist of flare-dip pairs. Flares are also detected in the UV range but not the dips. Aperiodic dipping behaviour is also observed during X-ray quiescence but not in the UV. The 0.2-100keV spectrum is featureless and described by a power law model with Gamma=1.7. The 100MeV-10GeV spectrum is instead represented by a ...

2010-01-01

Stimulated radiation of high - current relativistic electron beams

A parametric experimental investigation of an 85 millimeter bore angular contact ball bearing running in RP-1 fuel was performed at speeds of 10,000 to 24,000 rpm. Thrust loads were varied from 4450 to 17,800 Newtons (1000 to 4000 lbs.). Radial loads were varied from 1335 to 13,350 Newtons (300 to 3000 lbs.). RP-1 lubrication for the bearing was provided through a stationary jet ring located adjacent to the test bearing outer ring. Increases in both the thrust and radial loads resulted in increased bearing temperature, while increases in shaft speed resulted in much more dramatic increases in bearing temperature. These trends are typical for ball bearings operating under these types of conditions. Results are given for outer ring temperatures of the test bearing at the various test conditions employed. In addition, the heat energy removed from the bearing by the RP-1 was determined by measuring the increase in temperature as the RP-1 passed ...

1993-01-01

Gauge equivalence of representations of symmetry groups in quantum mechanics

The most propagated mechanisms of stimulated radiation of electron beam such as Cherenkov one-particle and collective effects, ondulator and magnetic bremsshrahlung radiations, Doppler anomalous effect, Thompson and Raman scattering and radiation are discussed. Relation of spontaneous radiation mechanisms of individual electron and stimulated radiation effects in electron beams has been elucidated, grounds of linear electrodynamics of radiative beam instabilities are stated, and main mechanisms of their nonlinear stabilization are elucidated as well. Various simulated processes in electron beams are considered from the unique point of view using a simple mathematical apparatus and such physical laws as conservation and Newton laws.

1987-01-01

Constraints on extra dimensions from cosmological and terrestrial measurements

The equivalence of representations of symmetry groups operating upon wave-functions in configuration space is studied with regard to the (intuitive) notion of physical equivalence. A refinement of the usual projective equivalence relation is introduced, called gauge equivalence, for which the allowed unitary equivalence transformations are gauge transformations. For a Euclidean as well as for Newton-Hooke symmetry group the gauge equivalence classes of unitary multiplier representations are determined. These examples support the assertion that equivalence from a physical viewpoint corresponds better to this new gauge equivalence concept than to the usual notion of projective equivalence. (author).

Constraints on Extra Dimensions from Cosmological and Terrestrial Measurements

If quantum fields exist in extra compact dimensions, they will give rise to a quantum vacuum or Casimir energy. That vacuum energy will manifest itself as a cosmological constant. The fact that supernova and cosmic microwave background data indicate that the cosmological constant is of the same order as the critical mass density to close the universe supplies a lower bound on the size of the extra dimensions. Recent laboratory constraints on deviations from Newton's law place an upper limit. The allowed region is so small as to suggest that either extra compact dimensions do not exist, or their number is about to be tightly constrained by experimental data.

2001-01-01

An efficient higher order family of root finders

If quantum fields exist in extra compact dimensions, they will give rise to a quantum vacuum or Casimir energy. That vacuum energy will manifest itself as a cosmological constant. The fact that supernova and cosmic microwave background data indicate that the cosmological constant is of the same order as the critical mass density to close the universe supplies a lower bound on the size of the extra dimensions. Recent laboratory constraints on deviations from Newton's law place an upper limit. The allowed region is so small as to suggest that either extra compact dimensions do not exist, or their properties are about to be tightly constrained by experimental data.

2000-01-01

British Library Electronic Table of Contents (United Kingdom)

A one parameter family of iterative methods for the simultaneous approximation of simple complex zeros of a polynomial, based on a cubically convergent Hansen-Patricks family, is studied. We show that the convergence of the basic family of the fourth order can be increased to five and six using Newtons and Halleys corrections, respectively. Since these corrections use the already calculated values, the computational efficiency of the accelerated methods is significantly increased. Further acceleration is achieved by applying the Gauss-Seidel approach (single-step mode). One of the most important problems in solving nonlinear equations, the construction of initial conditions which provide both the guaranteed and fast convergence, is considered for the proposed accelerated family. These cond...

2008-01-01

APSIS - an Artificial Planetary System in Space to probe extra-dimensional gravity and MOND

Three Dimensional Visualization for the Steam Injection into Water Pool using Electrical Resistance Tomography

A proposal is made to test Newton's inverse-square law using the perihelion shift of test masses (planets) in free fall within a spacecraft located at the Earth-Sun L2 point. Such an Artificial Planetary System In Space (APSIS) will operate in a drag-free environment with controlled experimental conditions and minimal interference from terrestrial sources of contamination. We demonstrate that such a space experiment can probe the presence of a "hidden" fifth dimension on the scale of a micron, if the perihelion shift of a "planet" can be measured to sub-arc-second accuracy. Some suggestions for spacecraft design are made.

2006-01-01

The X-ray Halo of G21.5-0.9

The direct injection of steam into a water pool is a method of heat transfer used in many process industries. The amount of research in this area however is limited to the nuclear industry, with applications relating to reactor cooling systems. Electrical resistance tomography (ERT), a low cost, non-invasive and which has high temporal resolution characteristics, can be used as a visualization tool for the resistivity distribution for the steam injection into water pool such as IRWST. In this paper, three dimensional resistivity distribution of the process is obtained through ERT using iterative Gauss-Newton method. Numerical experiments are performed by assuming different resistive objects in the water pool. Numerical results show that ERT is successful in estimating the resistivity distribution for the injection of steam in the water pool

2010-10-01

The Stefan problem solved via conjugate gradient-like iterative methods on a parallel vector machine

The emission of the plerion G21.5-0.9 appears more extended in X rays than in radio. This is an unexpected result because it would imply that short-lived X-ray electrons may reach distances even larger than radio electrons. Applying an empirical relationship between dust scattering optical depth and photoelectric column density, the measured column density leads to a large optical depth at 1 keV, of about 1. Therefore we investigate the hypothesis that the detected halo be an effect of dust scattering, re-analyzing an Cal/PV XMM-Newton observation of G21.5-0.9 and critically examining it in terms of a dust scattering model. We also present a spectral analysis of a prominent extended feature in the northern sector of the halo.

2003-01-01

1 2 3

Science.gov (United States)

The aim of this paper is to illustrate the validity and efficiency of iterative methods for solving large linear systems arising from the finite element discretization of the equation governing conduction-controlled solidification processes. Starting from the basic enthalpy equation, two alternative formulations are obtained and fixed-grid finite element discretizations are developed. These discretizations yield a set of nonlinear equations that are linearized using the Newton-Raphson scheme. The linearized equations are used as a basis for evaluating different iterative methods of the conjugated gradient type. Symmetric scaling and incomplete factorization preconditioning of the linear equations are used to improve the convergence properties of the iterative methods. Vectorization and parallelization are also employed to make full use of the CRAY-2 supercomputer. The results indicate that the implementation of currently available iterative solvers leads to ...

1991-01-01

The Cosmological Constant and Lorentz Invariance of the Vacuum State

Research on optimizing pass schedule of tandem cold mill

One hope to solve the cosmological constant problem is to identify a symmetry principle, based on which the cosmological constant can be reduced either to zero, or to a tiny value. Here, we note that requiring that the vacuum state is Lorentz invariant significantly reduces the theoretical value of the vacuum energy density. Hence, this also reduces the discrepancy between the observed value of the cosmological constant and its theoretical expectation, down from 123 orders of magnitude to 56 orders of magnitude. We find that, at one loop level, massless particles do not yield any contribution to the cosmological constant. Another important consequence of Lorentz symmetry is stabilization of the gravitational hierarchy: the cosmological constant (divided by Newton's constant) does not run as the quartic power of the renormalization group scale, but instead only logarithmically.

2011-01-01

Research on optimizing pass schedule of tandem cold mill

In this paper, research on pass schedule of tandem cold mill (TCM) is carried out. According to load (reduction, rolling force, motor power) balance, non-linear equations set with variables of inter-stand thickness is constructed. The pass schedule optimization is carried out by solving the non-linear equations set. As the traditional method, the Newton-Raphson method is used for solving the non-linear equations set. In this paper a new simple method is brought up. On basis of the monotone relations between thickness and load, the inter-stands thickness is adjusted dynamically. The solution of non-linear equations set can be converged by iterative calculation. This method can avoid the derivative calculation used by traditional method. So, this method is simple and calculation speed is high. It is suitable for on-line control. (author)

2000-07-01

Optimization and sensitivity analysis of controlling arrivals in the queueing system with single working vacation

2000-08-20

British Library Electronic Table of Contents (United Kingdom)

This paper analyzes the F-policy M/M/1/K queueing system with working vacation and an exponential startup time. The F-policy deals with the issue of controlling arrivals to a queueing system, and the server requires a startup time before allowing customers to enter the system. For the queueing systems with working vacation, the server can still provide service to customers rather than completely stop the service during a vacation period. The matrix-analytic method is applied to develop the steady-state probabilities, and then obtain several system characteristics. We construct the expected cost function and formulate an optimization problem to find the minimum cost. The direct search method and Quasi-Newton method are implemented to determine the optimal system capacity K, the optimal thre...

2010-01-01

Investigating the Ultraviolet Properties of Gravity with a Wilsonian Renormalization Group Equation

FUV and X-ray absorption in the Warm-Hot Intergalactic Medium

We review and extend in several directions recent results on the asymptotic safety approach to quantum gravity. The central issue in this approach is the search of a Fixed Point having suitable properties, and the tool that is used is a type of Wilsonian renormalization group equation. We begin by discussing various cutoff schemes, i.e. ways of implementing the Wilsonian cutoff procedure. We compare the beta functions of the gravitational couplings obtained with different schemes, studying first the contribution of matter fields and then the so-called Einstein-Hilbert truncation, where only the cosmological constant and Newton's constant are retained. In this context we make connection with old results, in particular we reproduce the results of the epsilon expansion and the perturbative one loop divergences. We then apply the Renormalization Group to higher derivative gravity. In the case of a general action quadratic in curvature we recover, within certain ...

2008-01-01

The Warm-Hot Intergalactic Medium (WHIM) arises from shock-heated gas collapsing in large-scale filaments and probably harbours a substantial fraction of the baryons in the local Universe. Absorption-line measurements in the ultraviolet (UV) and in the X-ray band currently represent the best method to study the WHIM at low redshifts. We here describe the physical properties of the WHIM and the concepts behind WHIM absorption line measurements of H I and high ions such as O VI, O VII, and O VIII in the far-ultraviolet and X-ray band. We review results of recent WHIM absorption line studies carried out with UV and X-ray satellites such as FUSE, HST, Chandra, and XMM-Newton and discuss their implications for our knowledge of the WHIM.

2008-01-01

Do Spinors Frame-Drag?

Using tree swallows to monitor impacts of aquatic contamination in Great Lakes areas of concern

We investigate the effect of the intrinsic spin of a fundamental spinor field on the surrounding spacetime geometry. We show that despite the lack of a rotating stress-energy source (and despite claims to the contrary) the intrinsic spin of a spin-half fermion gives rise to a frame-dragging effect analogous to that of orbital angular momentum, even in Einstein-Hilbert gravity where torsion is constrained to be zero. This resolves a paradox regarding the counter-force needed to restore Newton's third law in the well known spin-orbit interaction. In addition, the frame-dragging effect gives rise to a {\\it long-range} gravitationally mediated spin-spin dipole interaction coupling the {\\it internal} spins of two sources. We argue that despite the weakness of the interaction, the spin-spin interaction will dominate over the ordinary inverse square Newtonian interaction in any process of sufficiently high-energy for quantum field theoretical effects to be ...

2009-01-01

An X-ray source population study of the Andromeda galaxy M 31

Tree swallows were used to evaluate movement and potential impacts of contaminants from sediments in Newton Creek (diesel range organics: DROs) and Sheboygan River (PCBs), tributaries to Lake Superior and Lake Michigan, respectively. Contaminated sites occurred along the course of each river, while reference sites were located upstream or on a nearby river. Productivity was monitored and eggs, day 1 nestlings, and day 12 nestlings were collected from each nest. Whole body or egg homogenates were analyzed for PCBs or DROs. EROD activity in livers from day 12 nestlings is being determined for both PCB and DRO exposures. In the Newton Creek study, hatching success was similar for DRO and reference sites. DROs were detected in gastrointestinal tracts of 1 nestling from the reference and 1 from the contaminated site. DROs were not detected in any egg samples. In the Sheboygan River study, hatching success rates differed between 1 reference and 1 ...

1995-12-31

The jet-powered optical nebula of Cygnus X-1

XMM-Newton EPIC observations reveal the population of X-ray sources of the bright Local Group spiral galaxy M 31, a low-star-formation-rate galaxy like the Milky Way, down to a 0.2-4.5 keV luminosity of 4.4E34 erg/s. With the help of X-ray hardness ratios and optical and radio information different source classes can be distinguished. The survey detected 856 sources in an area of 1.24 square degrees. Sources within M 31 are 44 supernova remnants (SNR) and candidates, 18 super-soft sources (SSS), 16 X-ray binaries (XRBs) and candidates, as well as 37 globular cluster sources (GlC) and candidates, i.e. most likely low mass XRBs within the GlC. 567 hard sources may either be XRBs or Crab-like SNRs in M 31 or background AGN. 22 sources are new SNR candidates in M 31 based on X-ray selection criteria. Time variability information can be used to improve the source classification. Two GlC sources show type I X-ray bursts as known from Galactic neutron star low mass XRBs. ...

2005-01-01

The XMM-Newton wide-field survey in the COSMOS field. IV: X-ray spectral properties of Active Galactic Nuclei

We present H-alpha and [O III] (5007 Angstroms) images of the nebula powered by the jet of the black hole candidate and microquasar Cygnus X-1, observed with the 2.5m Isaac Newton Telescope (INT). The ring-like structure is luminous in [O III] and there exists a thin outer shell with a high [O III] / H-alpha flux ratio. This outer shell probably originates in the collisionally excited atoms close to the front of the bow shock. Its presence indicates that the gas is shock excited as opposed to photoionised, supporting the jet-powered scenario. The shock velocity was previously constrained at 20 = 100 km/s (1 sigma confidence) based on a comparison of the observed [O III] / H-alpha ratio in the bow shock with a number of radiative shock models. From this we further constrain the time-averaged power of the jet: P_Jet = (4 - 14)*10^36 erg/s. The H-alpha flux behind the shock front is typically 4*10^-15 erg/s/cm^2/arcsec^2, and we estimate an upper limit of ~8*10^-15 ...

2007-01-01

Some numerical problems in atomic physics

We present a detailed spectral analysis of point-like X-ray sources in the XMM-COSMOS field. Our sample of 135 sources only includes those that have more than 100 net counts in the 0.3-10 keV energy band and have been identified through optical spectroscopy. The majority of the sources are well described by a simple power-law model with either no absorption (76%) or a significant intrinsic, absorbing column (20%).As expected, the distribution of intrinsic absorbing column densities is markedly different between AGN with or without broad optical emission lines. We find within our sample four Type-2 QSOs candidates (L_X > 10^44 erg/s, N_H > 10^22 cm^-2), with a spectral energy distribution well reproduced by a composite Seyfert-2 spectrum, that demonstrates the strength of the wide field XMM/COSMOS survey to detect these rare and underrepresented sources.

2006-01-01

Probing the Origins of the CIV and Fe Kalpha Baldwin Effect

The numerical calculation of oscillatory wavefunctions occurring in the theory of electron scattering from positive ions is considered, the aim being to develop methods for evaluating the oscillatory solutions using a logarithmic grid. Thus, the permitted step-size for the solution depends only upon the potential, enabling calculations at different energies to be performed using the same discretisation. Two similar generalisations of the W.K.B. transformation for second order ordinary differential equations are made and applied to the radial Dirac equations. Iterative schemes for the solution of the resulting non-linear amplitude equations are discussed, one for use in the asymptotic region and the other a Newton scheme for continuing the solution closer to the origin. Using these techniques a reliable and efficient program for calculating relativistic scattering phase shifts from ionic potentials has been written. Relativistic Optical Model problems have been ...

1981-01-01

Particle Swarm Optimization: An efficient method for tracing periodic orbits in 3D galactic potentials

We use UV/optical and X-ray observations of 272 radio-quiet Type 1 AGNs and quasars to investigate the CIV Baldwin Effect (BEff). The UV/optical spectra are drawn from the Hubble Space Telescope, International Ultraviolet Explorer and Sloan Digital Sky Survey archives. The X-ray spectra are from the Chandra and XMM-Newton archives. We apply correlation and partial-correlation analyses to the equivalent widths, continuum monochromatic luminosities, and alpha_ox, which characterizes the relative X-ray to UV brightness. The equivalent width of the CIV 1549 emission line is correlated with both alpha_ox and luminosity. We find that by regressing l_UV with EW(CIV) and alpha_ox, we can obtain tighter correlations than by regressing l_UV with only EW(CIV). Both correlation and regression analyses imply that l_UV is not the only factor controlling the changes of EW(CIV); alpha_ox (or, equivalently, the soft X-ray emission) plays a fundamental role in the formation and ...

2009-01-01

Parameters determination of IEC-TDG by directly solving non-linear equation; Chokusetsuho ni yoru IEC-TDG hakei kaiseki

We propose the Particle Swarm Optimization (PSO) as an alternative method for locating periodic orbits in a three--dimensional (3D) model of barred galaxies. We develop an appropriate scheme that transforms the problem of finding periodic orbits into the problem of detecting global minimizers of a function, which is defined on the Poincar\\'{e} Surface of Section (PSS) of the Hamiltonian system. By combining the PSO method with deflection techniques, we succeeded in tracing systematically several periodic orbits of the system. The method succeeded in tracing the initial conditions of periodic orbits in cases where Newton iterative techniques had difficulties. In particular, we found families of 2D and 3D periodic orbits associated with the inner 8:1 to 12:1 resonances, between the radial 4:1 and corotation resonances of our 3D Ferrers bar model. The main advantages of the proposed algorithm is its simplicity, its ability to work using function values solely, as ...

2005-01-01

Line intensity enhancements in stellar coronal X-ray spectra due to opacity effects

An IEC-TDG waveform analysis program that can accurately calculate the waveform parameters in waveforms containing much noise as well as waveforms based on the IEC 1083-2 standard was developed in the impulse voltage and current test field. In the conventional curve fitting method, it is difficult to solve the non-linear equation directly when obtaining waveform parameters from the assigned discrete data. An approximate solution method is low in precision. Therefore, the residual between the standardized data and approximate function was calculated for high-speed Fourier transformation. Unrelated components were eliminated from the calculated frequency spectrum to synthesize waveforms. An L-M method is used for analysis. The solution is first tracked by a maximum diving method. If the convergence is judged to have been largely promoted, the method is switched into a Newton method to continue the calculation. The non-linear equation in a curve fitting part can be ...

1997-04-20

He-like ions as practical astrophysical plasma diagnostics: From stellar coronae to active galactic nuclei

Context. The I(15.01 A)/I(16.78 A) emission line intensity ratio in Fe XVII has been reported to deviate from its theoretical value in solar and stellar X-ray spectra. This is attributed to opacity in the 15.01 A line, leading to a reduction in its intensity, and was interpreted in terms of a geometry in which the emitters and absorbers are spatially distinct. Aims. We study the I(15.01 A)/I(16.78 A) intensity ratio for the active cool dwarf EV Lac, in both flare and quiescent spectra. Methods. The observations were obtained with the Reflection Grating Spectrometer on the XMM-Newton satellite. The emission measure distribution versus temperature reconstruction technique is used for our analysis. Results. We find that the 15.01 A line exhibits a significant enhancement in intensity over the optically thin value. To our knowledge, this is the first time that such an enhancement has been detected on such a sound statistical basis. We interpret this enhancement in ...

2008-01-01

Extension of QMD to the relativistic region and the calculation of proton induced reactions

We review X-ray plasma diagnostics based on the line ratios of He-like ions. Triplet/singlet line intensities can be used to determine electronic temperature and density, and were first developed for the study of the solar corona. Since the launches of the X-ray satellites Chandra and XMM-Newton, these diagnostics have been extended and used (from CV to Si XIII) for a wide variety of astrophysical plasmas such as stellar coronae, supernova remnants, solar system objects, active galactic nuclei, and X-ray binaries. Moreover, the intensities of He-like ions can be used to determine the ionization process(es) at work, as well as the distance between the X-ray plasma and the UV emission source for example in hot stars. In the near future thanks to the next generation of X-ray satellites (e.g., Astro-H and IXO), higher-Z He-like lines (e.g., iron) will be resolved, allowing plasmas with higher temperatures and densities to be probed. Moreover, the so-called satellite ...

2011-01-01

DETECTION OF AN X-RAY PULSATION FOR THE GAMMA-RAY PULSAR CENTERED IN CTA 1

When quantum molecular dynamics (QMD) is applied to the nuclear reactions in theory of relativity region, a number of problems arise, and in order to solve them, the prescription of the extension of ordinary nonrelativistic QMD is introduced, and the analysis of proton incidence reaction by using it is shown. By introducing the interaction corresponding to Lorentz transformation, the problems were solved. QMD is the semiclassical simulation that treats the motion of nucleons represented by Gauss wave packet. The motion of wave packet center is expressed by Newton equations and two-nucleon collision. The introduction of the interaction corresponding to Lorentz transformation is explained. As the result of the introduction, through the relative distance of two particles, the interaction becomes to depend on momentum. The phase distribution function of one body corresponding to Lorentz transformation is used for calculating the final state Pauli effect of collision ...

1995-06-01

Are Optically-Selected Quasars Universally X-Ray Luminous? X-Ray/UV Relations in Sloan Digital Sky Survey Quasars

We report the detection of X-ray pulsations with a period of ?315.87 ms from the 2009 XMM-Newton observation for the radio-quiet ?-ray pulsar, LAT PSR J0007+7303, centered in the supernova remnant CTA 1. The detected pulsed period is consistent with the ?-ray periodicity at the same epoch found with the Fermi Gamma-ray Space Telescope. The broader sinusoidal structure in the folded light curve of the X-ray emission is dissimilar to that of the ?-ray emission, and the phase of the peak is about 0.5 shifting from the peak in the ?-ray bands, indicating that the main component of the X-rays originates from different sites of the pulsar. We conclude that the main component of the X-ray pulsation is contributed by the thermal emission from the neutron star. Although with a significantly different characteristic age, PSR J0007+7303 is similar to Geminga in emission properties of X-rays and ?-rays; this makes PSR J0007+7303 the second radio-quiet ?-ray pulsar with ...

2010-12-10

1 2 3

A unified model of combined energy systems with different cycle modes and its optimum performance characteristics

We analyze archived Chandra and XMM-Newton X-ray observations of 536 Sloan Digital Sky Survey (SDSS) Data Release 5 (DR5) quasars (QSOs) at 1.7 <= z <= 2.7 in order to characterize the relative UV and X-ray spectral properties of QSOs that do not have broad UV absorption lines (BALs). We constrain the fraction of X-ray weak, non-BAL QSOs and find that such objects are rare; for example, sources underluminous by a factor of 10 comprise $\\la$2% of optically-selected SDSS QSOs. X-ray luminosities vary with respect to UV emission by a factor of $\\la$2 over several years for most sources. UV continuum reddening and the presence of narrow-line absorbing systems are not strongly associated with X-ray weakness in our sample. X-ray brightness is significantly correlated with UV emission line properties, so that relatively X-ray weak, non-BAL QSOs generally have weaker, blueshifted CIV$\\lambda$1549 emission and broader CIII]$\\lambda$1909 lines. The CIV emission ...

2008-01-01

Utilization of smooth models for seismic data tomographic inversion

A unified model is presented for a class of combined energy systems, in which the systems mainly consist of a heat engine, a combustor and a counter-flow heat exchanger and the heat engine in the systems may have different thermodynamic cycle modes such as the Brayton cycle, Carnot cycle, Stirling cycle, Ericsson cycle, and so on. Not only the irreversibilities of the heat leak and finite-rate heat transfer but also the different cycle modes of the heat engine are considered in the model. On the basis of Newton's law, expressions for the overall efficiency and power output of the combined energy system with an irreversible Brayton cycle are derived. The maximum overall efficiency and power output and other relevant parameters are calculated. The general characteristic curves of the system are presented for some given parameters. Several interesting cases are discussed in detail. The results obtained here are very general and significant and can be used to ...

2009-06-15

Traffic flow optimization on rural networks

Reflection tomography allows in theory the determination of the velocity distribution and the reflector geometries in the subsurface from the travel times of the seismic waves. The solution of the inverse problem is the model (velocity distribution and reflectors) that minimizes the misfits between observed travel times and travel times computed by raytracing (the forward problem). Two representations of the velocity field are generally used: blocky models (discontinuous) and smooth models. The object of this thesis is the study of tomography with smooth models which offer interesting properties at a practical level as well as at a mathematical and numerical level. We propose an original formulation of the inverse problem which allow the integration of a priori information adapted to smooth models. This formulation assures the existence, the uniqueness and the stability of the solution of the linearized inverse problem. The non-linear inverse problem is solved by a ...

1995-03-01

PROBING THE ORIGINS OF THE C IV AND Fe K? BALDWIN EFFECTS

Seasonal traffic on intercity corridors often creates severe congestion at various bottlenecks, such as land-drops and drawbridges. A Traffic Simulation And Optimization model (TSAO) has been developed in order to 1) evaluate traffic impacts (including travel time, vehicle operating costs, accidents, air pollution and fuel consumption) of a given routing pattern on a network, and 2) determine the optimal flow pattern. In its latest version, the TSAO model uses: 1) a macroscopic simulation approach to trace vehicle platoons through a network; 2) Greenshield's equation to represent volume-speed relationship on freeways and multi-lane highways, and Highway Capacity Manual method for two-lane rural highways; 3) Lighthill's shockwave function to simulate queue propagation; 4) event-scan time management for simulation; 5) variable link capacity for two-lane rural highways or to stimulate incidents; 6) AASHTO Guideline for determining vehicle operating costs and fuel ...

1983-01-01

Measuring q sub 0 using supernovae at z approx 0. 3

We use UV/optical and X-ray observations of 272 radio-quiet Type 1 active galactic nuclei and quasars to investigate the C IV Baldwin Effect (BEff). The UV/optical spectra are drawn from the Hubble Space Telescope, International Ultraviolet Explorer and Sloan Digital Sky Survey archives. The X-ray spectra are from the Chandra and XMM-Newton archives. We apply correlation and partial-correlation analyses to the equivalent widths (EWs), continuum monochromatic luminosities, and ?ox, which characterizes the relative X-ray to UV brightness. The EW of the C IV ?1549 emission line is correlated with both ?ox and luminosity. We find that by regressing l?(2500 A) with EW(C IV) and ?ox, we can obtain tighter correlations than by regressing l?(2500 A) with only EW(C IV). Both correlation and regression analyses imply that l?(2500 A) is not the only factor controlling the changes of EW(C IV); ?ox (or, equivalently, the soft X-ray emission) plays a fundamental role in the ...

2009-09-01

Gaussian pulse decomposition: an intuitive model of electrocardiogram waveforms.

The measurement of q{sub o} is extremely important for understanding the quantity of matter in our universe. The measurement of q{sub o} using supernovae of type Ia as standard candles is appealing because it requires less modeling than other methods using galaxies. The challenge with using supernovae to measure q{sub o} is in finding enough of them. In order to find supernovas, we have constructed a very popular f/1 camera for the 3.9m Anglo-Australian Telescope. The camera uses reducing optics that put a 17 in. {times} 17 in. field on a 1024 {times} 1024 pixel Thomson CCD. Using this system, we image to 23rd magnitude in five minutes. We have developed a software package that uses image subtraction to find supernovae that are approximately magnitude 22.4 or brighter in these images. One field can be processed every 6.6 minutes on a relatively unloaded VAX 6000-6510. We estimate that this system should find one supernova in every 105--139 images (about two nights of observation on the ...

1992-07-01

Science.gov (United States)

This study presents a novel approach to modeling the electrocardiogram (ECG): the Gaussian pulse decomposition. Constituent waves of the ECG are decomposed into and represented by Gaussian pulses using an iterative algorithm: the chip away decomposition (ChAD) algorithm. At each iteration, a nonlinear minimization method is used to fit a portion of the ECG waveform with a single Gaussian pulse, which is then subtracted from the ECG waveform. The process iterates on the resulting residual waveform until the normalized mean square error is below an acceptable level. Three different minimization methods were compared for their applicability to the ChAD algorithm; the Nelder-Mead simplex method was found to be more noise-tolerant than the Newton-Raphson method or the steepest descent method. Using morphologically different ECG waveforms from the MIT-BIH arrhythmia database, it was demonstrated that the ChAD algorithm is capable of modeling not only normal beats, but ...

Development and application of a third order scheme of finite differences centered in mesh; Desarrollo y aplicacion de un esquema de tercer orden de diferencias finitas centradas en malla

Demonstrating leadership inside and out : green building project trains employees, educates clients

In this work the development of a third order scheme of finite differences centered in mesh is presented and it is applied in the numerical solution of those diffusion equations in multi groups in stationary state and X Y geometry. Originally this scheme was developed by Hennart and del Valle for the monoenergetic diffusion equation with a well-known source and they show that the one scheme is of third order when comparing the numerical solution with the analytical solution of a model problem using several mesh refinements and boundary conditions. The scheme by them developed it also introduces the application of numeric quadratures to evaluate the rigidity matrices and of mass that its appear when making use of the finite elements method of Galerkin. One of the used quadratures is the open quadrature of 4 points, no-standard, of Newton-Cotes to evaluate in approximate form the elements of the rigidity matrices. The other quadrature is that of 3 points of Radau ...

2003-07-01

A 3-D hydrodynamic dispersion model for modeling tracer transport in Geothermal Reservoirs

This article described a full-scale green renovation project undertaken by a construction management firm in Newton, Massachusetts. Chapman Construction and Design began renovating its own office building in late 2007 and completed it in September 2008. The firm set out to achieve LEED Platinum certification from the United States Green Building Council. The project involved energy saving modifications to its mechanical systems; replacing plumbing fixtures with high efficiency alternatives to reduce water consumption by at least 40 per cent; installing skylights and additional windows to allow more natural light into the interior space; salvaging doors, lighting, metal studs and masonry whenever possible; and using sustainable products such as recyclable carpets, tiles and low-VOC paint. The main feature of the renovation was a new 47 kW (DC) photovoltaic (PV) solar power array installed on the rooftop. The PV array includes 208 panels and a SunPower SPR-225 system ...

2009-01-15

THE DOUBLE-DEGENERATE NUCLEUS OF THE PLANETARY NEBULA TS 01: A CLOSE BINARY EVOLUTION SHOWCASE

A 3-D hydrodynamic dispersion model for tracer transport is developed and implemented into the TOUGH2 EOS3 (T2R3D) module. The model formulation incorporates a full dispersion tensor, based on a 3-D velocity field with a 3-D, irregular grid in a heterogeneous geological system. Two different weighting schemes are proposed for spatial average of 3-D velocity fields and concentration gradients to evaluate the mass flux by dispersion and diffusion of a tracer or a radionuclide. This new module of the TOUGH2 code is designed to simulate processes of tracer/radionuclide transport using an irregular, 3-D integral finite difference grid in non-isothermal, three-dimensional, multiphase, porous/fractured subsurface systems. The numerical method for this transport module is based on the integral finite difference scheme, as in the TOUGH2 code. The major assumptions of the tracer transport module are: (a) a tracer or a radionuclide is present and transported only within the liquid phase, (b) ...