A smoothing technique for discrete delta functions with application to immersed boundary method in moving boundary simulations.

The effects of complex boundary conditions on flows are represented by a volume force in the immersed boundary methods. The problem with this representation is that the volume force exhibits non-physical oscillations in moving boundary simulations. A smoothing technique for discrete delta functions has been developed in this paper to suppress the non-physical oscillations in the volume forces. We have found that the non-physical oscillations are mainly due to the fact that the derivatives of the regular discrete delta functions do not satisfy certain moment conditions. It has been shown that the smoothed discrete delta functions constructed in this paper have one-order higher derivative than the regular ones. Moreover, not only the smoothed discrete delta functions satisfy the first two discrete moment conditions, but also their derivatives satisfy one-order higher moment condition than the regular ones. The smoothed discrete delta functions are tested by three test cases: a one-dimensional heat equation with a moving singular force, a two-dimensional flow past an oscillating cylinder, and the vortex-induced vibration of a cylinder. The numerical examples in these cases demonstrate that the smoothed discrete delta functions can effectively suppress the non-physical oscillations in the volume forces and improve the accuracy of the immersed boundary method with direct forcing in moving boundary simulations.

Numerical simulation of condensation of bubbles under microgravity conditions by moving mesh method in the double-staggered grid

A numerical 2D method for simulation of two-phase flows including phase change under microgravity conditions is presented in this paper, with a level set method being coupled with the moving mesh method in the double-staggered grid systems. When the grid lines bend very much in a curvilinear grid, great errors may be generated by using the collocated grid or the staggered grid. So the double-staggered grid was adopted in this paper. The level set method is used to track the liquid-vapor interface. The numerical analysis is fulfilled by solving the Navier-Stokes equations using the SIMPLER method, and the surface tension force is modeled by a continuum surface force approximation. A comparison of the numerical results obtained with different numerical strategies shows that the double-staggered grid moving-mesh method presented in this paper is more accurate than that used previously in the collocated grid system. Based on the method presented in this paper, the condensation of a single bubble in the cold water under different level of gravity is simulated. The results show that the condensation process under the normal gravity condition is different from the condensation process under microgravity conditions. The whole condensation time is much longer under the normal gravity than under the microgravity conditions.

Photographic database development of Nigerian fisheries: a complementary source of information documentation

The paper traces the history of the different documentation media used for information dissemination. Such early media are clay tablets, papyrus, and vellum or parchment codex. The invention of printing however revolutionized the information industry, enabling the production of books in multiple copies. Photography came into documentation mainly to preserve rare materials and those that easily deteriorate. This paper reports the efforts of National Institute for Freshwater Fisheries Research (NIFFR) and Kainji Lake Fisheries Promotion Project (KLFPPP), Nigeria, to develop an Object Oriented Database (OOD) using photographs. The photographs are stored in digitized form on commercial computers, using the program ACDSee 32 for classification, description and retrieval. Specifically the paper focuses on photographs in fisheries as visual communication and expression. Presently, the database contains photo documents about the following aspects of Kainji Lake fisheries: fishing gears and crafts, fish preservation methods

Studies on the nutritional value of cassava (Manihot utillisima) peels as energy source in the diet of Oreochromis niloticus fry fed in indoor glass aquaria

This study, though, has as its core objective cost reduction in aquaculture nutrition was equally designed to investigate the value of the peels of cassava (Manihot utillisima) as energy source in the diet of Oreochromis niloticus fry. Three levels of cassava peels diet and a control (100% yellow maize in the carbohydrate mixture) was prepared and tested on O. niloticus fry for ten (10) weeks. The fry with mean weight of 0.32g were grouped fifteen (15) in each of the glass aquaria measuring 60x30x30cm with a maximum capacity of 52 litres of water. The fry were fed twice daily at 10% biomass. Weekly, the fry were weighed to determine the weight increment or otherwise and the quality of feed adjusted accordingly. Water quality parameters like temperature, pH and dissolved oxygen (D.0) were monitored and found to be at desirable level. DT 3 (97 % cassava peels and 3% yellow maize) in the carbohydrate mixture gave the best growth performance. The fry fed, this diet gained mean weight of 1.18g for the period of the experiment. However, the poorest performance in terms of growth was from fry fed the control diet (100%yellow maize in the carbohydrate mixture) fry fed this diet gained mean weight of 0.80 for the duration of the experiment. Analysis of the various growth indices like SGR, PER, FCR and NPU shows that DT3 was the overall best diet with an SGR value of2.40 and FCR of 43.83. However, DT 1 (70% cassava peels and 30% yellow maize) gave the poorest SGR of 1.61 and FCR of 67.58. The difference in weight gain among the fry fed the three levels of cassava peels diet and the control was not statically significant (P>0.05)

