A study of a numerical solution of the steady two dimensions Navier-Stokes equations in a constricted channel problem by a compact fourth order method

We present a numerical solution for the steady 2D Navier-Stokes equations using a fourth order compact-type method. The geometry of the problem is a constricted symmetric channel, where the boundary can be varied, via a parameter, from a smooth constriction to one possessing a very sharp but smooth corner allowing us to analyse the behaviour of the errors when the solution is smooth or near singular. The set of non-linear equations is solved by the Newton method. Results have been obtained for Reynolds number up to 500. Estimates of the errors incurred have shown that the results are accurate and better than those of the corresponding second order method. (C) 2002 Elsevier B.V. All rights reserved.

Fourth-order method for solving the Navier-Stokes equations in a constricting channel

A fourth-order numerical method for solving the Navier-Stokes equations in streamfunction/vorticity formulation on a two-dimensional non-uniform orthogonal grid has been tested on the fluid flow in a constricted symmetric channel. The family of grids is generated algebraically using a conformal transformation followed by a non-uniform stretching of the mesh cells in which the shape of the channel boundary can vary from a smooth constriction to one which one possesses a very sharp but smooth corner. The generality of the grids allows the use of long channels upstream and downstream as well as having a refined grid near the sharp corner. Derivatives in the governing equations are replaced by fourth-order central differences and the vorticity is eliminated, either before or after the discretization, to form a wide difference molecule for the streamfunction. Extra boundary conditions, necessary for wide-molecule methods, are supplied by a procedure proposed by Henshaw et al. The ensuing set of non-linear equations is solved using Newton iteration. Results have been obtained for Reynolds numbers up to 250 for three constrictions, the first being smooth, the second having a moderately sharp corner and the third with a very sharp corner. Estimates of the error incurred show that the results are very accurate and substantially better than those of the corresponding second-order method. The observed order of the method has been shown to be close to four, demonstrating that the method is genuinely fourth-order. © 1977 John Wiley & Sons, Ltd.

Analysis of the accumulated energy by thermal stratification of liquid inside storage tanks with porous layers

A mathematical model was developed in order to study the behavior of thermal stratification of liquid in a typical storage tank with porous medium. The model employs a transient stream function-vorticity formulation to predict the development of stream function and temperature fields in a charging process. Parameters analyzed include Biot, Darcy, Reynolds and Richardson numbers, position, and the thickness of the porous medium. The results show the influence of these physical parameters that should be considered for a good design of storage tanks with thermal stratification.

A generalized alternating-direction implicit scheme for incompressible magnetohydrodynamic viscous flows at low magnetic Reynolds number

This paper presents numerical simulations of incompressible fluid flows in the presence of a magnetic field at low magnetic Reynolds number. The equations governing the flow are the Navier-Stokes equations of fluid motion coupled with Maxwell's equations of electromagnetics. The study of fluid flows under the influence of a magnetic field and with no free electric charges or electric fields is known as magnetohydrodynamics. The magnetohydrodynamics approximation is considered for the formulation of the non-dimensional problem and for the characterization of similarity parameters. A finite-difference technique is used to discretize the equations. In particular, an extension of the generalized Peaceman and Rachford alternating-direction implicit (ADI) scheme for simulating two-dimensional fluid flows is presented. The discretized conservation equations are solved in stream function-vorticity formulation. We compare the ADI and generalized ADI schemes, and show that the latter is more efficient in simulating low Reynolds number and magnetic Reynolds number problems. Numerical results demonstrating the applicability of this technique are also presented. The simulation of incompressible magneto hydrodynamic fluid flows is illustrated by numerical solution for two-dimensional cases. (c) 2007 Elsevier B.V. All rights reserved.

Aplicação da técnica de elementos finitos na estratificação de armazenadores térmicos com barreiras de meios porosos

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

A regime-dependent balanced control variable based on potential vorticity

