Comet 1P/Halley multifluid MHD model for the Giotto fly-by

The interaction of comets with the solar wind has been the focus of many studies including numerical modeling. We compare the results of our multifluid MHD simulation of comet 1P/Halley to data obtained during the flyby of the European Space Agency's Giotto spacecraft in 1986. The model solves the full set of MHD equations for the individual fluids representing the solar wind protons, the cometary light and heavy ions, and the electrons. The mass loading, charge-exchange, dissociative ion-electron recombination, and collisional interactions between the fluids are taken into account. The computational domain spans over several million kilometers, and the close vicinity of the comet is resolved to the details of the magnetic cavity. The model is validated by comparison to the corresponding Giotto observations obtained by the Ion Mass Spectrometer, the Neutral Mass Spectrometer, the Giotto magnetometer experiment, and the Johnstone Plasma Analyzer instrument. The model shows the formation of the bow shock, the ion pile-up, and the diamagnetic cavity and is able to reproduce the observed temperature differences between the pick-up ion populations and the solar wind protons. We give an overview of the global interaction of the comet with the solar wind and then show the effects of the Lorentz force interaction between the different plasma populations.

Solvation-Driven Charge Transfer and Localization in Metal Complexes

In any physicochemical process in liquids, the dynamical response of the solvent to the solutes out of equilibrium plays a crucial role in the rates and products: the solvent molecules react to the changes in volume and electron density of the solutes to minimize the free energy of the solution, thus modulating the activation barriers and stabilizing (or destabilizing) intermediate states. In charge transfer (CT) processes in polar solvents, the response of the solvent always assists the formation of charge separation states by stabilizing the energy of the localized charges. A deep understanding of the solvation mechanisms and time scales is therefore essential for a correct description of any photochemical process in dense phase and for designing molecular devices based on photosensitizers with CT excited states. In the last two decades, with the advent of ultrafast time-resolved spectroscopies, microscopic models describing the relevant case of polar solvation (where both the solvent and the solute molecules have a permanent electric dipole and the mutual interaction is mainly dipole−dipole) have dramatically progressed. Regardless of the details of each model, they all assume that the effect of the electrostatic fields of the solvent molecules on the internal electronic dynamics of the solute are perturbative and that the solvent−solute coupling is mainly an electrostatic interaction between the constant permanent dipoles of the solute and the solvent molecules. This well-established picture has proven to quantitatively rationalize spectroscopic effects of environmental and electric dynamics (time-resolved Stokes shifts, inhomogeneous broadening, etc.). However, recent computational and experimental studies, including ours, have shown that further improvement is required. Indeed, in the last years we investigated several molecular complexes exhibiting photoexcited CT states, and we found that the current description of the formation and stabilization of CT states in an important group of molecules such as transition metal complexes is inaccurate. In particular, we proved that the solvent molecules are not just spectators of intramolecular electron density redistribution but significantly modulate it. Our results solicit further development of quantum mechanics computational methods to treat the solute and (at least) the closest solvent molecules including the nonperturbative treatment of the effects of local electrostatics and direct solvent−solute interactions to describe the dynamical changes of the solute excited states during the solvent response.

Space requirement of a prefabricated bar on two interforaminal implants: a prospective clinical study

OBJECTIVES This clinical study measured the dimensional changes of existing lower complete dentures due to the integration of a prefabricated implant bar. Additionally, the impact of this dimensional change on patient satisfaction and oral function was analyzed. METHODS Twenty edentulous patients (10 men/10 women; aged 65.9 ± 11.8 years) received two interforaminal implants. Subsequent to surgery, a chair side adapted, prefabricated bar (SFI Bar(®), C+M, Biel, Switzerland) was inserted, and the matrix was polymerized into the existing lower denture. The change of the denture's lingual dimension was recorded by means of a bicolored, silicone denture duplicate that was sectioned in the oro-vestibular direction in the regions of the symphysis (S) and the implants (I-left, I-right). On the sections, the dimensional increase was measured using a light microscope. Six months after bar insertion, patients answered a standardized questionnaire. RESULTS All dentures exhibited increased lingual volume, more extensively at S than at I (P = 0.001). At S, the median diagonal size of the denture was doubled (+4.33 mm), and at I, the median increase was 50% (I-left/-right = +2.66/+2.62 mm). The original denture size influenced the volume increase (P = 0.024): smaller dentures led to a larger increase. The amount of denture increase did not have negative impact on either self-perceived oral function or patient satisfaction. Approximately, 95% of the patients were satisfied with the treatment results. CONCLUSIONS The lingual size of a lower denture was enlarged by the integration of a prefabricated bar without any negative side effects. Thus, this attachment system is suitable to convert an existing full denture into an implant-supported overdenture.

