An optimal control method for fluid structure interaction systems via adjoint boundary pressure

In recent year, in spite of the computational complexity, Fluid-structure interaction (FSI) problems have been widely studied due to their applicability in science and engineering. Fluid-structure interaction systems consist of one or more solid structures that deform by interacting with a surrounding fluid flow. FSI simulations evaluate the tensional state of the mechanical component and take into account the effects of the solid deformations on the motion of the interior fluids. The inverse FSI problem can be described as the achievement of a certain objective by changing some design parameters such as forces, boundary conditions and geometrical domain shapes. In this paper we would like to study the inverse FSI problem by using an optimal control approach. In particular we propose a pressure boundary optimal control method based on Lagrangian multipliers and adjoint variables. The objective is the minimization of a solid domain displacement matching functional obtained by finding the optimal pressure on the inlet boundary. The optimality system is derived from the first order necessary conditions by taking the Fréchet derivatives of the Lagrangian with respect to all the variables involved. The optimal solution is then obtained through a standard steepest descent algorithm applied to the optimality system. The approach presented in this work is general and could be used to assess other objective functionals and controls. In order to support the proposed approach we perform a few numerical tests where the fluid pressure on the domain inlet controls the displacement that occurs in a well defined region of the solid domain

Optimal pressure boundary control of steady multiscale fluid-structure interaction shell model derived from koiter equations

Fluid-structure interaction (FSI) problems are of great interest, due to their applicability in science and engineering. However, the coupling between large fluid domains and small moving solid walls presents numerous numerical difficulties and, in some configurations, where the thickness of the solid wall can be neglected, one can consider membrane models, which are derived from the Koiter shell equations with a reduction of the computational cost of the algorithm. With this assumption, the FSI simulation is reduced to the fluid equations on a moving mesh together with a Robin boundary condition that is imposed on the moving solid surface. In this manuscript, we are interested in the study of inverse FSI problems that aim to achieve an objective by changing some design parameters, such as forces, boundary conditions, or geometrical domain shapes. We study the inverse FSI membrane model by using an optimal control approach that is based on Lagrange multipliers and adjoint variables. In particular, we propose a pressure boundary optimal control with the purpose to control the solid deformation by changing the pressure on a fluid boundary. We report the results of some numerical tests for two-dimensional domains to demonstrate the feasibility and robustness of our method

A New Anisotropic Four-Parameter Turbulence Model for Low Prandtl Number Fluids

Due to their interesting thermal properties, liquid metals are widely studied for heat transfer applications where large heat fluxes occur. In the framework of the Reynolds-Averaged Navier– Stokes (RANS) approach, the Simple Gradient Diffusion Hypothesis (SGDH) and the Reynolds Analogy are almost universally invoked for the closure of the turbulent heat flux. Even though these assumptions can represent a reasonable compromise in a wide range of applications, they are not reliable when considering low Prandtl number fluids and/or buoyant flows. More advanced closure models for the turbulent heat flux are required to improve the accuracy of the RANS models dealing with low Prandtl number fluids. In this work, we propose an anisotropic four-parameter turbulence model. The closure of the Reynolds stress tensor and turbulent heat flux is gained through nonlinear models. Particular attention is given to the modeling of dynamical and thermal time scales. Numerical simulations of low Prandtl number fluids have been performed over the plane channel and backward-facing step configurations

Carbon monoxide poisoning in the home as an indicator of social integration of immigrant population

Verona is a town in Veneto Region (Italy) characterised by a high socio-economic status, but we register some health problems related to low income or social deprivation, such as unintentional carbon monoxide poisoning in private houses. Recently, indeed, there has been a rise in the number of accidents concerning immigrant populations. Since early 1990s, epidemiological informations in the matter weren?t complete and the number of cases was therefore probably underestimated. Presently no centralised system exists in Italy or in Veneto Region for recording the incidence of household CO poisoning. An ?Observatory on CO poisonings? was therefore activated in Verona Municipality in 1994 to overcome the problem of lacking data. The Observatory systematically collects data on cases of acute household CO poisoning occurring in the whole Verona area. These informations are used to assess public health interventions. During 13 years? activity 671 cases of CO poisoning and 22 deaths were recorded. The accidents were caused by the malfunctioning of home heating appliances. Statistics show a progressive increase in the number of not Italian involved in these episodes. Inspections revealed that immigrants often live in poor houses with unsafe systems or with extremely precarious heating systems. Data from this surveillance system on CO poisoning have shown that this is a public health problem involving immigrants much more often than local population. This system in real time supplies epidemiological and environmental data for improving public health intervention strategies. Such data collection systems could be useful to study other relevant problems of health disparities in lower socio-economical classes

SBO analysis of a generic PWR-900 with ASTEC and MELCOR codes

Abstract After the Fukushima accident, the interest of the public to nuclear safety has growth and the international technical nuclear community has increased his attention in the investigation and the characterization of Severe Accident (SA) scenarios. In order to simulate the different, complex and multi-physical phenomena involved in a SA, computational tools, known as SA codes, have been developed in the last decades. In order to give some insights on the modelling capabilities of these tools and the differences in the calculation results, also related to the user-effect, an analysis of an unmitigated Station Black Out (SBO) occurring in a generic Western three-loops PWR 900 MWe has been carried out by the authors in the framework of the NUGENIA TA-2 ASCOM project. The simulation results of ASTEC code (study carried out with ASTEC V2, IRSN all rights reserved, [2019]), developed by IRSN, and MELCOR 2.2 code, developed by SANDIA for USNRC, have been compared and analyzed. The SBO scenario considered takes into account the intervention of the accumulators as only accident mitigation strategy. Several figures of merits related to the thermal-hydraulic (e.g. primary pressure, cladding temperature, etc.) and to the core degradation (e.g. hydrogen production, etc.) have been considered to describe the accident evolution until the vessel failure, for the two codes comparison

