A characteristic based numerical method with tracking for nonlinear wave equations

Many physical problems can be modeled by scalar, first-order, nonlinear, hyperbolic, partial differential equations (PDEs). The solutions to these PDEs often contain shock and rarefaction waves, where the solution becomes discontinuous or has a discontinuous derivative. One can encounter difficulties using traditional finite difference methods to solve these equations. In this paper, we introduce a numerical method for solving first-order scalar wave equations. The method involves solving ordinary differential equations (ODEs) to advance the solution along the characteristics and to propagate the characteristics in time. Shocks are created when characteristics cross, and the shocks are then propagated by applying analytical jump conditions. New characteristics are inserted in spreading rarefaction fans. New characteristics are also inserted when values on adjacent characteristics lie on opposite sides of an inflection point of a nonconvex flux function, Solutions along characteristics are propagated using a standard fourth-order Runge-Kutta ODE solver. Shocks waves are kept perfectly sharp. In addition, shock locations and velocities are determined without analyzing smeared profiles or taking numerical derivatives. In order to test the numerical method, we study analytically a particular class of nonlinear hyperbolic PDEs, deriving closed form solutions for certain special initial data. We also find bounded, smooth, self-similar solutions using group theoretic methods. The numerical method is validated against these analytical results. In addition, we compare the errors in our method with those using the Lax-Wendroff method for both convex and nonconvex flux functions. Finally, we apply the method to solve a PDE with a convex flux function describing the development of a thin liquid film on a horizontally rotating disk and a PDE with a nonconvex flux function, arising in a problem concerning flow in an underground reservoir.

Optimal management of beaver population using a reduced-order distributed parameter model and single network adaptive critics

Beavers are often found to be in conflict with human interests by creating nuisances like building dams on flowing water (leading to flooding), blocking irrigation canals, cutting down timbers, etc. At the same time they contribute to raising water tables, increased vegetation, etc. Consequently, maintaining an optimal beaver population is beneficial. Because of their diffusion externality (due to migratory nature), strategies based on lumped parameter models are often ineffective. Using a distributed parameter model for beaver population that accounts for their spatial and temporal behavior, an optimal control (trapping) strategy is presented in this paper that leads to a desired distribution of the animal density in a region in the long run. The optimal control solution presented, imbeds the solution for a large number of initial conditions (i.e., it has a feedback form), which is otherwise nontrivial to obtain. The solution obtained can be used in real-time by a nonexpert in control theory since it involves only using the neural networks trained offline. Proper orthogonal decomposition-based basis function design followed by their use in a Galerkin projection has been incorporated in the solution process as a model reduction technique. Optimal solutions are obtained through a "single network adaptive critic" (SNAC) neural-network architecture.

An investigation of the first-order spin-state transition in Fe(Phen)(2)(NCS)(2) EXAFS and infrared spectroscopy

The structure of Fe(Phen)(2)(NCS)(2) has been examined across the first-order spin-state transition by EXAFS with full multiple scattering analysis. The EXAFS data at 298 K can be satisfactorily assigned to the high-spin state, but the analysis of the low-temperature data at 90 K is not entirely unequivocal, although consistent with the predominant presence of the low-spin state. That some proportion of the high-spin state remains at low temperatures, well below the first-order transition, is clearly evidenced in the infrared spectra, suggesting possible sublattice ordering.

A study of partial discharge characteristics in oil impregnated pressboard insulation

An attempt has been made to study the effect of time and test procedure on the behaviour of partial discharge (PD) pulses causing failure of oil-pressboard system under power frequency voltages using circular disc shaped samples and uniform field electrodes. Weibull statistics have been used to handle the large amount of PD data. The PD phenomena has been found to be stress and time dependent. On the basis of stress level, three different regions are identified and in one of the regions, the rate of deterioration of the sample is at a maximum. The work presents some interesting features of Weibull parameters as related to the condition of insulation studied in addition to its usual PD characteristics

Reduced Order Suboptimal Control Design for a Class of Nonlinear Distributed Parameter Systems Using POD and /spl thetas/-D Techniques