The Effects of Density, Spatial Pattern, and Competitive Symmetry on Size Variation in Simulated Plant Populations

Patterns of size inequality in crowded plant populations are often taken to be indicative of the degree of size asymmetry of competition, but recent research suggests that some of the patterns attributed to size‐asymmetric competition could be due to spatial structure. To investigate the theoretical relationships between plant density, spatial pattern, and competitive size asymmetry in determining size variation in crowded plant populations, we developed a spatially explicit, individual‐based plant competition model based on overlapping zones of influence. The zone of influence of each plant is modeled as a circle, growing in two dimensions, and is allometrically related to plant biomass. The area of the circle represents resources potentially available to the plant, and plants compete for resources in areas in which they overlap. The size asymmetry of competition is reflected in the rules for dividing up the overlapping areas. Theoretical plant populations were grown in random and in perfectly uniform spatial patterns at four densities under size‐asymmetric and size‐symmetric competition. Both spatial pattern and size asymmetry contributed to size variation, but their relative importance varied greatly over density and over time. Early in stand development, spatial pattern was more important than the symmetry of competition in determining the degree of size variation within the population, but after plants grew and competition intensified, the size asymmetry of competition became a much more important source of size variation. Size variability was slightly higher at higher densities when competition was symmetric and plants were distributed nonuniformly in space. In a uniform spatial pattern, size variation increased with density only when competition was size asymmetric. Our results suggest that when competition is size asymmetric and intense, it will be more important in generating size variation than is local variation in density. Our results and the available data are consistent with the hypothesis that high levels of size inequality commonly observed within crowded plant populations are largely due to size‐asymmetric competition, not to variation in local density.

Denoising Your Monte Carlo Renders: Recent Advances in Image Space-Adaptive Sampling and Reconstruction

With the ongoing shift in the computer graphics industry toward Monte Carlo rendering, there is a need for effective, practical noise-reduction techniques that are applicable to a wide range of rendering effects and easily integrated into existing production pipelines. This course surveys recent advances in image-space adaptive sampling and reconstruction algorithms for noise reduction, which have proven very effective at reducing the computational cost of Monte Carlo techniques in practice. These approaches leverage advanced image-filtering techniques with statistical methods for error estimation. They are attractive because they can be integrated easily into conventional Monte Carlo rendering frameworks, they are applicable to most rendering effects, and their computational overhead is modest.

Interstellar Flow and Temperature Determination with IBEX: Robustness and Sensitivity to Systematic Effects

The Interstellar Boundary Explorer (IBEX) samples the interstellar neutral (ISN) gas flow of several species every year from December through late March when the Earth moves into the incoming flow. The first quantitative analyses of these data resulted in a narrow tube in four-dimensional interstellar parameter space, which couples speed, flow latitude, flow longitude, and temperature, and center values with approximately 3° larger longitude and 3 km s⁻¹ lower speed, but with temperatures similar to those obtained from observations by the Ulysses spacecraft. IBEX has now recorded six years of ISN flow observations, providing a large database over increasing solar activity and using varying viewing strategies. In this paper, we evaluate systematic effects that are important for the ISN flow vector and temperature determination. We find that all models in use return ISN parameters well within the observational uncertainties and that the derived ISN flow direction is resilient against uncertainties in the ionization rate. We establish observationally an effective IBEX-Lo pointing uncertainty of ±0°18 in spin angle and confirm an uncertainty of ±0°1 in longitude. We also show that the IBEX viewing strategy with different spin-axis orientations minimizes the impact of several systematic uncertainties, and thus improves the robustness of the measurement. The Helium Warm Breeze has likely contributed substantially to the somewhat different center values of the ISN flow vector. By separating the flow vector and temperature determination, we can mitigate these effects on the analysis, which returns an ISN flow vector very close to the Ulysses results, but with a substantially higher temperature. Due to coupling with the ISN flow speed along the ISN parameter tube, we provide the temperature Tvisn∞=8710+440/-680 K for Visn∞=26 km s⁻¹ for comparison, where most of the uncertainty is systematic and likely due to the presence of the Warm Breeze.

