Finite Element Method For A Nonlocal Problem Of Kirchhoff Type

The nonlocal term in the nonlinear equations of Kirchhoff type causes difficulties when the equation is solved numerically by using the Newton-Raphson method. This is because the Jacobian of the Newton-Raphson method is full. In this article, the finite element system is replaced by an equivalent system for which the Jacobian is sparse. We derive quasi-optimal error estimates for the finite element method and demonstrate the results with numerical experiments.

Nonlinear viscoelasticity of freestanding and polymer-anchored vertically aligned carbon nanotube foams

Vertical arrays of carbon nanotubes (VACNTs) show unique mechanical behavior in compression, with a highly nonlinear response similar to that of open cell foams and the ability to recover large deformations. Here, we study the viscoelastic response of both freestanding VACNT arrays and sandwich structures composed of a VACNT array partially embedded between two layers of poly(dimethylsiloxane) (PDMS) and bucky paper. The VACNTs tested are similar to 2 mm thick foams grown via an injection chemical vapor deposition method. Both freestanding and sandwich structures exhibit a time-dependent behavior under compression. A power-law function of time is used to describe the main features observed in creep and stress-relaxation tests. The power-law exponents show nonlinear viscoelastic behavior in which the rate of creep is dependent upon the stress level and the rate of stress relaxation is dependent upon the strain level. The results show a marginal effect of the thin PDMS/bucky paper layers on the viscoelastic responses. At high strain levels (epsilon - 0.8), the peak stress for the anchored CNTs reaches similar to 45 MPa, whereas it is only similar to 15MPa for freestanding CNTs, suggesting a large effect of PDMS on the structural response of the sandwich structures. (C) 2012 American Institute of Physics. http://dx.doi.org/10.1063/1.3699184]

Wrinkled and slack membranes: nonlinear 3D elasticity solutions via smooth DMS-FEM and experiment

The smooth DMS-FEM, recently proposed by the authors, is extended and applied to the geometrically nonlinear and ill-posed problem of a deformed and wrinkled/slack membrane. A key feature of this work is that three-dimensional nonlinear elasticity equations corresponding to linear momentum balance, without any dimensional reduction and the associated approximations, directly serve as the membrane governing equations. Domain discretization is performed with triangular prism elements and the higher order (C1 or more) interelement continuity of the shape functions ensures that the errors arising from possible jumps in the first derivatives of the conventional C0 shape functions do not propagate because the ill-conditioned tangent stiffness matrices are iteratively inverted. The present scheme employs no regularization and exhibits little sensitivity to h-refinement. Although the numerically computed deformed membrane profiles do show some sensitivity to initial imperfections (nonplanarity) in the membrane profile needed to initiate transverse deformations, the overall patterns of the wrinkles and the deformed shapes appear to be less so. Finally, the deformed profiles, computed through the DMS FEM-based weak formulation, are compared with those obtained through an experiment on an ultrathin Kapton membrane, wherein wrinkles form because of the applied boundary displacement conditions. Comparisons with a reported experiment on a rectangular membrane are also provided. These exercises lend credence to the feasibility of the DMS FEM-based numerical route to computing post-wrinkled membrane shapes. Copyright (c) 2012 John Wiley & Sons, Ltd.

Linearization and Reconstruction of Non-linear Diffuse Optical Tomographic Image

Diffuse optical tomography (DOT) is one of the ways to probe highly scattering media such as tissue using low-energy near infra-red light (NIR) to reconstruct a map of the optical property distribution. The interaction of the photons in biological tissue is a non-linear process and the phton transport through the tissue is modelled using diffusion theory. The inversion problem is often solved through iterative methods based on nonlinear optimization for the minimization of a data-model misfit function. The solution of the non-linear problem can be improved by modeling and optimizing the cost functional. The cost functional is f(x) = x(T)Ax - b(T)x + c and after minimization, the cost functional reduces to Ax = b. The spatial distribution of optical parameter can be obtained by solving the above equation iteratively for x. As the problem is non-linear, ill-posed and ill-conditioned, there will be an error or correction term for x at each iteration. A linearization strategy is proposed for the solution of the nonlinear ill-posed inverse problem by linear combination of system matrix and error in solution. By propagating the error (e) information (obtained from previous iteration) to the minimization function f(x), we can rewrite the minimization function as f(x; e) = (x + e)(T) A(x + e) - b(T)(x + e) + c. The revised cost functional is f(x; e) = f(x) + e(T)Ae. The self guided spatial weighted prior (e(T)Ae) error (e, error in estimating x) information along the principal nodes facilitates a well resolved dominant solution over the region of interest. The local minimization reduces the spreading of inclusion and removes the side lobes, thereby improving the contrast, localization and resolution of reconstructed image which has not been possible with conventional linear and regularization algorithm.