In this paper it is argued that rotational wind is not the best choice of leading control variable for variational data assimilation, and an alternative is suggested and tested. A rotational wind parameter is used in most global variational assimilation systems as a pragmatic way of approximately representing the balanced component of the assimilation increments. In effect, rotational wind is treated as a proxy for potential vorticity, but one that it is potentially not a good choice in flow regimes characterised by small Burger number. This paper reports on an alternative set of control variables which are based around potential vorticity. This gives rise to a new formulation of the background error covariances for the Met Office's variational assimilation system, which leads to flow dependency. It uses similar balance relationships to traditional schemes, but recognises the existence of unbalanced rotational wind which is used with a new anti-balance relationship. The new scheme is described and its performance is evaluated and compared to a traditional scheme using a sample of diagnostics.

Comment on ``Piecewise potential vorticity inversion: Elementary tests''

Egger (2008) constructs some idealised experiments to test the usefulness of piecewise potential vorticity inversion (PPVI) in the diagnosis of Rossby wave dynamics and baroclinic development. He concludes that, ``PPVI does not help us to understand the dynamics of linear Rossby waves. It provides local tendencies of the streamfunction which are unrelated to the true ones. The same way, the motion of baroclinic waves in shear flow cannot be understood by using PPVI. Moreover, the effect of boundary temperatures as determined by PPVI is unrelated to the flow evolution.'' He goes further in arguing that we should not consider velocities as ``induced'' by PV anomalies defined by carving up the global domain. However, these conclusions partly reflect the limitations of his idealised experiments and the manner in which the PV components were partitioned from one another.

Steady periodic water waves with constant vorticity: regularity and local bifurcation

This paper studies periodic traveling gravity waves at the free surface of water in a flow of constant vorticity over a flat bed. Using conformal mappings the free-boundary problem is transformed into a quasilinear pseudodifferential equation for a periodic function of one variable. The new formulation leads to a regularity result and, by use of bifurcation theory, to the existence of waves of small amplitude even in the presence of stagnation points in the flow.

Some experiments with high order compact methods using a computer algebra software - Part II (non-uniform grid)

We generalize a procedure proposed by Mancera and Hunt [P.F.A. Mancera, R. Hunt, Some experiments with high order compact methods using a computer algebra software-Part 1, Appl. Math. Comput., in press, doi: 10.1016/j.amc.2005.05.015] for obtaining a compact fourth-order method to the steady 2D Navier-Stokes equations in the streamfunction formulation-vorticity using the computer algebra system Maple, which includes conformal mappings and non-uniform grids. To analyse the procedure we have solved a constricted stepped channel problem, where a fine grid is placed near the re-entrant corner by transformation of the independent variables. (c) 2006 Elsevier B.V. All rights reserved.

Effect of mesh distortion on the accuracy of high order vorticity-velocity CFD approaches

The numerical solution of the incompressible Navier-Stokes equations offers an alternative to experimental analysis of fluid-structure interaction (FSI). We would save a lot of time and effort and help cut back on costs, if we are able to accurately model systems by these numerical solutions. These advantages are even more obvious when considering huge structures like bridges, high rise buildings or even wind turbine blades with diameters as large as 200 meters. The modeling of such processes, however, involves complex multiphysics problems along with complex geometries. This thesis focuses on a novel vorticity-velocity formulation called the Kinematic Laplacian Equation (KLE) to solve the incompressible Navier-stokes equations for such FSI problems. This scheme allows for the implementation of robust adaptive ordinary differential equations (ODE) time integration schemes, allowing us to tackle each problem as a separate module. The current algortihm for the KLE uses an unstructured quadrilateral mesh, formed by dividing each triangle of an unstructured triangular mesh into three quadrilaterals for spatial discretization. This research deals with determining a suitable measure of mesh quality based on the physics of the problems being tackled. This is followed by exploring methods to improve the quality of quadrilateral elements obtained from the triangles and thereby improving the overall mesh quality. A series of numerical experiments were designed and conducted for this purpose and the results obtained were tested on different geometries with varying degrees of mesh density.

Effect of high order interpolation in the stability and efficiency of the time-integration process in vorticity-velocity CFD algorithms