Effects of exercise with oxygen prebreathe for decompression risk reduction

Astronauts performing extravehicular activities (EVA) are at risk for occupational hazards due to a hypobaric environment, in particular Decompression Sickness (DCS). DCS results from nitrogen gas bubble formation in body tissues and venous blood. Denitrogenation achieved through lengthy staged decompression protocols has been the mainstay of prevention of DCS in space. Due to the greater number and duration of EVAs scheduled for construction and maintenance of the International Space Station, more efficient alternatives to accomplish missions without compromising astronaut safety are desirable. ^ This multi-center, multi-phase study (NASA-Prebreathe Reduction Protocol study, or PRP) was designed to identify a shorter denitrogenation protocol that can be implemented before an EVA, based on the combination of adynamia and exercise enhanced oxygen prebreathe. Human volunteers recruited at three sites (Texas, North Carolina and Canada) underwent three different combinations (“PRP phases”) of intense and light exercise prior to decompression in an altitude chamber. The outcome variables were detection of venous gas embolism (VGE) by precordial Doppler ultrasound, and clinical manifestations of DCS. Independent variables included age, gender, body mass index, oxygen consumption peak, peak heart rate, and PRP phase. Data analysis was performed both by pooling results from all study sites, and by examining each site separately. ^ Ten percent of the subjects developed DCS and 20% showed evidence of high grade VGE. No cases of DCS occurred in one particular PRP phase with use of the combination of dual-cycle ergometry (10 minutes at 75% of VO2 peak) plus 24 minutes of light EVA exercise (p = 0.04). No significant effects were found for the remaining independent variables on the occurrence of DCS. High grade VGE showed a strong correlation with subsequent development of DCS (sensitivity, 88.2%; specificity, 87.2%). In the presence of high grade VGE, the relative risk for DCS ranged from 7.52 to 35.0. ^ In summary, a good safety level can be achieved with exercise-enhanced oxygen denitrogenation that can be generalized to the astronaut population. Exercise is beneficial in preventing DCS if a specific schedule is followed, with an individualized VO2 prescription that provides a safety level that can then be applied to space operations. Furthermore, VGE Doppler detection is a useful clinical tool for prediction of altitude DCS. Because of the small number of high grade VGE episodes, the identification of a high probability DCS situation based on the presence of high grade VGE seems justified in astronauts. ^

The Effects Of Reinforcement Inclusions On Wear Tolerance, Playing Quality And Physical Properties In A Silt Loam And Sand Rootzone Matrix

Reinforcement inclusions have been advocated to alleviate wear, compaction, and unstable surfaces in sports fields, but little research on the effects of these materials has been conducted in the USA. Experiments were established on a native silt loam and a sand rootzone matrix, seeded with a Kentucky bluegrass (Poa pratensis L.) blend, at the Joseph Troll Turf Research Center, University of Massachusetts, Amherst, USA to determine the effects of reinforcement inclusions on wear, surface hardness, traction, ball roll, ball bounce resilience, water infiltration rate, soil bulk density, air porosity, total porosity, and root weights. Three types of reinforcement inclusions (Sportgrass, Netlon, Turfgrids) were tested along with a non-reinforced control in a three year study. The treatments were set out in a randomized complete block design with four replications in both soils. No inclusion provided less wear or greater infiltration or air-filled porosity relative to the control. Reinforcement inclusions showed significant differences, however, in surface hardness, traction, and ball roll relative to the control, although this varied with the time of year. Infiltration rates, airfilled porosity, total pore space, bulk density, hardness, traction, ball roll, and ball rebound were greater on the sand rootzone than on the silt loam. Significant correlations were present between soil bulk density, surface hardness, traction, and ball roll. Based on our study, the use of reinforcement inclusions to provide better wear tolerance for sand or native soil athletic fields is not warranted. Certain playing surface characteristics, however, may be slightly improved with the use of reinforcement inclusions. The use of sands for sports surfaces is justified based upon the improvement in playing quality characteristics and soil physical properties important to a good playing surface.