Estimation of heavy rainfall during cyclonic storms from microwave observations using nonlinear approach over Indian Ocean

In this study, an effort has been made to study heavy rainfall events during cyclonic storms over Indian Ocean. This estimate is based on microwave observations from tropical rainfall measuring mission (TRMM) Microwave Imager (TMI). Regional scattering index (SI) developed for Indian region based on measurements at 19-, 21- and 85-GHz brightness temperature and polarization corrected temperature (PCT) at 85 GHz have been utilized in this study. These PCT and SI are collocated against Precipitation Radar (PR) onboard TRMM to establish a relationship between rainfall rate, PCT and SI. The retrieval technique using both linear and nonlinear regressions has been developed utilizing SI, PCT and the combination of SI and PCT. The results have been compared with the observations from PR. It was found that a nonlinear algorithm using combination of SI and PCT is more accurate than linear algorithm or nonlinear algorithm using either SI or PCT. Statistical comparison with PR exhibits the correlation coefficients (CC) of 0.68, 0.66 and 0.70, and root mean square error (RMSE) of 1.78, 1.96 and 1.68 mm/h from the observations of SI, PCT and combination of SI and PCT respectively using linear regressions. When nonlinear regression is used, the CC of 0.73, 0.71, 0.79 and RMSE of 1.64, 1.95, 1.54 mm/h are observed from the observations of SI, PCT and combination of SI and PCT, respectively. The error statistics for high rain events (above 10 mm/h) shows the CC of 0.58, 0.59, 0.60 and RMSE of 5.07, 5.47, 5.03 mm/h from the observations of SI, PCT and combination of SI and PCT, respectively, using linear regression, and on the other hand, use of nonlinear regression yields the CC of 0.66, 0.64, 0.71 and RMSE of 4.68, 5.78 and 4.02 mm/h from the observations of SI, PCT and combined SI and PCT, respectively.

Nonlinear Refraction and Two-Photon Absorption in Dense 2Bi(2)O(3)-B2O3 Glasses

High density transparent glasses (7.86 g/cc) were fabricated in the 2Bi(2)O(3)-B2O3 (BBO) system. Optical band gap of the obtained glasses was found to be 2.6eV. The refractive index measured for these glasses was 2.25 +/- 0.05 at lambda=543 nm. Nonlinear refraction and absorption studies were carried out on the BBO glasses using z-scan technique a lambda=532 nm of 10 ns pulse width. The nonlinear refractive index obtained was n(2)=12.1x10(-14) cm(2)/W and nonlinear absorption coefficient was beta=15.2 cm/GW. The n(2) and beta values of the BBO glasses were large compared to the other reported high index bismuth based oxide glass systems in the literature. These were attributed to the high density, high linear refractive index, low band gap and two photon absorption associated with these glasses. The electronic origin of large nonlinearities was discussed based on bond-orbital theory.

Nonlinear dynamic state estimation in instrumented structures with conditionally linear Gaussian substructures

Many problems of state estimation in structural dynamics permit a partitioning of system states into nonlinear and conditionally linear substructures. This enables a part of the problem to be solved exactly, using the Kalman filter, and the remainder using Monte Carlo simulations. The present study develops an algorithm that combines sequential importance sampling based particle filtering with Kalman filtering to a fairly general form of process equations and demonstrates the application of a substructuring scheme to problems of hidden state estimation in structures with local nonlinearities, response sensitivity model updating in nonlinear systems, and characterization of residual displacements in instrumented inelastic structures. The paper also theoretically demonstrates that the sampling variance associated with the substructuring scheme used does not exceed the sampling variance corresponding to the Monte Carlo filtering without substructuring. (C) 2012 Elsevier Ltd. All rights reserved.

