16 resultados para Numerical solutions of ODE’s
em Biblioteca Digital da Produção Intelectual da Universidade de São Paulo
Resumo:
Scaling methods allow a single solution to Richards' equation (RE) to suffice for numerous specific cases of water flow in unsaturated soils. During the past half-century, many such methods were developed for similar soils. In this paper, a new method is proposed for scaling RE for a wide range of dissimilar soils. Exponential-power (EP) functions are used to reduce the dependence of the scaled RE on the soil hydraulic properties. To evaluate the proposed method, the scaled RE was solved numerically considering two test cases: infiltration into relatively dry soils having initially uniform water content distributions, and gravity-dominant drainage occurring from initially wet soil profiles. Although the results for four texturally different soils ranging from sand to heavy clay (adopted from the UNSODA database) showed that the scaled solution were invariant for a wide range of flow conditions, slight deviations were observed when the soil profile was initially wet in the infiltration case or deeply wet in the drainage case. The invariance of the scaled RE makes it possible to generalize a single solution of RE to many dissimilar soils and conditions. Such a procedure reduces the numerical calculations and provides additional opportunities for solving the highly nonlinear RE for unsaturated water flow in soils.
Resumo:
In this work, different methods to estimate the value of thin film residual stresses using instrumented indentation data were analyzed. This study considered procedures proposed in the literature, as well as a modification on one of these methods and a new approach based on the effect of residual stress on the value of hardness calculated via the Oliver and Pharr method. The analysis of these methods was centered on an axisymmetric two-dimensional finite element model, which was developed to simulate instrumented indentation testing of thin ceramic films deposited onto hard steel substrates. Simulations were conducted varying the level of film residual stress, film strain hardening exponent, film yield strength, and film Poisson's ratio. Different ratios of maximum penetration depth h(max) over film thickness t were also considered, including h/t = 0.04, for which the contribution of the substrate in the mechanical response of the system is not significant. Residual stresses were then calculated following the procedures mentioned above and compared with the values used as input in the numerical simulations. In general, results indicate the difference that each method provides with respect to the input values depends on the conditions studied. The method by Suresh and Giannakopoulos consistently overestimated the values when stresses were compressive. The method provided by Wang et al. has shown less dependence on h/t than the others.
Resumo:
In the past few decades detailed observations of radio and X-ray emission from massive binary systems revealed a whole new physics present in such systems. Both thermal and non-thermal components of this emission indicate that most of the radiation at these bands originates in shocks. O and B-type stars and WolfRayet (WR) stars present supersonic and massive winds that, when colliding, emit largely due to the freefree radiation. The non-thermal radio and X-ray emissions are due to synchrotron and inverse Compton processes, respectively. In this case, magnetic fields are expected to play an important role in the emission distribution. In the past few years the modelling of the freefree and synchrotron emissions from massive binary systems have been based on purely hydrodynamical simulations, and ad hoc assumptions regarding the distribution of magnetic energy and the field geometry. In this work we provide the first full magnetohydrodynamic numerical simulations of windwind collision in massive binary systems. We study the freefree emission characterizing its dependence on the stellar and orbital parameters. We also study self-consistently the evolution of the magnetic field at the shock region, obtaining also the synchrotron energy distribution integrated along different lines of sight. We show that the magnetic field in the shocks is larger than that obtained when the proportionality between B and the plasma density is assumed. Also, we show that the role of the synchrotron emission relative to the total radio emission has been underestimated.
Resumo:
We study local rigidity and multiplicity of constant scalar curvature metrics in arbitrary products of compact manifolds. Using (equivariant) bifurcation theory we determine the existence of infinitely many metrics that are accumulation points of pairwise non-homothetic solutions of the Yamabe problem. Using local rigidity and some compactness results for solutions of the Yamabe problem, we also exhibit new examples of conformal classes (with positive Yamabe constant) for which uniqueness holds. (C) 2011 Elsevier Masson SAS. All rights reserved.
Resumo:
In this paper, we give sufficient conditions for the uniform boundedness and uniform ultimate boundedness of solutions of a class of retarded functional differential equations with impulse effects acting on variable times. We employ the theory of generalized ordinary differential equations to obtain our results. As an example, we investigate the boundedness of the solution of a circulating fuel nuclear reactor model.
Resumo:
We study an elliptic system of the form Lu = vertical bar v vertical bar(p-1) v and Lv = vertical bar u vertical bar(q-1) u in Omega with homogeneous Dirichlet boundary condition, where Lu := -Delta u in the case of a bounded domain and Lu := -Delta u + u in the cases of an exterior domain or the whole space R-N. We analyze the existence, uniqueness, sign and radial symmetry of ground state solutions and also look for sign changing solutions of the system. More general non-linearities are also considered.
Resumo:
This paper shows in detail the modelling of anisotropic polymeric foam under compression and tension loadings, including discussions on isotropic material models and the entire procedure to calibrate the parameters involved. First, specimens of poly(vinyl chloride) (PVC) foam were investigated through experimental analyses in order to understand the mechanical behavior of this anisotropic material. Then, isotropic material models available in the commercial software Abaqus (TM) were investigated in order to verify their ability to model anisotropic foams and how the parameters involved can influence the results. Due to anisotropy, it is possible to obtain different values for the same parameter in the calibration process. The obtained set of parameters are used to calibrate the model according to the application of the structure. The models investigated showed minor and major limitations to simulate the mechanical behavior of anisotropic PVC foams under compression, tension and multi-axial loadings. Results show that the calibration process and the choice of the material model applied to the polymeric foam can provide good quantitative results and save project time. Results also indicate what kind and order of error one will get if certain choices are made throughout the modelling process. Finally, even though the developed calibration procedure is applied to specific PVC foam, it still outlines a very broad drill to analyze other anisotropic cellular materials.
Resumo:
Mach number and thermal effects on the mechanisms of sound generation and propagation are investigated in spatially evolving two-dimensional isothermal and non-isothermal mixing layers at Mach number ranging from 0.2 to 0.4 and Reynolds number of 400. A characteristic-based formulation is used to solve by direct numerical simulation the compressible Navier-Stokes equations using high-order schemes. The radiated sound is directly computed in a domain that includes both the near-field aerodynamic source region and the far-field sound propagation. In the isothermal mixing layer, Mach number effects may be identified in the acoustic field through an increase of the directivity associated with the non-compactness of the acoustic sources. Baroclinic instability effects may be recognized in the non-isothermal mixing layer, as the presence of counter-rotating vorticity layers, the resulting acoustic sources being found less efficient. An analysis based on the acoustic analogy shows that the directivity increase with the Mach number can be associated with the emergence of density fluctuations of weak amplitude but very efficient in terms of noise generation at shallow angle. This influence, combined with convection and refraction effects, is found to shape the acoustic wavefront pattern depending on the Mach number.
Resumo:
This work presents numerical simulations of two fluid flow problems involving moving free surfaces: the impacting drop and fluid jet buckling. The viscoelastic model used in these simulations is the eXtended Pom-Pom (XPP) model. To validate the code, numerical predictions of the drop impact problem for Newtonian and Oldroyd-B fluids are presented and compared with other methods. In particular, a benchmark on numerical simulations for a XPP drop impacting on a rigid plate is performed for a wide range of the relevant parameters. Finally, to provide an additional application of free surface flows of XPP fluids, the viscous jet buckling problem is simulated and discussed. (C) 2011 Elsevier B.V. All rights reserved.
Resumo:
In this work we extend the first order formalism for cosmological models that present an interaction between a fermionic and a scalar field. Cosmological exact solutions describing universes filled with interacting dark energy and dark matter have been obtained. Viable cosmological solutions with an early period of decelerated expansion followed by late acceleration have been found, notably one which presents a dark matter component dominating in the past and a dark energy component dominating in the future. In another one, the dark energy alone is the responsible for both periods, similar to a Chaplygin gas case. Exclusively accelerating solutions have also been obtained.
Resumo:
We characterize the existence of periodic solutions of some abstract neutral functional differential equations with finite and infinite delay when the underlying space is a UMD space. (C) 2011 Elsevier Inc. All rights reserved.
Resumo:
STUDY BY MASS SPECTROMETRY OF SOLUTIONS OF [HYDROXY(TOSYLOXY)IODO]BENZENE: PROPOSED DISPROPORTIONATION MECHANISMS. Solutions of [hydroxy(tosyloxy)iodo]benzene (HTIB or Koser's reagent) in acetonitrile were analyzed using high resolution electrospray ionization mass spectrometry (ESI-MS) and electrospray ionization tandem mass spectrometry (ESI-MS/MS) under different conditions. Several species were characterized in these analyses. Based on these data, mechanisms were proposed for the disproportionation of the iodine(III) compounds in iodine(V) and iodine(I) species.
Resumo:
The stability of two recently developed pressure spaces has been assessed numerically: The space proposed by Ausas et al. [R.F. Ausas, F.S. Sousa, G.C. Buscaglia, An improved finite element space for discontinuous pressures, Comput. Methods Appl. Mech. Engrg. 199 (2010) 1019-1031], which is capable of representing discontinuous pressures, and the space proposed by Coppola-Owen and Codina [A.H. Coppola-Owen, R. Codina, Improving Eulerian two-phase flow finite element approximation with discontinuous gradient pressure shape functions, Int. J. Numer. Methods Fluids, 49 (2005) 1287-1304], which can represent discontinuities in pressure gradients. We assess the stability of these spaces by numerically computing the inf-sup constants of several meshes. The inf-sup constant results as the solution of a generalized eigenvalue problems. Both spaces are in this way confirmed to be stable in their original form. An application of the same numerical assessment tool to the stabilized equal-order P-1/P-1 formulation is then reported. An interesting finding is that the stabilization coefficient can be safely set to zero in an arbitrary band of elements without compromising the formulation's stability. An analogous result is also reported for the mini-element P-1(+)/P-1 when the velocity bubbles are removed in an arbitrary band of elements. (C) 2012 Elsevier B.V. All rights reserved.
Resumo:
The reduction of friction and wear in systems presenting metal-to-metal contacts, as in several mechanical components, represents a traditional challenge in tribology. In this context, this work presents a computational study based on the linear Archard's wear law and finite element modeling (FEM), in order to analyze unlubricated sliding wear observed in typical pin on disc tests. Such modeling was developed using finite element software Abaqus® with 3-D deformable geometries and elastic–plastic material behavior for the contact surfaces. Archard's wear model was implemented into a FORTRAN user subroutine (UMESHMOTION) in order to describe sliding wear. Modeling of debris and oxide formation mechanisms was taken into account by the use of a global wear coefficient obtained from experimental measurements. Such implementation considers an incremental computation for surface wear based on the nodal displacements by means of adaptive mesh tools that rearrange local nodal positions. In this way, the worn track was obtained and new surface profile is integrated for mass loss assessments. This work also presents experimental pin on disc tests with AISI 4140 pins on rotating AISI H13 discs with normal loads of 10, 35, 70 and 140 N, which represent, respectively, mild, transition and severe wear regimes, at sliding speed of 0.1 m/s. Numerical and experimental results were compared in terms of wear rate and friction coefficient. Furthermore, in the numerical simulation the stress field distribution and changes in the surface profile across the worn track of the disc were analyzed. The applied numerical formulation has shown to be more appropriate to predict mild wear regime than severe regime, especially due to the shorter running-in period observed in lower loads that characterizes this kind of regime.
Resumo:
Computational fluid dynamics, CFD, is becoming an essential tool in the prediction of the hydrodynamic efforts and flow characteristics of underwater vehicles for manoeuvring studies. However, when applied to the manoeuvrability of autonomous underwater vehicles, AUVs, most studies have focused on the de- termination of static coefficients without considering the effects of the vehicle control surface deflection. This paper analyses the hydrodynamic efforts generated on an AUV considering the combined effects of the control surface deflection and the angle of attack using CFD software based on the Reynolds-averaged Navier–Stokes formulations. The CFD simulations are also independently conducted for the AUV bare hull and control surface to better identify their individual and interference efforts and to validate the simulations by comparing the experimental results obtained in a towing tank. Several simulations of the bare hull case were conducted to select the k –ω SST turbulent model with the viscosity approach that best predicts its hydrodynamic efforts. Mesh sensitivity analyses were conducted for all simulations. For the flow around the control surfaces, the CFD results were analysed according to two different methodologies, standard and nonlinear. The nonlinear regression methodology provides better results than the standard methodology does for predicting the stall at the control surface. The flow simulations have shown that the occurrence of the control surface stall depends on a linear relationship between the angle of attack and the control surface deflection. This type of information can be used in designing the vehicle’s autopilot system.