Atmospheric, oceanic and hydrological polar motion excitation mechanisms based on a combination of geometric and gravimetric space observations

The goal of our study is to determine accurate time series of geophysical Earth rotation excitations to learn more about global dynamic processes in the Earth system. For this purpose, we developed an adjustment model which allows to combine precise observations from space geodetic observation systems, such as Satellite Laser Ranging (SLR), Global Navigation Satellite Systems (GNSS), Very Long Baseline Interferometry (VLBI), Doppler Orbit determination and Radiopositioning Integrated on Satellite (DORIS), satellite altimetry and satellite gravimetry in order to separate geophysical excitation mechanisms of Earth rotation. Three polar motion time series are applied to derive the polar motion excitation functions (integral effect). Furthermore we use five time variable gravity field solutions from Gravity Recovery and Climate Experiment (GRACE) to determine not only the integral mass effect but also the oceanic and hydrological mass effects by applying suitable filter techniques and a land-ocean mask. For comparison the integral mass effect is also derived from degree 2 potential coefficients that are estimated from SLR observations. The oceanic mass effect is also determined from sea level anomalies observed by satellite altimetry by reducing the steric sea level anomalies derived from temperature and salinity fields of the oceans. Due to the combination of all geodetic estimated excitations the weaknesses of the individual processing strategies can be reduced and the technique-specific strengths can be accounted for. The formal errors of the adjusted geodetic solutions are smaller than the RMS differences of the geophysical model solutions. The improved excitation time series can be used to improve the geophysical modeling.

Interactive effects of salinity and elevated CO2 levels on juvenile eastern oysters, Crassostrea virginica

Rising levels of atmospheric CO2 lead to acidification of the ocean and alter seawater carbonate chemistry, which can negatively impact calcifying organisms, including mollusks. In estuaries, exposure to elevated CO2 levels often co-occurs with other stressors, such as reduced salinity, which enhances the acidification trend, affects ion and acid-base regulation of estuarine calcifiers and modifies their response to ocean acidification. We studied the interactive effects of salinity and partial pressure of CO2 (PCO2) on biomineralization and energy homeostasis in juveniles of the eastern oyster, Crassostrea virginica, a common estuarine bivalve. Juveniles were exposed for 11 weeks to one of two environmentally relevant salinities (30 or 15 PSU) either at current atmospheric PCO2 (400 µatm, normocapnia) or PCO2 projected by moderate IPCC scenarios for the year 2100 (700-800 µatm, hypercapnia). Exposure of the juvenile oysters to elevated PCO2 and/or low salinity led to a significant increase in mortality, reduction of tissue energy stores (glycogen and lipid) and negative soft tissue growth, indicating energy deficiency. Interestingly, tissue ATP levels were not affected by exposure to changing salinity and PCO2, suggesting that juvenile oysters maintain their cellular energy status at the expense of lipid and glycogen stores. At the same time, no compensatory upregulation of carbonic anhydrase activity was found under the conditions of low salinity and high PCO2. Metabolic profiling using magnetic resonance spectroscopy revealed altered metabolite status following low salinity exposure; specifically, acetate levels were lower in hypercapnic than in normocapnic individuals at low salinity. Combined exposure to hypercapnia and low salinity negatively affected mechanical properties of shells of the juveniles, resulting in reduced hardness and fracture resistance. Thus, our data suggest that the combined effects of elevated PCO2 and fluctuating salinity may jeopardize the survival of eastern oysters because of weakening of their shells and increased energy consumption.

Modelización de costas rocosas acantiladas