Probing ultrafast carrier dynamics, nonlinear absorption and refraction in core-shell silicon nanowires

We investigate the relaxation dynamics of photogenerated carriers in silicon nanowires consisting of a crystalline core and a surrounding amorphous shell, using femtosecond time-resolved differential reflectivity and transmission spectroscopy at 3.15 eV and 1.57 eV photon energies. The complex behaviour of the differential transmission and reflectivity transients is the mixed contributions from the crystalline core and the amorphous silicon on the nanowire surface and the substrate where competing effects of state-filling and photoinduced absorption govern the carrier dynamics. Faster relaxation rates are observed on increasing the photogenerated carrier density. Independent experimental results on crystalline silicon-on-sapphire (SOS) help us in separating the contributions from the carrier dynamics in crystalline core and the amorphous regions in the nanowire samples. Further, single-beam z-scan nonlinear transmission experiments at 1.57 eV in both open- and close-aperture configurations yield two-photon absorption coefficient beta (similar to 3 cm/GW) and nonlinear refraction coefficient gamma (-2.5 x 10 (-aEuro parts per thousand 4) cm(2)/GW).

Large nonlinear refraction and two photon absorption in ferroelectric Bi2VO5.5 thin films

Ferroelectric c-oriented Bi2VO5.5 (BVO) thin films (thickness approximate to 300 nm) were fabricated by pulsed laser deposition on corning glass substrates. Nonlinear refractive index (n(2)) and two photon absorption coefficient (beta) were measured by Z-scan technique at 532 nm wavelength delivering pulses with 10 ns duration. Relatively large values of n(2) = 2.05 +/- 0.2 x 10(-10) cm(2)/W and beta = 9.36 +/- 0.3 cm/MW were obtained for BVO thin films. Origin of the large optical nonlinearities in BVO thin films was discussed based on bond-orbital theory of transition metal oxides. (c) 2012 Elsevier B.V. All rights reserved.

Urea-based constructs readily amplify and attenuate nonlinear optical activity in response to H-bonding and anion recognition

Urea-based molecular constructs are shown for the first time to be nonlinear optically (NLO) active in solution. We demonstrate self-assembly triggered large amplification and specific anion recognition driven attenuation of the NLO activity. This orthogonal modulation along with an excellent nonlinearity-transparency trade-off makes them attractive NLO probes for studies related to weak self-assembly and anion transportation by second harmonic microscopy.

Calculation of Heat of Adsorption of Gases and Refrigerants on Activated Carbons from Direct Measurements Fitted to the Dubinin-Astakhov Equation

The heat of adsorption of methane, ethane, carbon dioxide, R-507a and R-134a on several specimens of microporous activated carbons is derived from experimental adsorption data fitted to the Dubinin-Astakhov equation. These adsorption results are compared with literature data obtained from calorimetric measurements and from the pressure-temperature relation during isosteric heating/cooling. Because the adsorbed phase volume plays an important role, its dependence on temperature and pressure needs to be correctly assessed. In addition, for super-critical gas adsorption, the evaluation of the pseudo-saturation pressure also needs a judicious treatment. Based on the evaluation of carbon dioxide adsorption, it can be seen that sub-critical and super-critical adsorption show different temperature dependences of the isosteric heat of adsorption. The temperature and loading dependence of this property needs to be taken into account while designing practical systems. Some practical implications of these findings are enumerated.

SOLVING VOLUME INTEGRAL EQUATION USING ScaLAPACK

In this article, we investigate the performance of a volume integral equation code on BlueGene/L system. Volume integral equation (VIE) is solved for homogeneous and inhomogeneous dielectric objects for radar cross section (RCS) calculation in a highly parallel environment. Pulse basis functions and point matching technique is used to convert the volume integral equation into a set of simultaneous linear equations and is solved using parallel numerical library ScaLAPACK on IBM's distributed-memory supercomputer BlueGene/L by different number of processors to compare the speed-up and test the scalability of the code.

Lithium D-Isoascorbate Monohydrate, a New Nonlinear Optical Material