A new computational tool is presented in this paper for suboptimal control design of a class of nonlinear distributed parameter systems. First proper orthogonal decomposition based problem-oriented basis functions are designed, which are then used in a Galerkin projection to come up with a low-order lumped parameter approximation. Next, a suboptimal controller is designed using the emerging /spl thetas/-D technique for lumped parameter systems. This time domain sub-optimal control solution is then mapped back to the distributed domain using the same basis functions, which essentially leads to a closed form solution for the controller in a state feedback form. Numerical results for a real-life nonlinear temperature control problem indicate that the proposed method holds promise as a good suboptimal control design technique for distributed parameter systems.

Integral treatment for the representation of thermodynamic properties in multicomponent systems using interaction parameters

The partial thermodynamic functions of the solvent component of a ternary system have been deduced in terms of the interaction parameters by integration of several series which emerge from the Maclaurin infinite series based on the integral property of the system and subjected to appropriate boundary conditions. The series integration shows that the resulting partial functions are suitable for interpreting the thermodynamic properties of the system and are independent of compositional paths. In the present analysis, the higher order terms of these series are found to make insignificant contributions.

FIR system identification based on subspaces of a higher order cumulant matrix

A new linear algebraic approach for identification of a nonminimum phase FIR system of known order using only higher order (>2) cumulants of the output process is proposed. It is first shown that a matrix formed from a set of cumulants of arbitrary order can be expressed as a product of structured matrices. The subspaces of this matrix are then used to obtain the parameters of the FIR system using a set of linear equations. Theoretical analysis and numerical simulation studies are presented to characterize the performance of the proposed methods.

Brownian dynamics simulation of dense binary colloidal mixtures. II. Translational and bond-orientational order

We report the Brownian dynamics simulation results on the translational and bond-angle-orientational correlations for charged colloidal binary suspensions as the interparticle interactions are increased to form a crystalline (for a volume fraction phi = 0.2) or a glassy (phi = 0.3) state. The translational order is quantified in terms of the two- and four-point density autocorrelation functions whose comparisons show that there is no growing correlation length near the glass transition. The nearest-neighbor orientational order is determined in terms of the quadratic rotational invariant Q(l) and the bond-orientational correlation functions g(l)(t). The l dependence of Q(l) indicates that icosahedral (l = 6) order predominates at the cost of the cubic order (l = 4) near the glass as well as the crystal transition. The density and orientational correlation functions for a supercooled liquid freezing towards a glass fit well to the streched-exponential form exp[-(t/tau)(beta)]. The average relaxation times extracted from the fitted stretched-exponential functions as a function of effective temperatures T* obey the Arrhenius law for liquids freezing to a crystal whereas these obey the Vogel-Tamman-Fulcher law exp[AT(0)*/(T* - T-0*)] for supercooled Liquids tending towards a glassy state. The value of the parameter A suggests that the colloidal suspensions are ''fragile'' glass formers like the organic and molecular liquids.

First-order hyperpolarizabilities and pKa of weak organic acids in protic solvents are lineraly related

Second-order nonlinearities (beta) of five weak organic acids in protic solvents have been measured by the double-quantum Rayleigh scattering (DRS) technique. beta is found to bear a linear relationship to the pK(a) of these compounds in those solvents. A direct implication of this observation is that the DRS technique can be used to determine the pK(a) of weak organic acids in any solvent.

First-order hyperpolarizabilities of octupolar aromatic molecules: symmetrically substituted triazines

The first hyperpolarizabilities of some symmetrically substituted triazines have been measured and compared with those of the corresponding symmetrically substituted benzenes. The octupolar triazines have higher quadratic polarizabilities than the corresponding octupolar benzenes. The triazine ring seems to be a better central acceptor than the benzene ring, but if it acts as a donor as in sym-triphenyl triazine, the nonlinearity improves further.

First-Order Hyperpolarizabilities of Sulfophthalein Dyes

In this paper we report the first hyperpolarizabilities (beta) of 12, sulfophthalein dyes. Since these dyes are ionic in nature, their second-order nonlinearities were measured by the hyper-Rayleigh scattering technique in solution. The measured beta values are large and highly solvent dependent. Inclusion of solvent polarity in ab initio estimates of static second-order polarizability does not fully account for the experimental beta values. Contributions from the dissociated forms of the dye in different solvents seem to play an important role in enhancing beta in these systems.

