Equivalent differential equations for nonlinear dynamical systems

A technique for obtaining approximate periodic solutions to nonlinear ordinary differential equations is investigated. The approach is based on defining an equivalent differential equation whose exact periodic solution is known. Emphasis is placed on the mathematical justification of the approach. The relationship between the differential equation error and the solution error is investigated, and, under certain conditions, bounds are obtained on the latter. The technique employed is to consider the equation governing the exact solution error as a two point boundary value problem. Among other things, the analysis indicates that if an exact periodic solution to the original system exists, it is always possible to bound the error by selecting an appropriate equivalent system.

Three equivalence criteria for minimizing the differential equation error are compared, namely, minimum mean square error, minimum mean absolute value error, and minimum maximum absolute value error. The problem is analyzed by way of example, and it is concluded that, on the average, the minimum mean square error is the most appropriate criterion to use.

A comparison is made between the use of linear and cubic auxiliary systems for obtaining approximate solutions. In the examples considered, the cubic system provides noticeable improvement over the linear system in describing periodic response.

A comparison of the present approach to some of the more classical techniques is included. It is shown that certain of the standard approaches where a solution form is assumed can yield erroneous qualitative results.

Experimental study of nonlinear switching characteristics of conventional 2×2 fused tapered couplers

The nonlinear switching characteristics of fused fiber directional couplers were studied experimentally. By using femtosecond laser pulses with pulse width of 100 fs and wavelength of about 1550 nm from a system of Ti:sapphire laser and optical parametric amplifier (OPA), the nonlinear switching properties of a null coupler and a 100% coupler were measured. The experimental results were coincident with the simulations based on nonlinear propagation equations in fiber by using super-mode theory. Nonlinear loss in fiber was also measured to get the injected power at the coupler. After deducting the nonlinear loss and input efficiency, the nonlinear switching critical peak powers for a 100% and a null fused couplers were calculated to be 9410 and 9440 W, respectively. The nonlinear loss parameter P_(N) in an expression of α_(NL)=αP/P_(N) was obtained to be P_(N)=0.23 W.

Influence of recording conditions on crossed-beam photorefractive gratings in doubly doped LiNbO3 crystals

The influence of the recording conditions, including the widths of the recording beams, the width ratio of the recording beams, and the recording angles, on the properties of crossed-beam photorefractive gratings in doubly doped LiNbO3 crystals is studied. A theoretical model that combines the band transport model with two-dimensional coupled-wave theory is proposed. The numerical calculations of the space-charge field, the intensity profiles of the diffracted beam, and the diffraction efficiency are presented. (C) 2006 Optical Society of America.

Aberration measurement of projection optics in lithographic tools based on two-beam interference theory

The degradation of image quality caused by aberrations of projection optics in lithographic tools is a serious problem in optical lithography. We propose what we believe to be a novel technique for measuring aberrations of projection optics based on two-beam interference theory. By utilizing the partial coherent imaging theory, a novel model that accurately characterizes the relative image displacement of a fine grating pattern to a large pattern induced by aberrations is derived. Both even and odd aberrations are extracted independently from the relative image displacements of the printed patterns by two-beam interference imaging of the zeroth and positive first orders. The simulation results show that by using this technique we can measure the aberrations present in the lithographic tool with higher accuracy. (c) 2006 Optical Society of America.

A study of nonlinear phenomena in the propagation of electromagnetic waves in a weakly ionized gas

This thesis is a study of nonlinear phenomena in the propagation of electromagnetic waves in a weakly ionized gas externally biased with a magnetostatic field. The present study is restricted to the nonlinear phenomena rising from the interaction of electromagnetic waves in the ionized gas. The important effects of nonlinearity are wave-form distortion leads to cross modulation of one wave by a second amplitude-modulated wave.

The nonlinear effects are assumed to be small so that a perturbation method can be used. Boltzmann’s kinetic equation with an appropriate expression for the collision term is solved by expanding the electron distribution function into spherical harmonics in velocity space. In turn, the electron convection current density and the conductivity tensors of the nonlinear ionized gas are found from the distribution function. Finally, the expression for the current density and Maxwell’s equations are employed to investigate the effects of nonlinearity on the propagation of electromagnetic waves in the ionized gas, and also on the reflection of waves from an ionized gas of semi-infinite extent.

SteelConverter and Caltech VirtualShaker: Rapid Nonlinear Cloud-Based Structural Model Conversion and Analysis

