Convective Terms Research Articles

Reducing spatial discretization error on coarse CFD simulations using an openFOAM-embedded deep learning framework

AbstractWe propose a method for reducing the spatial discretization error of coarse computational fluid dynamics (CFD) problems by enhancing the quality of low-resolution simulations using deep learning. We feed the model with fine-grid data after projecting it to the coarse-grid discretization. We substitute the default differencing scheme for the convection term by a feed-forward neural network that interpolates velocities from cell centers to face values to produce velocities that approximate the down-sampled fine-grid data well. The deep learning framework incorporates the open-source CFD code OpenFOAM, resulting in an end-to-end differentiable model. We automatically differentiate the CFD physics using a discrete adjoint code version. We present a fast communication method between TensorFlow (Python) and OpenFOAM (c++) that accelerates the training process. We applied the model to the flow past a square cylinder problem, reducing the error from 120% to 25% in the velocity for simulations inside the training distribution compared to the traditional solver using an x8 coarser mesh. For simulations outside the training distribution, the error reduction in the velocities was about 50%. The training is affordable in terms of time and data samples since the architecture exploits the local features of the physics.

Adaptation of the Low Dissipation Low Dispersion Scheme for Reactive Multicomponent Flows on Unstructured Grids Using Density-Based Solvers

Abstract The low dissipation low dispersion (LD2) second-order accurate scheme is effective for scale-resolving simulations in finite volume computational fluid dynamics (CFD) solvers, thanks to its combination of a skew-symmetric split form for convective terms and matrix-valued artificial dissipation fluxes. However, reactive flow simulations face challenges due to steep gradients and varying gas properties. This study replaces the skew-symmetric scheme with the kinetic energy and entropy preserving (KEEP) scheme, utilizing quadratic and cubic split forms for convective terms, enhancing stability. The nonsmooth fluid interfaces in reactive flow simulations necessitate upwind fluxes for reactive scalars to limit total variation (TV), also requiring upwind fluxes for the mixture-dependent internal energy fluxes. Other convective terms use central discretizations from the KEEP scheme, leveraging LD2’s spatial reconstruction to minimize dispersive errors. Numerical assessments show this approach reduces spurious pressure oscillations in single and multicomponent flows. The absolute flux Jacobian for dissipation flux calculation is efficiently computed using an expanded Turkel’s approach for thermally perfect gas mixtures. Partial pressure derivatives are approximated when using the flamelet generated manifolds (FGM) combustion model. The proposed scheme is evaluated through scale-resolving simulations of the Cambridge burner flame SWB1 on an unstructured grid using the density-based solver TRACE, employing both finite rate chemistry (FRC) and FGM combustion models. Comparative analysis with the all-speed SLAU2 scheme shows the superior performance of the proposed scheme in handling turbulent reactive multicomponent flows.

On 24th-Order Multioperator Approximations in Schemes for Equations with Convective Terms

On 24th-Order Multioperator Approximations in Schemes for Equations with Convective Terms

Squeeze force of a Maxwell fluid between circular smooth surfaces with simple harmonic motion

The force and mechanical power required to maintain the simple harmonic motion (SHM) of the upper circular surface squeezing a viscoelastic fluid film is analyzed. The amplitude of the displacement of the upper surface is very small compared to the gap width as a function of time. The smoothness of the upper and lower surfaces is characterized by the slip model with two constant parameters, a slip length and a critical surface shear stress. The nonlinear convection terms in the momentum equation are neglected since the viscous forces dominate the inertial forces. The acceleration and deceleration terms are retained since the upper plate oscillates harmonically and the velocity in the fluid is strictly periodic. An exact solution of the governing equations is found as a function of the Deborah number, the Womersley number, the slip length, and the critical surface shear stress. A circular region without slip condition, bounded by a time-dependent radius, appears when the shear stress of the fluid does not exceed a critical surface shear stress. In addition, an annular region with slip up to the radius of the disk appears when the critical surface shear stress is exceeded. Our results show that viscoelastic and hydrophobic effects together with the Womersley number and a critical surface stress cause changes in the amplitude and phase lag of the waveform of the time-dependent radius and the force acting on the wall surface to maintain the SHM of the upper disk.

An efficient eigenvalue bounding method: CFL condition revisited