Numerical validation of a κ-ω-κ θ -ω θ heat transfer turbulence model for heavy liquid metals

The correct prediction of heat transfer in turbulent flows is relevant in almost all industrial applications but many of the heat transfer models available in literature are validated only for ordinary fluids with Pr ≃ 1. In commercial Computational Fluid Dynamics codes only turbulence models with a constant turbulent Prandtl number of 0.85 — 0.9 are usually implemented but in heavy liquid metals with low Prandtl numbers it is well known that these models fail to reproduce correlations based on experimental data. In these fluids heat transfer is mainly due to molecular diffusion and the time scales of temperature and velocity fields are rather different, so simple turbulence models based on similarity between temperature and velocity cannot reproduce experimental correlations. In order to reproduce experimental results and Direct Numerical Simulation data obtained for fluids with Pr ≃ 0.025 we introduce a κ-ε-κ θ -ε θ turbulence model. This model, however, shows some numerical instabilities mainly due to the strong coupling between κ and ε on the walls. In order to fix this problem we reformulate the model into a new four parameter κ-ω-κ θ -ω θ where the dissipation rate on the wall is completely independent on the fluctuations. The model improves numerical stability and convergence. Numerical simulations in plane and channel geometries are reported and compared with experimental, Direct Numerical Simulation results and with results obtained with the κ-ε formulation, in order to show the model capabilities and validate the improved κ-ω model

The Ramazzini Institute 13-week study on glyphosate-based herbicides at human-equivalent dose in Sprague Dawley rats: Study design and first in-life endpoints evaluation

Background: Glyphosate-based herbicides (GBHs) are the most widely used pesticides worldwide, and glyphosate is the active ingredient of such herbicides, including the formulation known as Roundup. The massive and increasing use of GBHs results in not only the global burden of occupational exposures, but also increased exposure to the general population. The current pilot study represents the first phase of a long-term investigation of GBHs that we are conducting over the next 5 years. In this paper, we present the study design, the first evaluation of in vivo parameters and the determination of glyphosate and its major metabolite aminomethylphosphonic acid (AMPA) in urine. Methods: We exposed Sprague-Dawley (SD) rats orally via drinking water to a dose of glyphosate equivalent to the United States Acceptable Daily Intake (US ADI) of 1.75 mg/kg bw/day, defined as the chronic Reference Dose (cRfD) determined by the US EPA, starting from prenatal life, i.e. gestational day (GD) 6 of their mothers. One cohort was continuously dosed until sexual maturity (6-week cohort) and another cohort was continuously dosed until adulthood (13-week cohort). Here we present data on general toxicity and urinary concentrations of glyphosate and its major metabolite AMPA. Results: Survival, body weight, food and water consumption of the animals were not affected by the treatment with either glyphosate or Roundup. The concentration of both glyphosate and AMPA detected in the urine of SD rats treated with glyphosate were comparable to that observed in animals treated with Roundup, with an increase in relation to the duration of treatment. The majority of glyphosate was excreted unchanged. Urinary levels of the parent compound, glyphosate, were around 100-fold higher than the level of its metabolite, AMPA. Conclusions: Glyphosate concentrations in urine showed that most part of the administered dose was excreted as unchanged parent compound upon glyphosate and Roundup exposure, with an increasing pattern of glyphosate excreted in urine in relation to the duration of treatment. The adjuvants and the other substances present in Roundup did not seem to exert a major effect on the absorption and excretion of glyphosate. Our results demonstrate that urinary glyphosate is a more relevant marker of exposure than AMPA in the rodent model

L'accertamento della qualità di coltivatore diretto per l'esercizio della prelazione, Nota a Cass., 26 febbraio 1988, n. 2052

--

A FEM variational approach to the drop spreading over dry surfaces.

A numerical study on contact angle dynamics during spreading and recoil of droplets impacting orthogonally on various surfaces is presented. The Navier-Stokes equations are used where inertia, viscous, capillary and contact angle hysteresis phenomena act together to influence the contact motion. This numerical study is based on a finite element variational formulation of the problem and, since the problem is clearly singular, a limiting approach is discussed. The boundary singular condition is imposed as a penalty term and the problem is reformulated through an optimality system. The equation that relates the apparent contact angle and the contact line velocity is imposed as an external constraint. The evolution of the interface is solved over axisymmetric domains by using a new numerical surface tension representation and a new tracking marker technique, which conserves almost exactly the mass. The surface tension term is computed by using the fem test functions and therefore only the first derivative is required. This mathematical and numerical modeling of the wetting phenomena circumvents numerous difficulties due to the representation of surface tension, the non-integrable stress singularity at the moving contact line and the inability to describe in close form the velocity dependence of the dynamic contact line. In this work in order to see the capability of the numerical model the wetting properties of the target surfaces range from wettable to non-wettable and several numerical simulations with partial and complete rebound are investigated. A simple model for the apparent contact angle is used in order to provide an insight to the dynamic behavior of the apparent contact angle and its dependence on contact line velocity