In search of slow-moving ionizing massive particles

Many particles proposed by theories, such as GUT monopoles, nuclearites and 1/5 charge superstring particles, can be categorized as Slow-moving, Ionizing, Massive Particles (SIMPs).

Detailed calculations of the signal-to-noise ratios in vanous acoustic and mechanical methods for detecting such SIMPs are presented. It is shown that the previous belief that such methods are intrinsically prohibited by the thermal noise is incorrect, and that ways to solve the thermal noise problem are already within the reach of today's technology. In fact, many running and finished gravitational wave detection ( GWD) experiments are already sensitive to certain SIMPs. As an example, a published GWD result is used to obtain a flux limit for nuclearites.

The result of a search using a scintillator array on Earth's surface is reported. A flux limit of 4.7 x 10^(-12) cm^(-2)sr^(-1)s^(-1) (90% c.l.) is set for any SIMP with 2.7 x 10^(-4) less than β less than 5 x 10^(-3) and ionization greater than 1/3 of minimum ionizing muons. Although this limit is above the limits from underground experiments for typical supermassive particles (10^(16)GeV), it is a new limit in certain β and ionization regions for less massive ones (~10^9 GeV) not able to penetrate deep underground, and implies a stringent limit on the fraction of the dark matter that can be composed of massive electrically and/ or magnetically charged particles.

The prospect of the future SIMP search in the MACRO detector is discussed. The special problem of SIMP trigger is examined and a circuit proposed, which may solve most of the problems of the previous ones proposed or used by others and may even enable MACRO to detect certain SIMP species with β as low as the orbital velocity around the earth.

Determination of source finiteness and depth of large earthquakes

In this thesis, a method to retrieve the source finiteness, depth of faulting, and the mechanisms of large earthquakes from long-period surface waves is developed and applied to several recent large events.

In Chapter 1, source finiteness parameters of eleven large earthquakes were determined from long-period Rayleigh waves recorded at IDA and GDSN stations. The basic data set is the seismic spectra of periods from 150 to 300 sec. Two simple models of source finiteness are studied. The first model is a point source with finite duration. In the determination of the duration or source-process times, we used Furumoto's phase method and a linear inversion method, in which we simultaneously inverted the spectra and determined the source-process time that minimizes the error in the inversion. These two methods yielded consistent results. The second model is the finite fault model. Source finiteness of large shallow earthquakes with rupture on a fault plane with a large aspect ratio was modeled with the source-finiteness function introduced by Ben-Menahem. The spectra were inverted to find the extent and direction of the rupture of the earthquake that minimize the error in the inversion. This method is applied to the 1977 Sumbawa, Indonesia, 1979 Colombia-Ecuador, 1983 Akita-Oki, Japan, 1985 Valparaiso, Chile, and 1985 Michoacan, Mexico earthquakes. The method yielded results consistent with the rupture extent inferred from the aftershock area of these earthquakes.

In Chapter 2, the depths and source mechanisms of nine large shallow earthquakes were determined. We inverted the data set of complex source spectra for a moment tensor (linear) or a double couple (nonlinear). By solving a least-squares problem, we obtained the centroid depth or the extent of the distributed source for each earthquake. The depths and source mechanisms of large shallow earthquakes determined from long-period Rayleigh waves depend on the models of source finiteness, wave propagation, and the excitation. We tested various models of the source finiteness, Q, the group velocity, and the excitation in the determination of earthquake depths.

The depth estimates obtained using the Q model of Dziewonski and Steim (1982) and the excitation functions computed for the average ocean model of Regan and Anderson (1984) are considered most reasonable. Dziewonski and Steim's Q model represents a good global average of Q determined over a period range of the Rayleigh waves used in this study. Since most of the earthquakes studied here occurred in subduction zones Regan and Anderson's average ocean model is considered most appropriate.

Our depth estimates are in general consistent with the Harvard CMT solutions. The centroid depths and their 90 % confidence intervals (numbers in the parentheses) determined by the Student's t test are: Colombia-Ecuador earthquake (12 December 1979), d = 11 km, (9, 24) km; Santa Cruz Is. earthquake (17 July 1980), d = 36 km, (18, 46) km; Samoa earthquake (1 September 1981), d = 15 km, (9, 26) km; Playa Azul, Mexico earthquake (25 October 1981), d = 41 km, (28, 49) km; El Salvador earthquake (19 June 1982), d = 49 km, (41, 55) km; New Ireland earthquake (18 March 1983), d = 75 km, (72, 79) km; Chagos Bank earthquake (30 November 1983), d = 31 km, (16, 41) km; Valparaiso, Chile earthquake (3 March 1985), d = 44 km, (15, 54) km; Michoacan, Mexico earthquake (19 September 1985), d = 24 km, (12, 34) km.