Direct and large-eddy simulations of turbulence are often solved using explicit temporal schemes. However, this imposes very small time-steps because the eigenvalues of the (linearized) dynamical system, re-scaled by the time-step, must lie inside the stability region. In practice, fast and accurate estimations of the spectral radii of both the discrete convective and diffusive terms are therefore needed. This is virtually always done using the so-called CFL condition. On the other hand, the large heterogeneity and complexity of modern supercomputing systems are nowadays hindering the efficient cross-platform portability of CFD codes. In this regard, our leitmotiv reads: relying on a minimal set of (algebraic) kernels is crucial for code portability and maintenance! In this context, this work focuses on the computation of eigenbounds for the above-mentioned convective and diffusive matrices which are needed to determine the time-step à la CFL. To do so, a new inexpensive method, that does not require to re-construct these time-dependent matrices, is proposed and tested. It just relies on a sparse-matrix vector product where only vectors change on time. Hence, both implementation in existing codes and cross-platform portability are straightforward. The effectiveness and robustness of the method are demonstrated for different test cases on both structured Cartesian and unstructured meshes. Finally, the method is combined with a self-adaptive temporal scheme, leading to significantly larger time-steps compared with other more conventional CFL-based approaches.

An artificial compressibility SAV finite element method for the time-dependent natural convection problem

In this article, we present, analyze, and test a first-order implicit-explicit type scheme based on the artificial compressibility (AC) method and the scalar auxiliary variable (SAV) approach for solving the time-dependent natural convection problem. The proposed method is a combination of a mixed finite element approximation for spatial discretization, the first-order backward Euler scheme for temporal discretization, and explicit treatment for nonlinear terms. It effectively decouples the velocity and pressure via artificial compressibility, thereby reducing computational complexity and execution time. Moreover, we use the SAV approach for the convective terms, it leads to the scheme is linear, and only require solving a sequence of linear differential equation with constant coefficients at each time step. It is proved that the solution of this method is bounded and satisfying the law of energy dissipation. The rigorous error estimate of velocity, pressure, and temperature are derived. Finally, some numerical tests are implemented to verify the theoretical analysis and illustrate the efficiency of the presented method.

Influence of thermal radiation and magnetic field on convective transfer within an inclined square enclosure filled with Cu-water nanofluid

The present work investigates the effects of thermal radiation and magnetic field on the laminar convective transfer process on copper-based nanofluid within an inclined square enclosure. The studied configuration is a rigid-walled cavity subjected to a horizontal temperature gradient where the left wall is maintained at a hot temperature Th and subjected to a magnetic field B0. In contrast, the right wall is kept at a cold temperature Tc. The horizontal walls are adiabatic. The dimensionless governing equations are solved using the finite difference method and the UPWIND scheme to solve the convective terms. Formulated using the stream function, vorticity, and temperature. Effects on mean Nusselt numbers are investigated at different Rayleigh numbers 104,105,106, Hartmann numbers HA=0to100, tilt angles γ=0,π/6,π/4andπ/3, nanofluid volume fractions 0%≤φ≤6% and different radiation parameters ( Rd=0,1,2,3 ). The findings show that the HTR increased 3.6 times by augmenting the Ra from 104 to 106. Also, increasing the Ha number from 0 to 100 causes an 83.16% reduction in the HTR. Considering the radiation mode of heat transfer causes an increase in HTR. The average Nusselt values increase by 59.72% for the concentration of 6% while increasing the radiation parameter from 0 to 3. On the other hand, an increase in volume fraction leads to a deterioration in mean Nusselt numbers (Numean).

Positive solutions of nonlinear elliptic equations involving unbounded variable exponents and convection term

Positive solutions of nonlinear elliptic equations involving unbounded variable exponents and convection term

A nonlocal convection–diffusion model with Gaussian‐type kernels and meshfree discretization

AbstractNonlocal models have demonstrated their indispensability in numerical simulations across a spectrum of critical domains, ranging from analyzing crack and fracture behavior in structural engineering to modeling anomalous diffusion phenomena in materials science and simulating convection processes in heterogeneous environments. In this study, we present a novel framework for constructing nonlocal convection–diffusion models using Gaussian‐type kernels. Our framework uniquely formulates the diffusion term by correlating the constant diffusion coefficient with the variance of the Gaussian kernel. Simultaneously, the convection term is defined by integrating the variable velocity field into the kernel as the expectation of a multivariate Gaussian distribution, facilitating a comprehensive representation of convective transport phenomena. We rigorously establish the well‐posedness of the proposed nonlocal model and derive a maximum principle to ensure its stability and reliability in numerical simulations. Furthermore, we develop a meshfree discretization scheme tailored for numerically simulating our model, designed to uphold both the discrete maximum principle and asymptotic compatibility. Through extensive numerical experiments, we validate the efficacy and versatility of our framework, demonstrating its superior performance compared to existing approaches.

