Composite scheme using localized relaxation with non-standard finite difference method for hyperbolic conservation laws

Non-standard finite difference methods (NSFDM) introduced by Mickens [Non-standard Finite Difference Models of Differential Equations, World Scientific, Singapore, 1994] are interesting alternatives to the traditional finite difference and finite volume methods. When applied to linear hyperbolic conservation laws, these methods reproduce exact solutions. In this paper, the NSFDM is first extended to hyperbolic systems of conservation laws, by a novel utilization of the decoupled equations using characteristic variables. In the second part of this paper, the NSFDM is studied for its efficacy in application to nonlinear scalar hyperbolic conservation laws. The original NSFDMs introduced by Mickens (1994) were not in conservation form, which is an important feature in capturing discontinuities at the right locations. Mickens [Construction and analysis of a non-standard finite difference scheme for the Burgers–Fisher equations, Journal of Sound and Vibration 257 (4) (2002) 791–797] recently introduced a NSFDM in conservative form. This method captures the shock waves exactly, without any numerical dissipation. In this paper, this algorithm is tested for the case of expansion waves with sonic points and is found to generate unphysical expansion shocks. As a remedy to this defect, we use the strategy of composite schemes [R. Liska, B. Wendroff, Composite schemes for conservation laws, SIAM Journal of Numerical Analysis 35 (6) (1998) 2250–2271] in which the accurate NSFDM is used as the basic scheme and localized relaxation NSFDM is used as the supporting scheme which acts like a filter. Relaxation schemes introduced by Jin and Xin [The relaxation schemes for systems of conservation laws in arbitrary space dimensions, Communications in Pure and Applied Mathematics 48 (1995) 235–276] are based on relaxation systems which replace the nonlinear hyperbolic conservation laws by a semi-linear system with a stiff relaxation term. The relaxation parameter (λ) is chosen locally on the three point stencil of grid which makes the proposed method more efficient. This composite scheme overcomes the problem of unphysical expansion shocks and captures the shock waves with an accuracy better than the upwind relaxation scheme, as demonstrated by the test cases, together with comparisons with popular numerical methods like Roe scheme and ENO schemes.

Dielectric anisotropy and relaxation studies of the homologous series of N-(4-n-alkyloxy-2-hydroxy-benzylidene)-4-carbethoxyaniline

Dielectric properties of the homologous series of newly synthesized nonchiral compounds N-(4-n-alkyloxy-2-hydroxy-benzylidene)-4-carbethoxyaniline, (n = 6, 8, 10, 12) having wide temperature range (∼60°C) smectic A (SmA) phase, have been studied by the impedance spectroscopy in the frequency range of 100 Hz to 1 MHz. Measurements have been carried out for two principal alignments (planar as well as homeotropic) of the SmA phase. Dielectric anisotropy (Δε' = ε'∥ - ε'⊥) for all the members of the series has been found to be negative for the whole temperature range of SmA phase. Magnitude of the dielectric anisotropy (|Δε'|) has been found to decrease with the number of alkyl chains. Relaxation frequencies corresponding to the rotation of the individual molecules about their short axes, lie below 1 MHz and obey the Arrhenius law by which activation energies have been determined. However, the relaxation frequencies corresponding to the rotation of the molecules about their short axes apparently lie above 10 MHz.

Influence of sintering conditions and doping on the dielectric relaxation originating from the surface layer effects in CaCu3Ti4O12 ceramics

Detailed investigations into the dielectric dispersion phenomenon in the giant dielectric constant material CaCu3Ti4O12 (CCTO) around room temperature revealed the existence of two successive dielectric relaxations. In the temperature domain, a new dielectric relaxation was clearly observed around 250K, in addition to the well-investigated dielectric relaxation close to 100K. The effect of sintering and doping (La3+) on the strength of these dielectric relaxations were studied in detail. The sintering temperature as well as its duration was found to have tremendous influence on the dielectric relaxation that was encountered around 250 K. This Maxwell-Wagner (M-W) type of relaxation was found to be originating from the surface layer containing the Cu-rich phase, which was ascribed to the difference in the oxygen content between the surface and the interior of the sample. Interestingly, this particular additional relaxation was not observed in La2/3Cu3Ti4O12, a low dielectric constant member of the CCTO family, in which the segregation of Cu-rich phase on the surface was absent. Indeed the correlation between the new relaxation and the presence of Cu-rich phase in CCTO ceramics was further corroborated by the absence of the same after removing the top and bottom layers. (C) 2007 Elsevier Ltd. All rights reserved.

Influence of Cooling Rate on the Enthalpy Relaxation and Fragility of a Metallic Glass