Single crystals of lithium D-isoascorbate monohydrate (LDAM), (C6H7O6Li center dot H2O), are grown by a solution growth method. The crystal structure of LDAM is solved using single crystal X-ray diffraction. The space group is orthorhombic P2(1)2(1)2(1) with four formula units per unit cell and lattice parameters a = 7.7836(3) angstrom, b = 8.7456(3) angstrom, and c = 11.0368(4) angstrom. Solubility of the material in water is determined thermogravimetrically and found to have a positive temperature coefficient of solubility. Large optical quality single crystals are subsequently grown from aqueous solution by a slow cooling method. The crystal has a bulky prismatic habit and among the prominent faces the c face appears as the only principal morphological face. The crystal exhibits a (010) cleavage. Dielectric spectroscopy reveals a nearly Debye type Cole-Cole behavior with anisotropy in relaxation. Optical transmission range is found to be from 300 to 1400 nm. The principal refractive indices of this biaxial crystal, measured using Brewster's angle method, at wavelengths 405, 543, and 632.8 nm, show high dispersion. The crystal is negative biaxial with 2V(z) = 107.8 degrees (405 nm) and belongs to the Hobden class 3. Theoretically generated type 1 and type 2 second order phase matching curves match very well with the experimental results. The second-order nonlinear coefficient d(14) was determined to be 7 x 10(-13) m/V. For the optimum phase matching direction (type 2), the second-order effective nonlinear coefficient and the walk off angle are determined to be 0.84 times d(14) and 3.5 degrees respectively. The crystal possesses high multiple surface damage thresholds of 18 GW/cm(2) and 8 GW/cm(2) at laser wavelengths 1064 and 532 nm, respectively.

Studies on lithium L-ascorbate dihydrate: An interesting chiral nonlinear optical crystal

Lithium L-Ascorbate dihydrate (LLA) is a new metal organic nonlinear optical crystal belonging to the saccharide family. Single crystals of LLA were grown from aqueous solution. Solubility of the crystal has a positive temperature coefficient facilitating growth by slow cooling. Rietveld refinement was used to confirm the phase formation. The crystal has prismatic habit with (010), (001) and (10-1) prominent faces. Thermal analysis shows that the crystal is stable up to 102 degrees C. Transmission spectrum of the crystal extends from 302 nm to 1600 nm. Dielectric spectroscopic analysis revealed Cole Cole behaviour and prominent piezoelectric resonance peaks were observed in the range of 100-200 kHz. Second harmonic generation (SHG) conversion efficiency of up to 2.56 times that of a phase matched KDP crystal was achieved when the (010) plate of LLA single crystal was rotated about the +ve c axis, by 9.4 degrees in the clockwise direction. We also observed SHG conical sections which were attributed to noncollinear phase matching. The observation of the third conical section suggests very high birefringence and large nonlinear coefficients. A detailed study of surface laser damage showed that the crystal has high multiple damage thresholds of 9.7 GW cm(-2) and 42 GW cm(-2) at 1064 nm and 532 nm radiation respectively. (C) 2012 Elsevier B.V. All rights reserved.

Nonlinear dynamic analysis of dragonfly inspired piezoelectrically driven flapping and pitching wing

The nonlinear equations for coupled elastic flapping-twisting motion of a dragonfly in- spired smart flapping wing are used for a flapping wing actuated from the root by a PZT unimorph in the piezofan configuration. Excitation by the piezoelectric harmonic force generates only the flap bending motion, which in turn, induces the elastic twist motion due to interaction between flexural and torsional vibrations modes. An unsteady aerodynamic model is used to obtain the aerodynamic forces. Numerical simulations are performed using a wing whose size is the same as the dragonfly Sympetrum Frequens wing. It is found that the value of average lift reaches to its maximum when the smart flapping wing is excited at a frequency closer to the natural frequency in torsion. Moreover, consideration of the elastic twisting of flapping wing leads to an increase in the lift force. It is also found that the flapping wing generates sufficient lift to support its own weight and carry a small pay- load. Therefore, the piezoelectrically actuated smart flapping wing based on the geometry of Sympetrum Frequens wing and undergoing flapping-twisting motions may be considered as a potential candidate for use in MAV applications.