On the Sub and Supersolution Method for Nonlinear Elliptic Equations with a Convective Term, in Orlicz Spaces

In this note we provide an overview of some existence (with sign information) and regularity results for differential equations, in which the method of sub and supersolutions plays an important role. We list some classical results and then we focus on the Dirichlet problem, for problems driven by a general differential operator, depending on (x,u,∇u), and with a convective term f. Our framework is that of Orlicz–Sobolev spaces. We also present several examples.

The localized RBF interpolation with its modifications for solving the incompressible two-phase fluid flows: A conservative Allen–Cahn–Navier–Stokes system

In this research work, we apply a numerical scheme based on the first-order time integration approach combined with the modifications of the meshless approximation for solving the conservative Allen–Cahn–Navier–Stokes equations. More precisely, we first utilize a first-order time discretization for the Navier–Stokes equations and the time-splitting technique of order one for the dynamics of the phase-field variable. Besides, we use the local interpolation based on the Matérn radial function for spatial discretization. We should solve a Poisson equation with the proper boundary conditions to have the divergence-free property during the numerical algorithm. Accordingly, the applied numerical procedure could not give a stable and accurate solution. Instead, we solve a regularization system in a discrete form. To prevent the instability of the numerical solution concerning the convection term, a biharmonic term with a small coefficient based on the high-order hyperviscosity formulation has been added, which has been approximated by a scalable interpolation based on the combination of polyharmonic spline with polynomials (known as the PHS+poly). The obtained full-discrete problem is solved using the biconjugate gradient stabilized method considering a proper preconditioner. We investigate the potency of the numerical scheme by presenting some simulations via uniform, hexagonal, and quasi-uniform nodes on rectangular and irregular domains. Besides, we have compared the proposed meshless method with the standard finite element method due to the used CPU time.

Toward large eddy simulation of shear-thinning liquid jets: A priori analysis of subgrid scale closures for multiphase flows

While direct numerical simulation (DNS) of multiphase flows has been the focus of many research investigations in recent years, large eddy simulation (LES) of multiphase flows remains a challenge. There is no standardized set of governing equations for multiphase LES. Different approaches and formulations have been discussed in the literature, each with its own advantages and disadvantages. In this paper, the conventional (non-weighted) filtering approach is compared with the density-weighted Favre filtering method by evaluating the subgrid scale (SGS) energy transfer for a simple test case of a shear-thinning droplet in air. The findings reveal that, unlike the Favre filtering approach, the conventional filtering method results in a notable amount of nonphysical backward scatter in the flow. Based on these results, the Favre filtering method appears preferable and is applied to the a priori analysis of shear-thinning liquid jets, where the viscosity has been modeled using the Carreau–Yasuda model. First, by explicitly filtering existing DNS data of shear-thinning jet breakup into stagnant air, the order of magnitude of different SGS terms is evaluated using the Favre filtering method. Consistent with earlier studies on Newtonian jets, the present study indicates that the diffusive term remains negligible, while the convective term plays a dominant role. Functional and structural models for the closure of the convective SGS term are assessed by means of a correlation analysis and an order of magnitude study. Existing structural models provide good results for both Newtonian and shear-thinning cases. Promising a posteriori model candidates are discussed.

An unconditionally stable scheme for the immersed boundary method with application in cardiac mechanics

The stability of the immersed boundary (IB) method is a challenge in simulating fluid–structure interaction problems, where time step constraints are significantly stricter than in pure fluid simulations. We propose a novel unconditionally stable scheme for the immersed boundary finite element (IBFE) method. The structure is handled implicitly and characterized by strain energy functions, rather than being modeled as fibers or membranes. Through energy estimate, we prove the unconditional stability of the fully discrete approximation in the absence of the convective term. In real simulations of cardiac mechanics problems, the time step is much larger, only limited by the Courant–Friedrichs–Lewy condition of the fluid. The novelty of this work lies in the combination of dual interpolation and distribution operators, the Jacobian-free Newton–Krylov method for solving nonlinear algebraic systems, and the semi-Lagrangian method for handling the convective term. To validate the effectiveness and accuracy of our approach, we present various benchmarks and conduct a quasi-static simulation of a three-dimensional real left ventricular model. We have shown that the numerical stability of our scheme is very robust even with much larger time step compared to conventional explicit IB methods. Our work paves the way for further works on efficient solvers of the IBFE method.