Structural relaxation behavior of a rapidly quenched (RQ) and a slowly cooled Pd40Cu30Ni10P20 metallic glass was investigated and compared. Differential scanning calorimetry was employed to monitor the relaxation enthalpies at the glass transition temperature, T-g , and the Kolrausch-Williams-Watts (KWW) stretched exponential function was used to describe its variation with annealing time. It was found that the rate of enthalpy recovery is higher in the ribbon, implying that the bulk is more resistant to relaxation at low temperatures of annealing. This was attributed to the possibility of cooling rate affecting the locations where the glasses get trapped within the potential energy landscape. The RQ process traps a larger amount of free volume, resulting in higher fragility, and in turn relaxes at the slightest thermal excitation (annealing). The slowly cooled bulk metallic glass (BMG), on the other hand, entraps lower free volume and has more short-range ordering, hence requiring a large amount of perturbation to access lower energy basins.

Vibrational phase relaxation of O–H stretch in bulk water: Role of large amplitude angular jumps and negative cross-correlations among the forces on the O–H bond

The ultrafast vibrational phase relaxation of O–H stretch in bulk water is investigated in molecular dynamics simulations. The dephasing time (T2) of the O–H stretch in bulk water calculated from the frequency fluctuation time correlation function (Cω(t)) is in the range of 70–80 femtosecond (fs), which is comparable to the characteristic timescale obtained from the vibrational echo peak shift measurements using infrared photon echo [W.P. de Boeij, M.S. Pshenichnikov, D.A. Wiersma, Ann. Rev. Phys. Chem. 49 (1998) 99]. The ultrafast decay of Cω(t) is found to be responsible for the ultrashort T2 in bulk water. Careful analysis reveals the following two interesting reasons for the ultrafast decay of Cω(t). (A) The large amplitude angular jumps of water molecules (within 30–40 fs time duration) provide a large scale contribution to the mean square vibrational frequency fluctuation and gives rise to the rapid spectral diffusion on 100 fs time scale. (B) The projected force, due to all the atoms of the solvent molecules on the oxygen (FO(t)) and hydrogen (FH(t)) atom of the O–H bond exhibit a large negative cross-correlation (NCC). We further find that this NCC is partly responsible for a weak, non-Arrhenius temperature dependence of the dephasing rate.

Subquadratic wavenumber dependence of the structural relaxation of supercooled liquid in the crossover regime

As a liquid is progressively supercooled toward its glass transition temperature, an intriguing weakening of the wavenumber (q) dependence of the structural relaxation time tau(q) in the intermediate-to-large q limit is observed both in experiments and simulation studies. Neither continuous Brownian diffusive dynamics nor discontinuous activated events can alone explain the anomalous wavenumber dependence. Here we use our recently developed theory that unifies the mode coupling theory for continuous dynamics, with the random first order transition theory treatment of activated discontinuous motion as a nucleationlike instanton process to understand the wavenumber dependence of density relaxation. The predicted smooth change in mechanism of relaxation from diffusive to activated, in the crossover regime, is wavevector dependent and appears to be responsible for the observed subquadratic,nalmost linear, q dependence of the relaxation time.

Dielectric relaxation and electrical conductivity in Bi5NbO10 oxygen ion conductors prepared by a modified sol–gel process

Crystalline Bi5NbO10 nanoparticles have been achieved through a modified sol–gel process using a mixture of ethylenediamine and ethanolamine as a solvent. The Bi5NbO10 nanoparticles were characterized by X-ray diffraction (XRD), differential scanning calorimetry/thermogravimetry (DSC/TG), Fourier transform infrared spectroscopy (FT-IR), transmission electron microscopy (TEM) and Raman spectroscopy. The results showed that well-dispersed 5–60 nm Bi5NbO10 nanoparticles were prepared through heat-treating the precursor at 650 °C and the high density pellets were obtained at temperatures lower than those commonly employed. The frequency and temperature dependence of the dielectric constant and the electrical conductivity of the Bi5NbO10 solid solutions were investigated in the 0.1 Hz to 1 MHz frequency range. Two distinct relaxation mechanisms were observed in the plots of dielectric loss and the imaginary part of impedance (Z″) versus frequency in the temperature range of 200–350 °C. The dielectric constant and the loss in the low frequency regime were electrode dependent. The ionic conductivity of Bi5NbO10 solid solutions at 700 °C is 2.86 Ω−1 m−1 which is in same order of magnitude for Y2O3-stabilized ZrO2 ceramics at same temperature. These results suggest that Bi5NbO10 is a promising material for an oxygen ion conductor.

Vibrational force constant and electron relaxation in H2 and Li2

The force constants of H2 and Li2 are evaluated employing their extended Hartree-Fock wavefunctions by a polynomial fit of their force curves. It is suggested that, based on incomplete multiconfiguration Hartree-Fock wavefunctions, force constants calculated from the energy derivatives are numerically more accurate than those obtained from the derivatives of the Hellmann-Feynman forces. It is observed that electrons relax during the nuclear vibrations in such a fashion as to facilitate the nuclear motions.

NMR relaxation study of disorder in condensed matter: solid solutions of betaine phosphate and phosphite

Proton NMR relaxation measurements have been carried out in anti-ferroelectric Betaine phosphate (BP), ferroelectric Betaine phosphite (BPI) and the mixed system BPI(1-x)BPx, at 11.4MHz and 23.3MHz from 300K to 80K for x=0.0, 0.25, 0.45, 0.85, and 1.0. The temperature dependence of spin lattice relaxation time T, exhibits two minima as expected from the BPP model in BP and BPI. The Larmor frequency dependence of T, in the mixed system is rather unusual and exhibits different slopes for the low temperature wings at the two frequencies, which is a clear experimental evidence of the presence of different methyl groups with different activation energies (E-a) indicating disorder.