In Chapter 3, the vertical extent of faulting of the 1983 Akita-Oki, and 1977 Sumbawa, Indonesia earthquakes are determined from fundamental and overtone Rayleigh waves. Using fundamental Rayleigh waves, the depths are determined from the moment tensor inversion and fault inversion. The observed overtone Rayleigh waves are compared to the synthetic overtone seismograms to estimate the depth of faulting of these earthquakes. The depths obtained from overtone Rayleigh waves are consistent with the depths determined from fundamental Rayleigh waves for the two earthquakes. Appendix B gives the observed seismograms of fundamental and overtone Rayleigh waves for eleven large earthquakes.

Two new integral transforms and their applications

This thesis is in two parts. In Part I the independent variable θ in the trigonometric form of Legendre's equation is extended to the range ( -∞, ∞). The associated spectral representation is an infinite integral transform whose kernel is the analytic continuation of the associated Legendre function of the second kind into the complex θ-plane. This new transform is applied to the problems of waves on a spherical shell, heat flow on a spherical shell, and the gravitational potential of a sphere. In each case the resulting alternative representation of the solution is more suited to direct physical interpretation than the standard forms.

In Part II separation of variables is applied to the initial-value problem of the propagation of acoustic waves in an underwater sound channel. The Epstein symmetric profile is taken to describe the variation of sound with depth. The spectral representation associated with the separated depth equation is found to contain an integral and a series. A point source is assumed to be located in the channel. The nature of the disturbance at a point in the vicinity of the channel far removed from the source is investigated.

Wave propagation in an inhomogeneous medium

This paper is in two parts. In the first part we give a qualitative study of wave propagation in an inhomogeneous medium principally by geometrical optics and ray theory. The inhomogeneity is represented by a sound-speed profile which is dependent upon one coordinate, namely the depth; and we discuss the general characteristics of wave propagation which result from a source placed on the sound channel axis. We show that our mathematical model of the sound- speed in the ocean actually predicts some of the behavior of the observed physical phenomena in the underwater sound channel. Using ray theoretic techniques we investigate the implications of our profile on the following characteristics of SOFAR propagation: (i) the sound energy traveling further away from the axis takes less time to travel from source to receiver than sound energy traveling closer to the axis, (ii) the focusing of sound energy in the sound channel at certain ranges, (iii) the overall ray picture in the sound channel.

In the second part a more penetrating quantitative study is done by means of analytical techniques on the governing equations. We study the transient problem for the Epstein profile by employing a double transform to formally derive an integral representation for the acoustic pressure amplitude, and from this representation we obtain several alternative representations. We study the case where both source and receiver are on the channel axis and greatly separated. In particular we verify some of the earlier results derived by ray theory and obtain asymptotic results for the acoustic pressure in the far-field.

An investigation into the value of yam (Dioscorea rotundata) peels as energy source in the diet of Oreochromis niloticus fry fed in indoor glass aquaria

The study was designed to investigate the value of the peels of yam (Dioscorea rotundata) as energy source in the diet of Oreochromis niloticus fry and to investigate the level of inclusion of this peels that will give optimum growth performance. Four diets, three levels of yam peels and a control, was prepared and tested on O. niloticus fry (mean weight of 0.27g) for ten weeks. Fifteen (15) O. niloticus fry were grouped in each of the glass aquaria, measuring 60x30x3Ocm and with a maximum capacity of 52 liters of water. The fry were fed twice daily at 10% biomass. The fry were weighed weekly to determine weight increment or otherwise and the quality of feed was adjusted accordingly. DTl (70% yam peels and 30% yellow maize) in the carbohydrate mixture gave the best performance. The fry fed this diet, gained a mean weight of 1.20g for the period of the experiment. The poorest performance in terms of growth was from fry fed the control diet (100% yellow maize in the carbohydrate mixture). Fry fed this diet gained mean weight of 0.80g for the duration of the experiment. Analysis of the various growth indices like SGR, PER, FCR and NPU shows that DTl was the overall best diet with an SGR value of I. 92 and FCR of 54.10. The difference in weight gain by fry fed the three levels of yam peels diet and the control diet (100% yellow maize) was not statistically significant (P>0.05)

Policy analysis of shoreline restoration options on private shorelines of Puget Sound