An artificial neural network model for recovering small-scale velocity in large-eddy simulation of isotropic turbulent flows

Small-scale motions in turbulent flows play a significant role in various small-scale processes, such as particle relative dispersion and collision, bubble or droplet deformation, and orientation dynamics of non-sphere particles. Recovering the small-scale flows that cannot be resolved in large eddy simulation (LES) is of great importance for such processes sensitive to the small-scale motions in turbulent flows. This study proposes a subgrid-scale model for recovering the small-scale turbulent velocity field based on the artificial neural network (ANN). The governing equations of small-scale turbulent velocity are linearized, and the pressure gradient and the nonlinear convection term are modeled with the aid of the ANN. Direct numerical simulation (DNS) and filtered direct numerical simulation (FDNS) provide the data required for training and validating the ANN. The large-scale velocity and velocity gradient tensor are selected as inputs for the ANN model. The linearized governing equations of small-scale turbulent velocity are numerically solved by coupling the large-scale flow field information. The results indicate that the model established by the ANN can accurately recover the small-scale velocity lost in FDNS due to filtering operation. With the ANN model, the flow fields at different Reynolds numbers agree well with the DNS results regarding velocity field statistics, flow field structures, turbulent energy spectra, and two-point, two-time Lagrangian correlation functions. This study demonstrates that the proposed ANN model can be applied to recovering the small-scale velocity field in the LES of isotropic turbulent flows at different Reynolds numbers.

Periodic wave solutions for a generalized reaction–convection–diffusion equation of high-order

In this paper, we investigate the uniqueness of periodic wave solutions for a generalized reaction–convection–diffusion equation with arbitrarily high-order reaction term or convection term. The main technique is to prove the monotonicity of the ratio of two Abelian integrals by a new criterion and Descartes’ rule of signs. A positive answer to a conjecture stated in Wei and Chen (2023) is given and some numeric simulations are carried out to illustrate the obtained theoretical results.

A linear second-order maximum bound principle preserving finite difference scheme for the generalized Allen–Cahn equation

In this paper, we consider the numerical method for generalized Allen–Cahn equation with nonlinear mobility and convection term. We propose a linear second-order finite difference scheme which preserves the discrete maximum bound principle (MBP). The scheme is discretized by stabilized Crank–Nicolson in time, upwind scheme for convection term and central-difference scheme for diffusion term. We show that the proposed scheme preserves the discrete MBP under some constraints on temporal step size and stabilizing parameter. Optimal L∞-error estimate is obtained for our scheme. Several numerical experiments are performed to validate our theoretical results.

CFD study of heat transfer in power‐law fluids over multiple corrugated circular cylinders in a heat exchanger

AbstractThe heat transfer in power‐law fluids across three corrugated circular cylinders placed in a triangular pitch arrangement is studied computationally in a confined channel. Continuity, momentum, and energy balance equations were solved using ANSYS FLUENT (Version 18.0). The flow is assumed to be steady, incompressible, two‐dimensional, and laminar. A square domain of side 300Dh is selected after a detailed domain study. An optimized grid with 98,187 cells is used in the study. The convergence criteria of 10−7 for the continuity, x‐momentum, and y‐momentum balances and 10−12 for the energy equation were used. Constant density and non‐Newtonian power‐law viscosity modules were used. The diffusive term is discretized using a central difference scheme. Convective terms are discretized using the Second‐Order Upwind scheme. Pressure–velocity coupling between continuity and momentum equations was implemented using the semi‐implicit method for pressure‐linked equation scheme. Streamlines show wake development behind the cylinders, which is very dominant at large ReN and n. Isotherm contours are cramped at higher values of ReN and PrN, implying higher heat transfer. Global parameters, like, Cd and Nu, are computed for the wide ranges of controlling dimensionless parameters, such as power‐law index (0.3 ≤ n ≤ 1.5), Reynolds (0.1 ≤ ReN ≤ 40), and Prandtl (0.72 ≤ PrN ≤ 500) numbers. The NuLocal plot attains a pitch near the corrugation of the surface due to abrupt changes in velocity and temperature gradients. Nu increases with ReN and/or PrN and decreases with n under ot herwise identical situations. Nu is correlated with pertinent parameters, namely, ReN, PrN, and n.