STEEL, the Caltech created nonlinear large displacement analysis software, is currently used by a large number of researchers at Caltech. However, due to its complexity, lack of visualization tools (such as pre- and post-processing capabilities) rapid creation and analysis of models using this software was difficult. SteelConverter was created as a means to facilitate model creation through the use of the industry standard finite element solver ETABS. This software allows users to create models in ETABS and intelligently convert model information such as geometry, loading, releases, fixity, etc., into a format that STEEL understands. Models that would take several days to create and verify now take several hours or less. The productivity of the researcher as well as the level of confidence in the model being analyzed is greatly increased.

It has always been a major goal of Caltech to spread the knowledge created here to other universities. However, due to the complexity of STEEL it was difficult for researchers or engineers from other universities to conduct analyses. While SteelConverter did help researchers at Caltech improve their research, sending SteelConverter and its documentation to other universities was less than ideal. Issues of version control, individual computer requirements, and the difficulty of releasing updates made a more centralized solution preferred. This is where the idea for Caltech VirtualShaker was born. Through the creation of a centralized website where users could log in, submit, analyze, and process models in the cloud, all of the major concerns associated with the utilization of SteelConverter were eliminated. Caltech VirtualShaker allows users to create profiles where defaults associated with their most commonly run models are saved, and allows them to submit multiple jobs to an online virtual server to be analyzed and post-processed. The creation of this website not only allowed for more rapid distribution of this tool, but also created a means for engineers and researchers with no access to powerful computer clusters to run computationally intensive analyses without the excessive cost of building and maintaining a computer cluster.

In order to increase confidence in the use of STEEL as an analysis system, as well as verify the conversion tools, a series of comparisons were done between STEEL and ETABS. Six models of increasing complexity, ranging from a cantilever column to a twenty-story moment frame, were analyzed to determine the ability of STEEL to accurately calculate basic model properties such as elastic stiffness and damping through a free vibration analysis as well as more complex structural properties such as overall structural capacity through a pushover analysis. These analyses showed a very strong agreement between the two softwares on every aspect of each analysis. However, these analyses also showed the ability of the STEEL analysis algorithm to converge at significantly larger drifts than ETABS when using the more computationally expensive and structurally realistic fiber hinges. Following the ETABS analysis, it was decided to repeat the comparisons in a software more capable of conducting highly nonlinear analysis, called Perform. These analyses again showed a very strong agreement between the two softwares in every aspect of each analysis through instability. However, due to some limitations in Perform, free vibration analyses for the three story one bay chevron brace frame, two bay chevron brace frame, and twenty story moment frame could not be conducted. With the current trend towards ultimate capacity analysis, the ability to use fiber based models allows engineers to gain a better understanding of a building’s behavior under these extreme load scenarios.

Following this, a final study was done on Hall’s U20 structure [1] where the structure was analyzed in all three softwares and their results compared. The pushover curves from each software were compared and the differences caused by variations in software implementation explained. From this, conclusions can be drawn on the effectiveness of each analysis tool when attempting to analyze structures through the point of geometric instability. The analyses show that while ETABS was capable of accurately determining the elastic stiffness of the model, following the onset of inelastic behavior the analysis tool failed to converge. However, for the small number of time steps the ETABS analysis was converging, its results exactly matched those of STEEL, leading to the conclusion that ETABS is not an appropriate analysis package for analyzing a structure through the point of collapse when using fiber elements throughout the model. The analyses also showed that while Perform was capable of calculating the response of the structure accurately, restrictions in the material model resulted in a pushover curve that did not match that of STEEL exactly, particularly post collapse. However, such problems could be alleviated by choosing a more simplistic material model.

Diffraction properties of volume holographic gratings for an ultrashort pulsed beam with different polarization states readout

The diffraction properties of volume holographic gratings are studied when the gratings are illuminated by an ultrashort pulsed beam with different polarization states. The developed coupled wave theory of Kogelnik is used. Considering the dispersion effect of the grating media, solutions for the diffracted and transmitted intensities, diffraction efficiencies and the bandwidths of the gratings are given in transmission volume holographic gratings and reflection volume holographic gratings. The bandwidths of the gratings are reduced by the dispersion effect of the grating media. They also have different influences on the diffraction of an ultrashort pulsed beam with different polarization states. For different values of the ratio of the spectral bandwidth of the input pulse to that of the grating, the changes of the spectral and temporal distributions of the diffracted intensities, as well as the diffraction efficiencies of the gratings are shown.