Highly Efficient and Thermally Stable Second-Order Nonlinear Optical Chromophores and Electrooptic Polymers

Organic polymeric electro-optic (E-O) materials have attracted significant attention because of their potential use as fast and efficient components of integrated photonic devices (1,2). However, the practical application of these materials in optical devices is somewhat limited by the stringent material requirements imposed by the device design, fabrication processes and operating environments. Among the various material requirements, the most notable ones are large electro-optic coefficients (r(33)) and high thermal stability (3). The design of poled polymeric materials with high electro-optic activity (r(33)) involves the optimization of the percent incorporation of efficient (large beta mu) second order nonlinear optical (NLO) chromophores into the polymer matrices and the effective creation of poling-induced non-centrosymmetric structures. The factors that affect the material stability are a) the inherent thermal stability of the NLO chromophores, b) the chemical stability of the NLO chromophores during the polymer processing conditions, and c) the long-term dipolar alignment stability at high temperatures. Although considerable progress has been made in achieving these properties (4), organic polymeric materials suitable for practical E-O device applications are yet to be developed. This chapter highlights some of our approaches in the optimization of molecular and material nonlinear optical and thermal properties.

Nonsimilar mixed convection flow of a non-Newtonian fluid past a vertical wedge

The flow and heat transfer characteristics of a second-order fluid over a vertical wedge with buoyancy forces have been analysed. The coupled nonlinear partial differential equations governing the nonsimilar mixed convection flow have been solved numerically using Keller box method. The effects of the buoyancy parameter, viscoelastic parameter, mass transfer parameter, pressure gradient parameter, Prandtl number and viscous dissipation parameter on the skin friction and heat transfer have been examined in detail. Particular cases of the present results match exactly with those available in the literature.

Vibration and Stability of Fluid Conveying Pipes with Stochastic Parameters

Flexible cantilever pipes conveying fluids with high velocity are analysed for their dynamic response and stability behaviour. The Young's modulus and mass per unit length of the pipe material have a stochastic distribution. The stochastic fields, that model the fluctuations of Young's modulus and mass density are characterized through their respective means, variances and autocorrelation functions or their equivalent power spectral density functions. The stochastic non self-adjoint partial differential equation is solved for the moments of characteristic values, by treating the point fluctuations to be stochastic perturbations. The second-order statistics of vibration frequencies and mode shapes are obtained. The critical flow velocity is-first evaluated using the averaged eigenvalue equation. Through the eigenvalue equation, the statistics of vibration frequencies are transformed to yield critical flow velocity statistics. Expressions for the bounds of eigenvalues are obtained, which in turn yield the corresponding bounds for critical flow velocities.

Purification and partial characterization of acyl carrier proteins from developing oil seeds of pisa (Actinodaphne hookeri) and ground nut (Arachis hypogaea)

Acyl carrier proteins (ACP) were purified to homogeneity in the active form from developing seeds of pisa (Actinodaphne hookeri) which synthesizes exclusively trilaurin and from ground nut (Arachis hypogaea) which synthesizes triacylglycerols containing long chain fatty acids. Two major isoforms of ACPs were purified from developing pisa seeds using DEAE-cellulose, Superose-6 FPLC and C-4 reversed phase HPLC chromatographic methods. In contrast, only a single form of ACP was present in ground nut seeds which was purified by anion-exchange and activated thiol-Sepharose 4B affinity chromatography. The two isoforms of ACPs from pisa showed nearly the same specific activity of 6,706 and 7,175 pmol per min per mg protein while ground nut ACP showed a specific activity of 3,893 pmol per min per mg protein when assayed using E. coli acyl-ACP synthetase and [1-C-14]palmitic acid. When compared with E. coli ACP, the purified ACPs from both the seeds showed considerable difference in their mobility in native PAGE, but showed similar mobility in SDS-PAGE under reducing conditions. In the absence of reducing agents formation of dimers was quite prominent. The ACPs from both the seed sources were acid- and heat-stable. The major isoform of pisa seed ACP and the ground nut ACP contain 91 amino acids with M(r) 11,616 and 1,228 respectively. However, there is significant variation in their amino acid composition. A comparision of the amino acid sequence in the N-terminal region of pisa and ground nut seed ACPs showed considerable homology between themselves and with other plant ACPs but not with E. coli ACP.