Temperature- and pressure-dependent study of Cl-35 NQR frequency and spin lattice relaxation time in 2,3-dichloroanisole

The temperature and pressure dependence of Cl-35 NQR frequency and spin lattice relaxation time (T-1) were investigated in 2,3-dichloroanisole. Two NQR signals were observed throughout the temperature and pressure range studied. T-1 were measured in the temperature range from 77 to 300 K and from atmospheric pressure to 5 kbar. Relaxation was found to be due to the torsional motion of the molecule and also reorientation f motion of the CH3 group. T-1 versus temperature data were analyzed on the basis of Woessner and Gutowsky model, and the activation energy for the reorientation of the CH3 group was estimated. The temperature dependence of the average torsional lifetimes of the molecules and the transition probabilities were also obtained. NQR frequency shows a nonlinear behavior with pressure, indicating both dynamic and static effects of pressure. The pressure coefficients were observed to be positive for both the lines. A thermodynamic analysis of the data was carried out to determine the constant volume temperature coefficients of the NQR frequency. The variation of spin lattice time with pressure was very small, showing that the relaxation is mainly due to the torsional motions of the molecules. Copyright (C) 2010 John Wiley & Sons, Ltd.

A molecular theory of collective orientational relaxation in pure and binary dipolar liquids

A molecular theory of collective orientational relaxation of dipolar molecules in a dense liquid is presented. Our work is based on a generalized, nonlinear, Smoluchowski equation (GSE) that includes the effects of intermolecular interactions through a mean‐field force term. The effects of translational motion of the liquid molecules on the orientational relaxation is also included self‐consistently in the GSE. Analytic expressions for the wave‐vector‐dependent orientational correlation functions are obtained for one component, pure liquid and also for binary mixtures. We find that for a dipolar liquid of spherical molecules, the correlation function ϕ(k,t) for l=1, where l is the rank of the spherical harmonics, is biexponential. At zero wave‐vector, one time constant becomes identical with the dielectric relaxation time of the polar liquid. The second time constant is the longitudinal relaxation time, but the contribution of this second component is small. We find that polar forces do not affect the higher order correlation functions (l>1) of spherical dipolar molecules in a linearized theory. The expression of ϕ(k,t) for a binary liquid is a sum of four exponential terms. We also find that the wave‐vector‐dependent relaxation times depend strongly on the microscopic structure of the dense liquid. At intermediate wave vectors, the translational diffusion greatly accelerates the rate of orientational relaxation. The present study indicates that one must pay proper attention to the microscopic structure of the liquid while treating the translational effects. An analysis of the nonlinear terms of the GSE is also presented. An interesting coupling between the number density fluctuation and the orientational fluctuation is uncovered.

Polarization relaxation, dielectric dispersion, and solvation dynamics in dense dipolar liquid

A unified treatment of polarization relaxation, dielectric dispersion and solvation dynamics in a dense, dipolar liquid is presented. It is shown that the information of solvent polarization relaxation that is obtained by macroscopic dielectric dispersion experiments is not sufficient to understand dynamics of solvation of a newly created ion or dipole. In solvation, a significant contribution comes from intermediate wave vector processes which depend critically on the short range (nearest‐neighbor) spatial and orientational order that are present in a dense, dipolar liquid. An analytic expression is obtained for the time dependent solvation energy that depends, in addition to the translational and rotational diffusion coefficients of the liquid, on the ratio of solute–solvent molecular sizes and on the microscopic structure of the polar liquid. Mean spherical approximation (MSA) theory is used to obtain numerical results for polarization relaxation, for wave vector and frequency dependent dielectric function and for time dependent solvation energy. We find that in the absence of translational contribution, the solvation of an ion is, in general, nonexponential. In this case, the short time decay is dominated by the longitudinal relaxation time but the long time decay is dominated by much slower large wave vector processes involving nearest‐neighbor molecules. The presence of a significant translational contribution drastically alters the decay behavior. Now, the long‐time behavior is given by the longitudinal relaxation time constant and the short time dynamics is controlled by the large wave vector processes. Thus, although the continuum model itself is conceptually wrong, a continuum model like result is recovered in the presence of a sizeable translational contribution. The continuum model result is also recovered in the limit of large solute to solvent size ratio. In the opposite limit of small solute size, the decay is markedly nonexponential (if the translational contribution is not very large) and a complete breakdown of the continuum model takes place. The significance of these results is discussed.

A soluble random-matrix model for relaxation in quantum systems

We study the relaxation of a degenerate two-level system interacting with a heat bath, assuming a random-matrix model for the system-bath interaction. For times larger than the duration of a collision and smaller than the Poincaré recurrence time, the survival probability of still finding the system at timet in the same state in which it was prepared att=0 is exactly calculated.