Simulation of an 1857-like Mw 7.9 San Andreas Fault Earthquake and the Response of Tall Steel Moment Frame Buildings in Southern California – A Prototype Study

In 1857, an earthquake of magnitude 7.9 occurred on the San Andreas fault, starting at Parkfield and rupturing in a southeasterly direction for more than 360 km. Such a unilateral rupture produces significant directivity toward the San Fernando and Los Angeles basins. The strong shaking in the basins due to this earthquake would have had significant long-period content (2-8 s), and the objective of this study is to quantify the impact of such an earthquake on two 18-story steel moment frame building models, hypothetically located at 636 sites on a 3.5 km grid in southern California. End-to-end simulations include modeling the source and rupture of a fault at one end, numerically propagating the seismic waves through the earth structure, simulating the damage to engineered structures and estimating the economic impact at the other end using high-performance computing. In this prototype study, we use an inferred finite source model of the magnitude 7.9, 2002 Denali fault earthquake in Alaska, and map it onto the San Andreas fault with the rupture originating at Parkfield and propagating southward over a distance of 290 km. Using the spectral element seismic wave propagation code, SPECFEM3D, we simulate an 1857-like earthquake on the San Andreas fault and compute ground motions at the 636 analysis sites. Using the nonlinear structural analysis program, FRAME3D, we subsequently analyze 3-D structural models of an existing tall steel building designed using the 1982 Uniform Building Code (UBC), as well as one designed according to the 1997 UBC, subjected to the computed ground motion at each of these sites. We summarize the performance of these structural models on contour maps of peak interstory drift. We then perform an economic loss analysis for the two buildings at each site, using the Matlab Damage and Loss Analysis (MDLA) toolbox developed to implement the PEER loss-estimation methodology. The toolbox includes damage prediction and repair cost estimation for structural and non-structural components and allows for the computation of the mean and variance of building repair costs conditional on engineering demand parameters (i.e. inter-story drift ratios and peak floor accelerations). Here, we modify it to treat steel-frame high-rises, including aspects such as mechanical, electrical and plumbing systems, traction elevators, and the possibility of irreparable structural damage. We then generate contour plots of conditional mean losses for the San Fernando and the Los Angeles basins for the pre-Northridge and modern code-designed buildings, allowing for comparison of the economic effects of the updated code for the scenario event. In principle, by simulating multiple seismic events, consistent with the probabilistic seismic hazard for a building site, the same basic approach could be used to quantify the uncertain losses from future earthquakes

Production of Enhanced Beam Halos via Collective Modes and Colored Noise

We investigate how collective modes and colored noise conspire to produce a beam halo with much larger amplitude than could be generated by either phenomenon separately. The collective modes are lowest-order radial eigenmodes calculated self-consistently for a configuration corresponding to a direct-current, cylindrically symmetric, warm-fluid Kapchinskij-Vladimirskij equilibrium. The colored noise arises from unavoidable machine errors and influences the internal space-charge force. Its presence quickly launches statistically rare particles to ever-growing amplitudes by continually kicking them back into phase with the collective-mode oscillations. The halo amplitude is essentially the same for purely radial orbits as for orbits that are initially purely azimuthal; orbital angular momentum has no statistically significant impact. Factors that do have an impact include the amplitudes of the collective modes and the strength and autocorrelation time of the colored noise. The underlying dynamics ensues because the noise breaks the Kolmogorov-Arnol'd-Moser tori that otherwise would confine the beam. These tori are fragile; even very weak noise will eventually break them, though the time scale for their disintegration depends on the noise strength. Both collective modes and noise are therefore centrally important to the dynamics of halo formation in real beams.Comment: For full resolution pictures please go to http://www.nicadd.niu.edu/research/beams

Generation of angular-momentum-dominated electron beams from a photoinjector