Shape optimization for a nonlinear elliptic problem related to thermal insulation

In this paper, we consider a minimization problem of a nonlinear functional I_{\beta,p}(D,\Omega) related to a thermal insulation problem with a convection term, where \Omega is a bounded connected open set in \R^{n} and D\subset\overline{\Omega} is a compact set. The Euler–Lagrange equation relative to I_{\beta,p} is a p -Laplace equation, 1<p<\infty , with a Robin boundary condition with parameter \beta>0 . The main aim is to study extremum problems for I_{\beta,p}(D,\Omega) , among domains D with given geometrical constraints and \Omega\setminus D of fixed thickness. In the planar case, we show that under perimeter constraint the disk maximizes I_{\beta,p} . In the n -dimensional case we restrict our analysis to convex sets showing that the same is true for the ball but under different geometrical constraints.

Numerical solution of convected wave equation in free field using artificial boundary method

AbstractIn this article, we propose two procedures focusing on the computation of the time‐dependent convected wave equation in a free field with a uniform background flow. Both procedures are based on a framework, expended from Du et al. (SIAM J. Sci. Comput. 40 (2018), A1430–A1445.), of constructing the Dirichlet‐to‐Dirichlet (DtD)‐type discrete absorbing boundary conditions (ABCs). The first procedure is dedicated to reducing the infinite problem into a finite problem by a direct application of the framework on the finite difference discretization of the convected wave equation. However, the presence of convection terms makes the stability analysis hard to implement, which motivates us to develop the second procedure. First, the convected wave equation is transformed into a standard wave equation by using the Prandtl‐Glauert‐Lorentz transformation. After that, we obtain the DtD‐type ABC by using the above framework, and on this basis, derive an equivalent Dirichlet‐to‐Neumann‐type ABCs, which makes stability and convergence analysis easy to be obtained by the classical energy method. The effectiveness and comparison of these two procedures are investigated through numerical experiments.

Existence of positive solutions for a class of singular elliptic problems with convection term and critical exponential growth

This paper uses the Galerkin method to investigate the existence of positive solution to a class of singular elliptic problems given by {−Δu=λ0uβ0+Λ0|∇u|γ0+f0(u)|x|α0+h0(x),u>0inΩ,u=0on∂Ω,\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document} $$\\begin{aligned} \ extstyle\\begin{cases} -\\Delta u= \\displaystyle \\frac {\\lambda _{0}}{u^{\\beta _{0}}} + \\Lambda _{0} |\ abla u|^{\\gamma _{0}}+ \\frac{f_{0}(u)}{|x|^{\\alpha _{0}}}+ h_{0}(x), \\ \\ u>0 \\ \\ \ ext{in} \\ \\Omega , \\\\ u=0 \\ \ ext{on} \\ \\ \\partial \\Omega , \\end{cases}\\displaystyle \\end{aligned}$$ \\end{document} where Ω⊂R2\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$\\Omega \\subset \\mathbb{R}^{2}$\\end{document} is a bounded smooth domain, 0<β0\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$0<\\beta _{0}$\\end{document}, γ0≤1\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$\\gamma _{0} \\leq 1$\\end{document}, α0∈[0,2)\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$\\alpha _{0} \\in [0,2)$\\end{document}, h0(x)≥0\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$h_{0}(x)\\geq 0$\\end{document}, h0≠0\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$h_{0}\ eq 0$\\end{document}, h0∈L∞(Ω)\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$h_{0}\\in L^{\\infty}(\\Omega )$\\end{document}, 0<∥h0∥∞<λ0<Λ0\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$0<\\|h_{0}\\|_{\\infty} < \\lambda _{0} < \\Lambda _{0}$\\end{document}, and f0\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$f_{0}$\\end{document} are continuous functions. More precisely, f0\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$f_{0}$\\end{document} has a critical exponential growth, that is, the nonlinearity behaves like exp(ϒ‾s2)\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$\\exp (\\overline{\\Upsilon}s^{2})$\\end{document} as |s|→∞\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$|s| \ o \\infty $\\end{document}, for some ϒ‾>0\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$\\overline{\\Upsilon}>0$\\end{document}.