El retroceso de las costas acantiladas es un fenómeno muy extendido sobre los litorales rocosos expuestos a la incidencia combinada de los procesos marinos y meteorológicos que se dan en la franja costera. Este fenómeno se revela violentamente como movimientos gravitacionales del terreno esporádicos, pudiendo causar pérdidas materiales y/o humanas. Aunque el conocimiento de estos riesgos de erosión resulta de vital importancia para la correcta gestión de la costa, el desarrollo de modelos predictivos se encuentra limitado desde el punto de vista geomorfológico debido a la complejidad e interacción de los procesos de desarrollo espacio-temporal que tienen lugar en la zona costera. Los modelos de predicción publicados son escasos y con importantes inconvenientes: a) extrapolación, extienden la información de registros históricos; b) empíricos, sobre registros históricos estudian la respuesta al cambio de un parámetro; c) estocásticos, determinan la cadencia y magnitud de los eventos futuros extrapolando las distribuciones de probabilidad extraídas de catálogos históricos; d) proceso-respuesta, de estabilidad y propagación del error inexplorada; e) en Ecuaciones en Derivadas Parciales, computacionalmente costosos y poco exactos. La primera parte de esta tesis detalla las principales características de los modelos más recientes de cada tipo y, para los más habitualmente utilizados, se indican sus rangos de aplicación, ventajas e inconvenientes. Finalmente como síntesis de los procesos más relevantes que contemplan los modelos revisados, se presenta un diagrama conceptual de la recesión costera, donde se recogen los procesos más influyentes que deben ser tenidos en cuenta, a la hora de utilizar o crear un modelo de recesión costera con el objetivo de evaluar la peligrosidad (tiempo/frecuencia) del fenómeno a medio-corto plazo. En esta tesis se desarrolla un modelo de proceso-respuesta de retroceso de acantilados costeros que incorpora el comportamiento geomecánico de materiales cuya resistencia a compresión no supere los 5 MPa. El modelo simula la evolución espaciotemporal de un perfil-2D del acantilado que puede estar formado por materiales heterogéneos. Para ello, se acoplan la dinámica marina: nivel medio del mar, cambios en el nivel medio del lago, mareas y oleaje; con la evolución del terreno: erosión, desprendimiento rocoso y formación de talud de derrubios. El modelo en sus diferentes variantes es capaz de incluir el análisis de la estabilidad geomecánica de los materiales, el efecto de los derrubios presentes al pie del acantilado, el efecto del agua subterránea, la playa, el run-up, cambios en el nivel medio del mar o cambios (estacionales o interanuales) en el nivel medio de la masa de agua (lagos). Se ha estudiado el error de discretización del modelo y su propagación en el tiempo a partir de las soluciones exactas para los dos primeros periodos de marea para diferentes aproximaciones numéricas tanto en tiempo como en espacio. Los resultados obtenidos han permitido justificar las elecciones que minimizan el error y los métodos de aproximación más adecuados para su posterior uso en la modelización. El modelo ha sido validado frente a datos reales en la costa de Holderness, Yorkshire, Reino Unido; y en la costa norte del lago Erie, Ontario, Canadá. Los resultados obtenidos presentan un importante avance en los modelos de recesión costera, especialmente en su relación con las condiciones geomecánicas del medio, la influencia del agua subterránea, la verticalización de los perfiles rocosos y su respuesta ante condiciones variables producidas por el cambio climático (por ejemplo, nivel medio del mar, cambios en los niveles de lago, etc.). The recession of coastal cliffs is a widespread phenomenon on the rocky shores that are exposed to the combined incidence of marine and meteorological processes that occur in the shoreline. This phenomenon is revealed violently and occasionally, as gravitational movements of the ground and can cause material or human losses. Although knowledge of the risks of erosion is vital for the proper management of the coast, the development of cliff erosion predictive models is limited by the complex interactions between environmental processes and material properties over a range of temporal and spatial scales. Published prediction models are scarce and present important drawbacks: extrapolation, that extend historical records to the future; empirical, that based on historical records studies the system response against the change in one parameter; stochastic, that represent of cliff behaviour based on assumptions regarding the magnitude and frequency of events in a probabilistic framework based on historical records; process-response, stability and error propagation unexplored; PDE´s, highly computationally expensive and not very accurate. The first part of this thesis describes the main features of the latest models of each type and, for the most commonly used, their ranges of application, advantages and disadvantages are given. Finally as a synthesis of the most relevant processes that include the revised models, a conceptual diagram of coastal recession is presented. This conceptual model includes the most influential processes that must be taken into account when using or creating a model of coastal recession to evaluate the dangerousness (time/frequency) of the phenomenon to medium-short term. A new process-response coastal recession model developed in this thesis has been designed to incorporate the behavioural and mechanical characteristics of coastal cliffs which are composed of with materials whose compressive strength is less than 5 MPa. The model simulates the spatial and temporal evolution of a cliff-2D profile that can consist of heterogeneous materials. To do so, marine dynamics: mean sea level, waves, tides, lake seasonal changes; is coupled with the evolution of land recession: erosion, cliff face failure and associated protective colluvial wedge. The model in its different variants can include analysis of material geomechanical stability, the effect of debris present at the cliff foot, groundwater effects, beach and run-up effects, changes in the mean sea level or changes (seasonal or inter-annual) in the mean lake level. Computational implementation and study of different numerical resolution techniques, in both time and space approximations, and the produced errors are exposed and analysed for the first two tidal periods. The results obtained in the errors analysis allow us to operate the model with a configuration that minimizes the error of the approximation methods. The model is validated through profile evolution assessment at various locations of coastline retreat on the Holderness Coast, Yorkshire, UK and on the north coast of Lake Erie, Ontario, Canada. The results represent an important stepforward in linking material properties to the processes of cliff recession, in considering the effect of groundwater charge and the slope oversteeping and their response to changing conditions caused by climate change (i.e. sea level, changes in lakes levels, etc.).

