984 resultados para Numerical-solution
Exact solution to the exit-time problem for an undamped free particle driven by Gaussian white noise
Resumo:
In a recent paper [Phys. Rev. Lett. 75, 189 (1995)] we have presented the exact analytical expression for the mean exit time, T(x,v), of a free inertial process driven by Gaussian white noise out of a region (0,L) in space. In this paper we give a detailed account of the method employed and present results on asymptotic properties and averages of T(x,v).
Resumo:
We consider damage spreading transitions in the framework of mode-coupling theory. This theory describes relaxation processes in glasses in the mean-field approximation which are known to be characterized by the presence of an exponentially large number of metastable states. For systems evolving under identical but arbitrarily correlated noises, we demonstrate that there exists a critical temperature T0 which separates two different dynamical regimes depending on whether damage spreads or not in the asymptotic long-time limit. This transition exists for generic noise correlations such that the zero damage solution is stable at high temperatures, being minimal for maximal noise correlations. Although this dynamical transition depends on the type of noise correlations, we show that the asymptotic damage has the good properties of a dynamical order parameter, such as (i) independence of the initial damage; (ii) independence of the class of initial condition; and (iii) stability of the transition in the presence of asymmetric interactions which violate detailed balance. For maximally correlated noises we suggest that damage spreading occurs due to the presence of a divergent number of saddle points (as well as metastable states) in the thermodynamic limit consequence of the ruggedness of the free-energy landscape which characterizes the glassy state. These results are then compared to extensive numerical simulations of a mean-field glass model (the Bernasconi model) with Monte Carlo heat-bath dynamics. The freedom of choosing arbitrary noise correlations for Langevin dynamics makes damage spreading an interesting tool to probe the ruggedness of the configurational landscape.
Resumo:
In this paper we give some ideas that can be useful to solve Schrödinger equations in the case when the Hamiltonian contains a large term. We obtain an expansion of the solution in reciprocal powers of the large coupling constant. The procedure followed consists in considering that the small part of the Hamiltonian engenders a motion adiabatic to the motion generated by the large part of the same.
Resumo:
ABSTRACT Soil solution samplers may have the same working principle, but they differ in relation to chemical and physical characteristics, cost and handling, and these aspects exert influence on the chemical composition of the soil solution obtained. This study was carried out to evaluate, over time, the chemical composition of solutions extracted by Suolo Acqua, with the hydrophilic membrane (HM) as a standard, using soils with contrasting characteristics, and to determine the relationship between electrical conductivity (EC) and concentration of ions and pH of soil solution samples. This study was carried out under laboratory conditions, using three soils samples with different clay and organic matter (OM) contents. Soil solution contents of F−, Cl−, NO−3, Br−, SO42−, Na+, NH4+, K+, Mg2+, Ca2+, were analyzed, as well as inorganic, organic, and total C contents, pH, and EC, in four successive sampling times. Soil solution chemical composition extracted by the Suolo Acqua sampler is similar to that collected by the HM, but the Suolo Acqua extracted more Na+ and soluble organic C than the HM solution. Solution EC, cation and anion concentrations, and soluble C levels are higher in the soil with greater clay and OM contents (Latossolo and Cambissolo in this case). Soil solution composition varied over time, with considerable changes in pH, EC, and nutrient concentrations, especially associated with soil OM. Thus, single and isolated sampling of the soil solution must be avoided, otherwise composition of the soil solution may not be correctly evaluated. Soil solution EC was regulated by pH, as well as the sum of cation and anion concentrations, and the C contents determined in the soil liquid phase.
Resumo:
We characterize the Schatten class membership of the canonical solution operator to $\overline{\partial}$ acting on $L^2(e^{-2\phi})$, where $\phi$ is a subharmonic function with $\Delta\phi$ a doubling measure. The obtained characterization is in terms of $\Delta\phi$. As part of our approach, we study Hankel operators with anti-analytic symbols acting on the corresponding Fock space of entire functions in $L^2(e^{-2\phi})$
Resumo:
La regressió basada en distàncies és un mètode de predicció que consisteix en dos passos: a partir de les distàncies entre observacions obtenim les variables latents, les quals passen a ser els regressors en un model lineal de mínims quadrats ordinaris. Les distàncies les calculem a partir dels predictors originals fent us d'una funció de dissimilaritats adequada. Donat que, en general, els regressors estan relacionats de manera no lineal amb la resposta, la seva selecció amb el test F usual no és possible. En aquest treball proposem una solució a aquest problema de selecció de predictors definint tests estadístics generalitzats i adaptant un mètode de bootstrap no paramètric per a l'estimació dels p-valors. Incluim un exemple numèric amb dades de l'assegurança d'automòbils.
Resumo:
Interfacial hydrodynamic instabilities arise in a range of chemical systems. One mechanism for instability is the occurrence of unstable density gradients due to the accumulation of reaction products. In this paper we conduct two-dimensional nonlinear numerical simulations for a member of this class of system: the methylene-blue¿glucose reaction. The result of these reactions is the oxidation of glucose to a relatively, but marginally, dense product, gluconic acid, that accumulates at oxygen permeable interfaces, such as the surface open to the atmosphere. The reaction is catalyzed by methylene-blue. We show that simulations help to disassemble the mechanisms responsible for the onset of instability and evolution of patterns, and we demonstrate that some of the results are remarkably consistent with experiments. We probe the impact of the upper oxygen boundary condition, for fixed flux, fixed concentration, or mixed boundary conditions, and find significant qualitative differences in solution behavior; structures either attract or repel one another depending on the boundary condition imposed. We suggest that measurement of the form of the boundary condition is possible via observation of oxygen penetration, and improved product yields may be obtained via proper control of boundary conditions in an engineering setting. We also investigate the dependence on parameters such as the Rayleigh number and depth. Finally, we find that pseudo-steady linear and weakly nonlinear techniques described elsewhere are useful tools for predicting the behavior of instabilities beyond their formal range of validity, as good agreement is obtained with the simulations.
Resumo:
The Multiscale Finite Volume (MsFV) method has been developed to efficiently solve reservoir-scale problems while conserving fine-scale details. The method employs two grid levels: a fine grid and a coarse grid. The latter is used to calculate a coarse solution to the original problem, which is interpolated to the fine mesh. The coarse system is constructed from the fine-scale problem using restriction and prolongation operators that are obtained by introducing appropriate localization assumptions. Through a successive reconstruction step, the MsFV method is able to provide an approximate, but fully conservative fine-scale velocity field. For very large problems (e.g. one billion cell model), a two-level algorithm can remain computational expensive. Depending on the upscaling factor, the computational expense comes either from the costs associated with the solution of the coarse problem or from the construction of the local interpolators (basis functions). To ensure numerical efficiency in the former case, the MsFV concept can be reapplied to the coarse problem, leading to a new, coarser level of discretization. One challenge in the use of a multilevel MsFV technique is to find an efficient reconstruction step to obtain a conservative fine-scale velocity field. In this work, we introduce a three-level Multiscale Finite Volume method (MlMsFV) and give a detailed description of the reconstruction step. Complexity analyses of the original MsFV method and the new MlMsFV method are discussed, and their performances in terms of accuracy and efficiency are compared.
Resumo:
INTRODUCTION: The importance of the micromovements in the mechanism of aseptic loosening is clinically difficult to evaluate. To complete the analysis of a series of total knee arthroplasties (TKA), we used a tridimensional numerical model to study the micromovements of the tibial implant. MATERIAL AND METHODS: Fifty one patients (with 57 cemented Porous Coated Anatomic TKAs) were reviewed (mean follow-up 4.5 year). Radiolucency at the tibial bone-cement interface was sought on the AP radiographs and divided in 7 areas. The distribution of the radiolucency was then correlated with the axis of the lower limb as measured on the orthoradiograms. The tridimensional numerical model is based on the finite element method. It allowed the measurement of the cemented prosthetic tibial implant's displacements and the micromovements generated at bone-ciment interface. A total load (2000 Newton) was applied at first vertically and asymetrically on the tibial plateau, thereby simulating an axial deviation of the lower limbs. The vector's posterior inclination then permitted the addition of a tangential component to the axial load. This type of effort is generated by complex biomechanical phenomena such as knee flexion. RESULTS: 81 per cent of the 57 knees had a radiolucent line of at least 1 mm, at one or more of the tibial cement-epiphysis jonctional areas. The distribution of these lucent lines showed that they came out more frequently at the periphery of the implant. The lucent lines appeared most often under the unloaded margin of the tibial plateau, when axial deviation of lower limbs was present. Numerical simulations showed that asymetrical loading on the tibial plateau induced a subsidence of the loaded margin (0-100 microns) and lifting off at the opposite border (0-70 microns). The postero-anterior tangential component induced an anterior displacement of the tibial implant (160-220 microns), and horizontal micromovements with non homogenous distribution at the bone-ciment interface (28-54 microns). DISCUSSION: Comparison of clinical and numerical results showed a relation between the development of radiolucent lines and the unloading of the tibial implant's margin. The deleterious effect of lower limbs' axial deviation is thereby proven. The irregular distribution of lucent lines under the tibial plateau was similar of the micromovements' repartition at the bone-cement interface when tangential forces were present. A causative relation between the two phenomenaes could not however be established. Numerical simulation is a truly useful method of study; it permits to calculate micromovements which are relative, non homogenous and of very low amplitude. However, comparative clinical studies remain as essential to ensure the credibility of results.