Distortion of beam shape by a rotating double-prism wide-angle laser beam scanner

We describe the rigorous results of a wide-angle laser beam scanner, obtained with the help of the vector refraction theory. Using the rigorous results, the distortion of the beam shape is discussed. The distortion to the beam varies with the different relative angles of double prisms. The scanner expands the beam in some directions while it contracts the beam in other directions. According to the conservation of energy, the distribution of the laser intensity is changed as well. (c) 2006 Society of Photo-Optical Instrumentation Engineers.

Absolute oscillator strengths for iron, copper, cadmium and gold measured by the atomic beam method

Absolute f-values for 7 transitions in the first spectra of 4 elements have been measured using the atomic beam absorption technique. The equivalent widths of the absorption lines are measured with a photoelectric scanner and the atomic beam density is determined by continuously weighing a part of it with a sensitive automatic microbalance. The complete theory is presented and corrections are calculated to cope with gas absorption by the deposit on the microbalance pan and atoms which do not stick to the pan. An additional correction for the failure of the assumption of effusive flow in the formation of the atomic beam at large densities has been measured experimentally.

The following f-values were measured:

Fe: f_λ3720 = 0.0430 ± 8%

Cu: f_λ3247 = 0.427 ± 4.5%, f_λ3274 = 0.206 ± 4.7%, f_λ2492 = 0.0037 ± 9%

Cd: f_λ3261 = 0.00190 ± 7%, f_λ2288 = 1.38 ± 12%

Au: f_λ2428 = 0.283 ± 5.3%

Comparison with other accurately measured f-values, where they exist, shows agreement within experimental errors.

A hybrid-parallel framework for the nonlinear seismic analysis of very tall buildings

FRAME3D, a program for the nonlinear seismic analysis of steel structures, has previously been used to study the collapse mechanisms of steel buildings up to 20 stories tall. The present thesis is inspired by the need to conduct similar analysis for much taller structures. It improves FRAME3D in two primary ways.

First, FRAME3D is revised to address specific nonlinear situations involving large displacement/rotation increments, the backup-subdivide algorithm, element failure, and extremely narrow joint hysteresis. The revisions result in superior convergence capabilities when modeling earthquake-induced collapse. The material model of a steel fiber is also modified to allow for post-rupture compressive strength.

Second, a parallel FRAME3D (PFRAME3D) is developed. The serial code is optimized and then parallelized. A distributed-memory divide-and-conquer approach is used for both the global direct solver and element-state updates. The result is an implicit finite-element hybrid-parallel program that takes advantage of the narrow-band nature of very tall buildings and uses nearest-neighbor-only communication patterns.

Using three structures of varied sized, PFRAME3D is shown to compute reproducible results that agree with that of the optimized 1-core version (displacement time-history response root-mean-squared errors are ~〖10〗^(-5) m) with much less wall time (e.g., a dynamic time-history collapse simulation of a 60-story building is computed in 5.69 hrs with 128 cores—a speedup of 14.7 vs. the optimized 1-core version). The maximum speedups attained are shown to increase with building height (as the total number of cores used also increases), and the parallel framework can be expected to be suitable for buildings taller than the ones presented here.

PFRAME3D is used to analyze a hypothetical 60-story steel moment-frame tube building (fundamental period of 6.16 sec) designed according to the 1994 Uniform Building Code. Dynamic pushover and time-history analyses are conducted. Multi-story shear-band collapse mechanisms are observed around mid-height of the building. The use of closely-spaced columns and deep beams is found to contribute to the building's “somewhat brittle” behavior (ductility ratio ~2.0). Overall building strength is observed to be sensitive to whether a model is fracture-capable.

Diffraction of an ultrashort pulsed beam with arbitrary polarization state from a volume holographic grating in LiNbO3 crystals

Based on a modified coupled wave theory of Kogelnik, we have studied the diffraction of an ultrashort pulsed beam with an arbitrary polarization state from a volume holographic grating in photorefractive LiNbO3 crystals. The results indicate that the diffracted intensity distributions in the spectral and temporal domains and the diffraction efficiency of the grating are both changed by the polarization state and spectral bandwidth of the input pulsed beam. A method is given of choosing the grating parameters and input conditions to obtain a large variation range of the spectral bandwidth of the diffracted pulsed beam with an appropriate diffraction efficiency. Our study presents a possibility of using a volume holographic grating recorded in anisotropic materials to shape a broadband ultrashort pulsed beam by modulating its polarization state.