208 resultados para Time Integration
em Biblioteca Digital da Produção Intelectual da Universidade de São Paulo (BDPI/USP)
Resumo:
The fluid flow over bodies with complex geometry has been the subject of research of many scientists and widely explored experimentally and numerically. The present study proposes an Eulerian Immersed Boundary Method for flows simulations over stationary or moving rigid bodies. The proposed method allows the use of Cartesians Meshes. Here, two-dimensional simulations of fluid flow over stationary and oscillating circular cylinders were used for verification and validation. Four different cases were explored: the flow over a stationary cylinder, the flow over a cylinder oscillating in the flow direction, the flow over a cylinder oscillating in the normal flow direction, and a cylinder with angular oscillation. The time integration was carried out by a classical 4th order Runge-Kutta scheme, with a time step of the same order of distance between two consecutive points in x direction. High-order compact finite difference schemes were used to calculate spatial derivatives. The drag and lift coefficients, the lock-in phenomenon and vorticity contour plots were used for the verification and validation of the proposed method. The extension of the current method allowing the study of a body with different geometry and three-dimensional simulations is straightforward. The results obtained show a good agreement with both numerical and experimental results, encouraging the use of the proposed method.
Resumo:
The performance of laser-induced breakdown spectrometry (LIBS) for the determination of Ba, Cd, Cr and Pb in toys has been evaluated by using a Nd:YAG laser operating at 1064 nm and an Echelle spectrometer with intensified charge-coupled device detector. Samples were purchased in different cities of Sao Paulo State market and analyzed directly without sample preparation. Laser-induced breakdown spectrometry experimental conditions (number of pulses, delay time. integration time gate and pulse energy) were optimized by using a Doehlert design. Laser-induced breakdown spectrometry signals correlated reasonably well with inductively coupled plasma optical emission spectrometry (ICP OES) concentrations after microwave-assisted acid digestion of selected samples. Thermal analysis was used for polymer identification and scanning electron microscopy to Visualize differences in crater geometry of different polymers employed for toy fabrication. Results indicate that laser-induced breakdown spectrometry can be proposed as a rapid screening method for investigation of potentially toxic elements in toys. The unique application of laser-induced breakdown spectrometry for identification of contaminants in successive layers of ink and polymer is also demonstrated. (C) 2009 Elsevier B.V. All rights reserved.
Resumo:
Following the approach developed for rods in Part 1 of this paper (Pimenta et al. in Comput. Mech. 42:715-732, 2008), this work presents a fully conserving algorithm for the integration of the equations of motion in nonlinear shell dynamics. We begin with a re-parameterization of the rotation field in terms of the so-called Rodrigues rotation vector, allowing for an extremely simple update of the rotational variables within the scheme. The weak form is constructed via non-orthogonal projection, the time-collocation of which ensures exact conservation of momentum and total energy in the absence of external forces. Appealing is the fact that general hyperelastic materials (and not only materials with quadratic potentials) are permitted in a totally consistent way. Spatial discretization is performed using the finite element method and the robust performance of the scheme is demonstrated by means of numerical examples.
Resumo:
A fully conserving algorithm is developed in this paper for the integration of the equations of motion in nonlinear rod dynamics. The starting point is a re-parameterization of the rotation field in terms of the so-called Rodrigues rotation vector, which results in an extremely simple update of the rotational variables. The weak form is constructed with a non-orthogonal projection corresponding to the application of the virtual power theorem. Together with an appropriate time-collocation, it ensures exact conservation of momentum and total energy in the absence of external forces. Appealing is the fact that nonlinear hyperelastic materials (and not only materials with quadratic potentials) are permitted without any prejudice on the conservation properties. Spatial discretization is performed via the finite element method and the performance of the scheme is assessed by means of several numerical simulations.
Resumo:
This work presents a numerical method suitable for the study of the development of internal boundary layers (IBL) and their characteristics for flows over various types of coastal cliffs. The IBL is an important meteorological occurrence for flows with surface roughness and topographical step changes. A two-dimensional flow program was used for this study. The governing equations were written using the vorticity-velocity formulation. The spatial derivatives were discretized by high-order compact finite differences schemes. The time integration was performed with a low storage fourth-order Runge-Kutta scheme. The coastal cliff (step) was specified through an immersed boundary method. The validation of the code was done by comparison of the results with experimental and observational data. The numerical simulations were carried out for different coastal cliff heights and inclinations. The results show that the predominant factors for the height of the IBL and its characteristics are the upstream velocity, and the height and form (inclination) of the coastal cliff. Copyright (C) 2010 John Wiley & Sons, Ltd.
Resumo:
We present a variable time step, fully adaptive in space, hybrid method for the accurate simulation of incompressible two-phase flows in the presence of surface tension in two dimensions. The method is based on the hybrid level set/front-tracking approach proposed in [H. D. Ceniceros and A. M. Roma, J. Comput. Phys., 205, 391400, 2005]. Geometric, interfacial quantities are computed from front-tracking via the immersed-boundary setting while the signed distance (level set) function, which is evaluated fast and to machine precision, is used as a fluid indicator. The surface tension force is obtained by employing the mixed Eulerian/Lagrangian representation introduced in [S. Shin, S. I. Abdel-Khalik, V. Daru and D. Juric, J. Comput. Phys., 203, 493-516, 2005] whose success for greatly reducing parasitic currents has been demonstrated. The use of our accurate fluid indicator together with effective Lagrangian marker control enhance this parasitic current reduction by several orders of magnitude. To resolve accurately and efficiently sharp gradients and salient flow features we employ dynamic, adaptive mesh refinements. This spatial adaption is used in concert with a dynamic control of the distribution of the Lagrangian nodes along the fluid interface and a variable time step, linearly implicit time integration scheme. We present numerical examples designed to test the capabilities and performance of the proposed approach as well as three applications: the long-time evolution of a fluid interface undergoing Rayleigh-Taylor instability, an example of bubble ascending dynamics, and a drop impacting on a free interface whose dynamics we compare with both existing numerical and experimental data.
Resumo:
Here, we study the stable integration of real time optimization (RTO) with model predictive control (MPC) in a three layer structure. The intermediate layer is a quadratic programming whose objective is to compute reachable targets to the MPC layer that lie at the minimum distance to the optimum set points that are produced by the RTO layer. The lower layer is an infinite horizon MPC with guaranteed stability with additional constraints that force the feasibility and convergence of the target calculation layer. It is also considered the case in which there is polytopic uncertainty in the steady state model considered in the target calculation. The dynamic part of the MPC model is also considered unknown but it is assumed to be represented by one of the models of a discrete set of models. The efficiency of the methods presented here is illustrated with the simulation of a low order system. (C) 2010 Elsevier Ltd. All rights reserved.
Resumo:
In this study the hypothesis that interceptive movements are controlled on the basis of expectancy of time to target arrival was tested. The study was conducted through assessment of temporal errors and kinematics of interceptive movements to a moving virtual target. Initial target velocity was kept unchanged in part of the trials, and in the others it was decreased 300 ms before the due time of target arrival at the interception position, increasing in 100 ms time to target arrival. Different probabilities of velocity decrease ranging from 25 to 100% were compared. The results revealed that while there were increasing errors between probabilities of 25 and 75% for unchanged target velocity, the opposite relationship was observed for target velocity decrease. Kinematic analysis indicated that movement timing adjustments to target velocity decrease were made online. These results support the conception that visuomotor integration in the interception of moving targets is mediated by an internal forward model whose weights can be flexibly adjusted according to expectancy of time to target arrival.
Resumo:
The main purpose of this paper is to present architecture of automated system that allows monitoring and tracking in real time (online) the possible occurrence of faults and electromagnetic transients observed in primary power distribution networks. Through the interconnection of this automated system to the utility operation center, it will be possible to provide an efficient tool that will assist in decisionmaking by the Operation Center. In short, the desired purpose aims to have all tools necessary to identify, almost instantaneously, the occurrence of faults and transient disturbances in the primary power distribution system, as well as to determine its respective origin and probable location. The compilations of results from the application of this automated system show that the developed techniques provide accurate results, identifying and locating several occurrences of faults observed in the distribution system.
Resumo:
This paper studies a simplified methodology to integrate the real time optimization (RTO) of a continuous system into the model predictive controller in the one layer strategy. The gradient of the economic objective function is included in the cost function of the controller. Optimal conditions of the process at steady state are searched through the use of a rigorous non-linear process model, while the trajectory to be followed is predicted with the use of a linear dynamic model, obtained through a plant step test. The main advantage of the proposed strategy is that the resulting control/optimization problem can still be solved with a quadratic programming routine at each sampling step. Simulation results show that the approach proposed may be comparable to the strategy that solves the full economic optimization problem inside the MPC controller where the resulting control problem becomes a non-linear programming problem with a much higher computer load. (C) 2010 Elsevier Ltd. All rights reserved.
Resumo:
This article deals with the efficiency of fractional integration parameter estimators. This study was based on Monte Carlo experiments involving simulated stochastic processes with integration orders in the range]-1,1[. The evaluated estimation methods were classified into two groups: heuristics and semiparametric/maximum likelihood (ML). The study revealed that the comparative efficiency of the estimators, measured by the lesser mean squared error, depends on the stationary/non-stationary and persistency/anti-persistency conditions of the series. The ML estimator was shown to be superior for stationary persistent processes; the wavelet spectrum-based estimators were better for non-stationary mean reversible and invertible anti-persistent processes; the weighted periodogram-based estimator was shown to be superior for non-invertible anti-persistent processes.
Resumo:
Our objective is to verify the modulatory effects of bromazepam on EEG theta absolute power when subjects were submitted to a visuomotor task (i.e., car driver task). Sample was composed of 14 students (9 males and 5 females), right handed, with ages varying between 23 and 42 years (mean = 32.5 +/- 9.5), absence of mental or physical impairments, no psychoactive or psychotropic substance use and no neuromuscular disorders (screened by a clinical examination). The results showed an interaction between condition and electrodes (p=0.034) in favor of F8 electrode compared with F7 in both experimental conditions (t-test; p=0.001). Additionally, main effects were observed for condition (p=0.001), period (p=0.001) and electrodes (p=0.031) in favor of F4 electrode compared with F3. In conclusion, Br 6 mg of bromazepam may interfere in sensorimotor processes in the task performance in an unpredictable scenario allowing that certain visuospatial factors were predominant. Therefore, the results may reflect that bromazepam effects influence the performance of the involved areas because of the acquisition and integration of sensory stimuli processes until the development of a motor behavior based on the same stimuli. (C) 2009 Elsevier Ireland Ltd. All rights reserved.
Resumo:
The study aimed to elucidate electrophysiological and cortical mechanisms involved in anticipatory actions when healthy subjects had to catch balls in free drop; specifically through quantitative electroencephalography (qEEG) alpha absolute power changes. Our hypothesis is that during the preparation of motor action (i.e., 2 s before ball`s drop) occurred integration among left medial frontal, left primary somatomotor and left posterior parietal cortices, showing a differentiated activity involving expectation, planning and preparedness. This hypothesis supports a lateralization of motor function. Although we contend that in right-handers the left hemisphere takes on a dominant role for the regulation of motor behavior. The sample was composed of 23 healthy subjects (13 male and 10 female), right handed, with ages varying between 25 and 40 years old (32.5 +/- 7.5), absence of mental and physical illness, right handed, and do not make use of any psychoactive or psychotropic substance at the time of the study. The experiment consisted of a task of catching balls in free drop. The three-way ANOVA analysis demonstrated an interaction between moment and position in left medial frontal cortex (F3 electrode), somatomotor cortex (C3 electrode) and posterior parietal cortex (P3 electrode: p < 0.001). Summarizing, through experimental task employed, it was possible to observe integration among frontal, central and parietal regions. This integration appears to be more predominant in expectation, planning and motor preparation. In this way, it established an absolute predominance of this mechanism under the left hemisphere. (C) 2008 Elsevier Ireland Ltd. All rights reserved.
Resumo:
Introduction. A fundamental aspect of planning future actions is the performance and control of motor tasks. This behaviour is done through sensory-motor integration. Aim. To explain the electrophysiological mechanisms in the cortex (modifications to the alpha band) that are involved in anticipatory actions when individuals have to catch a free-falling object. Subjects and methods. The sample was made up of 20 healthy subjects of both sexes (11 males and 9 females) with ages ranging between 25 and 40 years (32.5 +/- 7.5) who were free of mental or physical diseases (previous medical history); the subjects were right-handed (Edinburgh Inventory) and were not taking any psychoactive or psychotropic substances at the time of the study. The experiment consisted in a task in which subjects had to catch freely falling objects. The experiment was made up of six blocks of 15 tests, each of which lasted 2 minutes and 30 seconds before and two seconds after each ball was dropped. Results. An interaction of the factors moment and position was only observed for the right parietooccipital cortex, in the combination of electrodes P4-O2. Conclusion. These findings suggest that the right parietooccipital cortex plays an important role in increasing expectation and swiftness in the process of preparing for a motor task.
Resumo:
This study investigated the effects of the cement type and the water storage time on the push-out bond strength of a glass fiber post. Glass fiber posts (Fibrekor, Jeneric Pentron) were luted to post spaces using a self-cured resin cement (C&B Cement [CB]), a glass ionomer cement (Ketac Cem [KC]) or a resin-modified glass ionomer cement (GC FujiCEM [FC]) according to the manufacturers’ instructions. For each luting agent, the specimens were exposed to one of the following water storage times (n=5): 1 day (T1), 7 days (T7), 90 days (T90) and 180 days (T180). Push-out tests were performed after the storage times. Control specimens were not exposed to water storage, but subjected to the push-out test 10 min after post cementation. Data (in MPa) were analyzed by Kruskal-Wallis and Dunn`s test (α=0.05). Cement type and water storage time had a significant effect (p<0.05) on the push-out bond strength. CB showed significantly higher values of retention (p<0.05) than KC and FC, irrespective of the water storage time. Water storage increased significantly the push-out bond strength in T7 and T90, regardless of the cement type (p<0.05). The results showed that fiber posts luted to post spaces with the self-cured resin cement exhibited the best bonding performance throughout the 180-day water storage period. All cements exhibited a tendency to increase the bond strength after 7 and 90 days of water storage, decreasing thereafter.