The direct boundary element method: 2D site effects assessment on laterally varying layered media (methodology)

The Direct Boundary Element Method (DBEM) is presented to solve the elastodynamic field equations in 2D, and a complete comprehensive implementation is given. The DBEM is a useful approach to obtain reliable numerical estimates of site effects on seismic ground motion due to irregular geological configurations, both of layering and topography. The method is based on the discretization of the classical Somigliana's elastodynamic representation equation which stems from the reciprocity theorem. This equation is given in terms of the Green's function which is the full-space harmonic steady-state fundamental solution. The formulation permits the treatment of viscoelastic media, therefore site models with intrinsic attenuation can be examined. By means of this approach, the calculation of 2D scattering of seismic waves, due to the incidence of P and SV waves on irregular topographical profiles is performed. Sites such as, canyons, mountains and valleys in irregular multilayered media are computed to test the technique. The obtained transfer functions show excellent agreement with already published results.

Dynamic stabilization of L2 periodic orbits using attitude-orbit coupling effects

Numerical explorations show how the known periodic solutions of the Hill problem are modified in the case of the attitude-orbit coupling that may occur for large satellite structures. We focus on the case in which the elongation is the dominant satellite’s characteristic and find that a rotating structure may remain with its largest dimension in a plane parallel to the plane of the primaries. In this case, the effect produced by the non-negligible physical length is dynamically equivalent to the perturbation produced by an oblate central body on a mass-point satellite. Based on this, it is demonstrated that the attitude-orbital coupling of a long enough body may change the dynamical characteristics of a periodic orbit about the collinear Lagrangian points.

Self-similar motion of laser half-space plasmas. I. Deflagration regime

The one-dimensional self-similar motion of an initially cold, half-space plasma of electron density 0,produced by the (anomalous) absorption of a laser pulse of irradiation = (j>0f/T(0< (< T) at the critical density nc(«c/«0=echarge number), and a = (9/c/4m,) (T/C 2n l/4>oKe)213, where k, m, are Boltzmann's constant and the ion mass, and Ke X (electron temperature)5'2 = heat conductivity. If a >e- 4 ' 3 , a deflagration wave separates an isentropic compression with a shock bounding the undisturbed plasma, and an isentropic expansion flow to the vacuum. The structures of these three regions are completely determined; in particular, the Chapman-Jouguet condition is proved and the density behind the deflagration is found. The deflagration-compression thickness ratio is large (small) for a^e- 5 ' 3(a>e- 5 ' 3 ) . The compression to expansion ratio for both energy and thickness is 0(e"2). For Z,- large, a deflagration exists even if a~e~413. Condition a>e~4'3 may be applied to pulses that are not linear.

Self-similar motion of laser half-space plasmas. II. Thermal waves and intermediate regimes

The one-dimensional self-similar motion of an initially cold, half-space plasma of electron density n,produced by the (anomalous) absorption of a laser pulse of irradiation

charge number), and a parameter a€~4'3, a qualitative discussion of how plasma behavior changes with a, is given.