The numerical solution of the incompressible Navier-Stokes Equations offers an effective alternative to the experimental analysis of Fluid-Structure interaction i.e. dynamical coupling between a fluid and a solid which otherwise is very complex, time consuming and very expensive. To have a method which can accurately model these types of mechanical systems by numerical solutions becomes a great option, since these advantages are even more obvious when considering huge structures like bridges, high rise buildings, or even wind turbine blades with diameters as large as 200 meters. The modeling of such processes, however, involves complex multiphysics problems along with complex geometries. This thesis focuses on a novel vorticity-velocity formulation called the KLE to solve the incompressible Navier-stokes equations for such FSI problems. This scheme allows for the implementation of robust adaptive ODE time integration schemes and thus allows us to tackle the various multiphysics problems as separate modules. The current algorithm for KLE employs a structured or unstructured mesh for spatial discretization and it allows the use of a self-adaptive or fixed time step ODE solver while dealing with unsteady problems. This research deals with the analysis of the effects of the Courant-Friedrichs-Lewy (CFL) condition for KLE when applied to unsteady Stoke’s problem. The objective is to conduct a numerical analysis for stability and, hence, for convergence. Our results confirmthat the time step ∆t is constrained by the CFL-like condition ∆t ≤ const. hα, where h denotes the variable that represents spatial discretization.

Lawsonia Inermis-mediated Synthesis Of Silver Nanoparticles: Activity Against Human Pathogenic Fungi And Bacteria With Special Reference To Formulation Of An Antimicrobial Nanogel.

Lawsonia inermis mediated synthesis of silver nanoparticles (Ag-NPs) and its efficacy against Candida albicans, Microsporum canis, Propioniabacterium acne and Trichophyton mentagrophytes is reported. A two-step mechanism has been proposed for bioreduction and formation of an intermediate complex leading to the synthesis of capped nanoparticles was developed. In addition, antimicrobial gel for M. canis and T. mentagrophytes was also formulated. Ag-NPs were synthesized by challenging the leaft extract of L. inermis with 1 mM AgNO₃. The Ag-NPs were characterized by Ultraviolet-Visible (UV-Vis) spectrophotometer and Fourier transform infrared spectroscopy (FTIR). Transmission electron microscopy (TEM), nanoparticle tracking and analysis sytem (NTA) and zeta potential was measured to detect the size of Ag-NPs. The antimicrobial activity of Ag-NPs was evaluated by disc diffusion method against the test organisms. Thus these Ag-NPs may prove as a better candidate drug due to their biogenic nature. Moreover, Ag-NPs may be an answer to the drug-resistant microorganisms.

Pharmacokinetic And Pharmacodynamic Evaluation Of A Nanotechnological Topical Formulation Of Lidocaine/prilocaine (nanorap) In Healthy Volunteers.

Nanorap is a new nanotechnological formulation for topical anesthesia composed of lidocaine (2.5%) and prilocaine (2.5%). The present study evaluated the pharmacokinetics (PK) of Nanorap. For the determination of lidocaine and prilocaine in human plasma a new method using high-performance liquid-chromatography coupled to tandem mass spectrometry was developed. Nanorap pharmacodynamic (PD) and its physical proprieties were also evaluated. Nanorap was administered by topical application of 2g to healthy volunteers and blood samples were collected for the PK analysis. The drugs were extracted from plasma by liquid-liquid extraction with ether/hexane (80/20, v/v). The chromatography separation was performed on a Genesis C18 analytical column 4 µm (100 x 2.1 mm i.d.) with a mobile phase of methanol/acetonitrile/water (40/30/30, for lidocaine, and 50/30/20, for prilocaine, v/v/v) + 2 mM of ammonium acetate and ropivacaine as internal standard. The drugs were quantified using a mass spectrometer with an electrospray source in the ESI positive mode (ES+) configured for multiple reaction monitoring. The PD of Nanorap was evaluated with the use of a visual analogue scale. Nanorap was characterized by cryofracture. The chromatography run time was 5.5 min for lidocaine and 3.3 min for prilocaine and the lower limit of quantification was 0.05 ng/mL for both drugs. Mean Cmax was 6.62 and 1.72 ng/mL for lidocaine and prilocaine, respectively. Median Tmax was 6.5 hours for both drugs. Nanocapsules had a mean size of 88nm and mean drug association of 92.5% and 89% for lidocaine and prilocaine, respectively. The PD study showed that Nanorap has a sufficient analgesic effect (>30% reduction in pain) after 10 minutes of application. A new simple, selective and sensitive method for determination of lidocaine and prilocaine in human plasma was developed. Nanorap generated safe plasma levels of the drugs and satisfactory analgesic effect.