Various projects under study require an angular-momentum-dominated electron beam generated by a photoinjector. Some of the proposals directly use the angular-momentum-dominated beams (e.g. electron cooling of heavy ions), while others require the beam to be transformed into a flat beam (e.g. possible electron injectors for light sources and linear colliders). In this paper, we report our experimental study of an angular-momentum-dominated beam produced in a photoinjector, addressing the dependencies of angular momentum on initial conditions. We also briefly discuss the removal of angular momentum. The results of the experiment, carried out at the Fermilab/NICADD Photoinjector Laboratory, are found to be in good agreement with theoretical and numerical models.Comment: 8 pages, 7 figures, submitted to Phys. Rev. ST Accel. Beam

Age-dependent differences in human brain activity using a face- and location-matching task: An fMRI study

Purpose: To evaluate the differences of cortical activation patterns in young and elderly healthy subjects for object and spatial visual processing using a face- and location-matching task. Materials and Methods: We performed a face- and a location-matching task in 15 young (mean age: 28 +/- 9 years) and 19 elderly (mean age: 71 +/- 6 years) subjects. Each experiment consisted of 7 blocks alternating between activation and control condition. For face matching, the subjects had to indicate whether two displayed faces were identical or different. For location matching, the subjects had to press a button whenever two objects had an identical position. For control condition, we used a perception task with abstract images. Functional imaging was performed on a 1.5-tesla scanner using an EPI sequence. Results: In the face-matching task, the young subjects showed bilateral (right 1 left) activation in the occipito-temporal pathway (occipital gyrus, inferior and middle temporal gyrus). Predominantly right hemispheric activations were found in the fusiform gyrus, the right dorsolateral prefrontal cortex (inferior and middle frontal gyrus) and the superior parietal gyrus. In the elderly subjects, the activated areas in the right fronto-lateral cortex increased. An additional activated area could be found in the medial frontal gyrus (right > left). In the location-matching task, young subjects presented increased bilateral (right > left) activation in the superior parietal lobe and precuneus compared with face matching. The activations in the occipito-temporal pathway, in the right fronto-lateral cortex and the fusiform gyrus were similar to the activations found in the face-matching task. In the elderly subjects, we detected similar activation patterns compared to the young subjects with additional activations in the medial frontal gyrus. Conclusion: Activation patterns for object-based and spatial visual processing were established in the young and elderly healthy subjects. Differences between the elderly and young subjects could be evaluated, especially by using a face-matching task. Copyright (c) 2007 S. Karger AG, Basel

Chaos and the continuum limit in nonneutral plasmas and charged particle beams

This paper examines discreteness effects in nearly collisionless N-body systems of charged particles interacting via an unscreened r^-2 force, allowing for bulk potentials admitting both regular and chaotic orbits. Both for ensembles and individual orbits, as N increases there is a smooth convergence towards a continuum limit. Discreteness effects are well modeled by Gaussian white noise with relaxation time t_R = const * (N/log L)t_D, with L the Coulomb logarithm and t_D the dynamical time scale. Discreteness effects accelerate emittance growth for initially localised clumps. However, even allowing for discreteness effects one can distinguish between orbits which, in the continuum limit, feel a regular potential, so that emittance grows as a power law in time, and chaotic orbits, where emittance grows exponentially. For sufficiently large N, one can distinguish two different `kinds' of chaos. Short range microchaos, associated with close encounters between charges, is a generic feature, yielding large positive Lyapunov exponents X_N which do not decrease with increasing N even if the bulk potential is integrable. Alternatively, there is the possibility of larger scale macrochaos, characterised by smaller Lyapunov exponents X_S, which is present only if the bulk potential is chaotic. Conventional computations of Lyapunov exponents probe X_N, leading to the oxymoronic conclusion that N-body orbits which look nearly regular and have sharply peaked Fourier spectra are `very chaotic.' However, the `range' of the microchaos, set by the typical interparticle spacing, decreases as N increases, so that, for large N, this microchaos, albeit very strong, is largely irrelevant macroscopically. A more careful numerical analysis allows one to estimate both X_N and X_S.Comment: 13 pages plus 17 figure