Puget Sound shorelines have historically provided a diversity of habitats that support a variety of aquatic resources throughout the region. These valued natural resources are iconic to the region and remain central to both the economic vitality and community appreciation of Puget Sound. Deterioration of upland and nearshore shoreline habitats, have placed severe stress on many aquatic resources within the region (PSAT, 2007). Since a majority of Washington State shorelines are privately owned, regulatory authority to legislate restoration on private property is limited in scope and frequency. Washington States’ Shoreline Management Act (RCW 90.58) requires local jurisdictions to plan for appropriate future shoreline uses. Under the Act, future development can be regulated to protect existing ecological functions, but lost functions cannot be restored without purchase or compensation of restored areas. Therefore, questions remains as to the ecological resilience of the region when considering cumulative effect of existing/ongoing shoreline development constrained by limited shoreline restoration opportunities. In light of these questions, this analysis will explore opportunities to promote restoration on privately owned shorelines within Puget Sound. These efforts are intended to promote more efficient ecosystem management and improve ecosystem-wide ecological functions. From an economics perspective, results of past shoreline management can generally be characterized as both market and government failure in effectively protecting the publics’ interest in maintaining healthy shoreline resources. Therefore coastal development has proceeded in spite of negative externalities and market imbalances resulting in inefficient resource management driven by the individual ambitions of private shoreline property owners to develop their property to their highest and best use. Federally derived property rights will protect continuation of existing uses along privately owned shorelines; therefore, a fundamental challenge remains in sustainable management of existing shoreline resources while also restoring ecological functions lost to past mistakes in an effort to increase the ecologic resiliency within the region. (PDF contains 5 pages)

Physical circulation measurements in Pamlico Sound, used to locate new oyster sanctuaries

The population of eastern oyster, C. virginica, has declined over the last century on most areas of the east and gulf coasts. North Carolina’s restoration efforts depend on the construction of subtidal oyster reefs to be used as broodstock sanctuaries in Pamlico Sound, NC. Successful restoration of the oyster population requires several thriving reefs connected as a meta-population. C. virginica has a 2-3 week larval stage, during which it gradually settles through the water column. Larvae that can travel from one reef to another during that stage form the basis of a meta-population. (PDF contains 3 pages)

Time and space-resolved measurement of a gas-puff laser-plasma x-ray source

The dynamic interaction processes between a nano-second laser pulse and a gas-puff target, such as those of plasma formation, laser heating, and x-ray emission, have been investigated quantitatively. Time and space-resolved x-ray and optical measurement techniques were used in order to investigate time-resolved laser absorption and subsequent x-ray generation. Efficient absorption of the incident laser energy into the gas-puff target of 17%, 12%, 38%, and 91% for neon, argon, krypton, and xenon, respectively, was shown experimentally. It was found that the laser absorption starts and, simultaneously, soft x-ray emission occurs. The soft x-ray lasts much longer than the laser pulse due to the recombination. Temporal evolution of the soft x-ray emission region was analyzed by comparing the experimental results to the results of the model calculation, in which the laser light propagation through a gas-puff plasma was taken into account. (C) 2003 American Institute of Physics.

Eliciting and using stakeholder input to guide Puget Sound nearshore restoration decisionmaking

The goal of the Puget Sound Nearshore Ecosystem Restoration Project (PSNERP) is to improve system-wide functionality of nearshorei ecosystem processes. To achieve that goal, PSNERP plans to strategically restore nearshore sites throughout Puget Sound. PSNERP scientists are assessing changes to the nearshore, and will recommend an environmentally strategic restoration portfolio. Yet, PSNERP also needs stakeholder input to design a socially strategic portfolio. This research investigates the values and preferences of stakeholders in the Whidbey Sub-Basin of Puget Sound to help PSNERP be both socially and environmentally strategic. This investigation may be repeated in the six other Puget Sound Sub-Basins. The results will guide restoration portfolio design and future stakeholder involvement activities. This study examines four areas of stakeholder values and preferences: 1) beliefs about the causes, solutions, and severity of nearshore problems; 2) priorities for nearshore features, shoreforms, developments, and restoration objectives; 3) thoughts about ecosystem servicesiii and trade-offs among them; and 4) visions of a future, restored Puget Sound nearshore and the role of science in attaining this vision. The study is framed by two hypotheses from the Advocacy Coalition Framework (ACF), which suggests that groups of policy advocates form around shared “policy core beliefs” which can transcend traditional categories of stakeholders.(PDF contains 3 pages)

Neutron source from thin foil confined by two laser pulses

Neutron production from a thin deuterium-tritium (D-T) foil irradiated by two intense femtosecond laser pulses from opposite sides with zero phase difference is studied analytically and numerically. For the interaction of a laser pulse of amplitude a = 7, focal area 100 mu m(2) and areal density 4.4 x 10(18) cm(-2) with a D-T plasma foil, about 1.17 x 10(21) neutron s(-1) can be obtained, much more than from other methods. The profiles of the ion and electron densities are also calculated.